Concerted and birth-and-death evolution of multigene families

Electroneutral Na+/Cl- cotransport activity of zebrafish Slc12a10.1 expressed in Xenopus oocytes

Na+/Cl- cotransporter 2 (Ncc2 or Slc12a10) is a membrane transport protein that belongs to the electroneutral cation-chloride cotransporter family. The Slc12a10 gene (slc12a10) is widely present in bony vertebrates but is deleted or pseudogenized in birds, some bony fishes, and most mammals. Slc12a10 is highly homologous to Ncc (Slc12a3 or Ncc1); however, there are only a few reports measuring the activity of Slc12a10. In this study, we focused on zebrafish Slc12a10.1 (zSlc12a10.1) and analyzed its activity using Xenopus oocyte electrophysiology. Analysis using Na+-selective microelectrodes showed that intracellular sodium activity (aNai) in zSlc12a10.1 oocytes was significantly decreased in Na+- or Cl--free medium and recovered when Na+ or Cl- was readded to the medium. Similar analysis using a Cl--selective microelectrode showed that intracellular chloride activity (aCli) in zSlc12a10.1 oocytes significantly decreased in Na+- or Cl--free medium and recovered when Na+ or Cl- was readded to the medium. When a similar experiment was performed with a voltage clamp, the membrane current did not change when aNai of zSlc12a10.1 oocytes was decreased in Na+-free medium. Molecular phylogenetic and synteny analyses suggest that gene duplication between slc12a10.2 and slc12a10.3 in zebrafish is a relatively recent event, whereas gene duplication between slc12a10.1 and the ancestral gene of slc12a10.2/slc12a10.3 occurred at least about 2 million years ago. slc12a10 deficiency was observed in species belonging to Ictaluridae, Salmoniformes, Osmeriformes, Batrachoididae, Syngnathiformes, Gobiesociformes, Labriformes, and Tetraodontiformes. These results indicate that zebrafish Slc12a10.1 is an electroneutral Na+/Cl-cotransporter and establish its evolutionary position among various teleost slc12a10 paralogs.NEW & NOTEWORTHY Na+/Cl- cotransporter 2 (Slc12a10; Ncc2) is a protein highly homologous to Ncc (Slc12a3; Ncc1); however, there are only a few reports measuring the activity of Slc12a10. Electrophysiological analysis of Xenopus oocytes expressing zebrafish Slc12a10.1 showed that Slc12a10.1 acts as an electroneutral Na+/Cl-cotransporter. This is the third report on the activity of Slc12a10, following previous reports on Slc12a10 in eels.

Selection Across the Three-Dimensional Structure of Venom Proteins from North American Scolopendromorph Centipedes

Gene duplication followed by nucleotide differentiation is one of the simplest mechanisms to develop new functions for genes. However, the evolutionary processes underlying the divergence of multigene families remain controversial. We used multigene families found within the diversity of toxic proteins in centipede venom to test two hypotheses related to venom evolution: the two-speed mode of venom evolution and the rapid accumulation of variation in exposed residues (RAVER) model. The two-speed mode of venom evolution proposes that different types of selection impact ancient and younger venomous lineages with negative selection being the predominant form in ancient lineages and positive selection being the dominant form in younger lineages. The RAVER hypothesis proposes that, instead of different types of selection acting on different ages of venomous lineages, the different types of selection will selectively contribute to amino acid variation based on whether the residue is exposed to the solvent where it can potentially interact directly with toxin targets. This hypothesis parallels the longstanding understanding of protein evolution that suggests that residues found within the structural or active regions of the protein will be under negative or purifying selection, and residues that do not form part of these areas will be more prone to positive selection. To test these two hypotheses, we compared the venom of 26 centipedes from the order Scolopendromorpha from six currently recognized species from across North America using both transcriptomics and proteomics. We first estimated their phylogenetic relationships and uncovered paraphyly among the genus Scolopendra and evidence for cryptic diversity among currently recognized species. Using our phylogeny, we then characterized the diverse venom components from across the identified clades using a combination of transcriptomics and proteomics. We conducted selection-based analyses in the context of predicted three-dimensional properties of the venom proteins and found support for both hypotheses. Consistent with the two-speed hypothesis, we found a prevalence of negative selection across all proteins. Consistent with the RAVER hypothesis, we found evidence of positive selection on solvent-exposed residues, with structural and less-exposed residues showing stronger signal for negative selection. Through the use of phylogenetics, transcriptomics, proteomics, and selection-based analyses, we were able to describe the evolution of venom from an ancient venomous lineage and support principles of protein evolution that directly relate to multigene family evolution.

Copy number and sequence variation in rDNA of Daphnia pulex from natural populations: insights from whole-genome sequencing

Ribosomal DNA (rDNA) has a vital role in ribosome biogenesis as it contains the genes that encode ribosomal RNA (rRNA) separated by intergenic spacers (IGSs). The rRNA genes occur in hundreds to tens of thousands of copies per haploid genome in eukaryotes and are generally highly conserved with low variation within species. Due to the repetitive nature and large size of rDNA arrays, detecting intraindividual variation can be difficult. In this study, we use whole-genome sequences of 169 Daphnia pulex individuals from 10 natural populations to measure the copy number and sequence variation in rDNA. This revealed that variation in rDNA copy number between individuals spans an order of magnitude. We further observed a substantial level of sequence variation within individual genomes. As expected, single-nucleotide polymorphisms occurred in regions of lower functional constraint such as the IGS and expansion segments of the rRNA genes. The presence of strong linkage disequilibrium among variants facilitated identification of haplotypes within each population. Although there was evidence of recombination among haplotypes from different populations, it is insufficient to eliminate linkage disequilibrium within populations. Estimating copy number and haplotype diversity within individuals revealed that the level of intraindividual sequence variation is not strongly correlated with copy number. The observed patterns of variation highlight a complex evolutionary history of rDNA in D. pulex. Future research should explore the functional implications of rDNA copy number and sequence variation on organismal phenotypes.

Variations in the expression of odorant binding and chemosensory proteins in the developmental stages of whitefly Bemisia tabaci Asia II-1

The cotton whitefly, Bemisia tabaci, is considered as a species complex with 46 cryptic species, with Asia II-1 being predominant in Asia. This study addresses a significant knowledge gap in the characterization of odorant-binding proteins (OBPs) and chemosensory proteins (CSPs) in Asia II-1. We explored the expression patterns of OBPs and CSPs throughout their developmental stages and compared the motif patterns of these proteins. Significant differences in expression patterns were observed for the 14 OBPs and 14 CSPs of B. tabaci Asia II-1, with OBP8 and CSP4 showing higher expression across the developmental stages. Phylogenetic analysis reveals that OBP8 and CSP4 form distinct clades, with OBP8 appearing to be an ancestral gene, giving rise to the evolution of other odorant-binding proteins in B. tabaci. The genomic distribution of OBPs and CSPs highlights gene clustering on the chromosomes, suggesting functional conservation and evolutionary events following the birth-and-death model. Molecular docking studies indicate strong binding affinities of OBP8 and CSP4 with various odour compounds like β-caryophyllene, α-pinene, β-pinene and limonene, reinforcing their roles in host recognition and reproductive functions. This study elaborates on our understanding of the putative roles of different OBPs and CSPs in B. tabaci Asia II-1, hitherto unexplored. The dynamics of the expression of OBPs and CSPs and their interactions with odour compounds offer scope for developing innovative methods for controlling this global invasive pest.

GALEON: a comprehensive bioinformatic tool to analyse and visualize gene clusters in complete genomes

MOTIVATION

Gene clusters, defined as a set of genes encoding functionally related proteins, are abundant in eukaryotic genomes. Despite the increasing availability of chromosome-level genomes, the comprehensive analysis of gene family evolution remains largely unexplored, particularly for large and highly dynamic gene families or those including very recent family members. These challenges stem from limitations in genome assembly contiguity, particularly in repetitive regions such as large gene clusters. Recent advancements in sequencing technology, such as long reads and chromatin contact mapping, hold promise in addressing these challenges.

RESULTS

To facilitate the identification, analysis, and visualization of physically clustered gene family members within chromosome-level genomes, we introduce GALEON, a user-friendly bioinformatic tool. GALEON identifies gene clusters by studying the spatial distribution of pairwise physical distances among gene family members along with the genome-wide gene density. The pipeline also enables the simultaneous analysis and comparison of two gene families and allows the exploration of the relationship between physical and evolutionary distances. This tool offers a novel approach for studying the origin and evolution of gene families.

AVAILABILITY AND IMPLEMENTATION

GALEON is freely available from https://www.ub.edu/softevol/galeon and https://github.com/molevol-ub/galeon.

Fish antifreeze protein origin in sculpins by frameshifting within a duplicated housekeeping gene

