Polymorphic pseudogenes in the human genome - a comprehensive assessment.

BACKGROUND: Over the past decade, variations of the coding portion of the human genome have become increasingly evident. In this study, we focus on polymorphic pseudogenes, a unique and relatively unexplored type of pseudogene whose inactivating mutations have not yet been fixed in the human genome at the global population level. Thus, polymorphic pseudogenes are characterized by the presence in the population of both coding alleles and non-coding alleles originating from Loss-of-Function (LoF) mutations. These alleles can be found both in heterozygosity and in homozygosity in different human populations and thus represent pseudogenes that have not yet been fixed in the population. RESULTS: A methodical cross-population analysis of 232 polymorphic pseudogenes, including 35 new examples, reveals that human olfactory signalling, drug metabolism and immunity are among the systems most impacted by the variable presence of LoF variants at high frequencies. Within this dataset, a total of 179 genes presented polymorphic LoF variants in all analysed populations. Transcriptome and proteome analysis confirmed that although these genes may harbour LoF alleles, when the coding allele is present, the gene remains active and can play a functional role in various metabolic pathways, including drug/xenobiotic metabolism and immunity. The observation that many polymorphic pseudogenes are members of multigene families argues that genetic redundancy may play a key role in compensating for the inactivation of one paralogue. CONCLUSIONS: The distribution, expression and integration of cellular/biological networks in relation to human polymorphic pseudogenes, provide novel insights into the architecture of the human genome and the dynamics of gene gain and loss with likely functional impact.

Unveiling intraspecific diversity and evolutionary dynamics of the foodborne pathogen Bacillus paranthracis through high-quality pan-genome analysis.

Understanding the evolutionary dynamics of foodborne pathogens throughout host-associated habitats is of utmost importance. Bacterial pan-genomes, as dynamic entities, are strongly influenced by ecological lifestyles. As a phenotypically diverse species in the Bacillus cereus group, Bacillus paranthracis is recognized as an emerging foodborne pathogen and a probiotic simultaneously. This poorly understood species is a suitable study model for adaptive pan-genome evolution. In this study, we determined the biogeographic distribution, abundance, genetic diversity, and genotypic profiles of key genetic elements of B. paranthracis. Metagenomic read recruitment analyses demonstrated that B. paranthracis members are globally distributed and abundant in host-associated habitats. A high-quality pan-genome of B. paranthracis was subsequently constructed to analyze the evolutionary dynamics involved in ecological adaptation comprehensively. The open pan-genome indicated a flexible gene repertoire with extensive genetic diversity. Significant divergences in the phylogenetic relationships, functional enrichment, and degree of selective pressure between the different components demonstrated different evolutionary dynamics between the core and accessory genomes driven by ecological forces. Purifying selection and gene loss are the main signatures of evolutionary dynamics in B. paranthracis pan-genome. The plasticity of the accessory genome is characterized by horizontal gene transfer (HGT), massive gene losses, and weak purifying or positive selection, which might contribute to niche-specific adaptation. In contrast, although the core genome dominantly undergoes purifying selection, its association with HGT and positively selected mutations indicates its potential role in ecological diversification. Furthermore, host fitness-related dynamics are characterized by the loss of secondary metabolite biosynthesis gene clusters (BGCs) and CAZyme-encoding genes and the acquisition of antimicrobial resistance (AMR) and virulence genes via HGT. This study offers a case study of pan-genome evolution to investigate the ecological adaptations reflected by biogeographical characteristics, thereby advancing the understanding of intraspecific diversity and evolutionary dynamics of foodborne pathogens.

Simulating Weak Attacks in a New Duplication-Divergence Model with Node Loss.

A better understanding of protein-protein interaction (PPI) networks representing physical interactions between proteins could be beneficial for evolutionary insights as well as for practical applications such as drug development. As a statistical model for PPI networks, duplication-divergence models have been proposed, but they suffer from resulting in either very sparse networks in which most of the proteins are isolated, or in networks which are much denser than what is usually observed, having almost no isolated proteins. Moreover, in real networks, where a gene codes a protein, gene loss may occur. The loss of nodes has not been captured in duplication-divergence models to date. Here, we introduce a new duplication-divergence model which includes node loss. This mechanism results in networks in which the proportion of isolated proteins can take on values which are strictly between 0 and 1. To understand this new model, we apply strong and weak attacks to networks from duplication-divergence models with and without node loss, and compare the results to those obtained when carrying out similar attacks on two real PPI networks of E. coli and of S. cerevisiae. We find that the new model more closely reflects the damage caused by strong and weak attacks found in the PPI networks.

