Concerted and birth-and-death evolution of multigene families

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

An evolutionary mechanism to assimilate new nutrient sensors into the mTORC1 pathway

Animals must sense and respond to nutrient availability in their local niche. This task is coordinated in part by the mTOR complex 1 (mTORC1) pathway, which regulates growth and metabolism in response to nutrients1-5. In mammals, mTORC1 senses specific amino acids through specialized sensors that act through the upstream GATOR1/2 signaling hub6-8. To reconcile the conserved architecture of the mTORC1 pathway with the diversity of environments that animals can occupy, we hypothesized that the pathway might maintain plasticity by evolving distinct nutrient sensors in different metazoan phyla1,9,10. Whether such customization occurs-and how the mTORC1 pathway might capture new nutrient inputs-is not known. Here, we identify the Drosophila melanogaster protein Unmet expectations (Unmet, formerly CG11596) as a species-restricted nutrient sensor and trace its incorporation into the mTORC1 pathway. Upon methionine starvation, Unmet binds to the fly GATOR2 complex to inhibit dTORC1. S-adenosylmethionine (SAM), a proxy for methionine availability, directly relieves this inhibition. Unmet expression is elevated in the ovary, a methionine-sensitive niche11, and flies lacking Unmet fail to maintain the integrity of the female germline under methionine restriction. By monitoring the evolutionary history of the Unmet-GATOR2 interaction, we show that the GATOR2 complex evolved rapidly in Dipterans to recruit and repurpose an independent methyltransferase as a SAM sensor. Thus, the modular architecture of the mTORC1 pathway allows it to co-opt preexisting enzymes and expand its nutrient sensing capabilities, revealing a mechanism for conferring evolvability on an otherwise highly conserved system.

Role of F-box E3-ubiquitin ligases in plant development and stress responses

KEY MESSAGE

F-box E3-ubiquitin ligases regulate critical biological processes in plant development and stress responses. Future research could elucidate why and how plants have acquired a large number of F-box genes. The ubiquitin-proteasome system (UPS) is a predominant regulatory mechanism employed by plants to maintain the protein turnover in the cells and involves the interplay of three classes of enzymes, E1 (ubiquitin-activating), E2 (ubiquitin-conjugating), and E3 ligases. The diverse and most prominent protein family among eukaryotes, F-box proteins, are a vital component of the multi-subunit SCF (Skp1-Cullin 1-F-box) complex among E3 ligases. Several F-box proteins with multifarious functions in different plant systems have evolved rapidly over time within closely related species, but only a small part has been characterized. We need to advance our understanding of substrate-recognition regulation and the involvement of F-box proteins in biological processes and environmental adaptation. This review presents a background of E3 ligases with particular emphasis on the F-box proteins, their structural assembly, and their mechanism of action during substrate recognition. We discuss how the F-box proteins regulate and participate in the signaling mechanisms of plant development and environmental responses. We highlight an urgent need for research on the molecular basis of the F-box E3-ubiquitin ligases in plant physiology, systems biology, and biotechnology. Further, the developments and outlooks of the potential technologies targeting the E3-ubiquitin ligases for developing crop improvement strategies have been discussed.

Completing the TRB family: newly characterized members show ancient evolutionary origins and distinct localization, yet similar interactions

Telomere repeat binding proteins (TRBs) belong to a family of proteins possessing a Myb-like domain which binds to telomeric repeats. Three members of this family (TRB1, TRB2, TRB3) from Arabidopsis thaliana have already been described as associated with terminal telomeric repeats (telomeres) or short interstitial telomeric repeats in gene promoters (telo-boxes). They are also known to interact with several protein complexes: telomerase, Polycomb repressive complex 2 (PRC2) E(z) subunits and the PEAT complex (PWOs-EPCRs-ARIDs-TRBs). Here we characterize two novel members of the TRB family (TRB4 and TRB5). Our wide phylogenetic analyses have shown that TRB proteins evolved in the plant kingdom after the transition to a terrestrial habitat in Streptophyta, and consequently TRBs diversified in seed plants. TRB4-5 share common TRB motifs while differing in several others and seem to have an earlier phylogenetic origin than TRB1-3. Their common Myb-like domains bind long arrays of telomeric repeats in vitro, and we have determined the minimal recognition motif of all TRBs as one telo-box. Our data indicate that despite the distinct localization patterns of TRB1-3 and TRB4-5 in situ, all members of TRB family mutually interact and also bind to telomerase/PRC2/PEAT complexes. Additionally, we have detected novel interactions between TRB4-5 and EMF2 and VRN2, which are Su(z)12 subunits of PRC2.

Genome-wide analyses of Glutathione S-transferase gene family and expression profiling under deltamethrin exposure in non-biting midge Propsilocerus akamusi

Glutathione S-transferases (GSTs) are major enzymes in detoxification phase II, and have been functioned in resistance to various insecticides or oxidative stress. Herein, we selected the non-biting midge, Propsilocerus akamusi, widespread in Asian aquatic ecosystems, to uncover the gene location, structure, and phylogenetics relationship of GSTs at genome scale first time. Thirty-three cytosolic and four microsomal GST genes were identified and located on the four chromosomes. The cytosolic GSTs involved in the eight subclasses and five GST genes were unclassified. The expansion of GST genes in P. akamusi experienced duplication events on the delta, theta, xi, iota, and unclassified subclasses. The RNA-Seq analyses and RT-qPCR validation showed that the expression of PaGSTt2 gene is significantly elevated, with deltamethrin concentration increasing. The tertiary structure of PaGSTt2 enzyme was reconstructed, which was different from the other theta gene in the active site. In addition, the GST genes of six chironomids were first described based on the assembled genomes to explore the difference of those in the adaptation to kinds of environments. The GST frame for P. akmusi and its expression profiles provide valuable resources to understand their role in insecticide resistance of this species, as well as those of other biting midges.

Four ciliate-specific expansion events occurred during actin gene family evolution of eukaryotes

Actin gene family is a divergent and ancient eukaryotic cellular cytoskeletal gene family, and participates in many essential cellular processes. Ciliated protists offer us an excellent opportunity to investigate gene family evolution, since their gene families evolved faster in ciliates than in other eukaryotes. Nonetheless, actin gene family is well studied in few model ciliate species but little is known about its evolutionary patterns in ciliates. Here, we analyzed the evolutionary pattern of eukaryotic actin gene family based on genomes/transcriptomes of 36 species covering ten ciliate classes, as well as those of nine non-ciliate eukaryotic species. Results showed: (1) Except for conventional actins and actin-related proteins (Arps) shared by various eukaryotes, at least four ciliate-specific subfamilies occurred during evolution of actin gene family. Expansions of Act2 and ArpC were supposed to have happen in the ciliate common ancestor, while expansions of ActI and ActII may have occurred in the ancestor of Armophorea, Muranotrichea, and Spirotrichea. (2) The number of actin isoforms varied greatly among ciliate species. Environmental adaptability, whole genome duplication (WGD) or segmental duplication events, distinct spatial and temporal patterns of expression might play driving forces for the increasement of isoform numbers. (3) The 'birth and death' model of evolution could explain the evolution of actin gene family in ciliates. And actin genes have been generally under strong negative selection to maintain protein structures and physiological functions. Collectively, we provided meaningful information for understanding the evolution of eukaryotic actin gene family.

Functional cross-species conservation of guanylate-binding proteins in innate immunity

Guanylate binding proteins (GBPs) represent an evolutionary ancient protein family widely distributed among eukaryotes. They are interferon (IFN)-inducible guanosine triphosphatases that belong to the dynamin superfamily. GBPs are known to have a major role in the cell-autonomous innate immune response against bacterial, parasitic and viral infections and are also involved in inflammasome activation. Evolutionary studies depicted that GBPs present a pattern of gain and loss of genes in each family with several genes pseudogenized and some genes more divergent, indicative for the birth-and-death evolution process. Most species harbor large GBP gene clusters encoding multiple paralogs. Previous functional studies mainly focused on mouse and human GBPs, but more data are becoming available, broadening the understanding of this multifunctional protein family. In this review, we will provide new insights and give a broad overview about GBP evolution, conservation and their roles in all studied species, including plants, invertebrates and vertebrates, revealing how far the described features of GBPs can be transferred to other species.

