Molecular evolution of Hokkaido virus, a genotype of Orthohantavirus puumalaense, among Myodes rodents.

Viruses in the genus Orthohantavirus within the family Hantaviridae cause human hantavirus infections and represent a threat to public health. Hokkaido virus (HOKV), a genotype of Orthohantavirus puumalaense (Puumala virus; PUUV), was first identified in Tobetsu, Hokkaido, Japan. Although it is genetically related to the prototype of PUUV, the evolutionary pathway of HOKV is unclear. We conducted a field survey in a forest in Tobetsu in 2022 and captured 44 rodents. Complete coding genome sequences of HOKVs were obtained from five viral-RNA-positive rodents (four Myodes rufocanus bedfordiae and one Apodemus speciosus). Phylogenetic analysis revealed a close relationship between the phylogenies and geographical origins of M. rufocanus-related orthohantaviruses. Comparison of the phylogenetic trees of the S segments of orthohantaviruses and the cytochrome b genes of Myodes species suggested that Myodes-related orthohantaviruses evolved in Myodes rodent species as a result of genetic isolation and host switching.

Evolution of threat response-related polymorphisms at the SLC6A4 locus in callitrichid primates.

Variation in an upstream repetitive region at the SLC6A4 locus, which encodes the serotonin transporter, is associated with anxiety-related behaviour in a few primate species, including humans and rhesus macaques, and has been suggested to be related to ecological adaptability among macaques. In this study, we investigate evolution of SLC6A4 polymorphisms associated with anxiety-related behaviour in common marmosets (Callithrix jacchus). Assaying variation in the SLC6A4 repeat region across 14 species in eight genera of callitrichid primates (marmosets and tamarins), we find large interspecific variation in the number of repeats present (24-43). The black tufted-ear marmoset (C. penicillata) has sequence polymorphisms similar to those found in the common marmoset, which is its sister species, and no other species has intraspecific variation at these sites. We conclude that, similar to humans and macaques, the functional polymorphism at SLC6A4 in common marmosets has a recent evolutionary origin, and that the anxiety-related allele is evolutionarily derived. Common/black tufted-ear marmosets and rhesus/bonnet macaques share high ecological adaptability and behavioural flexibility that we propose may be related to the maintenance of the polymorphism.

Divergent evolution of male-determining loci on proto-Y chromosomes of the housefly.

Houseflies provide a good experimental model to study the initial evolutionary stages of a primary sex-determining locus because they possess different recently evolved proto-Y chromosomes that contain male-determining loci (M) with the same male-determining gene, Mdmd. We investigate M-loci genomically and cytogenetically revealing distinct molecular architectures among M-loci. M on chromosome V (MV) has two intact Mdmd copies in a palindrome. M on chromosome III (MIII) has tandem duplications containing 88 Mdmd copies (only one intact) and various repeats, including repeats that are XY-prevalent. M on chromosome II (MII) and the Y (MY) share MIII-like architecture, but with fewer repeats. MY additionally shares MV-specific sequence arrangements. Based on these data and karyograms using two probes, one derives from MIII and one Mdmd-specific, we infer evolutionary histories of polymorphic M-loci, which have arisen from unique translocations of Mdmd, embedded in larger DNA fragments, and diverged independently into regions of varying complexity.

The loci of environmental adaptation in a model eukaryote.

While the underlying genetic changes have been uncovered in some cases of adaptive evolution, the lack of a systematic study prevents a general understanding of the genomic basis of adaptation. For example, it is unclear whether protein-coding or noncoding mutations are more important to adaptive evolution and whether adaptations to different environments are brought by genetic changes distributed in diverse genes and biological processes or concentrated in a core set. We here perform laboratory evolution of 3360 Saccharomyces cerevisiae populations in 252 environments of varying levels of stress. We find the yeast adaptations to be primarily fueled by large-effect coding mutations overrepresented in a relatively small gene set, despite prevalent antagonistic pleiotropy across environments. Populations generally adapt faster in more stressful environments, partly because of greater benefits of the same mutations in more stressful environments. These and other findings from this model eukaryote help unravel the genomic principles of environmental adaptation.