The Change Rate of the Fbxl21 Gene and the Amino Acid Composition of Its Protein Correlate with the Species-Specific Lifespan in Placental Mammals.

This article proposes a methodology for establishing a relationship between the change rate of a given gene (relative to a given taxon) together with the amino acid composition of the proteins encoded by this gene and the traits of the species containing this gene. The methodology is illustrated based on the mammalian genes responsible for regulating the circadian rhythms that underlie a number of human disorders, particularly those associated with aging. The methods used are statistical and bioinformatic ones. A systematic search for orthologues, pseudogenes, and gene losses was performed using our previously developed methods. It is demonstrated that the least conserved Fbxl21 gene in the Euarchontoglires superorder exhibits a statistically significant connection of genomic characteristics (the median of dN/dS for a gene relative to all the other orthologous genes of a taxon, as well as the preference or avoidance of certain amino acids in its protein) with species-specific lifespan and body weight. In contrast, no such connection is observed for Fbxl21 in the Laurasiatheria superorder. This study goes beyond the protein-coding genes, since the accumulation of amino acid substitutions in the course of evolution leads to pseudogenization and even gene loss, although the relationship between the genomic characteristics and the species traits is still preserved. The proposed methodology is illustrated using the examples of circadian rhythm genes and proteins in placental mammals, e.g., longevity is connected with the rate of Fbxl21 gene change, pseudogenization or gene loss, and specific amino acid substitutions (e.g., asparagine at the 19th position of the CRY-binding domain) in the protein encoded by this gene.

Cytonuclear evolution in fully heterotrophic plants: lifestyles and gene function determine scenarios.

BACKGROUND: Evidence shows that full mycoheterotrophs and holoparasites often have reduced plastid genomes with rampant gene loss, elevated substitution rates, and deeply altered to conventional evolution in mitochondrial genomes, but mechanisms of cytonuclear evolution is unknown. Endoparasitic Sapria himalayana and mycoheterotrophic Gastrodia and Platanthera guangdongensis represent different heterotrophic types, providing a basis to illustrate cytonuclear evolution. Here, we focused on nuclear-encoded plastid / mitochondrial (N-pt / mt) -targeting protein complexes, including caseinolytic protease (ClpP), ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCo), oxidative phosphorylation system (OXPHOS), DNA recombination, replication, and repair (DNA-RRR) system, and pentatricopeptide repeat (PPR) proteins, to identify evolutionary drivers for cytonuclear interaction. RESULTS: The severity of gene loss of N-pt PPR and pt-RRR genes was positively associated with increased degree of heterotrophy in full mycoheterotrophs and S. himalayana, while N-mt PPR and mt-RRR genes were retained. Substitution rates of organellar and nuclear genes encoding N-pt/mt subunits in protein complexes were evaluated, cytonuclear coevolution was identified in S. himalayana, whereas disproportionate rates of evolution were observed in the OXPHOS complex in full mycoheterotrophs, only slight accelerations in substitution rates were identified in N-mt genes of full mycoheterotrophs. CONCLUSIONS: Nuclear compensatory evolution was identified in protein complexes encoded by plastid and N-pt genes. Selection shaping codon preferences, functional constraint, mt-RRR gene regulation, and post-transcriptional regulation of PPR genes all facilitate mito-nuclear evolution. Our study enriches our understanding of genomic coevolution scenarios in fully heterotrophic plants.

Convergent reductive evolution in bee-associated lactic acid bacteria.