Functional characterization of sex pheromone receptors in Spodoptera frugiperda, S. exigua, and S. litura moths

Moths possess an extremely sensitive and diverse sex pheromone processing system, in which pheromone receptors (PRs) are essential to ensure communication between mating partners. Functional properties of some PRs are conserved among species, which is important for reproduction. However, functional differentiation has occurred in some homologous PR genes, which may drive species divergence. Here, using genome analysis, 17 PR genes were identified from Spodoptera frugiperda, S. exigua, and S. litura, which belong to 6 homologous groups (odorant receptor [OR]6, 11, 13, 16, 56, and 62); of which 6 PR genes (OR6, OR11, OR13, OR16, OR56, and OR62) were identified in S. frugiperda and S. exigua, and 5 PR genes were identified in S. litura, excluding OR62. Using heterologous expression in Xenopus oocytes, we characterized the functions of PR orthologs including OR6, OR56, and OR62, which have not been clarified in previous studies. OR6 orthologs were specifically tuned to (Z,E)-9,12-tetradecadienyl acetate (Z9,E12-14:OAc), and OR62 orthologs were robustly tuned to Z7-12:OAc in S. frugiperda and S. exigua. The optimal ligand for OR56 was Z7-12:OAc in S. frugiperda, but responses were minimal in S. exigua and S. litura. In addition, SfruOR6 was male antennae-specific, whereas SfruOR56 and SfruOR62 were male antennae-biased. Our study further clarified the functional properties of PRs in 3 Spodoptera moth species, providing a comprehensive understanding of the mechanisms of intraspecific communication and interspecific isolation in Spodoptera.

Odorant-Binding Proteins and Chemosensory Proteins in Spodoptera frugiperda: From Genome-Wide Identification and Developmental Stage-Related Expression Analysis to the Perception of Host Plant Odors, Sex Pheromones, and Insecticides

Spodoptera frugiperda is a worldwide generalist pest with remarkable adaptations to environments and stresses, including developmental stage-related behavioral and physiological adaptations, such as diverse feeding preferences, mate seeking, and pesticide resistance. Insects’ odorant-binding proteins (OBPs) and chemosensory proteins (CSPs) are essential for the chemical recognition during behavioral responses or other physiological processes. The genome-wide identification and the gene expression patterns of all these identified OBPs and CSPs across developmental stage-related S. frugiperda have not been reported. Here, we screened for genome-wide SfruOBPs and SfruCSPs, and analyzed the gene expression patterns of SfruOBPs and SfruCSPs repertoires across all developmental stages and sexes. We found 33 OBPs and 22 CSPs in the S. frugiperda genome. The majority of the SfruOBP genes were most highly expressed in the adult male or female stages, while more SfruCSP genes were highly expressed in the larval or egg stages, indicating their function complementation. The gene expression patterns of SfruOBPs and SfruCSPs revealed strong correlations with their respective phylogenic trees, indicating a correlation between function and evolution. In addition, we analyzed the chemical-competitive binding of a widely expressed protein, SfruOBP31, to host plant odorants, sex pheromones, and insecticides. Further ligands binding assay revealed a broad functional related binding spectrum of SfruOBP31 to host plant odorants, sex pheromones, and insecticides, suggesting its potential function in food, mate seeking, and pesticide resistance. These results provide guidance for future research on the development of behavioral regulators of S. frugiperda or other environmentally friendly pest-control strategies.

Analysis of sequence diversity in Plasmodium falciparum glutamic acid-rich protein (PfGARP), an asexual blood stage vaccine candidate

Glutamic acid-rich protein of Plasmodium falciparum (PfGARP) binds to erythrocyte band 3 and may enhance cytoadherence of infected erythrocytes. Naturally acquired anti-PfGARP antibodies could confer protection against high parasitemia and severe symptoms. While whole genome sequencing analysis has suggested high conservation in this locus, little is known about repeat polymorphism in this vaccine candidate antigen. Direct sequencing was performed from the PCR-amplified complete PfGARP gene of 80 clinical isolates from four malaria endemic provinces in Thailand and an isolate from a Guinean patient. Publicly available complete coding sequences of this locus were included for comparative analysis. Six complex repeat (RI-RVI) and two homopolymeric glutamic acid repeat (E1 and E2) domains were identified in PfGARP. The erythrocyte band 3-binding ligand in domain RIV and the epitope for mAB7899 antibody eliciting in vitro parasite killing property were perfectly conserved across isolates. Repeat lengths in domains RIII and E1-RVI-E2 seemed to be correlated with parasite density of the patients. Sequence variation in PfGARP exhibited genetic differentiation across most endemic areas of Thailand. Phylogenetic tree inferred from this locus has shown that most Thai isolates formed closely related lineages, suggesting local expansion/contractions of repeat-encoding regions. Positive selection was observed in non-repeat region preceding domain RII which corresponded to a helper T cell epitope predicted to be recognized by a common HLA class II among Thai population. Predicted linear B cell epitopes were identified in both repeat and non-repeat domains. Besides length variation in some repeat domains, sequence conservation in non-repeat regions and almost all predicted immunogenic epitopes have suggested that PfGARP-derived vaccine may largely elicit strain-transcending immunity.

Extensive MHC class IIβ diversity across multiple loci in the small-spotted catshark (Scyliorhinus canicula)

The major histocompatibility complex (MHC) is a multigene family responsible for pathogen detection, and initiation of adaptive immune responses. Duplication, natural selection, recombination, and their resulting high functional genetic diversity spread across several duplicated loci are the main hallmarks of the MHC. Although these features were described in several jawed vertebrate lineages, a detailed MHC IIβ characterization at the population level is still lacking for chondrichthyans (chimaeras, rays and sharks), i.e. the most basal lineage to possess an MHC-based adaptive immune system. We used the small-spotted catshark (Scyliorhinus canicula, Carcharhiniformes) as a case-study species to characterize MHC IIβ diversity using complementary molecular tools, including publicly available genome and transcriptome datasets, and a newly developed high-throughput Illumina sequencing protocol. We identified three MHC IIβ loci within the same genomic region, all of which are expressed in different tissues. Genetic screening of the exon 2 in 41 individuals of S. canicula from a single population revealed high levels of sequence diversity, evidence for positive selection, and footprints of recombination. Moreover, the results also suggest the presence of copy number variation in MHC IIβ genes. Thus, the small-spotted catshark exhibits characteristics of functional MHC IIβ genes typically observed in other jawed vertebrates.

Past, Present, and Future of Naturally Occurring Antimicrobials Related to Snake Venoms

This review focuses on proteins and peptides with antimicrobial activity because these biopolymers can be useful in the fight against infectious diseases and to overcome the critical problem of microbial resistance to antibiotics. In fact, snakes show the highest diversification among reptiles, surviving in various environments; their innate immunity is similar to mammals and the response of their plasma to bacteria and fungi has been explored mainly in ecological studies. Snake venoms are a rich source of components that have a variety of biological functions. Among them are proteins like lectins, metalloproteinases, serine proteinases, L-amino acid oxidases, phospholipases type A₂, cysteine-rich secretory proteins, as well as many oligopeptides, such as waprins, cardiotoxins, cathelicidins, and β-defensins. In vitro, these biomolecules were shown to be active against bacteria, fungi, parasites, and viruses that are pathogenic to humans. Not only cathelicidins, but all other proteins and oligopeptides from snake venom have been proteolyzed to provide short antimicrobial peptides, or for use as templates for developing a variety of short unnatural sequences based on their structures. In addition to organizing and discussing an expressive amount of information, this review also describes new β-defensin sequences of Sistrurus miliarius that can lead to novel peptide-based antimicrobial agents, using a multidisciplinary approach that includes sequence phylogeny.

Chromosomal Mapping of the Histone Multigene Family and U2 snRNA in Hypancistrus Species (Siluriformes, Loricariidae) from the Brazilian Amazon