Antifreeze proteins (AFPs) are found in a variety of marine cold-water fishes where they prevent freezing by binding to nascent ice crystals. Their diversity (types I, II, III and antifreeze glycoproteins), as well as their scattered taxonomic distribution hint at their complex evolutionary history. In particular, type I AFPs appear to have arisen in response to the Late Cenozoic Ice Age that began ~ 34 million years ago via convergence in four different groups of fish that diverged from lineages lacking this AFP. The progenitor of the alanine-rich α-helical type I AFPs of sculpins has now been identified as lunapark, an integral membrane protein of the endoplasmic reticulum. Following gene duplication and loss of all but three of the 15 exons, the final exon, which encoded a glutamate- and glutamine-rich segment, was converted to an alanine-rich sequence by a combination of frameshifting and mutation. Subsequent gene duplications produced numerous isoforms falling into four distinct groups. The origin of the flounder type I AFP is quite different. Here, a small segment from the original antiviral protein gene was amplified and the rest of the coding sequence was lost, while the gene structure was largely retained. The independent origins of type I AFPs with up to 83% sequence identity in flounder and sculpin demonstrate strong convergent selection at the level of protein sequence for alanine-rich single alpha helices that bind to ice. Recent acquisition of these AFPs has allowed sculpins to occupy icy seawater niches with reduced competition and predation from other teleost species.

The Genome Organization of 5S rRNA Genes in the Model Organism Tribolium castaneum and Its Sibling Species Tribolium freemani

5S ribosomal DNAs (rDNAs) are arranged in tandem and are often under-represented in genome assemblies. In the present study, we performed a global and in-depth analysis of the 5S rDNAs in the model insect Tribolium castaneum and its closely related species Tribolium freemani. To accomplish this goal, we used our recently published genome assemblies based on Nanopore and PacBio long-read sequencing. Although these closely related species share the 5S rRNA gene sequence with high homology, they show a different organization of the 5S rDNA locus. Analysis of 5S rDNA arrays in T. castaneum revealed a typical tandemly repeated organization characterized by repeat units consisting of the 121 bp long 5S rRNA gene and the 71 bp long nontranscribed spacer (NTS). In contrast, T. freemani showed a much more complex organization of 5S rDNA arrays characterized by two patterns. The first is based on the association of 5S rRNA gene with arrays of a satellite DNA, representing the NTS sequence of the 5S rDNA genes in T. freemani. The second, more complex type is characterized by a somewhat less frequent occurrence of the 5S rRNA gene and its association with longer satellite DNA arrays that are regularly interrupted by Jockey-like retrotransposons. This organization, in which the ribosomal gene is associated with two completely different repetitive elements such as satellite DNAs and retrotransposons, suggests that the 5S rRNA gene, regardless of its crucial function in the genome, could be a subject of extremely dynamic genomic rearrangements.

Diversity of ribosomes at the level of rRNA variation associated with human health and disease

Ribosomal DNA and RNA (rDNA and rRNA) sequences are usually discarded from sequencing analyses. But with hundreds of copies of rDNA genes it is unknown whether they possess sequence variations that form different types of ribosomes that affect human physiology and disease. Here, we developed an algorithm for variant-calling between paralog genes (termed RGA) and compared rDNA variations found in short- and long-read sequencing data from the 1,000 Genomes Project (1KGP) and Genome In A Bottle (GIAB). We additionally developed a novel protocol for long-read sequencing full-length rRNA (RIBO-RT) from actively translating ribosomes. Our analyses identified hundreds of rDNA variants, most of which, surprisingly, are short insertion-deletions (indels) and dozens of highly abundant rRNA variants that are incorporated into translationally active ribosomes. To visualize variant ribosomes at the single cell level, we developed an in-situ rRNA sequencing method (SWITCH-seq) which revealed that variants are co-expressed within individual cells. Strikingly, by analyzing rDNA, we found that variants assemble into distinct ribosome subtypes. We discovered that these subtypes acquire different rRNA structures by successfully employing dimethyl sulfate (DMS) probing of full length rRNA. With this atlas we investigated rRNA variation changes across human tissues and cancer types. This revealed tissue-specific rRNA subtype expression in endoderm/ectoderm-derived tissues. In cancer, low abundant rRNA variants can become highly expressed, which suggests the presence of cancer-specific ribosomes. Together, this study identifies and comprehensively characterizes the diversity of ribosomes at the level of rRNA variants which is dominated by indel variants, their chromosomal location and unique structure as well as the association of ribosome variation with tissue-specific biology and cancer.

Cis-regulatory evolution of the recently expanded Ly49 gene family

Comparative genomics has revealed the rapid expansion of multiple gene families involved in immunity. Members within each gene family often evolved distinct roles in immunity. However, less is known about the evolution of their epigenome and cis-regulation. Here we systematically profile the epigenome of the recently expanded murine Ly49 gene family that mainly encode either inhibitory or activating surface receptors on natural killer cells. We identify a set of cis-regulatory elements (CREs) for activating Ly49 genes. In addition, we show that in mice, inhibitory and activating Ly49 genes are regulated by two separate sets of proximal CREs, likely resulting from lineage-specific losses of CRE activity. Furthermore, we find that some Ly49 genes are cross-regulated by the CREs of other Ly49 genes, suggesting that the Ly49 family has begun to evolve a concerted cis-regulatory mechanism. Collectively, we demonstrate the different modes of cis-regulatory evolution for a rapidly expanding gene family.

Nematode histone H2A variant evolution reveals diverse histories of retention and loss and evidence for conserved core-like variant histone genes

Histone variants are paralogs that replace canonical histones in nucleosomes, often imparting novel functions. However, how histone variants arise and evolve is poorly understood. Reconstruction of histone protein evolution is challenging due to large differences in evolutionary rates across gene lineages and sites. Here we used intron position data from 108 nematode genomes in combination with amino acid sequence data to find disparate evolutionary histories of the three H2A variants found in Caenorhabditis elegans: the ancient H2A.ZHTZ-1, the sperm-specific HTAS-1, and HIS-35, which differs from the canonical S-phase H2A by a single glycine-to-alanine C-terminal change. Although the H2A.ZHTZ-1 protein sequence is highly conserved, its gene exhibits recurrent intron gain and loss. This pattern suggests that specific intron sequences or positions may not be important to H2A.Z functionality. For HTAS-1 and HIS-35, we find variant-specific intron positions that are conserved across species. Patterns of intron position conservation indicate that the sperm-specific variant HTAS-1 arose more recently in the ancestor of a subset of Caenorhabditis species, while HIS-35 arose in the ancestor of Caenorhabditis and its sister group, including the genus Diploscapter. HIS-35 exhibits gene retention in some descendent lineages but gene loss in others, suggesting that histone variant use or functionality can be highly flexible. Surprisingly, we find the single amino acid differentiating HIS-35 from core H2A is ancestral and common across canonical Caenorhabditis H2A sequences. Thus, we speculate that the role of HIS-35 lies not in encoding a functionally distinct protein, but instead in enabling H2A expression across the cell cycle or in distinct tissues. This work illustrates how genes encoding such partially-redundant functions may be advantageous yet relatively replaceable over evolutionary timescales, consistent with the patchwork pattern of retention and loss of both genes. Our study shows the utility of intron positions for reconstructing evolutionary histories of gene families, particularly those undergoing idiosyncratic sequence evolution.

Copy number variation and elevated genetic diversity at immune trait loci in Atlantic and Pacific herring

BACKGROUND

Genome-wide comparisons of populations are widely used to explore the patterns of nucleotide diversity and sequence divergence to provide knowledge on how natural selection and genetic drift affect the genome. In this study we have compared whole-genome sequencing data from Atlantic and Pacific herring, two sister species that diverged about 2 million years ago, to explore the pattern of genetic differentiation between the two species.

RESULTS

The genome comparison of the two species revealed high genome-wide differentiation but with islands of remarkably low genetic differentiation, as measured by an FST analysis. However, the low FST observed in these islands is not caused by low interspecies sequence divergence (dxy) but rather by exceptionally high estimated intraspecies nucleotide diversity (π). These regions of low differentiation and elevated nucleotide diversity, termed high-diversity regions in this study, are not enriched for repeats but are highly enriched for immune-related genes. This enrichment includes genes from both the adaptive immune system, such as immunoglobulin, T-cell receptor and major histocompatibility complex genes, as well as a substantial number of genes with a role in the innate immune system, e.g. novel immune-type receptor, tripartite motif and tumor necrosis factor receptor genes. Analysis of long-read based assemblies from two Atlantic herring individuals revealed extensive copy number variation in these genomic regions, indicating that the elevated intraspecies nucleotide diversities were partially due to the cross-mapping of short reads.