Evolutionary conservation of VSX2 super-enhancer modules in retinal development.

Super-enhancers (SEs) are expansive regions of genomic DNA that regulate the expression of genes involved in cell identity and cell fate. We recently identified developmental stage- and cell type-specific modules within the murine Vsx2 SE. Here, we show that the human VSX2 SE modules have similar developmental stage- and cell type-specific activity in reporter gene assays. By inserting the human sequence of one VSX2 SE module into a mouse with microphthalmia, eye size was rescued. To understand the function of these SE modules during human retinal development, we deleted individual modules in human embryonic stem cells and generated retinal organoids. Deleting one module results in small organoids, recapitulating the small-eyed phenotype of mice with microphthalmia, while deletion of the other module led to disruptions in bipolar neuron development. This prototypical SE serves as a model for understanding developmental stage- and cell type-specific effects of neurogenic transcription factors with complex expression patterns. Moreover, by elucidating the gene regulatory mechanisms, we can begin to examine how dysregulation of these mechanisms contributes to phenotypic diversity and disease.

Identification of homologs of the Chlamydia trachomatis effector CteG reveals a family of Chlamydiaceae type III secreted proteins that can be delivered into host cells.

Chlamydiae are a large group of obligate endosymbionts of eukaryotes that includes the Chlamydiaceae family, comprising several animal pathogens. Among Chlamydiaceae, Chlamydia trachomatis causes widespread ocular and urogenital infections in humans. Like many bacterial pathogens, all Chlamydiae manipulate host cells by injecting them with type III secretion effector proteins. We previously characterized the C. trachomatis effector CteG, which localizes at the host cell Golgi and plasma membrane during distinct phases of the chlamydial infectious cycle. Here, we show that CteG is a Chlamydiaceae-specific effector with over 60 homologs phylogenetically categorized into two distinct clades (CteG I and CteG II) and exhibiting several inparalogs and outparalogs. Notably, cteG I homologs are syntenic to C. trachomatis cteG, whereas cteG II homologs are syntenic among themselves but not with C. trachomatis cteG. This indicates a complex evolution of cteG homologs, which is unique among C. trachomatis effectors, marked by numerous events of gene duplication and loss. Despite relatively modest sequence conservation, nearly all tested CteG I and CteG II proteins were identified as type III secretion substrates using Yersinia as a heterologous bacterial host. Moreover, most of the type III secreted CteG I and CteG II homologs were delivered by C. trachomatis into host cells, where they localized at the Golgi region and cell periphery. Overall, this provided insights into the evolution of bacterial effectors and revealed a Chlamydiaceae family of type III secreted proteins that underwent substantial divergence during evolution while conserving the capacity to localize at specific host cell compartments.

Enrichment of rare codons at 5' ends of genes is a spandrel caused by evolutionary sequence turnover and does not improve translation.