Loricariidae (Siluriformes) comprises ∼1026 species of neotropical fish, being considered the most diverse among the Siluriformes. Studies on repetitive DNA sequences have provided important data on the evolution of the genomes of members of this family, especially of the Hypostominae subfamily. In this study, the chromosomal mapping of the histone multigene family and U2 snRNA was performed in two species belonging to the Hypancistrus genus, Hypancistrus sp. "pão" (2n = 52, 22m + 18sm +12st) and Hypancistrus zebra (2n = 52, 16m + 20sm +16st). The presence of dispersed signals of histones H2A, H2B, H3, and H4 in the karyotype of both species, with each sequence displaying a varied level of accumulation and dispersion of these sequences between them was observed; in addition, U2 snDNA probe only showed positive results in H. zebra, which present this multigene in the terminal region of three chromosomal pairs. The obtained results resemble data already analyzed in the literature, in which the action of transposable elements interfere in the organization of these multigene families, in addition to other evolutionary processes that shape the evolution of the genome, such as circular or ectopic recombination. This study also shows that the dispersion of the multigene histone family is quite complex, and from this, these data serve as a point of discussion for the evolutionary processes that occur in the Hypancistrus karyotype.

A highly diverse set of novel immunoglobulin-like transcript (NILT) genes in zebrafish indicates a wide range of functions with complex relationships to mammalian receptors

Multiple novel immunoglobulin-like transcripts (NILTs) have been identified from salmon, trout, and carp. NILTs typically encode activating or inhibitory transmembrane receptors with extracellular immunoglobulin (Ig) domains. Although predicted to provide immune recognition in ray-finned fish, we currently lack a definitive framework of NILT diversity, thereby limiting our predictions for their evolutionary origin and function. In order to better understand the diversity of NILTs and their possible roles in immune function, we identified five NILT loci in the Atlantic salmon (Salmo salar) genome, defined 86 NILT Ig domains within a 3-Mbp region of zebrafish (Danio rerio) chromosome 1, and described 41 NILT Ig domains as part of an alternative haplotype for this same genomic region. We then identified transcripts encoded by 43 different NILT genes which reflect an unprecedented diversity of Ig domain sequences and combinations for a family of non-recombining receptors within a single species. Zebrafish NILTs include a sole putative activating receptor but extensive inhibitory and secreted forms as well as membrane-bound forms with no known signaling motifs. These results reveal a higher level of genetic complexity, interindividual variation, and sequence diversity for NILTs than previously described, suggesting that this gene family likely plays multiple roles in host immunity.

Genetic characteristics and growth patterns of the hybrid grouper derived from the hybridization of Epinephelus fuscoguttatus (female) × Epinephelus polyphekadion (male)

Hybridization is one of the primary methods used to cultivate farmed grouper species. The hybrid grouper derived from crossing Epinephelus fuscoguttatus (♀) and E. polyphekadion (♂) exhibits growth superiority over its parents. The genetic characteristics and growth patterns of the hybrid grouper have not yet been defined. This study confirms the ploidy level of the hybrid grouper (2n = 48) using chromosome count analysis and flow cytometry. The 5S rDNA family was used to evaluate genetic diversity. Only one 5S class (~400 bp) was detected in the hybrid grouper, which could be used to distinguish between two different types based on nucleotide sequences, likely representing homologous unit classes from the female and male parental species. Growth patterns of 5-8-month-old hybrid groupers were also monitored. In this phase, a positive allometric growth pattern in body mass with total length was found. Body height and body mass were significantly correlated based on correlation and path coefficient, suggesting that body height could serve as an excellent index to increase body mass. These results aid our understanding of the genetic evolution of the hybrid grouper and inform the development of improved rearing techniques.

How Many Sirtuin Genes Are Out There? Evolution of Sirtuin Genes in Vertebrates With a Description of a New Family Member

Studying the evolutionary history of gene families is a challenging and exciting task with a wide range of implications. In addition to exploring fundamental questions about the origin and evolution of genes, disentangling their evolution is also critical to those who do functional/structural studies to allow a deeper and more precise interpretation of their results in an evolutionary context. The sirtuin gene family is a group of genes that are involved in a variety of biological functions mostly related to aging. Their duplicative history is an open question, as well as the definition of the repertoire of sirtuin genes among vertebrates. Our results show a well-resolved phylogeny that represents an improvement in our understanding of the duplicative history of the sirtuin gene family. We identified a new sirtuin gene family member (SIRT3.2) that was apparently lost in the last common ancestor of amniotes but retained in all other groups of jawed vertebrates. According to our experimental analyses, elephant shark SIRT3.2 protein is located in mitochondria, the overexpression of which leads to an increase in cellular levels of ATP. Moreover, in vitro analysis demonstrated that it has deacetylase activity being modulated in a similar way to mammalian SIRT3. Our results indicate that there are at least eight sirtuin paralogs among vertebrates and that all of them can be traced back to the last common ancestor of the group that existed between 676 and 615 millions of years ago.

Polymorphisms of rDNA genes in Cyberlindnera yeast suggest birth-and-death evolution events

In eukaryotes, the ribosome machinery is encoded by repeats of the ribosomal RNA genes: 26/28S, 18S, 5.8S, and 5S, structured in tandem arrays and frequently homogenized within a genome. This homogenization is thought to be driven by concerted evolution, evolving as a unit, which contributes to its target as the species barcode in modern taxonomy. However, high heterogeneity of rDNA genes has been reported, including in Saccharomycotina yeasts. Here, we describe the polymorphisms and heterogeneity of D1/D2 domains (26S rRNA) and the intergenic transcribed spacer of a new yeast species with affinities to the genus Cyberlindnera and their evolution. Both regions are not homogenized, failing the prediction of concerted evolution. Phylogenetic network analysis of cloned sequences revealed that Cyberlindnera sp. rDNAs are diverse and evolved by reticulation rather than by bifurcating tree evolution model. Predicted rRNA secondary structures also confirmed structural differences, except for some conserved hairpin loops. We hypothesize that some rDNA is inactive within this species and evolves by birth-and-death rather than concerted evolution. Our findings propel further investigation into the evolution of rDNA genes in yeasts.

Pheromone Response and Mating Behavior in Fission Yeast

Most ascomycete fungi, including the fission yeast Schizosaccharomyces pombe, secrete two peptidyl mating pheromones: C-terminally modified and unmodified peptides. S. pombe has two mating types, plus and minus, which secrete two different pheromones, P-factor (unmodified) and M-factor (modified), respectively. These pheromones are specifically recognized by receptors on the cell surface of cells of opposite mating types, which trigger a pheromone response. Recognition between pheromones and their corresponding receptors is important for mate discrimination; therefore, genetic changes in pheromone or receptor genes affect mate recognition and cause reproductive isolation that limits gene flow between populations. Such genetic variation in recognition via the pheromone/receptor system may drive speciation. Our recent studies reported that two pheromone receptors in S. pombe might have different stringencies in pheromone recognition. In this review, we focus on the molecular mechanism of pheromone response and mating behavior, emphasizing pheromone diversification and its impact on reproductive isolation in S. pombe and closely related fission yeast species. We speculate that the "asymmetric" system might allow flexible adaptation to pheromone mutational changes while maintaining stringent recognition of mating partners. The loss of pheromone activity results in the extinction of an organism's lineage. Therefore, genetic changes in pheromones and their receptors may occur gradually and/or coincidently before speciation. Our findings suggest that the M-factor plays an important role in partner discrimination, whereas P-factor communication allows flexible adaptation to create variations in S. pombe. Our inferences provide new insights into the evolutionary mechanisms underlying pheromone diversification.

Uncovering a 500 million year old history and evidence of pseudogenization for TLR15

Introduction

Toll like receptors (TLRs) are at the front line of pathogen recognition and host immune response. Many TLR genes have been described to date with some being found across metazoans while others are restricted to specific lineages. A cryptic member of the TLR gene family, TLR15, has a unique phylogenetic distribution. Initially described in extant species of birds and reptiles, an ortholog has been reported for cartilaginous fish.