CONCLUSIONS

This study demonstrates that copy number variation is a characteristic feature of immune trait loci in herring. Another important implication is that these loci are blind spots in classical genome-wide screens for genetic differentiation using short-read data, not only in herring, likely also in other species harboring qualitatively similar variation at immune trait loci. These loci stood out in this study because of the relatively high genome-wide baseline for FST values between Atlantic and Pacific herring.

Contrasting patterns of 5S rDNA repeats in European and Asian ecotypes of greater duckweed, Spirodela polyrhiza (Lemnaceae)

Ribosomal DNA (rDNA) contains highly conserved, specifically organized sequences encoding ribosomal RNAs (rRNAs) separated by variable non-transcribed intergenic spacers (NTSs) and is abundant in eukaryotic genomes. These characteristics make the rDNA an informative molecular target to study genome organization, molecular evolution, and phylogenetics. In this study, we characterized the 5S rDNA repeats in the greater duckweed Spiroldela polyrhiza, a species known for its small size, rapid growth, highly conserved genome organization, and low mutation rate. Sequence analysis of at least 12 individually cloned PCR fragments containing the 5S rDNA units for each of six ecotypes that originated from Europe (Ukraine) and Asia (China) revealed two distinct types of 5S rDNA repeats containing NTSs of different lengths and nucleotide compositions. The shorter 5S rDNA repeat units had a highly homogeneous 400-bp NTS, with few ecotype- or region-specific single-nucleotide polymorphisms (SNPs). The longer 5S rDNA units had NTSs of 1056-1084 bp with characteristic intra- and inter-genomic variants due to specific SNPs and insertions/deletions of 4-15-bp DNA elements. We also detected significant variability in the ratio of short/long 5S rDNA variants between ecotypes of S. polyrhiza. The contrasting dynamics of the two types of 5S rDNA units, combined with the unusually low repeat copy number (for plants) in S. polyrhiza (46-220 copies per genome), shows that this species could serve as an excellent model for examining the mechanisms of concerted evolution and functional significance of rDNA variability.

Convergent gene losses and pseudogenizations in multiple lineages of stomachless fishes

The regressive evolution of independent lineages often results in convergent phenotypes. Several teleost groups display secondary loss of the stomach, and four gastric genes, atp4a, atp4b, pgc, and pga2 have been co-deleted in agastric (stomachless) fish. Analyses of genotypic convergence among agastric fishes showed that four genes, slc26a9, kcne2, cldn18a, and vsig1, were co-deleted or pseudogenized in most agastric fishes of the four major groups. kcne2 and vsig1 were also deleted or pseudogenized in the agastric monotreme echidna and platypus, respectively. In the stomachs of sticklebacks, these genes are expressed in gastric gland cells or surface epithelial cells. An ohnolog of cldn18 was retained in some agastric teleosts but exhibited an increased non-synonymous substitution when compared with gastric species. These results revealed novel convergent gene losses at multiple loci among the four major groups of agastric fish, as well as a single gene loss in the echidna and platypus.

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

Animals sense and respond to nutrient availability in their environments, a task coordinated in part by the mTOR complex 1 (mTORC1) pathway. mTORC1 regulates growth in response to nutrients and, in mammals, senses specific amino acids through specialized sensors that bind the GATOR1/2 signaling hub. Given that animals can occupy diverse niches, we hypothesized that the pathway might evolve distinct sensors in different metazoan phyla. Whether such customization occurs, and how the mTORC1 pathway might capture new inputs, is unknown. Here, we identify the Drosophila melanogaster protein Unmet expectations (CG11596) as a species-restricted methionine sensor that directly binds the fly GATOR2 complex in a fashion antagonized by S-adenosylmethionine (SAM). We find that in Dipterans GATOR2 rapidly evolved the capacity to bind Unmet and to thereby repurpose a previously independent methyltransferase as a SAM sensor. Thus, the modular architecture of the mTORC1 pathway allows it to co-opt preexisting enzymes to expand its nutrient sensing capabilities, revealing a mechanism for conferring evolvability on an otherwise conserved system.

Comparative genomic analysis of chemosensory-related gene families in gastropods

Chemoreception is critical for the survival and reproduction of animals. Except for a reduced group of insects and chelicerates, the molecular identity of chemosensory proteins is poorly understood in invertebrates. Gastropoda is the extant mollusk class with the greatest species richness, including marine, freshwater, and terrestrial lineages, and likely, highly diverse chemoreception systems. Here, we performed a comprehensive comparative genome analysis taking advantage of the chromosome-level information of two Gastropoda species, one of which belongs to a lineage that underwent a whole genome duplication event. We identified thousands of previously uncharacterized chemosensory-related genes, the majority of them encoding G protein-coupled receptors (GPCR), mostly organized into clusters distributed across all chromosomes. We also detected gene families encoding degenerin epithelial sodium channels (DEG-ENaC), ionotropic receptors (IR), sensory neuron membrane proteins (SNMP), Niemann-Pick type C2 (NPC2) proteins, and lipocalins, although with a lower number of members. Our phylogenetic analysis of the GPCR gene family across protostomes revealed: (i) remarkable gene family expansions in Gastropoda; (ii) clades including members from all protostomes; and (iii) species-specific clades with a substantial number of receptors. For the first time, we provide new and valuable knowledge into the evolution of the chemosensory gene families in invertebrates other than arthropods.

The Dynamic Interplay Between Ribosomal DNA and Transposable Elements: A Perspective From Genomics and Cytogenetics

Although both are salient features of genomes, at first glance ribosomal DNAs and transposable elements are genetic elements with not much in common: whereas ribosomal DNAs are mainly viewed as housekeeping genes that uphold all prime genome functions, transposable elements are generally portrayed as selfish and disruptive. These opposing characteristics are also mirrored in other attributes: organization in tandem (ribosomal DNAs) versus organization in a dispersed manner (transposable elements); evolution in a concerted manner (ribosomal DNAs) versus evolution by diversification (transposable elements); and activity that prolongs genomic stability (ribosomal DNAs) versus activity that shortens it (transposable elements). Re-visiting relevant instances in which ribosomal DNA-transposable element interactions have been reported, we note that both repeat types share at least four structural and functional hallmarks: (1) they are repetitive DNAs that shape genomes in evolutionary timescales, (2) they exchange structural motifs and can enter co-evolution processes, (3) they are tightly controlled genomic stress sensors playing key roles in senescence/aging, and (4) they share common epigenetic marks such as DNA methylation and histone modification. Here, we give an overview of the structural, functional, and evolutionary characteristics of both ribosomal DNAs and transposable elements, discuss their roles and interactions, and highlight trends and future directions as we move forward in understanding ribosomal DNA-transposable element associations.

Multiple amino acid changes are responsible for the shift of tuning breadth along the evolutionary trajectory of a moth pheromone receptor

Sex pheromone recognition is essential for mating in many insects and plays a major role in maintaining reproductive barriers. A previous study from our lab reported the evolutionary history of the pheromone receptor OR5 in Spodoptera moths. Using heterologous expression in Xenopus oocytes and site-directed mutagenesis, we found that eight amino acid substitutions were sufficient to recapitulate the evolution from an ancestral broadly-tuned to a highly specific receptor. Here, we confirmed this result using expression in Drosophila olfactory neurons. This further confirmed that multiple amino acid changes explain the shift in tuning breadth of Spodoptera OR5 during evolution.

Fossil-calibrated phylogenies of Southern cave wētā show dispersal and extinction confound biogeographic signal

The biota of continents and islands are commonly considered to have a source-sink relationship, but small islands can harbour distinctive taxa. The distribution of four monotypic genera of Orthoptera on young subantarctic islands indicates a role for long-distance dispersal and extinction. Phylogenetic relationships were inferred from whole mtDNA genomes and nuclear sequences (45S cassette; four histones). We used a fossil and one palaeogeographic event to calibrate molecular clock analysis. We confirm that neither the Australian nor Aotearoa-New Zealand Rhaphidophoridae faunas are monophyletic. The radiation of Macropathinae may have begun in the late Jurassic, but trans-oceanic dispersal is required to explain the current distribution of some lineages within this subfamily. Dating the most recent common ancestor of seven island endemic species with their nearest mainland relative suggests that each existed long before their island home was available. Time estimates from our fossil-calibrated molecular clock analysis suggest several lineages have not been detected on mainland New Zealand, Australia, or elsewhere most probably due to their extinction, providing evidence that patterns of extinction, which are not consistently linked to range size or lineage age, confound biogeographic signal.

Phylogenomics of the Ecdysteroid Kinase-like (EcKL) Gene Family in Insects Highlights Roles in Both Steroid Hormone Metabolism and Detoxification