Distantly related organisms may evolve similar traits when exposed to similar environments or engaging in certain lifestyles. Several members of the Lactobacillaceae [lactic acid bacteria (LAB)] family are frequently isolated from the floral niche, mostly from bees and flowers. In some floral LAB species (henceforth referred to as bee-associated LAB), distinctive genomic (e.g., genome reduction) and phenotypic (e.g., preference for fructose over glucose or fructophily) features were recently documented. These features are found across distantly related species, raising the hypothesis that specific genomic and phenotypic traits evolved convergently during adaptation to the floral environment. To test this hypothesis, we examined representative genomes of 369 species of bee-associated and non-bee-associated LAB. Phylogenomic analysis unveiled seven independent ecological shifts toward the bee environment in LAB. In these species, we observed significant reductions of genome size, gene repertoire, and GC content. Using machine leaning, we could distinguish bee-associated from non-bee-associated species with 94% accuracy, based on the absence of genes involved in metabolism, osmotic stress, or DNA repair. Moreover, we found that the most important genes for the machine learning classifier were seemingly lost, independently, in multiple bee-associated lineages. One of these genes, acetaldehyde-alcohol dehydrogenase (adhE), encodes a bifunctional aldehyde-alcohol dehydrogenase which has been associated with the evolution of fructophily, a rare phenotypic trait that is pervasive across bee-associated LAB species. These results suggest that the independent evolution of distinctive phenotypes in bee-associated LAB has been largely driven by independent losses of the same sets of genes.IMPORTANCESeveral LAB species are intimately associated with bees and exhibit unique biochemical properties with potential for food applications and honeybee health. Using a machine learning-based approach, our study shows that adaptation of LAB to the bee environment was accompanied by a distinctive genomic trajectory deeply shaped by gene loss. Several of these gene losses occurred independently in distantly related species and are linked to some of their unique biotechnologically relevant traits, such as the preference for fructose over glucose (fructophily). This study underscores the potential of machine learning in identifying fingerprints of adaptation and detecting instances of convergent evolution. Furthermore, it sheds light onto the genomic and phenotypic particularities of bee-associated bacteria, thereby deepening the understanding of their positive impact on honeybee health.

Loss of Lateral suppressor gene is associated with evolution of root nodule symbiosis in Leguminosae.

BACKGROUND: Root nodule symbiosis (RNS) is a fascinating evolutionary event. Given that limited genes conferring the evolution of RNS in Leguminosae have been functionally validated, the genetic basis of the evolution of RNS remains largely unknown. Identifying the genes involved in the evolution of RNS will help to reveal the mystery. RESULTS: Here, we investigate the gene loss event during the evolution of RNS in Leguminosae through phylogenomic and synteny analyses in 48 species including 16 Leguminosae species. We reveal that loss of the Lateral suppressor gene, a member of the GRAS-domain protein family, is associated with the evolution of RNS in Leguminosae. Ectopic expression of the Lateral suppressor (Ls) gene from tomato and its homolog MONOCULM 1 (MOC1) and Os7 from rice in soybean and Medicago truncatula result in almost completely lost nodulation capability. Further investigation shows that Lateral suppressor protein, Ls, MOC1, and Os7 might function through an interaction with NODULATION SIGNALING PATHWAY 2 (NSP2) and CYCLOPS to repress the transcription of NODULE INCEPTION (NIN) to inhibit the nodulation in Leguminosae. Additionally, we find that the cathepsin H (CTSH), a conserved protein, could interact with Lateral suppressor protein, Ls, MOC1, and Os7 and affect the nodulation. CONCLUSIONS: This study sheds light on uncovering the genetic basis of the evolution of RNS in Leguminosae and suggests that gene loss plays an essential role.

The highly repetitive genome of Myxobolus rasmusseni, an emerging myxozoan parasite of fathead minnows.

Myxozoans are a monophyletic taxon of approximately 2,400 described species of parasites from the phylum Cnidaria. The recent focus on their negative impacts on fisheries, on their evolution from free-living ancestors, and on their emergence into new fish host populations has stressed the critical need for genomic resources for this parasitic group. Here, we describe the genome assembly and annotation of Myxobolus rasmusseni, an emerging parasite of fathead minnows in Alberta, Canada. The assembly is 174.6 Mb in size, 68% of which is made up of repetitive elements, making it one of the most repetitive animal genomes sequenced to date. Through comparisons to other myxozoans, we show that widespread gene loss, a known phenomenon of this group of parasites, is consistent with closely related species. Additionally, we assembled the M. rasmusseni mitochondrial genome which is nearly twice the size of the typical animal mitochondrial genome yet contains only five of the canonical mitochondrial protein-coding genes and open reading frames not found in other myxozoans. These results add to our understanding of the gene- and genome-level diversity observed in myxozoans.

Nematode genome announcement: A chromosome-scale genome assembly for the Pristionchus pacificus reference mapping strain PS1843.