Methods

Here, we significantly expanded the evolutionary analysis of TLR15 gene evolution, taking advantage of large genomic and transcriptomic resources available from different lineages of vertebrates. Additionally, we objectively search for TLR15 in lobe-finned and ray-finned fish, as well as in cartilaginous fish and jawless vertebrates.

Results and discussion

We confirm the presence of TLR15 in early branching jawed vertebrates - the cartilaginous fish, as well as in basal Sarcopterygii - in lungfish. However, within cartilaginous fish, the gene is present in Holocephalans (all three families) but not in Elasmobranchs (its sister-lineage). Holocephalans have long TLR15 protein sequences that disrupt the typical TLR structure, and some species display a pseudogene sequence due to the presence of frameshift mutations and early stop codons. Additionally, TLR15 has low expression levels in holocephalans when compared with other TLR genes. In turn, lungfish also have long TLR15 protein sequences but the protein structure is not compromised. Finally, TLR15 presents several sites under negative selection. Overall, these results suggest that TLR15 is an ancient TLR gene and is experiencing ongoing pseudogenization in early-branching vertebrates.

Hybridization and diversity of aquatic macrophyte Sparganium L. (Typhaceae) as revealed by high-throughput nrDNA sequencing

Sparganium is an emergent aquatic macrophyte widely spread in temperate and subtropical zones. Taxa of this genus feature high phenotypic plasticity and can produce interspecific hybrids. By means of high-throughput sequencing of the internal transcribed spacer (ITS1) of 35S rDNA, the status of 15 Eurasian Sparganium species and subspecies was clarified and the role of hybridization events in the recent evolution of the genus was investigated. It has been shown that a number of species such as S. angustifolium, S. fallax and S. subglobosum have homogenized rDNA represented by one major ribotype. The rDNA of other taxa is represented by two or more major ribotypes. Species with high rDNA heterogeneity are apparently of hybrid origin. Based on the differences in rDNA patterns, intraspecific diversity was identified in S. probatovae and S. emersum. Thus, we have concluded that Sparganium has extensive interspecific hybridization at the subgenus level, and there may also be occasional hybridization between species from different subgenera.

Rapid molecular diversification and homogenization of clustered major ampullate silk genes in Argiope garden spiders

The evolutionary diversification of orb-web weaving spiders is closely tied to the mechanical performance of dragline silk. This proteinaceous fiber provides the primary structural framework of orb web architecture, and its extraordinary toughness allows these structures to absorb the high energy of aerial prey impact. The dominant model of dragline silk molecular structure involves the combined function of two highly repetitive, spider-specific, silk genes (spidroins)-MaSp1 and MaSp2. Recent genomic studies, however, have suggested this framework is overly simplistic, and our understanding of how MaSp genes evolve is limited. Here we present a comprehensive analysis of MaSp structural and evolutionary diversity across species of Argiope (garden spiders). This genomic analysis reveals the largest catalog of MaSp genes found in any spider, driven largely by an expansion of MaSp2 genes. The rapid diversification of Argiope MaSp genes, located primarily in a single genomic cluster, is associated with profound changes in silk gene structure. MaSp2 genes, in particular, have evolved complex hierarchically organized repeat units (ensemble repeats) delineated by novel introns that exhibit remarkable evolutionary dynamics. These repetitive introns have arisen independently within the genus, are highly homogenized within a gene, but diverge rapidly between genes. In some cases, these iterated introns are organized in an alternating structure in which every other intron is nearly identical in sequence. We hypothesize that this intron structure has evolved to facilitate homogenization of the coding sequence. We also find evidence of intergenic gene conversion and identify a more diverse array of stereotypical amino acid repeats than previously recognized. Overall, the extreme diversification found among MaSp genes requires changes in the structure-function model of dragline silk performance that focuses on the differential use and interaction among various MaSp paralogs as well as the impact of ensemble repeat structure and different amino acid motifs on mechanical behavior.

Novel selectively amplified DNA sequences in the germline genome of the Japanese hagfish, Eptatretus burgeri

In the Japanese hagfish Eptatretus burgeri, 16 chromosomes (eliminated [E]-chromosomes) have been lost in somatic cells (2n = 36), which is equivalent to approx. 21% of the genomic DNA in germ cells (2n = 52). At least seven of the 12 eliminated repetitive DNA families isolated in eight hagfish species were selectively amplified in the germline genome of this species. One of them, EEEb1 (eliminated element of E. burgeri 1) is exclusively localized on all E-chromosomes. Herein, we identified four novel eliminated repetitive DNA families (named EEEb3-6) through PCR amplification and suppressive subtractive hybridization (SSH) combined with Southern-blot hybridization. EEEb3 was mosaic for 5S rDNA and SINE elements. EEEb4 was GC-rich repeats and has one pair of direct and inverted repeats, whereas EEEb5 and EEEb6 were AT-rich repeats with one pair and two pairs of sub-repeats, respectively. Interestingly, all repeat classes except EEEb3 were transcribed in the testes, although no open reading frames (ORF) were identified. We conducted fluorescence in situ hybridization (FISH) to examine the chromosomal localizations of EEEb3-6 and EEEb2, which was previously isolated from the germline genome of E. burgeri. All sequences were only found on all EEEb1-positive E-chromosomes. Copy number estimation of the repeated elements by slot-blot hybridization revealed that (i) the EEEb1-6 family members occupied 39.9% of the total eliminated DNA, and (ii) a small number of repeats were retained in somatic cells, suggesting that there is incomplete elimination of the repeated elements. These results provide new insights into the mechanisms involved in the chromosome elimination and the evolution of E-chromosomes.

Meta-analysis of major histocompatibility complex (MHC) class IIA reveals polymorphism and positive selection in many vertebrate species

Pathogen-mediated selection and sexual selection are important drivers of evolution. Both processes are known to target genes of the major histocompatibility complex (MHC), a gene family encoding cell-surface proteins that display pathogen peptides to the immune system. The MHC is also a model for understanding processes such as gene duplication and trans-species allele sharing. The class II MHC protein is a heterodimer whose peptide-binding groove is encoded by an MHC-IIA gene and an MHC-IIB gene. However, our literature review found that class II MHC papers on infectious disease or sexual selection included IIA data only 18% and 9% of the time, respectively. To assess whether greater emphasis on MHC-IIA is warranted, we analysed MHC-IIA sequence data from 50 species of vertebrates (fish, amphibians, birds, mammals) to test for polymorphism and positive selection. We found that the number of MHC-IIA alleles within a species was often high, and covaried with sample size and number of MHC-IIA genes assayed. While MHC-IIA variability tended to be lower than that of MHC-IIB, the difference was only ~25%, with ~3 fewer IIA alleles than IIB. Furthermore, the unexpectedly high MHC-IIA variability showed clear signatures of positive selection in most species, and positive selection on MHC-IIA was stronger in fish than in other surveyed vertebrate groups. Our findings underscore that MHC-IIA can be an important target of selection. Future studies should therefore expand the characterization of MHC-IIA at both allelic and genomic scales, and incorporate MHC-IIA into models of fitness consequences of MHC variation.

Varying strength of selection contributes to the intragenomic diversity of rRNA genes

Ribosome biogenesis in eukaryotes is supported by hundreds of ribosomal RNA (rRNA) gene copies that are encoded in the ribosomal DNA (rDNA). The multiple copies of rRNA genes are thought to have low sequence diversity within one species. Here, we present species-wide rDNA sequence analysis in Saccharomyces cerevisiae that challenges this view. We show that rDNA copies in this yeast are heterogeneous, both among and within isolates, and that many variants avoided fixation or elimination over evolutionary time. The sequence diversity landscape across the rDNA shows clear functional stratification, suggesting different copy-number thresholds for selection that contribute to rDNA diversity. Notably, nucleotide variants in the most conserved rDNA regions are sufficiently deleterious to exhibit signatures of purifying selection even when present in only a small fraction of rRNA gene copies. Our results portray a complex evolutionary landscape that shapes rDNA sequence diversity within a single species and reveal unexpectedly strong purifying selection of multi-copy genes.