The evolutionary dynamics of large gene families can offer important insights into the functions of their individual members. While the ecdysteroid kinase-like (EcKL) gene family has previously been linked to the metabolism of both steroid molting hormones and xenobiotic toxins, the functions of nearly all EcKL genes are unknown, and there is little information on their evolution across all insects. Here, we perform comprehensive phylogenetic analyses on a manually annotated set of EcKL genes from 140 insect genomes, revealing the gene family is comprised of at least 13 subfamilies that differ in retention and stability. Our results show the only two genes known to encode ecdysteroid kinases belong to different subfamilies and therefore ecdysteroid metabolism functions must be spread throughout the EcKL family. We provide comparative phylogenomic evidence that EcKLs are involved in detoxification across insects, with positive associations between family size and dietary chemical complexity, and we also find similar evidence for the cytochrome P450 and glutathione S-transferase gene families. Unexpectedly, we find that the size of the clade containing a known ecdysteroid kinase is positively associated with host plant taxonomic diversity in Lepidoptera, possibly suggesting multiple functional shifts between hormone and xenobiotic metabolism. Our evolutionary analyses provide hypotheses of function and a robust framework for future experimental studies of the EcKL gene family. They also open promising new avenues for exploring the genomic basis of dietary adaptation in insects, including the classically studied coevolution of butterflies with their host plants.

Evolving Together: Cassandra Retrotransposons Gradually Mirror Promoter Mutations of the 5S rRNA Genes

The 5S rRNA genes are among the most conserved nucleotide sequences across all species. Similar to the 5S preservation we observe the occurrence of 5S-related nonautonomous retrotransposons, so-called Cassandras. Cassandras harbor highly conserved 5S rDNA-related sequences within their long terminal repeats, advantageously providing them with the 5S internal promoter. However, the dynamics of Cassandra retrotransposon evolution in the context of 5S rRNA gene sequence information and structural arrangement are still unclear, especially: (1) do we observe repeated or gradual domestication of the highly conserved 5S promoter by Cassandras and (2) do changes in 5S organization such as in the linked 35S-5S rDNA arrangements impact Cassandra evolution? Here, we show evidence for gradual co-evolution of Cassandra sequences with their corresponding 5S rDNAs. To follow the impact of 5S rDNA variability on Cassandra TEs, we investigate the Asteraceae family where highly variable 5S rDNAs, including 5S promoter shifts and both linked and separated 35S-5S rDNA arrangements have been reported. Cassandras within the Asteraceae mirror 5S rDNA promoter mutations of their host genome, likely as an adaptation to the host's specific 5S transcription factors and hence compensating for evolutionary changes in the 5S rDNA sequence. Changes in the 5S rDNA sequence and in Cassandras seem uncorrelated with linked/separated rDNA arrangements. We place all these observations into the context of angiosperm 5S rDNA-Cassandra evolution, discuss Cassandra's origin hypotheses (single or multiple) and Cassandra's possible impact on rDNA and plant genome organization, giving new insights into the interplay of ribosomal genes and transposable elements.

Evolutionary and functional characterization of lagomorph guanylate-binding proteins: a story of gain and loss and shedding light on expression, localization and innate immunity-related functions

Guanylate binding proteins (GBPs) are an evolutionarily ancient family of proteins that are widely distributed among eukaryotes. They belong to the dynamin superfamily of GTPases, and their expression can be partially induced by interferons (IFNs). GBPs are involved in the cell-autonomous innate immune response against bacterial, parasitic and viral infections. Evolutionary studies have shown that GBPs exhibit a pattern of gene gain and loss events, indicative for the birth-and-death model of evolution. Most species harbor large GBP gene clusters that encode multiple paralogs. Previous functional and in-depth evolutionary studies have mainly focused on murine and human GBPs. Since rabbits are another important model system for studying human diseases, we focus here on lagomorphs to broaden our understanding of the multifunctional GBP protein family by conducting evolutionary analyses and performing a molecular and functional characterization of rabbit GBPs. We observed that lagomorphs lack GBP3, 6 and 7. Furthermore, Leporidae experienced a loss of GBP2, a unique duplication of GBP5 and a massive expansion of GBP4. Gene expression analysis by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) and transcriptome data revealed that leporid GBP expression varied across tissues. Overexpressed rabbit GBPs localized either uniformly and/or discretely to the cytoplasm and/or to the nucleus. Oryctolagus cuniculus (oc)GBP5L1 and rarely ocGBP5L2 were an exception, colocalizing with the trans-Golgi network (TGN). In addition, four ocGBPs were IFN-inducible and only ocGBP5L2 inhibited furin activity. In conclusion, from an evolutionary perspective, lagomorph GBPs experienced multiple gain and loss events, and the molecular and functional characteristics of ocGBP suggest a role in innate immunity.

rDNA and mtDNA analysis for the identification of genetic characters in the hybrid grouper derived from hybridization of Cromileptes altivelis (female) × Epinephelus lanceolatus (male)

BACKGROUND

Hybridization is a useful strategy to produce offspring with more desirable phenotypic characteristics than those of parents. The hybrid grouper derived from the cross of Cromileptes altivelis (♀, 2n = 48) with Epinephelus lanceolatus (♂, 2n = 48) exhibits improved growth compared with its female parent, which makes it valuable to aquaculture. However, the genetic traits of the hybrid grouper are poorly understood.

RESULTS

The observations showed that the hybrid grouper was diploid (2n = 48) and displayed intermediate morphology with the parent's measurable characteristics. The ribosomal DNA (rDNA) and mitochondria DNA (mtDNA) were characterized at molecular and phylogenetic level. High similarity and low genetic distance of 5S rDNA and mtDNA sequences between the hybrid grouper and C. altivelis showed that the hybrid grouper had a closer genetic relationship with female parents. The reconstructed phylogenetic tree based on COI gene and D-loop region of mtDNA recovered that mtDNA was maternally inherited in the hybrid grouper. Additionally, the DNA methylation level of 5S rDNA intergenic spacers (IGS) sequence was tested in here. The results showed that the DNA methylation status of the hybrid grouper was significantly lower than that of C. altivelis.

CONCLUSION

Results of this study provide important data on the genetic characteristics of the hybrid derived from the cross of C. altivelis and E. lanceolatus, and contribute the knowledge of both evolution and marine fish breeding.

MHCtools 1.5: Analysis of MHC Sequencing Data in R

The genes of the major histocompatibility complex (MHC) play a vital role in the vertebrate immune system and have attracted considerable interest in evolutionary biology. While the MHC has been characterized in detail in humans (human leukocyte antigen, HLA) and in model organisms such as the mouse, studies in non-model organisms often lack prior knowledge about structure, genetic variability, and evolutionary properties of this locus. MHC genotyping in non-model species commonly relies on PCR-based amplicon sequencing, and while several published protocols facilitate generation of MHC sequence data, there is a lack of transparent and standardized tools for downstream data analysis.Here, I present the R package MHCtools version 1.5, which contains 15 tools that (i) assist accurate MHC genotyping from high-throughput amplicon sequencing data, and provide standardized methods to analyze (ii) MHC diversity, (iii) MHC supertypes, and (iv) MHC haplotypes.I hope that MHCtools will be helpful in future studies of the MHC in non-model species and that it may help to advance our understanding of the important roles of the MHC in ecology and evolution.

A vertebrate-wide catalogue of T1R receptors reveals diversity in taste perception

Taste is a vital chemical sense for feeding behaviour. In mammals, the umami and sweet taste receptors comprise three members of the taste receptor type 1 (T1R/TAS1R) family: T1R1, T1R2 and T1R3. Because their functional homologues exist in teleosts, only three TAS1R genes generated by gene duplication are believed to have been inherited from the common ancestor of bony vertebrates. Here, we report five previously uncharacterized TAS1R members in vertebrates, TAS1R4, TAS1R5, TAS1R6, TAS1R7 and TAS1R8, based on genome-wide survey of diverse taxa. We show that mammalian and teleost fish TAS1R2 and TAS1R3 genes are paralogues. Our phylogenetic analysis suggests that the bony vertebrate ancestor had nine TAS1Rs resulting from multiple gene duplications. Some TAS1Rs were lost independently in descendent lineages resulting in retention of only three TAS1Rs in mammals and teleosts. Combining functional assays and expression analysis of non-teleost fishes we show that the novel T1Rs form heterodimers in taste-receptor cells and recognize a broad range of ligands such as essential amino acids, including branched-chain amino acids, which have not been previously considered as T1R ligands. This study reveals diversity of taste sensations in both modern vertebrates and their ancestors, which might have enabled vertebrates to adapt to diverse habitats on Earth.

Regulation potential of transcribed simple repeated sequences in developing neurons