Pristionchus pacificus is a free-living nematode that shares many features with Caenorhabditis elegans, such as its short generation time and hermaphroditism, but also exhibits novel traits, i.e., a mouth-form dimorphism that enables predation. The availability of various genetic tools and genomic resources make it a powerful model organism for comparative studies. Here, we present an updated genome of the P. pacificus strain PS1843 (Washington) that is most widely used for genetic analysis. Assembly of PacBio reads together with reference-guided scaffolding resulted in a chromosome-scale genome spanning 171Mb for the PS1843 strain. Whole genome alignments between the P. pacificus PS1843 genome and the genome of the P. pacificus reference strain PS312 (California) revealed megabase-sized regions on chromosomes III, IV, and X that explain the majority of genome size difference between both strains. The improved PS1843 genome will be useful for future forward genetic studies and evolutionary genomic comparisons at the intra-species level.

Bacterial cross-feeding can promote gene retention by lowering gene expression costs.

Gene loss is expected in microbial communities when the benefit of obtaining a biosynthetic precursor from a neighbor via cross-feeding outweighs the cost of retaining a biosynthetic gene. However, gene cost primarily comes from expression, and many biosynthetic genes are only expressed when needed. Thus, one can conversely expect cross-feeding to repress biosynthetic gene expression and promote gene retention by lowering gene cost. Here we examined long-term bacterial cocultures pairing Escherichia coli and Rhodopseudomonas palustris for evidence of gene loss or retention in response to cross-feeding of non-essential adenine. Although R. palustris continued to externalize adenine in long-term cultures, E. coli did not accumulate mutations in purine synthesis genes, even after 700 generations. E. coli purine synthesis gene expression was low in coculture, suggesting that gene repression removed selective pressure for gene loss. In support of this explanation, R. palustris also had low transcript levels for iron-scavenging siderophore genes in coculture, likely because E. coli facilitated iron acquisition by R. palustris. R. palustris siderophore gene mutations were correspondingly rare in long-term cocultures but were prevalent in monocultures where transcript levels were high. Our data suggests that cross-feeding does not always drive gene loss, but can instead promote gene retention by repressing costly expression.

Azooxanthellate Palythoa (Cnidaria: Anthozoa) Genomes Reveal Toxin-related Gene Clusters and Loss of Neuronal Genes in Hexacorals.

Zoantharia is an order among the Hexacorallia (Anthozoa: Cnidaria), and includes at least 300 species. Previously reported genomes from scleractinian corals and actiniarian sea anemones have illuminated part of the hexacorallian diversification. However, little is known about zoantharian genomes and the early evolution of hexacorals. To explore genome evolution in this group of hexacorals, here, we report de novo genome assemblies of the zoantharians Palythoa mizigama (Pmiz) and Palythoa umbrosa (Pumb), both of which are members of the family Sphenopidae, and uniquely live in comparatively dark coral reef caves without symbiotic Symbiodiniaceae dinoflagellates. Draft genomes generated from ultra-low input PacBio sequencing totaled 373 and 319 Mbp for Pmiz and Pumb, respectively. Protein-coding genes were predicted in each genome, totaling 30,394 in Pmiz and 24,800 in Pumb, with each set having ∼93% BUSCO completeness. Comparative genomic analyses identified 3,036 conserved gene families, which were found in all analyzed hexacoral genomes. Some of the genes related to toxins, chitin degradation, and prostaglandin biosynthesis were expanded in these two Palythoa genomes and many of which aligned tandemly. Extensive gene family loss was not detected in the Palythoa lineage and five of ten putatively lost gene families likely had neuronal function, suggesting biased gene loss in Palythoa. In conclusion, our comparative analyses demonstrate evolutionary conservation of gene families in the Palythoa lineage from the common ancestor of hexacorals. Restricted loss of gene families may imply that lost neuronal functions were effective for environmental adaptation in these two Palythoa species.

Alterations of pleiotropic neuropeptide-receptor gene couples in Cetacea.