Comparative analysis of detoxification-related gene superfamilies across five hemipteran species

Background

Hemiptera is one of the most speciose orders of insects, and the most speciose considering Hemimetabola. Through their evolutive history, hemipterans with different feeding habits have adapted to deal with different chemical challenges. Three major gene families are involved in xenobiotic detoxification in insects: the cytochromes P450 (CYPs), carboxyl/cholinesterases (CCEs), and glutathione transferases (GSTs). Here we perform a comparative analysis on the complement of these gene superfamilies across five hemipteran species; four heteropterans (the pentatomid plant feeders Nezara viridula and Halyomorpha halys; the hematophagous Cimex lectularius, Cimicidae, and Rhodnius prolixus, Reduviidae), and one Auchenorrhyncha plant feeder (Nilaparvata lugens).

Results

Our results point to an expansion of several enzyme families associated with xenobiotic detoxification in heteropterans with respect to other species and the existence of a dynamic evolution pattern including CYP3 clan, hormone and pheromone processing class in the CCE superfamily, and sigma class in GST superfamily. Other detoxification-related families are reduced in the hemipteran species analyzed here: reduction or even absence of epsilon class and reduced delta class in GST superfamily; absence of mitochondrial CYP12 family; absence of CYP9 family in CYP3 clan; and reduction or even absence of some dietary/detoxification groups of CCEs. Interestingly, the most polyphagous species analyzed here (H. halys) is also the one that presents the largest repertoire of detoxification enzymes. Gene cluster analysis suggests that this could be due to gene duplication events.

Conclusions

The evolutionary analysis performed here reveals characteristics that are both common and particular for heteropterans. The composition and organization of detoxification-related gene families could shed light on evolutionary forces that shaped their divergence. These families are important for both the detoxification of diet products and for conferring tolerance or resistance to synthetic insecticides. Furthermore, we present the first comprehensive analysis of detoxification gene superfamilies in N. viridula, an understudied species in spite of its economic relevance as a crop pest. The information obtained is of interest for basic insect science as well as for the control of harmful species and the management of insecticide resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08974-y.

Placing human gene families into their evolutionary context

Following the draft sequence of the first human genome over 20 years ago, we have achieved unprecedented insights into the rules governing its evolution, often with direct translational relevance to specific diseases. However, staggering sequence complexity has also challenged the development of a more comprehensive understanding of human genome biology. In this context, interspecific genomic studies between humans and other animals have played a critical role in our efforts to decode human gene families. In this review, we focus on how the rapid surge of genome sequencing of both model and non-model organisms now provides a broader comparative framework poised to empower novel discoveries. We begin with a general overview of how comparative approaches are essential for understanding gene family evolution in the human genome, followed by a discussion of analyses of gene expression. We show how homology can provide insights into the genes and gene families associated with immune response, cancer biology, vision, chemosensation, and metabolism, by revealing similarity in processes among distant species. We then explain methodological tools that provide critical advances and show the limitations of common approaches. We conclude with a discussion of how these investigations position us to gain fundamental insights into the evolution of gene families among living organisms in general. We hope that our review catalyzes additional excitement and research on the emerging field of comparative genomics, while aiding the placement of the human genome into its existentially evolutionary context.

Ecological constraints on highly evolvable olfactory receptor genes and morphology in neotropical bats

Although evolvability of genes and traits may promote specialization during species diversification, how ecology subsequently restricts such variation remains unclear. Chemosensation requires animals to decipher a complex chemical background to locate fitness-related resources, and thus the underlying genomic architecture and morphology must cope with constant exposure to a changing odorant landscape; detecting adaptation amidst extensive chemosensory diversity is an open challenge. In phyllostomid bats, an ecologically diverse clade that evolved plant visiting from a presumed insectivorous ancestor, the evolution of novel food detection mechanisms is suggested to be a key innovation, as plant-visiting species rely strongly on olfaction, supplementarily using echolocation. If this is true, exceptional variation in underlying olfactory genes and phenotypes may have preceded dietary diversification. We compared olfactory receptor (OR) genes sequenced from olfactory epithelium transcriptomes and olfactory epithelium surface area of bats with differing diets. Surprisingly, although OR evolution rates were quite variable and generally high, they are largely independent of diet. Olfactory epithelial surface area, however, is relatively larger in plant-visiting bats and there is an inverse relationship between OR evolution rates and surface area. Relatively larger surface areas suggest greater reliance on olfactory detection and stronger constraint on maintaining an already diverse OR repertoire. Instead of the typical case in which specialization and elaboration are coupled with rapid diversification of associated genes, here the relevant genes are already evolving so quickly that increased reliance on smell has led to stabilizing selection, presumably to maintain the ability to consistently discriminate among specific odorants-a potential ecological constraint on sensory evolution.

Metabolization and sequestration of plant specialized metabolites in insect herbivores: Current and emerging approaches

Herbivorous insects encounter diverse plant specialized metabolites (PSMs) in their diet, that have deterrent, anti-nutritional, or toxic properties. Understanding how they cope with PSMs is crucial to understand their biology, population dynamics, and evolution. This review summarizes current and emerging cutting-edge methods that can be used to characterize the metabolic fate of PSMs, from ingestion to excretion or sequestration. It further emphasizes a workflow that enables not only to study PSM metabolism at different scales, but also to tackle and validate the genetic and biochemical mechanisms involved in PSM resistance by herbivores. This review thus aims at facilitating research on PSM-mediated plant-herbivore interactions.

Functional Divergence in a Multi-gene Family Is a Key Evolutionary Innovation for Anaerobic Growth in Saccharomyces cerevisiae

The amplification and diversification of genes into large multi-gene families often mark key evolutionary innovations, but this process often creates genetic redundancy that hinders functional investigations. When the model budding yeast Saccharomyces cerevisiae transitions to anaerobic growth conditions, the cell massively induces the expression of seven serine/threonine-rich anaerobically-induced cell wall mannoproteins (anCWMPs): TIP1, TIR1, TIR2, TIR3, TIR4, DAN1, and DAN4. Here, we show that these genes likely derive evolutionarily from a single ancestral anCWMP locus, which was duplicated and translocated to new genomic contexts several times both prior to and following the budding yeast whole genome duplication (WGD) event. Based on synteny and their phylogeny, we separate the anCWMPs into four gene subfamilies. To resolve prior inconclusive genetic investigations of these genes, we constructed a set of combinatorial deletion mutants to determine their contributions toward anaerobic growth in S. cerevisiae. We found that two genes, TIR1 and TIR3, were together necessary and sufficient for the anCWMP contribution to anaerobic growth. Overexpressing either gene alone was insufficient for anaerobic growth, implying that they encode non-overlapping functional roles in the cell during anaerobic growth. We infer from the phylogeny of the anCWMP genes that these two important genes derive from an ancient duplication that predates the WGD event, whereas the TIR1 subfamily experienced gene family amplification after the WGD event. Taken together, the genetic and molecular evidence suggests that one key anCWMP gene duplication event, several auxiliary gene duplication events, and functional divergence underpin the evolution of anaerobic growth in budding yeasts.

BGFD: an integrated multi-omics database of barley gene families

BACKGROUND

A gene family comprises a group of genes with similar functional domains that play various roles in plant growth, development, and responses to environmental stimuli. Barley (Hordeum vulgare L.) is the fourth most cultivated cereal crop worldwide, and it is an important model species for genetic studies. Systematic identification and annotation of gene families are key for studies of molecular function and evolutionary history.

RESULTS