Simple repeated sequences (SRSs), defined as tandem iterations of microsatellite- to satellite-sized DNA units, occupy a substantial part of the human genome. Some of these elements are known to be transcribed in the context of repeat expansion disorders. Mounting evidence suggests that the transcription of SRSs may also contribute to normal cellular functions. Here, we used genome-wide bioinformatics approaches to systematically examine SRS transcriptional activity in cells undergoing neuronal differentiation. We identified thousands of long noncoding RNAs containing >200-nucleotide-long SRSs (SRS-lncRNAs), with hundreds of these transcripts significantly upregulated in the neural lineage. We show that SRS-lncRNAs often originate from telomere-proximal regions and that they have a strong potential to form multivalent contacts with a wide range of RNA-binding proteins. Our analyses also uncovered a cluster of neurally upregulated SRS-lncRNAs encoded in a centromere-proximal part of chromosome 9, which underwent an evolutionarily recent segmental duplication. Using a newly established in vitro system for rapid neuronal differentiation of induced pluripotent stem cells, we demonstrate that at least some of the bioinformatically predicted SRS-lncRNAs, including those encoded in the segmentally duplicated part of chromosome 9, indeed increase their expression in developing neurons to readily detectable levels. These and other lines of evidence suggest that many SRSs may be expressed in a cell type and developmental stage-specific manner, providing a valuable resource for further studies focused on the functional consequences of SRS-lncRNAs in the normal development of the human brain, as well as in the context of neurodevelopmental disorders.

Mwhitiwhiti Aotearoa: Phylogeny and synonymy of the silent alpine grasshopper radiation of New Zealand (Orthoptera: Acrididae)

Aotearoa New Zealand has a fauna of endemic alpine grasshoppers, consisting of thirteen species distributed among four genera. The many re-classifications of species within this group and the presence of species complexes highlight the uncertainty that surrounds relationships within and between these genera. High-throughput Next Generation Sequencing was used to assemble the complete mitochondrial genomes, 45S ribosomal cassettes and histone sequences of New Zealands four endemic alpine genera: Alpinacris, Brachaspis, Paprides and Sigaus. Phylogenetic analysis of these molecular datasets, as individual genes, partitions and combinations returned a consistent topology that is incompatible with the current classification. The genera Sigaus, Alpinacris, and Paprides all exhibit paraphyly. A consideration of the pronotum, epiphallus and terminalia of adult specimens reveals species-specific differences, but fails to provide compelling evidence for species groups justifying distinct genera. In combination with phylogenetic, morphological and spatial evidence we propose a simplified taxonomy consisting of a single genus for the mwhitiwhiti Aotearoa species radiation.

Tracking the Diversity and Chromosomal Distribution of the Olfactory Receptor Gene Repertoires of Three Anurans Species

Olfaction is a crucial capability for most vertebrates and is realized through olfactory receptors in the nasal cavity. The enormous diversity of olfactory receptors has been created by gene duplication, following a birth-and-death model of evolution. The olfactory receptor genes of the amphibians have received relatively little attention up to now, although recent studies have increased the number of species for which data are available. This study analyzed the diversity and chromosomal distribution of the OR genes of three anuran species (Engystomops pustulosus, Bufo bufo and Hymenochirus boettgeri). The OR genes were identified through searches for homologies, and sequence filtering and alignment using bioinformatic tools and scripts. A high diversity of OR genes was found in all three species, ranging from 917 in B. bufo to 1194 in H. boettgeri, and a total of 2076 OR genes in E. pustulosus. Six OR groups were recognized using an evolutionary gene tree analysis. While E. pustulosus has one of the highest numbers of genes of the gamma group (which detect airborne odorants) yet recorded in an anuran, B. bufo presented the smallest number of pseudogene sequences ever identified, with no pseudogenes in either the beta or epsilon groups. Although H. boettgeri shares many morphological adaptations for an aquatic lifestyle with Xenopus, and presented a similar number of genes related to the detection of water-soluble odorants, it had comparatively far fewer genes related to the detection of airborne odorants. This study is the first to describe the complete OR repertoire of the three study species and represents an important contribution to the understanding of the evolution and function of the sense of smell in vertebrates.

The Hypervariable Tpr Multigene Family of Theileria Parasites, Defined by a Conserved, Membrane-Associated, C-Terminal Domain, Includes Several Copies with Defined Orthology Between Species

Multigene families often play an important role in host-parasite interactions. One of the largest multigene families in Theileria parva, the causative agent of East Coast fever, is the T. parva repeat (Tpr) gene family. The function of the putative Tpr proteins remains unknown. The initial publication of the T. parva reference genome identified 39 Tpr family open reading frames (ORFs) sharing a conserved C-terminal domain. Twenty-eight of these are clustered in a central region of chromosome 3, termed the "Tpr locus", while others are dispersed throughout all four nuclear chromosomes. The Tpr locus contains three of the four assembly gaps remaining in the genome, suggesting the presence of additional, as yet uncharacterized, Tpr gene copies. Here, we describe the use of long-read sequencing to attempt to close the gaps in the reference assembly of T. parva (located among multigene families clusters), characterize the full complement of Tpr family ORFs in the T. parva reference genome, and evaluate their evolutionary relationship with Tpr homologs in other Theileria species. We identify three new Tpr family genes in the T. parva reference genome and show that sequence similarity among paralogs in the Tpr locus is significantly higher than between genes outside the Tpr locus. We also identify sequences homologous to the conserved C-terminal domain in five additional Theileria species. Using these sequences, we show that the evolution of this gene family involves conservation of a few orthologs across species, combined with gene gains/losses, and species-specific expansions.

Odorant receptor orthologues from moths display conserved responses to cis-jasmone

In insects, the odorant receptor (OR) multigene family evolves by the birth-and-death evolutionary model, according to which the OR repertoire of each species has undergone specific gene gains and losses depending on their chemical environment, resulting in taxon-specific OR lineage radiations with different sizes in the phylogenetic trees. Despite the general divergence in the gene family across different insect orders, the ORs in moths seem to be genetically conserved across species, clustered into 23 major clades containing multiple orthologous groups with single-copy gene from each species. We hypothesized that ORs in these orthologous groups are tuned to ecologically important compounds and functionally conserved. cis-Jasmone is one of the compounds that not only primes the plant defense of neighboring receiver plants, but also functions as a behavior regulator to various insects. To test our hypothesis, using Xenopus oocyte recordings, we functionally assayed the orthologues of BmorOR56, which has been characterized as a specific receptor for cis-jasmone. Our results showed highly conserved response specificity of the BmorOR56 orthologues, with all receptors within this group exclusively responding to cis-jasmone. This is supported by the dN/dS analysis, showing that strong purifying selection is acting on this group. Moreover, molecular docking showed that the ligand binding pockets of BmorOR56 orthologues to cis-jasmone are similar. Taken together, our results suggest the high conservation of OR for ecologically important compounds across Heterocera.

Molecular Identification and Antioxidant Activity Determination among Coffee Varieties Cultivated in Nepal

Coffee is the most popular beverage containing numerous phytochemical components that have antioxidant activity capable of scavenging free radicals. Antioxidant and phenolic contents have considerable benefits for human health. The aim of this study was the molecular identification of 9 coffee samples from the Nepal Agricultural Research Council, Lalitpur, Nepal, and the determination of the antioxidant activity and total phenolic content of green and roasted coffee beans. Molecular identification was performed using ITS-specific PCR followed by sequencing and phylogenetic tree construction using the maximum parsimony method. The DPPH assay was used to determine the antioxidant activity, and the Folin-Ciocalteu (F-C) assay was used to determine the total phenolic content. All the samples belonged to the taxa Coffea arabica. The antioxidant activity in roasted beans varied from 2.49 to 4.62 AAE mg/g and from 1.4 to 3.9 AAE mg/g in green beans. The total phenolic content varied from 2.58 to 3.38 GAE mg/g and from 4.16 to 5.36 GAE mg/g for the roasted beans and green beans, respectively. The data revealed that the highest antioxidant content (4.62 AAE mg/g) was found in roasted coffee and that the highest phenolic content (5.36 GAE mg/g) was found in green coffee. The study concludes that roasting increases the antioxidant activity but decreases the phenolic content of coffee.

Single-haplotype comparative genomics provides insights into lineage-specific structural variation during cat evolution