BACKGROUND: Habitat transitions have considerable consequences in organism homeostasis, as they require the adjustment of several concurrent physiological compartments to maintain stability and adapt to a changing environment. Within the range of molecules with a crucial role in the regulation of different physiological processes, neuropeptides are key agents. Here, we examined the coding status of several neuropeptides and their receptors with pleiotropic activity in Cetacea. RESULTS: Analysis of 202 mammalian genomes, including 41 species of Cetacea, exposed an intricate mutational landscape compatible with gene sequence modification and loss. Specifically for Cetacea, in the 12 genes analysed we have determined patterns of loss ranging from species-specific disruptive mutations (e.g. neuropeptide FF-amide peptide precursor; NPFF) to complete erosion of the gene across the cetacean stem lineage (e.g. somatostatin receptor 4; SSTR4). CONCLUSIONS: Impairment of some of these neuromodulators may have contributed to the unique energetic metabolism, circadian rhythmicity and diving response displayed by this group of iconic mammals.

Horizontal gene transfer from chloroplast to mitochondria of seagrasses in the yellow-Bohai seas.

Seagrasses are ideal for studying plant adaptation to marine environments. In this study, the mitochondrial (mt) and chloroplast (cp) genomes of Ruppia sinensis were sequenced. The results showed an extensive gene loss in seagrasses, including a complete loss of cp-rpl19 genes in Zosteraceae, most cp-ndh genes in Hydrocharitaceae, and mt-rpl and mt-rps genes in all seagrasses, except for the mt-rpl16 gene in Phyllospadix iwatensis. Notably, most ribosomal protein genes were lost in the mt and cp genomes. The deleted cp genes were not transferred to the mt genomes through horizontal gene transfer. Additionally, a significant DNA transfer between seagrass organelles was found, with the mt genomes of Zostera containing numerous sequences from the cp genome. Rearrangement analyses revealed an unreported inversion of the cp genome in R. sinensis. Moreover, four positively selected genes (atp8, nad5, atp4, and ccmFn) and five variable regions (matR, atp4, atp8, rps7, and ccmFn) were identified.

Plastome evolution in Santalales involves relaxed selection prior to loss of ndh genes and major boundary shifts of the inverted repeat.

BACKGROUND AND AIMS: Biological aspects of haustorial parasitism have significant effects on the configuration of the plastid genome. Approximately half the diversity of haustorial parasites belongs to the order Santalales, where a clearer picture of plastome evolution in relation to parasitism is starting to emerge. However, in previous studies of plastome evolution there is still a notable under-representation of members from non-parasitic and deep-branching hemiparasitic lineages, limiting evolutionary inference around the time of transition to a parasitic lifestyle. To expand taxon sampling relevant to this transition we therefore targeted three families of non-parasites (Erythropalaceae, Strombosiaceae, and Coulaceae), two families of root-feeding hemiparasites (Ximeniaceae and Olacaceae), and two families of uncertain parasitic status (Aptandraceae and Octoknemaceae). With data from these lineages we aimed to explore plastome evolution in relation to evolution of parasitism. METHODS: From 29 new samples we sequenced and annotated plastomes and the nuclear ribosomal cistron. We examined phylogenetic patterns, plastome evolution, and patterns of relaxed or intensified selection in plastid genes. Available transcriptome data were analyzed to investigate potential transfer of infA to the nuclear genome. RESULTS: Phylogenetic relationships indicate a single functional loss of all plastid ndh genes (ndhA-K) in a clade formed by confirmed parasites and Aptandraceae, and the loss coincides with major size and boundary shifts of the inverted repeat (IR) region. Depending on an autotrophic or heterotrophic lifestyle in Aptandraceae, plastome changes are either correlated with or predate evolution of parasitism. Phylogenetic patterns also indicate repeated loss of infA from the plastome, and based on presence of transcribed sequences with presequences corresponding to thylakoid luminal transit peptides, we infer that the genes were transferred to the nuclear genome. CONCLUSIONS: Except for the loss of the ndh complex, relatively few genes have been lost from the plastome in deep-branching root parasites in Santalales. Prior to loss of the ndh genes, they show signs of relaxed selection indicative of their dispensability. To firmly establish a potential correlation between ndh gene loss, plastome instability and evolution of parasitism, it is pertinent to refute or confirm a parasitic lifestyle all Santalales clades.

Convergent reductive evolution in bee-associated lactic acid bacteria.