We constructed a multi-omics database containing 5593 genes of 77 gene families called the Barley Gene Family Database (BGFD: http://barleygfdb.com ). BGFD is a free, user-friendly, and web-accessible platform that provides data on barley family genes. BGFD provides intuitive visual displays to facilitate studies of the physicochemical properties, gene structure, phylogenetic relationships, and motif organization of genes. Massive multi-omics datasets have been acquired and processed to generate an atlas of expression pattern profiles and genetic variation in BGFD. The platform offers several practical toolkits to conduct searches, browse, and employ BLAST functions, and the data are downloadable.

CONCLUSIONS

BGFD will aid research on the domestication and adaptive evolution of barley; it will also facilitate the screening of candidate genes and exploration of important agronomic traits in barley.

CC and CXC chemokines in turbot (Scophthalmus maximus L.): Identification, evolutionary analyses, and expression profiling after Aeromonas salmonicida infection

Chemokines are a superfamily of structurally related cytokines, which exert essential roles in guiding cell migration in development, homeostasis, and immunity. CC and CXC chemokines are the two major subfamilies in teleost species. In this study, a total of seventeen CC and CXC chemokines, with inclusion of twelve CC and five CXC chemokines, were systematically identified from the turbot genome, making turbot the teleost harboring the least number of CC and CXC chemokines among all teleost species ever reported. Phylogeny, synteny, and genomic organization analyses were performed to annotate these genes, and multiple chemokine genes were identified in the turbot genome, due to the tandem duplications (CCL19 and CCL20), the whole genome duplications (CCL20, CCL25, and CXCL12), and the teleost-specific members (CCL34-36, CCL44, and CXCL18). In addition, chemokines were ubiquitously expressed in nine examined healthy tissues, with high expression levels observed in liver, gill, and spleen. Moreover, most chemokines were significantly differentially expressed in gill and spleen after Aeromonas salmonicida infection, and exhibited tissue-specific and time-dependent manner. Finally, protein-protein interaction network (PPI) analysis indicated that turbot chemokines interacted with a few immune-related genes such as interleukins, cathepsins, stats, and TLRs. These results should be valuable for comparative immunological studies and provide insights for further functional characterization of chemokines in teleost.

Spodoptera littoralis genome mining brings insights on the dynamic of expansion of gustatory receptors in polyphagous noctuidae

The bitter taste, triggered via gustatory receptors, serves as an important natural defense against the ingestion of poisonous foods in animals, and the increased host breadth is usually linked to an increase in the number of gustatory receptor genes. This has been especially observed in polyphagous insect species, such as noctuid species from the Spodoptera genus. However, the dynamic and physical mechanisms leading to these gene expansions and the evolutionary pressures behind them remain elusive. Among major drivers of genome dynamics are the transposable elements but, surprisingly, their potential role in insect gustatory receptor expansion has not been considered yet. In this work, we hypothesized that transposable elements and possibly positive selection would be involved in the highly dynamic evolution of gustatory receptor in Spodoptera spp. We first sequenced de novo the full 465 Mb genome of S. littoralis, and manually annotated the main chemosensory genes, including a large repertoire of 373 gustatory receptor genes (including 19 pseudogenes). We also improved the completeness of S. frugiperda and S. litura gustatory receptor gene repertoires. Then, we annotated transposable elements and revealed that a particular category of class I retrotransposons, the SINE transposons, was significantly enriched in the vicinity of gustatory receptor gene clusters, suggesting a transposon-mediated mechanism for the formation of these clusters. Selection pressure analyses indicated that positive selection within the gustatory receptor gene family is cryptic, only 7 receptors being identified as positively selected. Altogether, our data provide a new good quality Spodoptera genome, pinpoint interesting gustatory receptor candidates for further functional studies and bring valuable genomic information on the mechanisms of gustatory receptor expansions in polyphagous insect species.

Comparative Genomic Analysis of Seven Vibrio alginolyticus Strains Isolated From Shrimp Larviculture Water With Emphasis on Chitin Utilization

The opportunistic pathogen Vibrio alginolyticus is gaining attention because of its disease-causing risks to aquatic animals and humans. In this study, seven Vibrio strains isolated from different shrimp hatcheries in Southeast China were subjected to genome sequencing and subsequent comparative analysis to explore their intricate relationships with shrimp aquaculture. The seven isolates had an average nucleotide identity of ≥ 98.3% with other known V. alginolyticus strains. The species V. alginolyticus had an open pan-genome, with the addition of ≥ 161 novel genes following each new genome for seven isolates and 14 publicly available V. alginolyticus strains. The percentages of core genes of the seven strains were up to 83.1–87.5%, indicating highly conserved functions, such as chitin utilization. Further, a total of 14 core genes involved in the chitin degradation pathway were detected on the seven genomes with a single copy, 12 of which had undergone significant purifying selection (dN/dS < 1). Moreover, the seven strains could utilize chitin as the sole carbon-nitrogen source. In contrast, mobile genetic elements (MGEs) were identified in seven strains, including plasmids, prophages, and genomic islands, which mainly encoded accessory genes annotated as hypothetical proteins. The infection experiment showed that four of the seven strains might be pathogenic because the survival rates of Litopenaeus vannamei postlarvae were significantly reduced (P < 0.05) when compared to the control. However, no obvious correlation was noted between the number of putative virulence factors and toxic effects of the seven strains. Collectively, the persistence of V. alginolyticus in various aquatic environments may be attributed to its high genomic plasticity via the acquisition of novel genes by various MGEs. In view of the strong capability of chitin utilization by diverse vibrios, the timely removal of massive chitin-rich materials thoroughly in shrimp culture systems may be a key strategy to inhibit proliferation of vibrios and subsequent infection of shrimp. In addition, transcontinental transfer of potentially pathogenic V. alginolyticus strains should receive great attention to avoid vibriosis.

Copy number changes in co-expressed odorant receptor genes enable selection for sensory differences in drosophilid species

Despite numerous examples of chemoreceptor gene family expansions and contractions, how these relate to modifications in the sensory neuron populations in which they are expressed remains unclear. Drosophila melanogaster's odorant receptor (Or) family is ideal for addressing this question because most Ors are expressed in distinct olfactory sensory neuron (OSN) types. Between-species changes in Or copy number may therefore indicate increases or reductions in the number of OSN populations. Here we investigated the Or67a subfamily, which exhibits copy number variation in D. melanogaster and its closest relatives: D. simulans, D. sechellia and D. mauritiana. These species' common ancestor had three Or67a paralogues that had already diverged adaptively. Following speciation, two Or67a paralogues were lost independently in D. melanogaster and D. sechellia, with ongoing positive selection shaping the intact genes. Unexpectedly, the functionally diverged Or67a paralogues in D. simulans are co-expressed in a single neuron population, which projects to a glomerulus homologous to that innervated by Or67a neurons in D. melanogaster. Thus, while sensory pathway neuroanatomy is conserved, independent selection on co-expressed receptors has contributed to species-specific peripheral coding. This work reveals a type of adaptive change largely overlooked for olfactory evolution, raising the possibility that similar processes influence other cases of insect Or co-expression.

Karyotypes of Manatees: New Insights into Hybrid Formation (Trichechus inunguis × Trichechus m. manatus) in the Amazon Estuary

Great efforts have been made to preserve manatees. Recently, a hybrid zone was described between Trichechus inunguis (TIN) and the Trichechus manatus manatus (TMM) in the Amazon estuary. Cytogenetic data on these sirenians are limited, despite being fundamental to understanding the hybridization/introgression dynamics and genomic organization in Trichechus. We analyzed the karyotype of TMM, TIN, and two hybrid specimens ("Poque" and "Vitor") by classical and molecular cytogenetics. G-band analysis revealed that TMM (2n = 48) and TIN (2n = 56) diverge by at least six Robertsonian translocations and a pericentric inversion. Hybrids had 2n = 50, however, with Autosomal Fundamental Number (FNA) = 88 in "Poque" and FNA = 74 in "Vitor", and chromosomal distinct pairs in heterozygous; additionally, "Vitor" exhibited heteromorphisms and chromosomes whose pairs could not be determined. The U2 snDNA and Histone H3 multi genes are distributed in small clusters along TIN and TMM chromosomes and have transposable Keno and Helitron elements (TEs) in their sequences. The different karyotypes observed among manatee hybrids may indicate that they represent different generations formed by crossing between fertile hybrids and TIN. On the other hand, it is also possible that all hybrids recorded represent F1 and the observed karyotype differences must result from mechanisms of elimination.

Comparisons of chemosensory gene repertoires in human and non-human feeding Anopheles mosquitoes link olfactory genes to anthropophily

We investigate the genetic basis of anthropophily (human host use) in a non-model mosquito species group, the Anopheles farauti complex from the southwest Pacific. This complex has experienced multiple transitions from anthropophily to zoophily, contrasting with well-studied systems (the global species Aedes aegypti and the African Anopheles gambiae complex) that have evolved to be specialist anthropophiles. By performing tests of selection and assessing evolutionary patterns for >200 olfactory genes from nine genomes, we identify several candidate genes associated with differences in anthropophily in this complex. Based on evolutionary patterns (phylogenetic relationships, fixed amino acid differences, and structural differences) as well as results from selection analyses, we identify numerous genes that are likely to play an important role in mosquitoes’ ability to detect humans as hosts. Our findings contribute to the understanding of the evolution of insect olfactory gene families and mosquito host preference as well as having potential applied outcomes.

•

Genomes of Anopheles mosquitoes with differing host preferences were sequenced

•

Evolutionary comparisons were performed on >200 insect chemosensory genes

•

These comparisons revealed candidate genes involved in human feeding

•

Two of the main candidates identified were co-receptor Ir8a and Or75

Single amino acid variations drive functional divergence of cytochrome P450s in Helicoverpa species

Divergence of gene function is a hallmark of evolution, but assessing such divergence in one species or between species requires information on functional alterations of the alleles and homologs. Here, we explore the functional divergence of two paralogs, CYP6AE19 and CYP6AE20, from Helicoverpa armigera, and two close orthologs, CYP6B8 and CYP6B7, from two related species (Helicoverpa zea and H. armigera); although there is high sequence identity within each pair of enzymes, the latter P450 of each pair has lost metabolic competence towards the plant allelochemical xanthotoxin. Multiple chimeric and single/double site mutants were created by exchanging the diverse substrate recognition sites (SRSs) and amino acids within each pair of P450s. Heterologous expression in Sf9 cells and in vitro metabolism studies showed that the exchange of SRS4 swapped the activity of CYP6AE19 and CYP6AE20, and subsequent site-directed mutagenesis demonstrated that the CYP6AE20 V318M substitution causes a gain-of-function towards xanthotoxin. Meanwhile, a single amino acid substitution (L489P) in SRS6 was found to swap activity between the CYP6B orthologs. Sequence alignments of CYP6AE paralogs and all reported insect xanthotoxin-metabolizing P450s suggest M318 and P489 are essential for the catalytic activities of CYP6AE paralogs and CYP6B orthologs, respectively, but P450s in different subfamilies may have different mechanisms towards the same substrate. Our findings demonstrate that a single amino acid substitution can suffice to alter substrate metabolism and this functional divergence resulting from natural mutations will help to further our understanding of the process of natural selection of P450 genes and their role in insect-host plant interactions.

Endangered Schizothoracin Fish in the Tarim River Basin Are Threatened by Introgressive Hybridization

Simple Summary

Interspecific hybridization and introgression may threaten the survival of endangered species. Big-head Schizothoracin (Aspiorhynchus laticeps) and Tarim Schizothoracin (Schizothorax biddulphi) are two recognized types of fish of Schizothoracinae which belong to the family Cyprinidae. Big-head Schizothoracin and Tarim Schizothoracin are sympatrically distributed in the Tarim River basin in Xinjiang, China, and have been listed as first-class and second-class national key protected animals, respectively. Based on morphological characteristics, Schizothorax esocinus (another Schizothoracin that occurs in the Tarim River basin) is speculated to be hybrid offspring of Big-head Schizothoracin and Tarim Schizothoracin, but there is no direct genetic evidence for this point. In this study, the hybridization status of Schizothorax esocinus was confirmed by comparing genetic constitutions of mitochondrial DNA (mtDNA) and inter transcribed spacer (ITS2) encoded by the nuclear genome. Extensive hybridization and introgression between Big-head Schizothoracin and Tarim Schizothoracin was detected. Since both Big-head Schizothoracin and Tarim Schizothoracin are critically endangered fishes, risk management due to hybridization is recommended to be considered as a part of current conservation programs.

Abstract

Big-head Schizothoracin (Aspiorhynchus laticeps) and Tarim Schizothoracin (Schizothorax biddulphi) are locally sympatric in the Tarim River Basin. Although another Schizothoracin (Schizothorax esocinus) in Tarim River basin has been speculated to be hybrid offspring of Big-head Schizothoracin and Tarim Schizothoracin, there was no genetic evidence. Previous studies on the genetics and evolution of Schizothoracins in Xinjiang Province were mostly based on mitochondrial DNA (mtDNA), whose characteristics of maternal inheritance made it hard to answer the question of whether there was hybridization and introgression between Big-head Schizothoracin and Tarim Schizothoracin. In this study, cytochrome b (cytb) gene of mtDNA and internal transcribed spacer 2 (ITS2) that is encoded by the nuclear genome were genotyped within the entire samples at the same time. Our results confirmed that Schizothorax esocinus was the hybrid offspring of Big-head Schizothoracin and Tarim Schizothoracin. The heterozygous ITS2 genotypes and/or Aspiorhynchus laticeps-like mtDNA were also detected in a subset of samples that should have been identified as pure Schizothorax biddulphi based on morphology. The ITS2 is characterized by multi-copy, concert evolution, and biparental inheritance. Thus, by comparing with mtDNA data, broad-scale bidirectional hybridization and introgression between Big-head Schizothoracin and Tarim Schizothoracin were revealed. Although interspecific hybridization may play a positive role in ecology and evolution, interspecific hybrids could threaten their parental species by the swamping of genetics and demography. As both parents of hybridization are critically endangered fishes, in this case, it is urgently necessary to strengthen the scientific assessment of the risks of the hybrids and the control of the hybridization and introgression between Aspiorhynchus laticeps and Schizothorax biddulphi in the Tarim River Basin.

Paradox lost: Concerted evolution and centromeric instability: Centromeres are hospitable habitats for repeats that evolve adaptations for proliferation within the nucleus sometimes at organismal cost.: Centromeres are hospitable habitats for repeats that evolve adaptations for proliferation within the nucleus sometimes at organismal cost

Homologous centromeres compete for segregation to the secondary oocyte nucleus at female meiosis I. Centromeric repeats also compete with each other to populate centromeres in mitotic cells of the germline and have become adapted to use the recombinational machinery present at centromeres to promote their own propagation. Repeats are not needed at centromeres, rather centromeres appear to be hospitable habitats for the colonization and proliferation of repeats. This is probably an indirect consequence of two distinctive features of centromeric DNA. Centromeres are subject to breakage by the mechanical forces exerted by microtubules and meiotic crossing-over is suppressed. Centromeric proteins acting in trans are under selection to mitigate the costs of centromeric repeats acting in cis. Collateral costs of mitotic competition at centromeres may help to explain the high rates of aneuploidy observed in early human embryos.

Major and minor U small nuclear RNAs genes characterization in a neotropical fish genome: Chromosomal remodeling and repeat units dispersion in Parodontidae

In association with many proteins, small nuclear RNAs (snRNAs) organize the spliceosomes that play a significant role in processing precursor mRNAs during gene expression. According to snRNAs genic arrangements, two kinds of spliceosomes (major and minor) can be organized into eukaryotic cells. Although in situ localization of U1 and U2 snDNAs have been performed in fish karyotypes, studies with genomic characterization and functionality of U snRNAs integrated into chromosomal changes on Teleostei are still scarce. This study aimed to achieve a genomic characterization of the U snRNAs genes in Apareiodon sp. (2n = 54, ZZ/ZW), apply these data to recognize functional/defective copies, and map chromosomal changes involving snDNAs in Parodontidae species karyotype diversification. Nine snRNA multigene families (U1, U2, U4, U5, U6, U11, U12, U4atac and U6atac) arranged in putatively functional copies in the genome were analyzed. Proximal Sequence Elements (PSE) and TATA-box promoters occurrence, besides an entire transcribed region and conserved secondary structures, qualify them for spliceosome activity. In addition, several defective copies or pseudogenes were identified for the snRNAs that make up the major spliceosome. In situ localization of snDNAs in five species of Parodontidae demonstrated that U1, U2, and U4 snDNAs were involved in chromosomal location changes or units dispersion. The U snRNAs defective/pseudogenes units dispersion could be favored by the probable occurrence of active retrotransposition enzymes in the Apareiodon genome. The U2 and U4 snDNAs sites were involved in independent events in the differentiation of sex chromosomes among Parodontidae lineages. The study characterized U snRNA genes that compose major and minor spliceosomes in the Apareiodon sp. genome and proposes that their defective copies trigger chromosome differentiation and diversification events in Parodontidae.

MHCtools - an R package for MHC high-throughput sequencing data: genotyping, haplotype and supertype inference, and downstream genetic analyses in non-model organisms

The major histocompatibility complex (MHC) plays a central role in the vertebrate adaptive immune system and has been of long-term interest in evolutionary biology. While several protocols have been developed for MHC genotyping, there is a lack of transparent and standardized tools for downstream analysis of MHC data. Here, we present the R package MHCtools and demonstrate the use of its functions to (i) assist accurate MHC genotyping from high-throughput amplicon-sequencing data, (ii) infer functional MHC supertypes using bootstrapped clustering analysis, (iii) identify segregating MHC haplotypes from family data, and (iv) analyse functional and genetic distances between MHC sequences. We employed MHCtools to analyse MHC class I (MHC-I) amplicon data of 559 great reed warblers (Acrocephalus arundinaceus). We identified 390 MHC-I alleles which clustered into 14 functional supertypes. A phylogenetic analysis and analyses of positive selection suggested that the MHC-I alleles belonged to several distinct functional groups. We furthermore identified 107 segregating haplotypes among 116 families, and found substantial variation in diversity with 4-21 MHC-I alleles and 3-13 MHC-I supertypes per haplotype. Finally, we show that the great reed warbler haplotypes harboured combinations of MHC-I supertypes with greater functional divergence than observed in simulated populations of possible haplotypes, a result that is in accordance with the divergent allele advantage hypothesis. Our study demonstrates the power of MHCtools to support genotyping and analysis of MHC in non-model species, which we hope will encourage broad implementation among researchers in MHC genetics and evolution.

Histone renegades: Unusual H2A histone variants in plants and animals

H2A variants are histones that carry out specialized nucleosome function during the eukaryote genome packaging. Most genes encoding H2A histone variants arose in the distant past, and have highly conserved domains and structures. Yet, novel H2A variants have continued to arise throughout the radiation of eukaryotes and disturbed the apparent tranquility of nucleosomes. These species-specific H2A variants contributed to the functional diversification of nucleosomes through changes in both their structure and expression patterns. In this short review, we discuss the evolutionary trajectories of these histone renegades in plants and animal genomes.

TLR7 and TLR8 evolution in lagomorphs: different patterns in the different lineages

Toll-like receptors (TLRs) are one of the most ancient and widely studied innate immune receptors responsible for host defense against invading pathogens. Among the known TLRs, TLR7 and TLR8 sense and recognize single-stranded (ss) RNAs with a dynamic evolutionary history. While TLR8 was lost in birds and duplicated in turtles and crocodiles, TLR7 is duplicated in some birds, but in other tetrapods, there is only one copy. In mammals, with the exception of lagomorphs, TLR7 and TLR8 are highly conserved. Here, we aim to study the evolution of TLR7 and TLR8 in mammals, with a special focus in the order Lagomorpha. By searching public sequence databases, conducting evolutionary analysis, and evaluating gene expression, we were able to confirm that TLR8 is absent in hares but widely expressed in the European rabbit. In contrast, TLR7 is absent in the European rabbit and quite divergent in hares. Our results suggest that, in lagomorphs, more in particular in leporids, TLR7 and TLR8 genes have evolved faster than in any other mammalian group. The long history of interaction with viruses and their location in highly dynamic telomeric regions might explain the pattern observed.

Understanding Diversity, Evolution, and Structure of Small Heat Shock Proteins in Annelida Through in Silico Analyses

Small heat shock proteins (sHsps) are oligomeric stress proteins characterized by an α-crystallin domain (ACD). These proteins are localized in different subcellular compartments and play critical roles in the stress physiology of tissues, organs, and whole multicellular eukaryotes. They are ubiquitous proteins found in all living organisms, from bacteria to mammals, but they have never been studied in annelids. Here, a data set of 23 species spanning the annelid tree of life, including mostly transcriptomes but also two genomes, was interrogated and 228 novel putative sHsps were identified and manually curated. The analysis revealed very high protein diversity and showed that a significant number of sHsps have a particular dimeric architecture consisting of two tandemly repeated ACDs. The phylogenetic analysis distinguished three main clusters, two of them containing both monomeric sHsps, and ACDs located downstream in the dimeric sHsps, and the other one comprising the upstream ACDs from those dimeric forms. Our results support an evolutionary history of these proteins based on duplication events prior to the Spiralia split. Monomeric sHsps 76) were further divided into five subclusters. Physicochemical properties, subcellular location predictions, and sequence conservation analyses provided insights into the differentiating elements of these putative functional groups. Strikingly, three of those subclusters included sHsps with features typical of metazoans, while the other two presented characteristics resembling non-metazoan proteins. This study provides a solid background for further research on the diversity, evolution, and function in the family of the sHsps. The characterized annelid sHsps are disclosed as essential for improving our understanding of this important family of proteins and their pleotropic functions. The features and the great diversity of annelid sHsps position them as potential powerful molecular biomarkers of environmental stress for acting as prognostic tool in a diverse range of environments.