Previously, Tuller et al. found that the first 30-50 codons of the genes of yeast and other eukaryotes are slightly enriched for rare codons. They argued that this slowed translation, and was adaptive because it queued ribosomes to prevent collisions. Today, the translational speeds of different codons are known, and indeed rare codons are translated slowly. We re-examined this 5' slow translation 'ramp.' We confirm that 5' regions are slightly enriched for rare codons; in addition, they are depleted for downstream Start codons (which are fast), with both effects contributing to slow 5' translation. However, we also find that the 5' (and 3') ends of yeast genes are poorly conserved in evolution, suggesting that they are unstable and turnover relatively rapidly. When a new 5' end forms de novo, it is likely to include codons that would otherwise be rare. Because evolution has had a relatively short time to select against these codons, 5' ends are typically slightly enriched for rare, slow codons. Opposite to the expectation of Tuller et al., we show by direct experiment that genes with slowly translated codons at the 5' end are expressed relatively poorly, and that substituting faster synonymous codons improves expression. Direct experiment shows that slow codons do not prevent downstream ribosome collisions. Further informatic studies suggest that for natural genes, slow 5' ends are correlated with poor gene expression, opposite to the expectation of Tuller et al. Thus, we conclude that slow 5' translation is a 'spandrel'--a non-adaptive consequence of something else, in this case, the turnover of 5' ends in evolution, and it does not improve translation.

The immune modules conserved across the tree of life: Towards a definition of ancestral immunity.

Immune defence mechanisms exist across the tree of life in such diversity that prokaryotic antiviral responses have historically been considered unrelated to eukaryotic immunity. Mechanisms of defence in divergent eukaryotes were similarly believed to be largely clade specific. However, recent data indicate that a subset of modules (domains and proteins) from prokaryote defence systems are conserved in eukaryotes and populate many stages of innate immune pathways. In this Essay, we propose the notion of ancestral immunity, which corresponds to the set of immune modules conserved between prokaryotes and eukaryotes. After offering a typology of ancestral immunity, we speculate on the selective pressures that could have led to the differential conservation of specific immune modules across domains of life. The exploration of ancestral immunity is in its infancy and appears full of promises to illuminate immune evolution, and also to identify and decipher immune mechanisms of economic, ecological, and therapeutic importance.

Comparative mitochondrial genome analysis provides new insights into the classification of Modiolinae.

BACKGROUND: Mitochondrial genomes have become a powerful tool for studying molecular genetics and phylogeny of mollusks. Currently, the position of Modiolinae within Mytilidae and the taxonomic and phylogenetic relationships within Modiolinae were still controversial. This study focuses on the complete mitochondrial genomes of two species: Modiolus modulaides (Röding, 1798) and Modiolus auriculatus Krauss, 1848, which have not been sequenced before. METHODS AND RESULTS: We assembled and characterized the mitochondrial genomes of M. modulaides and M. auriculatus and then analyzed the phylogenetic relationships. The mitochondrial genomes of M. modulaides and M. auriculatus were 15,422 bp and 16,027 bp, respectively. Both of them were composed of 36 functional genes, including 12 protein-coding genes, 22 transfer RNAs, and 2 ribosomal RNAs. All protein-coding genes showed A + T bias, positive GC skews, and negative AT skews in nucleotide composition. Phylogenetic analysis based on the mitochondrial genomes showed that Modiolinae and Bathymodiolinae clustered together to form a sister relationship. Seven Modiolinae species were divided into two clades: L1 (M. modulaides, M. auriculatus and Modiolus philippinarum Hanley, 1843) and L2 [Modiolus modiolus (Linnaeus, 1758), Modiolus kurilensis Bernard, 1983, Modiolus nipponicus (Oyama, 1950), and Modiolus comptus (Sowerby III, 1915)]. The divergence time of the two clades was approximately 105.75 Ma. Furthermore, the transfer RNA gene rearrangement, longer genetic distance, and greater genetic differentiation were confirmed between the L1 and L2 clades, as well as differences in the external characteristics of the shells of the two clades. CONCLUSIONS: Based on the molecular data, it was speculated that species from the L1 clade might belong to other genera or new genera. This study provides molecular information for further taxonomic and phylogenetic studies of Mytilidae.

70 Years - from DNA Double Helix via Approaching Systems Genomics to a Generalized Unified Evolution Theory.

70 Years - from DNA Double Helix via Approaching Systems Genomics to a Generalized Unified Evolution Theory.

A broad survey of choanoflagellates revises the evolutionary history of the Shaker family of voltage-gated K+ channels in animals.

The Shaker family of voltage-gated K+ channels has been thought of as an animal-specific ion channel family that diversified in concert with nervous systems. It comprises four functionally independent gene subfamilies (Kv1-4) that encode diverse neuronal K+ currents. Comparison of animal genomes predicts that only the Kv1 subfamily was present in the animal common ancestor. Here, we show that some choanoflagellates, the closest protozoan sister lineage to animals, also have Shaker family K+ channels. Choanoflagellate Shaker family channels are surprisingly most closely related to the animal Kv2-4 subfamilies which were believed to have evolved only after the divergence of ctenophores and sponges from cnidarians and bilaterians. Structural modeling predicts that the choanoflagellate channels share a T1 Zn2+ binding site with Kv2-4 channels that is absent in Kv1 channels. We functionally expressed three Shakers from Salpingoeca helianthica (SheliKvT1.1-3) in Xenopus oocytes. SheliKvT1.1-3 function only in two heteromultimeric combinations (SheliKvT1.1/1.2 and SheliKvT1.1/1.3) and encode fast N-type inactivating K+ channels with distinct voltage dependence that are most similar to the widespread animal Kv1-encoded A-type Shakers. Structural modeling of the T1 assembly domain supports a preference for heteromeric assembly in a 2:2 stoichiometry. These results push the origin of the Shaker family back into a common ancestor of metazoans and choanoflagellates. They also suggest that the animal common ancestor had at least two distinct molecular lineages of Shaker channels, a Kv1 subfamily lineage predicted from comparison of animal genomes and a Kv2-4 lineage predicted from comparison of animals and choanoflagellates.

Chromosome-Scale Genome Assembly for Soft-Stem Bulrush (Schoenoplectus tabernaemontani) Confirms a Clade-Specific Whole-Genome Duplication in Cyperaceae.

Schoenoplectus tabernaemontani (C. C. Gmelin) Palla is a typical macrophyte in diverse wetland ecosystems. This species holds great potential in decontamination applications and carbon sequestration. Previous studies have shown that this species may have experienced recent polyploidization. This would make S. tabernaemontani a unique model to study the processes and consequences of whole-genome duplications in the context of the well-documented holocentric chromosomes and dysploidy events in Cyperaceae. However, the inference was not completely solid because it lacked homology information that is essential to ascertain polyploidy. We present here the first chromosome-level genome assembly for S. tabernaemontani. By combining Oxford Nanopore Technologies (ONT) long reads and Illumina short reads, plus chromatin conformation via the Hi-C method, we assembled a genome spanning 507.96 Mb, with 99.43% of Hi-C data accurately mapped to the assembly. The assembly contig N50 value was 3.62 Mb. The overall BUSCO score was 94.40%. About 68.94% of the genome was comprised of repetitive elements. A total of 36,994 protein-coding genes were predicted and annotated. Long terminal repeat retrotransposons accounted for ∼26.99% of the genome, surpassing the content observed in most sequenced Cyperid genomes. Our well-supported haploid assembly comprised 21 pseudochromosomes, each harboring putative holocentric centromeres. Our findings corroborated a karyotype of 2n = 2X = 42. We also confirmed a recent whole-genome duplication occurring after the divergence between Schoenoplecteae and Bolboschoeneae. Our genome assembly expands the scope of sequenced genomes within the Cyperaceae family, encompassing the fifth genus. It also provides research resources on Cyperid evolution and wetland conservation.

Pseudogenization of NK3 homeobox 2 (Nkx3.2) in monotremes provides insight into unique gastric anatomy and physiology.

The enzymatic breakdown and regulation of food passage through the vertebrate antral stomach and pyloric sphincter (antropyloric region) is a trait conserved over 450 million years. Development of the structures involved is underpinned by a highly conserved signalling pathway involving the hedgehog, bone morphogenetic protein and Wingless/Int-1 (Wnt) protein families. Monotremes are one of the few vertebrate lineages where acid-based digestion has been lost, and this is consistent with the lack of genes for hydrochloric acid secretion and gastric enzymes in the genomes of the platypus (Ornithorhynchus anatinus) and short-beaked echidna (Tachyglossus aculeatus) . Furthermore, these species feature unique gastric phenotypes, both with truncated and aglandular antral stomachs and the platypus with no pylorus. Here, we explore the genetic underpinning of monotreme gastric phenotypes, investigating genes important in antropyloric development using the newest monotreme genomes (mOrnAna1.pri.v4 and mTacAcu1) together with RNA-seq data. We found that the pathway constituents are generally conserved, but surprisingly, NK3 homeobox 2 (Nkx3.2) was pseudogenized in both platypus and echidna. We speculate that the unique sequence evolution of Grem1 and Bmp4 sequences in the echidna lineage may correlate with their pyloric-like restriction and that the convergent loss of gastric acid and stomach size genotypes and phenotypes in teleost and monotreme lineages may be a result of eco-evolutionary dynamics. These findings reflect the effects of gene loss on phenotypic evolution and further elucidate the genetic control of monotreme stomach anatomy and physiology.

The Structural Biology of Catalase Evolution.

Catalases are essential enzymes for removal of hydrogen peroxide, enabling aerobic and anaerobic metabolism in an oxygenated atmosphere. Monofunctional heme catalases, catalase-peroxidases, and manganese catalases, evolved independently more than two billion years ago, constituting a classic example of convergent evolution. Herein, the diversity of catalase sequences is analyzed through sequence similarity networks, providing the context for sequence distribution of major catalase families, and showing that many divergent catalase families remain to be experimentally studied.

Time-calibrated phylogenetic and chromosomal mobilome analyses of Staphylococcus aureus CC398 reveal geographical and host-related evolution.

An international collection of Staphylococcus aureus of clonal complex (CC) 398 from diverse hosts spanning all continents and a 30 year-period is studied based on whole-genome sequencing (WGS) data. The collection consists of publicly available genomic data from 2994 strains and 134 recently sequenced Swiss methicillin-resistant S. aureus (MRSA) CC398 strains. A time-calibrated phylogeny reveals the presence of distinct phylogroups present in Asia, North and South America and Europe. European MRSA diverged from methicillin-susceptible S. aureus (MSSA) at the beginning of the 1950s. Two major European phylogroups (EP4 and EP5), which diverged approximately 1974, are the main drivers of MRSA CC398 spread in Europe. Within EP5, an emergent MRSA lineage spreading among the European horse population (EP5-Leq) diverged approximately 1996 from the pig lineage (EP5-Lpg), and also contains human-related strains. EP5-Leq is characterized by staphylococcal cassette chromosome mec (SCCmec) IVa and spa type t011 (CC398-IVa-t011), and EP5-Lpg by CC398-SCCmecVc-t011. The lineage-specific antibiotic resistance and virulence gene patterns are mostly mediated by the acquisition of mobile genetic elements like SCCmec, S. aureus Genomic Islands (SaGIs), prophages and transposons. Different combinations of virulence factors are present on S. aureus pathogenicity islands (SaPIs), and novel antimicrobial resistance gene containing elements are associated with certain lineages expanding in Europe. This WGS-based analysis reveals the actual evolutionary trajectory and epidemiological trend of the international MRSA CC398 population considering host, temporal, geographical and molecular factors. It provides a baseline for global WGS-based One-Health studies of adaptive evolution of MRSA CC398 as well as for local outbreak investigations.

Revisiting the evolution of bow-tie architecture in signaling networks.

Bow-tie architecture is a layered network structure that has a narrow middle layer with multiple inputs and outputs. Such structures are widely seen in the molecular networks in cells, suggesting that a universal evolutionary mechanism underlies the emergence of bow-tie architecture. The previous theoretical studies have implemented evolutionary simulations of the feedforward network to satisfy a given input-output goal and proposed that the bow-tie architecture emerges when the ideal input-output relation is given as a rank-deficient matrix with mutations in network link intensities in a multiplicative manner. Here, we report that the bow-tie network inevitably appears when the link intensities representing molecular interactions are small at the initial condition of the evolutionary simulation, regardless of the rank of the goal matrix. Our dynamical system analysis clarifies the mechanisms underlying the emergence of the bow-tie structure. Further, we demonstrate that the increase in the input-output matrix reduces the width of the middle layer, resulting in the emergence of bow-tie architecture, even when evolution starts from large link intensities. Our data suggest that bow-tie architecture emerges as a side effect of evolution rather than as a result of evolutionary adaptation.

Rewinding the Ratchet: Rare Recombination Locally Rescues Neo-W Degeneration and Generates Plateaus of Sex-Chromosome Divergence.

Natural selection is less efficient in the absence of recombination. As a result, nonrecombining sequences, such as sex chromosomes, tend to degenerate over time. Although the outcomes of recombination arrest are typically observed after many millions of generations, recent neo-sex chromosomes can give insight into the early stages of this process. Here, we investigate the evolution of neo-sex chromosomes in the Spanish marbled white butterfly, Melanargia ines, where a Z-autosome fusion has turned the homologous autosome into a nonrecombining neo-W chromosome. We show that these neo-sex chromosomes are likely limited to the Iberian population of M. ines, and that they arose around the time when this population split from North-African populations, around 1.5 million years ago. Recombination arrest of the neo-W chromosome has led to an excess of premature stop-codons and frame-shift mutations, and reduced gene expression compared to the neo-Z chromosome. Surprisingly, we identified two regions of ∼1 Mb at one end of the neo-W that are both less diverged from the neo-Z and less degraded than the rest of the chromosome, suggesting a history of rare but repeated genetic exchange between the two neo-sex chromosomes. These plateaus of neo-sex chromosome divergence suggest that neo-W degradation can be locally reversed by rare recombination between neo-W and neo-Z chromosomes.

Irreducible Complexity of Hox Gene: Path to the Canonical Function of the Hox Cluster.

The evolution of major taxa is often associated with the emergence of new gene families. In all multicellular animals except sponges and comb jellies, the genomes contain Hox genes, which are crucial regulators of development. The canonical function of Hox genes involves colinear patterning of body parts in bilateral animals. This general function is implemented through complex, precisely coordinated mechanisms, not all of which are evolutionarily conserved and fully understood. We suggest that the emergence of this regulatory complexity was preceded by a stage of cooperation between more ancient morphogenetic programs or their individual elements. Footprints of these programs may be present in modern animals to execute non-canonical Hox functions. Non-canonical functions of Hox genes are involved in maintaining terminal nerve cell specificity, autophagy, oogenesis, pre-gastrulation embryogenesis, vertical signaling, and a number of general biological processes. These functions are realized by the basic properties of homeodomain protein and could have triggered the evolution of ParaHoxozoa and Nephrozoa subsequently. Some of these non-canonical Hox functions are discussed in our review.

Recurrent evolution of adhesive defence systems in amphibians by parallel shifts in gene expression.

Natural selection can drive organisms to strikingly similar adaptive solutions, but the underlying molecular mechanisms often remain unknown. Several amphibians have independently evolved highly adhesive skin secretions (glues) that support a highly effective antipredator defence mechanism. Here we demonstrate that the glue of the Madagascan tomato frog, Dyscophus guineti, relies on two interacting proteins: a highly derived member of a widespread glycoprotein family and a galectin. Identification of homologous proteins in other amphibians reveals that these proteins attained a function in skin long before glues evolved. Yet, major elevations in their expression, besides structural changes in the glycoprotein (increasing its structural disorder and glycosylation), caused the independent rise of glues in at least two frog lineages. Besides providing a model for the chemical functioning of animal adhesive secretions, our findings highlight how recruiting ancient molecular templates may facilitate the recurrent evolution of functional innovations.

Conservation of cis-Regulatory Syntax Underlying Deuterostome Gastrulation.

Throughout embryonic development, the shaping of the functional and morphological characteristics of embryos is orchestrated by an intricate interaction between transcription factors and cis-regulatory elements. In this study, we conducted a comprehensive analysis of deuterostome cis-regulatory landscapes during gastrulation, focusing on four paradigmatic species: the echinoderm Strongylocentrotus purpuratus, the cephalochordate Branchiostoma lanceolatum, the urochordate Ciona intestinalis, and the vertebrate Danio rerio. Our approach involved comparative computational analysis of ATAC-seq datasets to explore the genome-wide blueprint of conserved transcription factor binding motifs underlying gastrulation. We identified a core set of conserved DNA binding motifs associated with 62 known transcription factors, indicating the remarkable conservation of the gastrulation regulatory landscape across deuterostomes. Our findings offer valuable insights into the evolutionary molecular dynamics of embryonic development, shedding light on conserved regulatory subprograms and providing a comprehensive perspective on the conservation and divergence of gene regulation underlying the gastrulation process.