Distantly related organisms may evolve similar traits when exposed to similar environments or engaging in certain lifestyles. Several members of the Lactobacillaceae (LAB) family are frequently isolated from the floral niche, mostly from bees and flowers. In some floral LAB species (henceforth referred to as bee-associated), distinctive genomic (e.g., genome reduction) and phenotypic (e.g., preference for fructose over glucose or fructophily) features were recently documented. These features are found across distantly related species, raising the hypothesis that specific genomic and phenotypic traits evolved convergently during adaptation to the floral environment. To test this hypothesis, we examined representative genomes of 369 species of bee-associated and non-bee-associated LAB. Phylogenomic analysis unveiled seven independent ecological shifts towards the floral niche in LAB. In these bee-associated LAB, we observed pervasive, significant reductions of genome size, gene repertoire, and GC content. Using machine leaning, we could distinguish bee-associated from non-bee-associated species with 94% accuracy, based on the absence of genes involved in metabolism, osmotic stress, or DNA repair. Moreover, we found that the most important genes for the machine learning classifier were seemingly lost, independently, in multiple bee-associated lineages. One of these genes, adhE, encodes a bifunctional aldehyde-alcohol dehydrogenase associated with the evolution of fructophily, a rare phenotypic trait that was recently identified in many floral LAB species. These results suggest that the independent evolution of distinctive phenotypes in bee-associated LAB has been largely driven by independent loss of the same set of genes. Importance: Several lactic acid bacteria (LAB) species are intimately associated with bees and exhibit unique biochemical properties with potential for food applications and honeybee health. Using a machine-learning based approach, our study shows that adaptation of LAB to the bee environment was accompanied by a distinctive genomic trajectory deeply shaped by gene loss. Several of these gene losses occurred independently in distantly related species and are linked to some of their unique biotechnologically relevant traits, such as the preference of fructose over glucose (fructophily). This study underscores the potential of machine learning in identifying fingerprints of adaptation and detecting instances of convergent evolution. Furthermore, it sheds light onto the genomic and phenotypic particularities of bee-associated bacteria, thereby deepening the understanding of their positive impact on honeybee health.

The rate of W chromosome degeneration across multiple avian neo-sex chromosomes.

When sex chromosomes evolve recombination suppression, the sex-limited chromosome (Y/W) commonly degenerate by losing functional genes. The rate of Y/W degeneration is believed to slow down over time as the most essential genes are maintained by purifying selection, but supporting data are scarce especially for ZW systems. Here, we study W degeneration in Sylvioidea songbirds where multiple autosomal translocations to the sex chromosomes, and multiple recombination suppression events causing separate evolutionary strata, have occurred during the last ~ 28.1-4.5 million years (Myr). We show that the translocated regions have maintained 68.3-97.7% of their original gene content, compared to only 4.2% on the much older ancestral W chromosome. By mapping W gene losses onto a dated phylogeny, we estimate an average gene loss rate of 1.0% per Myr, with only moderate variation between four independent lineages. Consistent with previous studies, evolutionarily constrained and haploinsufficient genes were preferentially maintained on W. However, the gene loss rate did not show any consistent association with strata age or with the number of W genes at strata formation. Our study provides a unique account on the pace of W gene loss and reinforces the significance of purifying selection in maintaining essential genes on sex chromosomes.

Genomic decoding of Theobroma grandiflorum (cupuassu) at chromosomal scale: evolutionary insights for horticultural innovation.

BACKGROUND: Theobroma grandiflorum (Malvaceae), known as cupuassu, is a tree indigenous to the Amazon basin, valued for its large fruits and seed pulp, contributing notably to the Amazonian bioeconomy. The seed pulp is utilized in desserts and beverages, and its seed butter is used in cosmetics. Here, we present the sequenced telomere-to-telomere genome of cupuassu, disclosing its genomic structure, evolutionary features, and phylogenetic relationships within the Malvaceae family. FINDINGS: The cupuassu genome spans 423 Mb, encodes 31,381 genes distributed in 10 chromosomes, and exhibits approximately 65% gene synteny with the Theobroma cacao genome, reflecting a conserved evolutionary history, albeit punctuated with unique genomic variations. The main changes are pronounced by bursts of long-terminal repeat retrotransposons at postspecies divergence, retrocopied and singleton genes, and gene families displaying distinctive patterns of expansion and contraction. Furthermore, positively selected genes are evident, particularly among retained and dispersed tandem and proximal duplicated genes associated with general fruit and seed traits and defense mechanisms, supporting the hypothesis of potential episodes of subfunctionalization and neofunctionalization following duplication, as well as impact from distinct domestication process. These genomic variations may underpin the differences observed in fruit and seed morphology, ripening, and disease resistance between cupuassu and the other Malvaceae species. CONCLUSIONS: The cupuassu genome offers a foundational resource for both breeding improvement and conservation biology, yielding insights into the evolution and diversity within the genus Theobroma.