Evolution of Guanylate Binding Protein (GBP) Genes in Muroid Rodents (Muridae and Cricetidae) Reveals an Outstanding Pattern of Gain and Loss

Guanylate binding proteins (GBPs) are paramount in the host immunity by providing defense against invading pathogens. Multigene families related to the immune system usually show that the duplicated genes can either undergo deletion, gain new functions, or become non-functional. Here, we show that in muroids, the Gbp genes followed an unusual pattern of gain and loss of genes. Muroids present a high diversity and plasticity regarding Gbp synteny, with most species presenting two Gbp gene clusters. The phylogenetic analyses revealed seven different Gbps groups. Three of them clustered with GBP2, GBP5 and GBP6 of primates. Four new Gbp genes that appear to be exclusive to muroids were identified as Gbpa, b, c and d. A duplication event occurred in the Gbpa group in the common ancestor of Muridae and Cricetidae (~20 Mya), but both copies were deleted from the genome of Mus musculus, M. caroli and Cricetulus griseus. The Gbpb gene emerged in the ancestor of Muridae and Cricetidae and evolved independently originating Gbpb1 in Muridae, Gbpb2 and Gbpb3 in Cricetidae. Since Gbpc appears only in three species, we hypothesize that it was present in the common ancestor and deleted from most muroid genomes. The second Gbp gene cluster, Gbp6, is widespread across all muroids, indicating that this cluster emerged before the Muridae and Cricetidae radiation. An expansion of Gbp6 occurred in M. musculus and M. caroli probably to compensate the loss of Gbpa and b. Gbpd is divided in three groups and is present in most muroids suggesting that a duplication event occurred in the common ancestor of Muridae and Cricetidae. However, in Grammomys surdaster and Mus caroli, Gbpd2 is absent, and in Arvicanthis niloticus, Gbpd1 appears to have been deleted. Our results further demonstrated that primate GBP1, GBP3 and GBP7 are absent from the genome of muroids and showed that the Gbp gene annotations in muroids were incorrect. We propose a new classification based on the phylogenetic analyses and the divergence between the groups. Extrapolations to humans based on functional studies of muroid Gbps should be re-evaluated. The evolutionary analyses of muroid Gbp genes provided new insights about the evolution and function of these genes.