The role of structurally dynamic genomic regions in speciation is poorly understood due to challenges inherent in diploid genome assembly. Here we reconstructed the evolutionary dynamics of structural variation in five cat species by phasing the genomes of three interspecies F1 hybrids to generate near-gapless single-haplotype assemblies. We discerned that cat genomes have a paucity of segmental duplications relative to great apes, explaining their remarkable karyotypic stability. X chromosomes were hotspots of structural variation, including enrichment with inversions in a large recombination desert with characteristics of a supergene. The X-linked macrosatellite DXZ4 evolves more rapidly than 99.5% of the genome clarifying its role in felid hybrid incompatibility. Resolved sensory gene repertoires revealed functional copy number changes associated with ecomorphological adaptations, sociality and domestication. This study highlights the value of gapless genomes to reveal structural mechanisms underpinning karyotypic evolution, reproductive isolation and ecological niche adaptation.

Phylogeny, Ecology, and Gene Families Covariation Shaped the Olfactory Subgenome of Rodents

Olfactory receptor (OR) genes represent the largest multigenic family in mammalian genomes and encode proteins that bind environmental odorant molecules. The OR repertoire is extremely variable among species and is subject to many gene duplications and losses, which have been linked to ecological adaptations in mammals. Although they have been studied on a broad taxonomic scale (i.e., placental), finer sampling has rarely been explored in order to better capture the mechanisms that drove the evolution of the OR repertoire. Among placental mammals, rodents are well-suited for this task, as they exhibit diverse life history traits, and genomic data are available for most major families and a diverse array of lifestyles. In this study, 53 rodent published genomes were mined for their OR subgenomes. We retrieved more than 85,000 functional and pseudogene OR sequences that were subsequently classified into phylogenetic clusters. Copy number variation among rodents is similar to that of other mammals. Using our OR counts along with comparative phylogenetic approaches, we demonstrated that ecological niches such as diet, period of activity, and a fossorial lifestyle strongly impacted the proportion of OR pseudogenes. Within the OR subgenome, phylogenetic inertia was the main factor explaining the relative variations of the 13 OR gene families. However, a striking exception was a convergent 10-fold expansion of the OR family 14 among the phylogenetically divergent subterranean mole-rat lineages belonging to Bathyergidae and Spalacidae families. This study illustrates how the diversity of the OR repertoire has evolved among rodents, both shaped by selective forces stemming from species life history traits and neutral evolution along the rodent phylogeny.

An updated synthesis of and outstanding questions in the olfactory and vomeronasal systems in bats: Genetics asks questions only anatomy can answer

The extensive diversity observed in bat nasal chemosensory systems has been well-documented at the histological level. Understanding how this diversity evolved and developing hypotheses as to why particular patterns exist require a phylogenetic perspective, which was first outlined in the work of anatomist Kunwar Bhatnagar. With the onset of genetics and genomics, it might be assumed that the puzzling patterns observed in the morphological data have been clarified. However, there is still a widespread mismatch of genetic and morphological correlations among bat chemosensory systems. Novel genomic evidence has set up new avenues to explore that demand more evidence from anatomical structures. Here, we outline the progress that has been made in both morphological and molecular studies on the olfactory and vomeronasal systems in bats since the work of Bhatnagar. Genomic data of olfactory and vomeronasal receptors demonstrate the strong need for further morphological sampling, with a particular focus on receiving brain regions, glands, and ducts.

Ancient duplication, coevolution, and selection at the MHC class IIA and IIB genes of birds

Introduction

The Major Histocompatibility Complex (MHC) of vertebrates is a dynamically evolving multigene family primarily responsible for recognizing non-self peptide antigens and triggering a pathogen-specific adaptive immune response. In birds, the MHC was previously thought to evolve via concerted evolution with high degree of gene homogenization and the rapid loss of orthology. However, the discovery of two ancient avian MHC-IIB gene lineages (DAB1 and DAB2) originating before the radiation of extant birds indicated that despite the action of concerted evolution, orthology may be detectable for long evolutionary periods.

Methods

Here, we take advantage of rapidly accumulating digital genomic resources to search for the signal of an ancient duplication at the avian MHC-IIA genes, as well as to compare phylogenetic distribution and selection between MHC-IIA and IIB gene lineages.

Results

The analysis of MHC sequences from over 230 species representing ca. 70 bird families revealed the presence of two ancient MHC-IIA gene lineages (DAA1 and DAA2) and showed that their phylogenetic distribution matches exactly the distribution of DAB1 and DAB2 lineages, suggesting tight coevolution. The early post-duplication divergence of DAA1 and DAA2 was driven by positive selection fixing radical amino acid differences within the membrane-proximal domain and, most probably, being functionally related to the interactions between α2 and β2 chains of the MHC-II heterodimer. We detected no evidence for an overall (gene-wide) relaxation or intensification of selection at either DAA1/DAB1 or DAA2/DAB2, but codon-specific differences in selection signature were found at the peptide-binding sites between the two gene lineages, perhaps implying specialization to different pathogen regimes.

Discussion

Our results suggest that specific pairing of MHC-II α and β chains may have an adaptive significance, a conclusion that advances knowledge on the macroevolution of the avian MHC-II and opens exciting novel directions for future research.

Genomic Features and Comparative Genomic Analysis of Streptococcus sp. v1. nov., Isolated from an Endophthalmitis Patient

Endophthalmitis is an acute inflammatory intraocular condition that can cause permanent vision loss. The treatment strategy and visual outcome partly depend on the identification of the agents of pathogens. In this study, metagenomic sequencing was conducted to investigate the microbial and antibiotic resistance genes (ARGs) composition in the vitreous (intraocular body fluid) of an endophthalmitis patient, who progressed rapidly and accompanied by severe pain. Metagenomic sequencing data revealed that the vitreous sample was predominated by Streptococcus, with a low-diversity microbiome in the vitreous. This strain harbor's the ARGs mainly against beta-lactam, macrolide-lincosamide-streptogramin, and multidrug. Additionally, metagenome-assembled genome sequence of Streptococcus sp. v1. nov. was identified. The Tetra Correlation Search (TCS) analysis uncovered that the closest relative of the Streptococcus sp. v1. nov. was Streptococcus mitis SK321. Pan/core genome analysis for Streptococcus sp. v1. nov. and TCS top 25 hits strains revealed that most unique genes of Streptococcus sp. v1. nov. were linked to ATP-binding cassette transport system, which could indicate unique virulence and pathogenic potentials of Streptococcus sp. v1. nov. In addition, a total of 7 virulence factors were identified, and the overwhelming of them were classified into "offensive virulence factors". The high pathogenicity of Streptococcus sp. v1. nov. could be a reason for the patient's rapid disease progression. Our study was first isolated an ocular pathogen with highly virulent based on metagenomic sequencing and bioinformatics analysis, which has important reference value for revealing the composition and genome characteristics of pathogens in endophthalmitis patient in the future.

A supergene in seaweed flies modulates male traits and female perception

Supergenes, tightly linked sets of alleles, offer some of the most spectacular examples of polymorphism persisting under long-term balancing selection. However, we still do not understand their evolution and persistence, especially in the face of accumulation of deleterious elements. Here, we show that an overdominant supergene in seaweed flies, Coelopa frigida, modulates male traits, potentially facilitating disassortative mating and promoting intraspecific polymorphism. Across two continents, the Cf-Inv(1) supergene strongly affected the composition of male cuticular hydrocarbons (CHCs) but only weakly affected CHC composition in females. Using gas chromatography-electroantennographic detection, we show that females can sense male CHCs and that there may be differential perception between genotypes. Combining our phenotypic results with RNA-seq data, we show that candidate genes for CHC biosynthesis primarily show differential expression for Cf-Inv(1) in males but not females. Conversely, candidate genes for odorant detection were differentially expressed in both sexes but showed high levels of divergence between supergene haplotypes. We suggest that the reduced recombination between supergene haplotypes may have led to rapid divergence in mate preferences as well as increasing linkage between male traits, and overdominant loci. Together this probably helped to maintain the polymorphism despite deleterious effects in homozygotes.

U2 and U4 snDNA Comparative Chromosomal Mapping in the Neotropical Fish Genera Apareiodon and Parodon (Characiformes: Parodontidae)

Small nuclear DNA (snDNA) are valuable cytogenetic markers for comparative studies in chromosome evolution because different distribution patterns were found among species. Parodontidae, a Neotropical fish family, is known to have female heterogametic sex chromosome systems in some species. The U2 and U4 snDNA sites have been found to be involved in Z and W chromosome differentiation in Apareiodon sp., Apareiodon affinis, and Parodon hilarii. However, few studies have evaluated snDNA sites as propulsors of chromosome diversification among closely related fish species. In this study, we investigated the distribution of U2 and U4 snDNA clusters in the chromosomes of 10 populations/species belonging to Apareiodon and Parodon, aiming to identify chromosomal homeologies or diversification. In situ localization data revealed a submetacentric pair carrying the U2 snDNA site among the populations/species analyzed. Furthermore, all studied species demonstrated homeology in the location of U4 snDNA cluster in the proximal region of metacentric pair 1, besides an additional signal showing up with a divergence in Apareiodon. Comparative chromosomal mapping of U4 snDNA also helped to reinforce the proposal of the ZZ/ZW1W2 sex chromosome system origin in an A. affinis population. According to cytogenetic data, the study corroborates the diversification in Parodontidae paired species with uncertain taxonomy.