Horizontally transferred glycoside hydrolase 26 may aid hemipteran insects in plant tissue digestion.

Glycoside hydrolases are enzymes that break down complex carbohydrates into simple sugars by catalyzing the hydrolysis of glycosidic bonds. There have been multiple instances of adaptive horizontal gene transfer of genes belonging to various glycoside hydrolase families from microbes to insects, as glycoside hydrolases can metabolize constituents of the carbohydrate-rich plant cell wall. In this study, we characterize the horizontal transfer of a gene from the glycoside hydrolase family 26 (GH26) from bacteria to insects of the order Hemiptera. Our phylogenies trace the horizontal gene transfer to the common ancestor of the superfamilies Pentatomoidea and Lygaeoidea, which include stink bugs and seed bugs. After horizontal transfer, the gene was assimilated into the insect genome as indicated by the gain of an intron, and a eukaryotic signal peptide. Subsequently, the gene has undergone independent losses and expansions in copy number in multiple lineages, suggesting an adaptive role of GH26s in some insects. Finally, we measured tissue-level gene expression of multiple stink bugs and the large milkweed bug using publicly available RNA-seq datasets. We found that the GH26 genes are highly expressed in tissues associated with plant digestion, especially in the principal salivary glands of the stink bugs. Our results are consistent with the hypothesis that this horizontally transferred GH26 was co-opted by the insect to aid in plant tissue digestion and that this HGT event was likely adaptive.

Gene loss and relaxed selection of plaat1 in vertebrates adapted to low-light environments.

Gene loss is an important mechanism for evolution in low-light or cave environments where visual adaptations often involve a reduction or loss of eyesight. The plaat gene family encodes phospholipases essential for the degradation of organelles in the lens of the eye. These phospholipases translocate to damaged organelle membranes, inducing them to rupture. This rupture is required for lens transparency and is essential for developing a functioning eye. Plaat3 is thought to be responsible for this role in mammals, while plaat1 is thought to be responsible in other vertebrates. We used a macroevolutionary approach and comparative genomics to examine the origin, loss, synteny and selection of plaat1 across bony fishes and tetrapods. We showed that plaat1 (probably ancestral to all bony fish + tetrapods) has been lost in squamates and is significantly degraded in lineages of low-visual-acuity and blind mammals and fishes. Our findings suggest that plaat1 is important for visual acuity across bony vertebrates, and that its loss through relaxed selection and pseudogenization may have played a role in the repeated evolution of visual systems in low-light environments. Our study sheds light on the importance of gene-loss in trait evolution and provides insights into the mechanisms underlying visual acuity in low-light environments.

The genomic evolution of visual opsin genes in amphibians.

Among tetrapod (terrestrial) vertebrates, amphibians remain more closely tied to an amphibious lifestyle than amniotes, and their visual opsin genes may be adapted to this lifestyle. Previous studies have discussed physiological, morphological, and molecular changes in the evolution of amphibian vision. We predicted the locations of the visual opsin genes, their neighboring genes, and the tuning sites of the visual opsins, in 39 amphibian genomes. We found that all of the examined genomes lacked the Rh2 gene. The caecilian genomes have further lost the SWS1 and SWS2 genes; only the Rh1 and LWS genes were retained. The loss of the SWS1 and SWS2 genes in caecilians may be correlated with their cryptic lifestyles. The opsin gene syntenies were predicted to be highly similar to those of other bony vertebrates. Moreover, dual syntenies were identified in allotetraploid Xenopus laevis and X. borealis. Tuning site analysis showed that only some Caudata species might have UV vision. In addition, the S164A that occurred several times in LWS evolution might either functionally compensate for the Rh2 gene loss or fine-tuning visual adaptation. Our study provides the first genomic evidence for a caecilian LWS gene and a genomic viewpoint of visual opsin genes by reviewing the gains and losses of visual opsin genes, the rearrangement of syntenies, and the alteration of spectral tuning in the course of amphibians' evolution.