Ancient fish lineages illuminate toll-like receptor diversification in early vertebrate evolution

Since its initial discovery over 50 years ago, understanding the evolution of the vertebrate RAG- mediated adaptive immune response has been a major area of research focus for comparative geneticists. However, how the evolutionary novelty of an adaptive immune response impacted the diversity of receptors associated with the innate immune response has received considerably less attention until recently. Here, we investigate the diversification of vertebrate toll-like receptors (TLRs), one of the most ancient and well conserved innate immune receptor families found across the Tree of Life, integrating genomic data that represent all major vertebrate lineages with new transcriptomic data from Polypteriformes, the earliest diverging ray-finned fish lineage. Our analyses reveal TLR sequences that reflect the 6 major TLR subfamilies, TLR1, TLR3, TLR4, TLR5, TLR7, and TLR11, and also currently unnamed, yet phylogenetically distinct TLR clades. We additionally recover evidence for a pulse of gene gain coincident with the rise of the RAG-mediated adaptive immune response in jawed vertebrates, followed by a period of rapid gene loss during the Cretaceous. These gene losses are primarily concentrated in marine teleost fish and synchronous with the mid Cretaceous anoxic event, a period of rapid extinction for marine species. Finally, we reveal a mismatch between phylogenetic placement and gene nomenclature for up to 50% of TLRs found in clades such as ray-finned fishes, cyclostomes, amphibians, and elasmobranchs. Collectively, these results provide an unparalleled perspective of TLR diversity and offer a ready framework for testing gene annotations in non-model species.

Evolution of five environmentally responsive gene families in a pine-feeding sawfly, Neodiprion lecontei (Hymenoptera: Diprionidae)

A central goal in evolutionary biology is to determine the predictability of adaptive genetic changes. Despite many documented cases of convergent evolution at individual loci, little is known about the repeatability of gene family expansions and contractions. To address this void, we examined gene family evolution in the redheaded pine sawfly Neodiprion lecontei, a noneusocial hymenopteran and exemplar of a pine-specialized lineage evolved from angiosperm-feeding ancestors. After assembling and annotating a draft genome, we manually annotated multiple gene families with chemosensory, detoxification, or immunity functions before characterizing their genomic distributions and molecular evolution. We find evidence of recent expansions of bitter gustatory receptor, clan 3 cytochrome P450, olfactory receptor, and antimicrobial peptide subfamilies, with strong evidence of positive selection among paralogs in a clade of gustatory receptors possibly involved in the detection of bitter compounds. In contrast, these gene families had little evidence of recent contraction via pseudogenization. Overall, our results are consistent with the hypothesis that in response to novel selection pressures, gene families that mediate ecological interactions may expand and contract predictably. Testing this hypothesis will require the comparative analysis of high-quality annotation data from phylogenetically and ecologically diverse insect species and functionally diverse gene families. To this end, increasing sampling in under-sampled hymenopteran lineages and environmentally responsive gene families and standardizing manual annotation methods should be prioritized.

Genome-wide identification of the AOMT gene family in wax apple and functional characterization of SsAOMTs to anthocyanin methylation

Introduction

Anthocyanins are major pigments in the peels of red-series wax apple fruits, and two principal components of them, namely, the cyanin and the peonidin, are non-methoxylated and methoxylated anthocyanins, respectively. Anthocyanin O-methyltransferases (AOMTs) are an important group of enzymes that have the ability to catalyze anthocyanins methylation to promote the solubility, stability, and bioactivity of anthocyanins. Although AOMT genes have been studied in a variety of plants, the function of them in wax apple is generally not well understood.

Methods

The anthocyanin composition in peels of two wax apple cultivars was determined by High Performance Liquid Chromatography Tandem Mass Spectrometry (HPLS-MS). The genome-wide analysis of the AOMT genes was performed with bioinformatics technology, and the expression patterns of different plant tissues, cultivars, fruit ripening stages, and exogenous abscisic acid (ABA) treatments were analyzed by transcriptome sequencing analysis and real-time quantitative PCR verification. An initial functional evaluation was carried out in vitro using recombinant the Anthocyanin O-methyltransferase Gene 5 of S. samarangense (SsAOMT5) protein.

Results

Only two main compositions of anthocyanin were found in peels of two wax apple cultivars, and it was worth noting that Tub Ting Jiang cultivar contained non-methoxylated anthocyanin (Cy3G) only, whereas Daye cultivar contained both non-methoxylated and methoxylated (Pn3G) anthocyanins. A total of six SsAOMT genes were identified in the whole genome of wax apple, randomly distributing on three chromosomes. A phylogenic analysis of the protein sequences divided the SsAOMT gene family into three subgroups, and all SsAOMTs had highly conserved domains of AOMT family. In total, four types of stress- related and five types of hormone- related cis-elements were discovered in the promoter region of the SsAOMTs. Expression pattern analysis showed that SsAOMT5 and SsAOMT6 were expressed in all tissues to varying degrees; notably, the expression of SsAOMT5 was high in the flower and fruit and significantly higher in Daye peels than those of other cultivars in the fruit ripening period. Exogenous ABA treatment significantly increased anthocyanin accumulation, but the increase of methoxylated anthocyanin content did not reach significant level compared with those without ABA treatment, whereas the expression of SsAOMT5 upregulated under ABA treatment. We identified two homologous SsAOMT5 genes from Daye cultivar (DSsAOMT5) and Tub Ting Jiang cultivar (TSsAOMT5); the results of functional analyses to two SsAOMT5 recombinant proteins in vitro demonstrated that DSsAOMT5 showed methylation modification activity, but TSsAOMT5 did not.

Conclusion

In conclusion, SsAOMT5 was responsible for methylated anthocyanin accumulation in the peels of wax apple and played an important role in red coloration in wax apple peels.

Novel genes and alleles of the BTB/POZ protein family in Oryza rufipogon

The BTB/POZ family of proteins is widespread in plants and animals, playing important roles in development, growth, metabolism, and environmental responses. Although members of the expanded BTB/POZ gene family (OsBTB) have been identified in cultivated rice (Oryza sativa), their conservation, novelty, and potential applications for allele mining in O. rufipogon, the direct progenitor of O. sativa ssp. japonica and potential wide-introgression donor, are yet to be explored. This study describes an analysis of 110 BTB/POZ encoding gene loci (OrBTB) across the genome of O. rufipogon as outcomes of tandem duplication events. Phylogenetic grouping of duplicated OrBTB genes was supported by the analysis of gene sequences and protein domain architecture, shedding some light on their evolution and functional divergence. The O. rufipogon genome encodes nine novel BTB/POZ genes with orthologs in its distant cousins in the family Poaceae (Sorghum bicolor, Brachypodium distachyon), but such orthologs appeared to have been lost in its domesticated descendant, O. sativa ssp. japonica. Comparative sequence analysis and structure comparisons of novel OrBTB genes revealed that diverged upstream regulatory sequences and regulon restructuring are the key features of the evolution of this large gene family. Novel genes from the wild progenitor serve as a reservoir of potential new alleles that can bring novel functions to cultivars when introgressed by wide hybridization. This study establishes a foundation for hypothesis-driven functional genomic studies and their applications for widening the genetic base of rice cultivars through the introgression of novel genes or alleles from the exotic gene pool.

Intragenomic rDNA variation - the product of concerted evolution, mutation, or something in between?

The classical model of concerted evolution states that hundreds to thousands of ribosomal DNA (rDNA) units undergo homogenization, making the multiple copies of the individual units more uniform across the genome than would be expected given mutation frequencies and gene redundancy. While the universality of this over 50-year-old model has been confirmed in a range of organisms, advanced high throughput sequencing techniques have also revealed that rDNA homogenization in many organisms is partial and, in rare cases, even apparently failing. The potential underpinning processes leading to unexpected intragenomic variation have been discussed in a number of studies, but a comprehensive understanding remains to be determined. In this work, we summarize information on variation or polymorphisms in rDNAs across a wide range of taxa amongst animals, fungi, plants, and protists. We discuss the definition and description of concerted evolution and describe whether incomplete concerted evolution of rDNAs predominantly affects coding or non-coding regions of rDNA units and if it leads to the formation of pseudogenes or not. We also discuss the factors contributing to rDNA variation, such as interspecific hybridization, meiotic cycles, rDNA expression status, genome size, and the activity of effector genes involved in genetic recombination, epigenetic modifications, and DNA editing. Finally, we argue that a combination of approaches is needed to target genetic and epigenetic phenomena influencing incomplete concerted evolution, to give a comprehensive understanding of the evolution and functional consequences of intragenomic variation in rDNA.

Microbial gasdermins: More than a billion years of pyroptotic-like cell death

In the recent past, the concept of immunity has been extended to eukaryotic and prokaryotic microorganisms, like fungi and bacteria. The latest findings have drawn remarkable evolutionary parallels between metazoan and microbial defense-related genes, unveiling a growing number of shared transkingdom components of immune systems. One such component is the gasdermin family of pore-forming proteins - executioners of a highly inflammatory immune cell death program in mammals, termed pyroptosis. Pyroptotic cell death limits the spread of intracellular pathogens by eliminating infected cells and coordinates the broader inflammatory response to infection. The microbial gasdermins have similarly been implicated in defense-related cell death reactions in fungi, bacteria and archaea. Moreover, the discovery of the molecular regulators of gasdermin cytotoxicity in fungi and bacteria, has established additional evolutionary links to mammalian pyroptotic pathways. Here, we focus on the gasdermin proteins in microorganisms and their role in organismal defense and provide perspective on this remarkable case study in comparative immunology.

Unusually Divergent Ubiquitin Genes and Proteins in Plasmodium Species

Ubiquitin is an extraordinarily highly conserved 76 amino acid protein encoded by three different types of gene, where the primary translation products are fusions either of ubiquitin with one of two ribosomal proteins (RPs) or of multiple ubiquitin monomers from head to tail. Here, we investigate the evolution of ubiquitin genes in mammalian malaria parasites (Plasmodium species). The ubiquitin encoded by the RPS27a fusion gene is highly divergent, as previously found in a variety of protists. However, we also find that two other forms of divergent ubiquitin sequence, each previously thought to be extremely rare, have arisen recently during the divergence of Plasmodium subgenera. On two occasions, in two distinct lineages, the ubiquitin encoded by the RPL40 fusion gene has rapidly diverged. In addition, in one of these lineages, the polyubiquitin genes have undergone a single codon insertion, previously considered a unique feature of Rhizaria. There has been disagreement whether the multiple ubiquitin coding repeats within a genome exhibit concerted evolution or undergo a birth-and-death process; the Plasmodium ubiquitin genes show clear signs of concerted evolution, including the spread of this codon insertion to multiple repeats within the polyubiquitin gene.

Evolution of primate interferon-induced transmembrane proteins (IFITMs): a story of gain and loss with a differentiation into a canonical cluster and IFITM retrogenes

Interferon-inducible transmembrane proteins (IFITMs) are a family of transmembrane proteins. The subgroup of immunity-related (IR-)IFITMs is involved in adaptive and innate immune responses, being especially active against viruses. Here, we suggest that IFITMs should be classified as (1) a canonical IFITM gene cluster, which is located on the same chromosome, and (2) IFITM retrogenes, with a random and unique location at different positions within the genome. Phylogenetic analyses of the canonical cluster revealed the existence of three novel groups of primate IFITMs (pIFITM) in the IR-IFITM clade: the prosimian pIFITMs(pro), the new world monkey pIFITMs(nwm) and the old world monkey pIFITMs(owm). Therefore, we propose a new nomenclature: IR-pIFITM1, IR-pIFITM2, IR-pIFITM3, IR-pIFITMnwm, IR-pIFITMowm, and IR-pIFITMpro. We observed divergent evolution for pIFITM5 and pIFITM10, and evidence for concerted evolution and a mechanism of birth-and-death evolution model for the IR-pIFITMs. In contrast, the IFITMs scattered throughout the genomes possessed features of retrogenes retrotransposed by class 1 transposable elements. The origin of the IFITM retrogenes correspond to more recent events. We hypothesize that the transcript of a canonical IFITM3 has been constantly retrotransposed using class 1 transposable elements resulting in the IFITM retro(pseudo)genes. The unique pattern of each species has most likely been caused by constant pseudogenization and loss of the retro(pseudo)genes. This suggests a third mechanism of evolution for the IR-IFITMs in primates, similar to the birth-and-death model of evolution, but via a transposable element mechanism, which resulted in retro(pseudo)genes.

Improved haplotype resolution of highly duplicated MHC genes in a long-read genome assembly using MiSeq amplicons

Long-read sequencing offers a great improvement in the assembly of complex genomic regions, such as the major histocompatibility complex (MHC) region, which can contain both tandemly duplicated MHC genes (paralogs) and high repeat content. The MHC genes have expanded in passerine birds, resulting in numerous MHC paralogs, with relatively high sequence similarity, making the assembly of the MHC region challenging even with long-read sequencing. In addition, MHC genes show rather high sequence divergence between alleles, making diploid-aware assemblers incorrectly classify haplotypes from the same locus as sequences originating from different genomic regions. Consequently, the number of MHC paralogs can easily be over- or underestimated in long-read assemblies. We therefore set out to verify the MHC diversity in an original and a haplotype-purged long-read assembly of one great reed warbler Acrocephalus arundinaceus individual (the focal individual) by using Illumina MiSeq amplicon sequencing. Single exons, representing MHC class I (MHC-I) and class IIB (MHC-IIB) alleles, were sequenced in the focal individual and mapped to the annotated MHC alleles in the original long-read genome assembly. Eighty-four percent of the annotated MHC-I alleles in the original long-read genome assembly were detected using 55% of the amplicon alleles and likewise, 78% of the annotated MHC-IIB alleles were detected using 61% of the amplicon alleles, indicating an incomplete annotation of MHC genes. In the haploid genome assembly, each MHC-IIB gene should be represented by one allele. The parental origin of the MHC-IIB amplicon alleles in the focal individual was determined by sequencing MHC-IIB in its parents. Two of five larger scaffolds, containing 6-19 MHC-IIB paralogs, had a maternal and paternal origin, respectively, as well as a high nucleotide similarity, which suggests that these scaffolds had been incorrectly assigned as belonging to different loci in the genome rather than as alternate haplotypes of the same locus. Therefore, the number of MHC-IIB paralogs was overestimated in the haploid genome assembly. Based on our findings we propose amplicon sequencing as a suitable complement to long-read sequencing for independent validation of the number of paralogs in general and for haplotype inference in multigene families in particular.

CropGF: a comprehensive visual platform for crop gene family mining and analysis

A gene family refers to a group of genes that share a common ancestry and encode proteins or RNA molecules with similar functions or structural features. Gene families play a crucial role in determining the traits of plants and can be utilized to develop new crop varieties. Therefore, a comprehensive database of gene family is significant for gaining deep insight into crops. To address this need, we have developed CropGF (https://bis.zju.edu.cn/cropgf), a comprehensive visual platform that encompasses six important crops (rice, wheat, maize, barley, sorghum and foxtail millet) and one model plant (Arabidopsis), as well as genomics, transcriptomics and proteomics data for gene family mining and analysis, covering a total of 314 611 genes and 4399 types of domains. CropGF provides a versatile search system that allows for the identification of gene families and their members in a single crop or multiple crops. Users can customize their search based on gene family domains and/or homology using keywords or BLAST. To enhance usability, we have collected the corresponding ID information from various public databases for both genes and domains. Furthermore, CropGF comprises numerous downstream analysis modules, such as ka/ks analysis, phylogenetic tree construction, subcellular localization analysis and more. These visually-displayed modules provide intuitive insights into gene expression patterns, gene family expansion and functional relationships across different molecular levels and different species. We believe that CropGF will be a valuable resource for deep mining and analysis in future studies of crop gene families. Database URL https://bis.zju.edu.cn/cropgf.

Pervasive tandem duplications and convergent evolution shape coral genomes

BACKGROUND

Over the last decade, several coral genomes have been sequenced allowing a better understanding of these symbiotic organisms threatened by climate change. Scleractinian corals are reef builders and are central to coral reef ecosystems, providing habitat to a great diversity of species.

RESULTS

In the frame of the Tara Pacific expedition, we assemble two coral genomes, Porites lobata and Pocillopora cf. effusa, with vastly improved contiguity that allows us to study the functional organization of these genomes. We annotate their gene catalog and report a relatively higher gene number than that found in other public coral genome sequences, 43,000 and 32,000 genes, respectively. This finding is explained by a high number of tandemly duplicated genes, accounting for almost a third of the predicted genes. We show that these duplicated genes originate from multiple and distinct duplication events throughout the coral lineage. They contribute to the amplification of gene families, mostly related to the immune system and disease resistance, which we suggest to be functionally linked to coral host resilience.

CONCLUSIONS

At large, we show the importance of duplicated genes to inform the biology of reef-building corals and provide novel avenues to understand and screen for differences in stress resilience.

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.