Search

Argonaute (AGO) proteins are central players in RNA interference in eukaryotes. They associate with small RNAs (sRNA) and lead to transcriptional or posttranscriptional silencing of targets, thereby regulating diverse biological processes. The molecular and biological functions of AGO proteins have been extensively characterized, particularly in a few angiosperm species, leading to the recognition that the AGO family has expanded to accommodate diverse sRNAs thereby performing diverse biological functions. However, understanding of the expansion of AGO proteins in plants is still limited, due to a dearth of knowledge of AGO proteins in green algal groups. Here, we identified more than 2900 AGO proteins from 244 plant species, including green algae, and performed a large-scale phylogenetic analysis. The phylogeny shows that the plant AGO family gave rise to four clades after the emergence of hydrobiontic algae and prior to the emergence of land plants. Subsequent parallel expansion in ferns and angiosperms resulted in eight main clades in angiosperms: AGO2, AGO7, AGO6, AGO4, AGO1, AGO10a, AGO10b and AGO5. On the basis of this phylogeny, we identified two novel AGO4 orthologs that Arabidopsis does not have, and redefined AGO10, which is composed of AGO10a and AGO10b. Finally, we propose a hypothetical evolutionary model of AGO proteins in plants. Our studies provide a deeper understanding of the phylogenetic relationships of AGO family members in the green lineage, which would help to further reveal their roles as RNAi effectors.

International audience; The deployment of plant varieties carrying resistance genes (R) exerts strong selection pressure on pathogen populations. Rapidly evolving avirulence genes (Avr) allow pathogens to escape R-mediated plant immunity through a variety of mechanisms, leading to virulence. The poplar rust fungus Melampsora larici-populina is a damaging pathogen of poplars in Europe. It underwent a major adaptive event in 1994, with the breakdown of the poplar RMlp7 resistance gene. Population genomics studies identified a locus in the genome of M. larici-populina that probably corresponds to the candidate avirulence gene AvrMlp7. Here, to further characterize this effector, we used a population genetics approach on a comprehensive set of 281 individuals recovered throughout a 28-year period encompassing the resistance breakdown event. Using two dedicated molecular tools, genotyping at the candidate locus highlighted two different alterations of a predominant allele found mainly before the resistance breakdown: a nonsynonymous mutation and a complete deletion of this locus. This results in six diploid genotypes: three genotypes related to the avirulent phenotype and three related to the virulent phenotype. The temporal survey of the candidate locus revealed that both alterations were found in association during the resistance breakdown event. They pre-existed before the breakdown in a heterozygous state with the predominant allele cited above. Altogether, these results suggest that the association of both alterations at the candidate locus AvrMlp7 drove the poplar rust adaptation to RMlp7-mediated immunity. This study demonstrates for the first time a case of adaptation from standing genetic variation in rust fungi during a qualitative resistance breakdown.

Triterpenes are among the most diverse plant natural products, and their diversity is closely related to various triterpene skeletons catalyzed by different 2,3-oxidosqualene cyclases (OSCs). Celastrol, a friedelane-type triterpene with significant bioactivities, is specifically distributed in higher plants, such as Celastraceae species. Friedelin is an important precursor for the biosynthesis of celastrol, and it is synthesized through the cyclization of 2,3-oxidosqualene, with the highest number of rearrangements being catalyzed by friedelane-type triterpene cyclases. However, the molecular mechanisms underlying the catalysis of friedelin production by friedelane-type triterpene cyclases have not yet been fully elucidated. In this study, transcriptome data of four celastrol-producing plants from Celastraceae were used to identify a total of 21 putative OSCs. Through functional characterization, the friedelane-type triterpene cyclases were separately verified in the four plants. Analysis of the selection pressure showed that purifying selection acted on these OSCs, and the friedelane-type triterpene cyclases may undergo weaker selective restriction during evolution. Molecular docking and site-directed mutagenesis revealed that changes in some amino acids that are unique to friedelane-type triterpene cyclases may lead to variations in catalytic specificity or efficiency, thereby affecting the synthesis of friedelin. Our research explored the functional diversity of triterpene synthases from a multispecies perspective. It also provides some references for further research on the relative mechanisms of friedelin biosynthesis.

Measles is one of the most infectious diseases of humans. It is caused by measles virus (MeV) and can lead to serious illness, lifelong complications and even death. Whole-genome sequencing (WGS) is now available to study molecular epidemiology and identify MeV transmission pathways. In the present study, WGS of 23 MeV strains of genotype H1, collected in Mainland China between 2006 and 2018, were generated and compared to 31 WGSs from the public domain in order to analyze genomic characteristics, evolutionary rates and date of emergence of H1 genotype. The non-coding region between M and F protein genes (M/F NCR) was the most variable region throughout the genome. Whereas the nucleotide substitution rate of H1 WGS was around 0.75×10-3 substitution per site per year, the M/F NCR had an evolutionary rate 3 times higher, with 2.44×10-3 substitution per site per year. Phylogenetic analysis identified 3 distinct genetic groups. The Time of the Most Recent Common Ancestor (TMRCA) of H1 genotype was estimated at approximately 1988, while the first genetic group appeared around 1995 followed by 2 other genetic groups in 1999-2002. Bayesian skyline plot showed that the genetic diversity of H1 genotype remained stable even though the number of MeV cases decreased 50 times between 2014 (52,628) and 2020 (993). The current COVID-19 pandemic might have some effect on measles epidemic and further studies will be necessary to assess genetic diversity of H1 genotype in a post COVID area. This article is protected by copyright. All rights reserved.

Adaptive evolution of vision-related genes has been frequently observed in the process of invasion of new environments in a wide range of animal taxa. The typical example is that of the molecular evolution of rhodopsin associated with habitat changes in aquatic animals. However, few studies have investigated rhodopsin evolution during adaptive radiation across various habitats. In the present study, we examined the link between molecular evolutionary patterns in the rhodopsin gene and macroscopic habitat changes in Gymnogobius species (Gobiidae), which have adaptively radiated to diverse aquatic habitats including the sea, brackish waters, rivers and lakes. Analysis of amino acid substitutions in rhodopsin in the phylogenetic framework revealed convergent substitutions in 4-5 amino acids in three groups (four species), including two spectral tuning amino acid sites known to change rhodopsin's absorption wavelength. Positive selection was detected in the basal branches of each of these three groups, suggesting adaptive molecular convergence of rhodopsin. However, no significant correlation was observed between amino acid substitutions and the species' habitat changes, suggesting molecular adaptation to some unidentified micro-ecological environments. Taken together, these results emphasize the importance of considering not only macroscopic habitats but also micro-ecological environments when elucidating the driving forces of adaptive evolution of the visual system.

Plant hormone cytokinins are important regulators of plant development, response to environmental stresses and interplay with other plant hormones. Cytokinin dehydrogenases (CKXs) are proteins responsible for the irreversible break-down of cytokinins to the adenine and aldehyde. Even though plant CKXs have been extensively studied, homologous proteins from other taxa remain mainly uncharacterised. Here we present our study on the molecular evolution and divergence of the CKX from bacteria, fungi, amoeba and viridiplantae. Although CKXs are present in eukaryotes and prokaryotes, they are missing in algae and metazoan taxa. The prevalent domain architecture consists of the FAD-binding and cytokinin binding domains, whereas some bacteria appear to have only cytokinin binding domain proteins. The CKXs play important role in the various aspects of plant life including control of plant development, response to biotic and abiotic stress, influence nutrition. Results of our study suggested that CKX originates from the FAD-linked C-terminal oxidase and has a defence-oriented function. The obtained results significantly extend the current understanding of the cytokinin dehydrogenases structure–function from the relationship to homologues from other taxa and provide a starting point baseline for their future functional characterization.

Partial mtDNA cytochrome oxidase subunit I (COI) fragments and near entire stretch of 12S rDNA (12S) and control region (Dloop) of the Japanese spiny lobster (Panulirus japonicus) (n = 3) were amplified by PCR and used for direct nucleotide sequencing and for clone library-based nucleotide sequence analysis. Nucleotide sequences of a total of 75 clones in COI, 77 in 12S and 92 in Dloop were determined. Haplotypes of the clones matched with those obtained by direct sequencing were determined to be genuine mtDNA sequence of the individual. Phylogenetic analysis revealed several distinct groups of haplotypes in all three regions. Genuine mtDNA sequences were observed to form a group with their closely related variables, and most of these variables may be due to amplification error but a few to be heteroplasmy. Haplotypes determined as nuclear mitochondrial pseudogenes (NUMTs) formed distinct groups. Nucleotide sequence divergence (K2P distance) between genuine haplotypes and NUMTs were substantial (7.169–23.880% for COI, 1.336–23.434% for 12S, and 7.897–71.862% for Dloop). These values were comparable to or smaller than those between species of the genus Panulirus, indicating that integration of mtDNA into the nuclear genome is a continuous and dynamic process throughout pre- and post-speciation events. Double peaks in electropherograms obtained by direct nucleotide sequencing were attributed to common nucleotides shared by multiple NUMTs. Information on the heteroplasmy and NUMTs would be very important for addressing their impact on direct nucleotide sequencing and for quality control of nucleotide sequences obtained.

Background The cereal family Poaceae is one of the largest and most diverse angiosperm families. The central component of centromere specification and function is the centromere-specific histone H3 (CENH3). Some cereal species (maize, rice) have one copy of the gene encoding this protein, while some (wheat, barley, rye) have two. We applied a homology-based approach to sequenced cereal genomes, in order to finally trace the mutual evolution of the structure of the CENH3 genes and the nearby regions in various tribes. Results We have established that the syntenic group or the CENH3 locus with the CENH3 gene and the boundaries defined by the CDPK2 and bZIP genes first appeared around 50 Mya in a common ancestor of the subfamilies Bambusoideae, Oryzoideae and Pooideae. This locus came to Pooideae with one copy of CENH3 in the most ancient tribes Nardeae and Meliceae. The βCENH3 gene as a part of the locus appeared in the tribes Stipeae and Brachypodieae around 35–40 Mya. The duplication was accompanied by changes in the exon-intron structure. Purifying selection acts mostly on αCENH3s, while βCENH3s form more heterogeneous structures, in which clade-specific amino acid motifs are present. In barley species, the βCENH3 gene assumed an inverted orientation relative to αCENH3 and the CDPK2 gene was substituted with LHCB-l. As the evolution and domestication of plant species went on, the locus was growing in size due to an increasing distance between αCENH3 and βCENH3 because of a massive insertion of the main LTR-containing retrotransposon superfamilies, gypsy and copia, without any evolutionary preference on either of them. A comparison of the molecular structure of the locus in the A, B and D subgenomes of the hexaploid wheat T. aestivum showed that invasion by mobile elements and concomitant rearrangements took place in an independent way even in evolutionarily close species. Conclusions The CENH3 duplication in cereals was accompanied by changes in the exon-intron structure of the βCENH3 paralog. The observed general tendency towards the expansion of the CENH3 locus reveals an amazing diversity of ways in which different species implement the scenario described in this paper.

The F420 deazaflavin cofactor is an intriguing molecule as it structurally resembles the canonical flavin cofactor, although behaves as a nicotinamide cofactor due to its obligate hydride-transfer reactivity and similar low redox potential. Since its discovery, numerous enzymes relying on it have been described. The known deazaflavoproteins are taxonomically restricted to Archaea and Bacteria. The biochemistry of the deazaflavoenzymes is diverse and they exhibit great structural variability. In this study a thorough sequence and structural homology evolutionary analysis was performed in order to generate an overarching classification of the F420 -dependent oxidoreductases. Five different deazaflavoenzyme Classes (I-V) are described according to their structural folds as follows: Class I encompassing the TIM-barrel F420 -dependent enzymes; Class II including the Rossmann fold F420 -dependent enzymes; Class III comprising the β-roll F420 -dependent enzymes; Class IV which exclusively gathers the SH3 barrel F420 -dependent enzymes and Class V including the 3 layer ββα sandwich F420 -dependent enzymes. This classification provides a framework for the identification and biochemical characterization of novel deazaflavoenzymes. This article is protected by copyright. All rights reserved.

DNA polymerase template switching between short, non-identical inverted repeats (IRs) is a genetic mechanism that leads to the homogenization of IR arms and to IR spacer inversion, which cause multinucleotide mutations (MNMs). It is unknown if and how template switching affects gene evolution. In this study, we performed a phylogenetic analysis to determine the effect of template switching between IR arms on coding DNA of Saccharomyces cerevisiae. To achieve this, perfect IRs that co-occurred with MNMs between a strain and its parental node were identified in S. cerevisiae strains. We determined that template switching introduced MNMs into 39 protein-coding genes through S. cerevisiae evolution, resulting in both arm homogenization and inversion of the IR spacer. These events in turn resulted in nonsynonymous substitutions and up to five neighboring amino acid replacements in a single gene. The study demonstrates that template switching is a powerful generator of multiple substitutions within codons. Additionally, some template switching events occurred more than once during S. cerevisiae evolution. Our findings suggest that template switching constitutes a general mutagenic mechanism that results in both nonsynonymous substitutions and parallel evolution, which are traditionally considered as evidence for positive selection, without the need for adaptive explanations.

The leucine-rich repeat receptor-like protein kinase (LRR-RLK) gene family is the largest family of the receptor-like protein kinases (RLKs) superfamily in higher plants, which is involved in regulating the plant growth and development, stress responses, signal transduction and so on. However, no comprehensive analyses of LRR-RLKs have been reported in sugarcane. Here, we performed a comprehensive analysis of the LRR-RLK gene family in sugarcane ancestor species Saccharum spontaneum. A total of 437 LRR-RLK genes were identified and categorized into 14 groups based on a maximum likelihood phylogenetic tree. The chromosome location showed an uneven distribution on all 32 chromosomes in sugarcane. Subsequently, the exon–intron organization structure and conserved motif arrangement were relatively conserved among the same groups or subgroups and between Arabidopsis and S. spontaneum genomes. Furthermore, the promoter sequences analyses showed that sugarcane LRR-RLK genes (SsLRR-RLKs) were strongly regulated by various environmental stimuli, phytohormonal factors and transcription factors (TFs). Eventually, the expression profiles of SsLRR-RLK genes at different stresses were analyzed based on RNA-seq data, suggesting their potential roles in the regulation of sugarcane responses to diverse abiotic and biotic stress. Overall, the findings provide insight into the potential functional roles and lay the foundation for further functional study.

Alpha satellite DNA is a major DNA component of primate centromeres. We previously reported that Azara's owl monkey has two types of alpha satellite DNA, OwlAlp1 and OwlAlp2. OwlAlp2 (344 bp) exhibits a sequence similarity throughout its entire length with alpha satellite DNA of closely related species. OwlAlp1 (185 bp) corresponds to the part of OwlAlp2. Based on the observation that the CENP-A protein binds to OwlAlp1, we proposed that OwlAlp1 is a relatively new repetitive DNA that replaced OwlAlp2 as the centromeric satellite DNA. However, a detailed picture of the evolutionary process of this centromere DNA replacement remains largely unknown. Here, we performed a phylogenetic analysis of OwlAlp1 and OwlAlp2 sequences, and also compared our results to alpha satellite DNA sequences of other primate species. We found that: (i) OwlAlp1 exhibits a higher similarity to OwlAlp2 than to alpha satellite DNA of other species, (ii) OwlAlp1 has a single origin, and (iii) sequence variation is lower in OwlAlp1 than in OwlAlp2. We conclude that OwlAlp1 underwent a recent and rapid expansion in the owl monkey lineage. This centromere DNA replacement could have been facilitated by the heterochromatin reorganization that is associated with the adaptation of owl monkeys to a nocturnal lifestyle.

The complete mitochondrial genome (mitogenome) of Saturnia japonica (Lepidoptera: Saturniidae) was sequenced and annotated. It is a circular molecule of 15, 376 bp, composed of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNA), and an adenine (A) + thymine (T)-rich region. All protein-coding genes (PCGs) are initiated by the ATN codon except for cytochrome c oxidase subunit 1 (cox1) gene that is seemingly initiated by the CGA codon. Except for cox2 and nad4, which were terminated by incomplete stop codon T or TA, the rest were terminated by canonical stop codon TAA. The A + T-rich region is high conservative, including ‘ATAGA’ motif followed by a 19 bp poly-T stretch, a microsatellite-like element (AT)9 and also a poly-A element, with a total length of 332 bp. The Asn codon was the most frequently used codon, followed by Ile, Leu2, Lys, Met, Phe, and Tyr, while Cys was the least frequently used codon. Phylogenetic relationships analysis based on the 13 PCGs by using maximum likelihood (ML) and neighbor Joining (NJ) revealed that S. japonica belongs to the Saturniidae family. In this study, the annotation and characteristics of the mitogenome of S. japonica were resolved for the first time, which laid a foundation for species classification and the molecular evolution of Lepidoptera: Saturniidae.

The ability to predict emerging variants of SARS-CoV-2 would be of enormous value, as it would enable proactive design of vaccines in advance of such emergence. We estimated diversity of each site on a multiple sequence alignment (MSA) of the Spike (S) proteins from close relatives of SARS-CoV-2 that infected bat and pangolin before the pandemic. Then we compared the locations of high diversity sites in this MSA and those of mutations found in multiple emerging lineages of human-infecting SARS-CoV-2. This comparison revealed a significant correspondence, which suggests that a limited number of sites in this protein are repeatedly substituted in different lineages of this group of viruses. It follows, therefore, that the sites of future emerging mutations in SARS-CoV-2 can be predicted by analyzing their relatives (outgroups) that have infected non-human hosts. We discuss a possible evolutionary basis for these substitutions and provide a list of frequently substituted sites that potentially include future emerging variants in SARS-CoV-2., In order to improve our ability to predict emerging variants of SARS-CoV-2, Katoh and Standley analysed the molecular evolution of the Spike protein. They found a significant correspondence in the location of mutations between recently emerging SARS-CoV-2 variants and their relatives that infected bat and pangolin before the pandemic, which could suggest that the sites of future mutations could be predicted by analyzing their relatives that have infected non-human hosts.

Which variables determine the constraints on gene sequence evolution is one of the most central questions in molecular evolution. In the fission yeast Schizosaccharomyces pombe, an important model organism, the variables influencing the rate of sequence evolution have yet to be determined. Previous studies in other single celled organisms have generally found gene expression levels to be most significant, with numerous other variables such as gene length and functional importance identified as having a smaller impact. Using publicly available data, we used partial least squares regression, principal components regression, and partial correlations to determine the variables most strongly associated with sequence evolution constraints. We identify centrality in the protein–protein interactions network, amino acid composition, and cellular location as the most important determinants of sequence conservation. However, each factor only explains a small amount of variance, and there are numerous variables having a significant or heterogeneous influence. Our models explain more than half of the variance in dN, raising the possibility that future refined models could quantify the role of stochastics in evolutionary rate variation.

Vertebrate photoreceptors are categorized into two broad classes, rods and cones, responsible for dim- and bright-light vision, respectively. While many molecular features that distinguish rods and cones are known, gene expression differences among cone subtypes remain poorly understood. Teleost fishes are renowned for the diversity of their photoreceptor systems. Here, we used single-cell RNA-seq to profile adult photoreceptors in zebrafish, a teleost. We found that in addition to the four canonical zebrafish cone types, there exist subpopulations of green and red cones (previously shown to be located in the ventral retina) that express red-shifted opsin paralogs (opn1mw4 or opn1lw1) as well as a unique combination of cone phototransduction genes. Furthermore, the expression of many paralogous phototransduction genes is partitioned among cone subtypes, analogous to the partitioning of the phototransduction paralogs between rods and cones seen across vertebrates. The partitioned cone-gene pairs arose via the teleost-specific whole-genome duplication or later clade-specific gene duplications. We also discovered that cone subtypes express distinct transcriptional regulators, including many factors not previously implicated in photoreceptor development or differentiation. Overall, our work suggests that partitioning of paralogous gene expression via the action of differentially expressed transcriptional regulators enables diversification of cone subtypes in teleosts.

Porcine circovirus 3 (PCV3) is a recently discovered virus that has been detected in the swine population worldwide. PCV3 infection has been associated with several signs, but its pathogenicity is currently uncertain. This review article aimed to analyse the PCV3 strains that circulate in different countries in North and South America. We demonstrated the main regions of polymorphisms in the capsid protein structure. Furthermore, we found that PCV3 has at least six different lineages circulating in the Americas. Additional studies are required to determine the role of PCV3 in different clinical syndromes and its epidemiology in swine herds in North and South American countries.

Mammals have evolved different tooth phenotypes that are hypothesized to be associated with feeding habits. However, the genetic basis for the linkage has not been well explored. In this study, we investigated 13 tooth-related genes, including seven enamel-related genes (AMELX, AMBN, ENAM, AMTN, ODAM, KLK4 and MMP20) and six dentin-related genes (DSPP, COL1A1, DMP1, IBSP, MEPE and SPP1), from 63 mammals to determine their evolutionary history. Our results showed that different evolutionary histories have evolved among divergent feeding habits in mammals. There was stronger positive selection for eight genes (ENAM, AMTN, ODAM, KLK4, DSPP, DMP1, COL1A1, MEPE) in herbivore lineages. In addition, AMELX, AMBN, ENAM, AMTN, MMP20 and COL1A1 underwent accelerated evolution in herbivores. While relatively strong positive selection was detected in IBSP, SPP1, and DSPP, accelerated evolution was only detected for MEPE and SPP1 genes among the carnivorous lineages. We found positive selection on AMBN and ENAM genes for omnivorous primates in the catarrhini clade. Interestingly, a significantly positive association between the evolutionary rate of ENAM, ODAM, KLK4, MMP20 and the average enamel thickness was found in primates. Additionally, we found molecular convergence in some amino acid sites of tooth-related genes among the lineages whose feeding habit are similar. The positive selection of related genes might promote the formation and bio-mineralization of tooth enamel and dentin, which would make the tooth structure stronger. Our results revealed that mammalian tooth-related genes have experienced variable evolutionary histories, which provide some new insights into the molecular basis of dietary adaptation in mammals.

Many genes responsible for melanin biosynthesis in fungi were physically linked together. The PKS gene clusters in most of the melanin-producing fungi were regulated by the Cmr1. It was found that a close rearrangement of the PKS gene clusters had evolved in most of the melanin-producing fungi and various functions of melanin in them were beneficial to their adaptation to the changing environments. The melanin-producing fungi had undergone at least five large-scale differentiations, making their PKS gene clusters be quickly evolved and the fungi be adapted to different harsh environments. The recent gene losses and expansion were remarkably frequent in the PKS gene clusters, leading to their rapid evolution and adaptation of their hosts to different environments. The PKS gene and the CMR1 gene in them were subject to a strong co-evolution, but the horizontal gene transfer events might have occurred in the genome-duplicated species, Aspergillus and Penicillium.

The forest tree family Aceraceae is widespread in the northern hemisphere and it has ecological and economic importance. However, the phylogenetic relationships and classifications within the family are still controversial due to transitional intraspecific morphological characteristics and introgression hybridization among species. In this study, we determined the evolutionary relationships and molecular evolution of Aceraceae based on plastid phylogenomics and two nuclear gene variations. Phylogenetic analysis based on the plastid genomes suggested that Aceraceae species can be divided into two larger sub-clades corresponding to the two genera Acer and Dipteronia. Conjoint analysis of the plastid and nuclear gene sequences supported the classification with two genera in the family. Molecular dating showed that the two genera diverged 60.2 million years ago, which is generally consistently with previously reported results. Divergence hotspots and positively selected genes identified in the plastid genomes could be useful genetic resources in Aceraceae.

RNA interference (RNAi) plays important roles in organism development through post-transcriptional regulation of specific target mRNAs. Target specificity is largely controlled by base-pair complementarity between micro-RNA (miRNA) regulatory elements and short regions of the target mRNA. The pattern of miRNA production in a cell interacts with the cell’s mRNA transcriptome to generate a specific network of post-transcriptional regulation that can play critical roles in cellular metabolism, differentiation, tissue/organ development and developmental timing. In plants, miRNA production is orchestrated in the nucleus by a suite of proteins that control transcription of the pri-miRNA gene, post-transcriptional processing and nuclear export of the mature miRNA. In the model plant, Arabidopsis thaliana, post-transcriptional processing of miRNAs is controlled by a pair of physically-interacting proteins, HYL1 and DCL1. However, the evolutionary history of the HYL1-DCL1 interaction is unknown, as is its structural basis. Here we use ancestral sequence reconstruction and functional characterization of ancestral HYL1 in vitro and in vivo to better understand the origin and evolution of the HYL1-DCL1 interaction and its impact on miRNA production and plant development. We found the ancestral plant HYL1 evolved high affinity for both double-stranded RNA (dsRNA) and its DCL1 partner very early in plant evolutionary history, before the divergence of mosses from seed plants (~500 Ma), and these high-affinity interactions remained largely conserved throughout plant evolutionary history. Structural modeling and molecular binding experiments suggest that the second of two double-stranded RNA-binding motifs (DSRMs) in HYL1 may interact tightly with the first of two C-terminal DCL1 DSRMs to mediate the HYL1-DCL1 physical interaction necessary for efficient miRNA production. Transgenic expression of the nearly 200 Ma-old ancestral flowering-plant HYL1 in A. thaliana was sufficient to rescue many key aspects of plant development disrupted by HYL1− knockout and restored near-native miRNA production, suggesting that the functional partnership of HYL1-DCL1 originated very early in and was strongly conserved throughout the evolutionary history of terrestrial plants. Overall, our results are consistent with a model in which miRNA-based gene regulation evolved as part of a conserved plant ‘developmental toolkit’; its role in generating developmental novelty is probably related to the relatively rapid evolution of miRNA genes.

Phosphoprotein phosphatase (PPP) enzymes are ubiquitous proteins involved in cellular signaling pathways and other functions. Here we have traced the origin of the PPP sequences of Eukaryotes and their radiation. Using a bacterial PPP Hidden Markov Model (HMM) we uncovered “BacterialPPP-Like” sequences in Archaea. A HMM derived from eukaryotic PPP enzymes revealed additional, unique sequences in Archaea and Bacteria that were more like the eukaryotic PPP enzymes then the bacterial PPPs. These sequences formed the basis of phylogenetic tree inference and sequence structural analysis allowing the history of these sequence types to be elucidated. Our phylogenetic tree data strongly suggest that eukaryotic PPPs ultimately arose from ancestors in the Asgard archaea. We have clarified the radiation of PPPs within Eukaryotes, substantially expanding the range of known organisms with PPP subtypes (Bsu1, PP7, PPEF/RdgC) previously thought to have a more restricted distribution. Surprisingly, sequences from the Methanosarcinaceae (Euryarchaeota) form a strongly supported sister group to eukaryotic PPPs in our phylogenetic analysis. This strongly suggests an intimate association between an Asgard ancestor and that of the Methanosarcinaceae. This is highly reminiscent of the syntrophic association recently demonstrated between the cultured Lokiarchaeal species Prometheoarchaeum and a methanogenic bacterial species.

We examined the evolutionary relationship of the ASPM (abnormal spindle-like microcephaly associated) and MCPH1 (microcephalin-1) genes with brain volume among humans and other primates. We obtained sequences of these genes from 14 simiiform species including hominins. Two phylogenetic analyses of ASPM exon 3 and MCPH1 exons 8 and 11 were performed to maximize taxon sampling or sequence extension to compare the nucleotide substitution and encephalization rates, and examine signals of selection. Further assessment of selection among humans was done through the analysis of non-synonymous and synonymous substitutions (dN/dS), and linkage disequilibrium (LD) patterns. We found that the accelerated evolution of brain size in hominids, is related to synchronic acceleration in the substitution rates of ASPM and MCPH1, and to signals of positive selection, especially in hominins. The dN/dS and LD analyses in Homo detected sites under positive selection and some regions with haplotype blocks at several candidate sites surrounded by blocks in LD-equilibrium. Accelerations and signals of positive selection in ASPM and MCPH1 occurred in different lineages and periods being ASPM more closely related with the brain evolution of hominins. MCPH1 evolved under positive selection in different lineages of the Catarrhini, suggesting independent evolutionary roles of this gene among primates.

Identifying our most distant animal relatives has emerged as one of the most challenging problems in phylogenetics. This debate has major implications for our understanding of the origin of multicellular animals and of the earliest events in animal evolution, including the origin of the nervous system. Some analyses identify sponges as our most distant animal relatives (Porifera-sister hypothesis), and others identify comb jellies (Ctenophora-sister hypothesis). These analyses vary in many respects, making it difficult to interpret previous tests of these hypotheses. To gain insight into why different studies yield different results, an important next step in the ongoing debate, we systematically test these hypotheses by synthesizing 15 previous phylogenomic studies and performing new standardized analyses under consistent conditions with additional models. We find that Ctenophora-sister is recovered across the full range of examined conditions, and Porifera-sister is recovered in some analyses under narrow conditions when most outgroups are excluded and site-heterogeneous CAT models are used. We additionally find that the number of categories in site-heterogeneous models is sufficient to explain the Porifera-sister results. Furthermore, our cross-validation analyses show CAT models that recover Porifera-sister have hundreds of additional categories and fail to fit significantly better than site-heterogenuous models with far fewer categories. Systematic and standardized testing of diverse phylogenetic models suggests that we should be skeptical of Porifera-sister results both because they are recovered under such narrow conditions and because the models in these conditions fit the data no better than other models that recover Ctenophora-sister.

Vitellogenin (Vg), a storage protein, has been significantly studied for its egg yolk precursor role in oviparous animals. Recent studies found that vitellogenin and its Vg‐like homologs were fundamentally involved in many other biological processes in social insects such as female caste differences and oxidative stress resilience. In this study, we conducted the first large‐scale molecular evolutionary analyses of vitellogenin coding genes (Vg) and Vg‐like genes of bumble bees, a primitively eusocial insect belonging to the genus Bombus. We obtained sequences for each of the four genes (Vg, Vg‐like‐A, Vg‐like‐B, and Vg‐like‐C) from 27 bumble bee genomes (nine were newly sequenced in this study), and sequences from the two closest clades of Bombus, including five Apis species and five Tetragonula species. Our molecular evolutionary analyses show that in bumble bee, the conventional Vg experienced strong positive selection, while the Vg‐like genes showed overall relaxation of purifying selection. In Apis and Tetragonula; however, all four genes were found under purifying selection. Furthermore, the conventional Vg showed signs of strong positive selection in most subgenera in Bombus, apart from the obligate parasitic subgenus Psithyrus which has no caste differentiation. Together, these results indicate that the conventional Vg, a key pleiotropic gene in social insects, is the most rapidly evolving copy, potentially due to its multiple known social functions for both worker and queen castes. This study shows that concerted evolution and purifying selection shaped the evolution of the Vg gene family following their ancient gene duplication and may be the leading forces behind the evolution of new potential protein function enabling functional social pleiotropy., In this study, we conducted the first large‐scale molecular evolutionary analyses of vitellogenin‐coding genes (Vg) and Vg‐like genes (Vg‐like‐A, Vg‐like‐B, and Vg‐like‐C) of bumble bees, a primitively eusocial insects belonging to the genus Bombus. The overall dN/dS ratio of the conventional Vg was 1.311, and this gene was under strong positive selection in most (9/10) subgenera in Bombus, apart from the obligate parasitic subgenus Psithyrus (dN/dS = 0.713) which has no caste differentiation.

Reconstructing the histories of complex adaptations and identifying the evolutionary mechanisms underlying their origins are two of the primary goals of evolutionary biology. Taricha newts, which contain high concentrations of the deadly toxin tetrodotoxin (TTX) as an antipredator defense, have evolved resistance to self-intoxication, which is a complex adaptation requiring changes in six paralogs of the voltage-gated sodium channel (Nav) gene family, the physiological target of TTX. Here, we reconstruct the origins of TTX self-resistance by sequencing the entire Nav gene family in newts and related salamanders. We show that moderate TTX resistance evolved early in the salamander lineage in three of the six Nav paralogs, preceding the proposed appearance of tetrodotoxic newts by ∼100 million years. TTX-bearing newts possess additional unique substitutions across the entire Nav gene family that provide physiological TTX resistance. These substitutions coincide with signatures of positive selection and relaxed purifying selection, as well as gene conversion events, that together likely facilitated their evolution. We also identify a novel exon duplication within Nav1.4 encoding an expressed TTX-binding site. Two resistance-conferring changes within newts appear to have spread via nonallelic gene conversion: in one case, one codon was copied between paralogs, and in the second, multiple substitutions were homogenized between the duplicate exons of Nav1.4. Our results demonstrate that gene conversion can accelerate the coordinated evolution of gene families in response to a common selection pressure.

Unmitigated spread in Brazil Despite an extensive network of primary care availability, Brazil has suffered profoundly during the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic. Using daily data from state health offices, Castro et al. analyzed the pattern of spread of COVID-19 cases and deaths in the country from February to October 2020. Clusters of deaths before cases became apparent indicated unmitigated spread. SARS-CoV-2 circulated undetected in Brazil for more than a month as it spread north from Sã o Paulo. In Manaus, transmission reached unprecedented levels after a momentary respite in mid-2020. Faria et al. tracked the evolution of a new, more aggressive lineage called P.1, which has 17 mutations, including three (K417T, E484K, and N501Y) in the spike protein. After a period of accelerated evolution, this variant emerged in Brazil during November 2020. Coupled with the emergence of P.1, disease spread was accelerated by stark local inequalities and political upheaval, which compromised a prompt federal response. Science , abh1558 and abh2644, this issue p. 821 and p. 815

Members of the genus Acidithiobacillus, now ranked within the class Acidithiobacillia, are model bacteria for the study of chemolithotrophic energy conversion under extreme conditions. Knowledge of the genomic and taxonomic diversity of Acidithiobacillia is still limited. Here, we present a systematic analysis of nearly 100 genomes from the class sampled from a wide range of habitats. Some of these genomes are new and others have been reclassified on the basis of advanced genomic analysis, thus defining 19 Acidithiobacillia lineages ranking at different taxonomic levels. This work provides the most comprehensive classification and pangenomic analysis of this deep-branching class of Proteobacteria to date. The phylogenomic framework obtained illuminates not only the evolutionary past of this lineage, but also the molecular evolution of relevant aerobic respiratory proteins, namely the cytochrome bo3 ubiquinol oxidases.

Obligate scavenging on the dead and decaying animal matter is a rare dietary specialization that in extant vertebrates is restricted to vultures. These birds perform essential ecological services, yet many vulture species have undergone recent steep population declines and are now endangered. To test for molecular adaptations underlying obligate scavenging in vultures, and to assess whether genomic features might have contributed to their population declines, we generated high-quality genomes of the Himalayan and bearded vultures, representing both independent origins of scavenging within the Accipitridae, alongside a sister taxon, the upland buzzard. By comparing our data to published sequences from other birds, we show that the evolution of obligate scavenging in vultures has been accompanied by widespread positive selection acting on genes underlying gastric acid production, and immunity. Moreover, we find evidence of parallel molecular evolution, with amino acid replacements shared among divergent lineages of these scavengers. Our genome-wide screens also reveal that both the Himalayan and bearded vultures exhibit low levels of genetic diversity, equating to around a half of the mean genetic diversity of other bird genomes examined. However, demographic reconstructions indicate that population declines began at around the Last Glacial Maximum, predating the well-documented dramatic declines of the past three decades. Taken together, our genomic analyses imply that vultures harbor unique adaptations for processing carrion, but that modern populations are genetically depauperate and thus especially vulnerable to further genetic erosion through anthropogenic activities.

The circadian system of cyanobacteria is built upon a central oscillator consisting of three genes, kaiA, kaiB, and kaiC. The KaiA protein plays a key role in phosphorylation/dephosphorylation cycles of KaiC, which occur over the 24-h period. We conducted a comprehensive evolutionary analysis of the kaiA genes across cyanobacteria. The results show that, in contrast to the previous reports, kaiA has an ancient origin and is as old as cyanobacteria. The kaiA homologs are present in nearly all analyzed cyanobacteria, except Gloeobacter, and have varying domain architecture. Some Prochlorococcales, which were previously reported to lack the kaiA gene, possess a drastically truncated homolog. The existence of the diverse kaiA homologs suggests significant variation of the circadian mechanism, which was described for the model cyanobacterium, Synechococcus elongatus PCC7942. The major structural modifications in the kaiA genes (duplications, acquisition and loss of domains) have apparently been induced by global environmental changes in the different geological periods.

How much does natural selection, as opposed to genetic drift, admixture, and gene flow, contribute to the evolution of invasive species following introduction to a new environment? Here we assess how evolution can shape biological invasions by examining population genomic variation in non‐native guppies (Poecilia reticulata) introduced to the Hawaiian Islands approximately a century ago. By examining 18 invasive populations from four Hawaiian islands and four populations from the native range in northern South America, we reconstructed the history of introductions and evaluated population structure as well as the extent of ongoing gene flow across watersheds and among islands. Patterns of differentiation indicate that guppies have developed significant population structure, with little natural or human‐mediated gene flow having occurred among populations following introduction. Demographic modeling and admixture graph analyses together suggest that guppies were initially introduced to O‘ahu and Maui and then translocated to Hawai‘i and Kaua‘i. We detected evidence for only one introduction event from the native range, implying that any adaptive evolution in introduced populations likely utilized the genetic variation present in the founding population. Environmental association tests accounting for population structure identified loci exhibiting signatures of adaptive variation related to predators and landscape characteristics but not nutrient regimes. When paired with high estimates of effective population sizes and detectable population structure, the presence of environment‐associated loci supports the role of natural selection in shaping contemporary evolution of Hawaiian guppy populations. Our findings indicate that local adaptation may engender invasion success, particularly in species with life histories that facilitate rapid evolution. Finally, evidence of low gene flow between populations suggests that removal could be an effective approach to control invasive guppies across the Hawaiian archipelago.

Background There has been a rapid increase in the brain size relative to body size during mammalian evolutionary history. In particular, the enlarged and globular brain is the most distinctive anatomical feature of modern humans that set us apart from other extinct and extant primate species. Genetic basis of large brain size in modern humans has largely remained enigmatic. Genes associated with the pathological reduction of brain size (primary microcephaly-MCPH) have the characteristics and functions to be considered ideal candidates to unravel the genetic basis of evolutionary enlargement of human brain size. For instance, the brain size of microcephaly patients is similar to the brain size of Pan troglodyte and the very early hominids like the Sahelanthropus tchadensis and Australopithecus afarensis. Results The present study investigates the molecular evolutionary history of subset of autosomal recessive primary microcephaly (MCPH) genes; CEP135, ZNF335, PHC1, SASS6, CDK6, MFSD2A, CIT, and KIF14 across 48 mammalian species. Codon based substitutions site analysis indicated that ZNF335, SASS6, CIT, and KIF14 have experienced positive selection in eutherian evolutionary history. Estimation of divergent selection pressure revealed that almost all of the MCPH genes analyzed in the present study have maintained their functions throughout the history of placental mammals. Contrary to our expectations, human-specific adoptive evolution was not detected for any of the MCPH genes analyzed in the present study. Conclusion Based on these data it can be inferred that protein-coding sequence of MCPH genes might not be the sole determinant of increase in relative brain size during primate evolutionary history.

Asgard is a recently discovered superphylum of archaea that appears to include the closest archaeal relatives of eukaryotes1–5. Debate continues as to whether the archaeal ancestor of eukaryotes belongs within the Asgard superphylum or whether this ancestor is a sister group to all other archaea (that is, a two-domain versus a three-domain tree of life)6–8. Here we present a comparative analysis of 162 complete or nearly complete genomes of Asgard archaea, including 75 metagenome-assembled genomes that—to our knowledge—have not previously been reported. Our results substantially expand the phylogenetic diversity of Asgard and lead us to propose six additional phyla that include a deep branch that we have provisionally named Wukongarchaeota. Our phylogenomic analysis does not resolve unequivocally the evolutionary relationship between eukaryotes and Asgard archaea, but instead—depending on the choice of species and conserved genes used to build the phylogeny—supports either the origin of eukaryotes from within Asgard (as a sister group to the expanded Heimdallarchaeota–Wukongarchaeota branch) or a deeper branch for the eukaryote ancestor within archaea. Our comprehensive protein domain analysis using the 162 Asgard genomes results in a major expansion of the set of eukaryotic signature proteins. The Asgard eukaryotic signature proteins show variable phyletic distributions and domain architectures, which is suggestive of dynamic evolution through horizontal gene transfer, gene loss, gene duplication and domain shuffling. The phylogenomics of the Asgard archaea points to the accumulation of the components of the mobile archaeal ‘eukaryome’ in the archaeal ancestor of eukaryotes (within or outside Asgard) through extensive horizontal gene transfer. Comparative analysis of 162 genomes of Asgard archaea results in six newly proposed phyla, including a deep branch that is provisionally named Wukongarchaeota, and sheds light on the evolutionary history of this clade.

Sulfate-reducing bacteria Candidatus Desulforudis audaxviator (CDA) were originally discovered in deep fracture fluids accessed via South African gold mines and have since been found in geographically widespread deep subsurface locations. In order to constrain models for subsurface microbial evolution, we compared CDA genomes from Africa, North America and Eurasia using single cell genomics. Unexpectedly, 126 partial single amplified genomes from the three continents, a complete genome from of an isolate from Eurasia, and metagenome-assembled genomes from Africa and Eurasia shared >99.2% average nucleotide identity, low frequency of SNP’s, and near-perfectly conserved prophages and CRISPRs. Our analyses reject sample cross-contamination, recent natural dispersal, and unusually strong purifying selection as likely explanations for these unexpected results. We therefore conclude that the analyzed CDA populations underwent only minimal evolution since their physical separation, potentially as far back as the breakup of Pangea between 165 and 55 Ma ago. High-fidelity DNA replication and repair mechanisms are the most plausible explanation for the highly conserved genome of CDA. CDA presents a stark contrast to the current model organisms in microbial evolutionary studies, which often develop adaptive traits over far shorter periods of time.

SFTS virus (SFTSV) is a novel bunyavirus, which was discovered as the etiological agent of severe fever with thrombocytopenia syndrome (SFTS) in China in 2009, and was now prevalent in at least 25 provinces in China. SFTS was subsequently identified in South Korea and Japan in 2012. To explore the molecular evolution and genetic characteristics of this newly identified pathogen, we reported 72 whole genome sequences of SFTSV, and built a dataset of SFTSV genome sequences containing 292 L-segment, 302 M-segment and 502 S-segment. We clearly divided SFTSV into six genotypes, Genotype A-F. It was found that genotype F was the dominant epidemic genotype of Japan, South Korea, and Zhejiang province of China. The coalescent analysis supported that SFTSV originated in the early 18th century from Zhejiang province, and Genotype F was the most primitive one. Henan, Hubei, and Anhui provinces which are located in Dabie Mountain area were mainly epidemic of Genotype A, which emerged relatively late but distributed widely. A total of 37 recombination events were identified, making SFTSV with a high recombination frequency (L segment 5.1%, M segment 3.6%, S segment 0.8%) among negative-strand segmented RNA viruses. It was identified that 19 reassortant strains belonged to 12 reassortment forms of SFTSV genome containing 6 newly identified forms. The reassortment virus and recombination in tick were both found for the first time. We also found many of genotype-specific mutation sites, 7 of which could be considered as potential molecular marker for genotype classification. This study promoted a more comprehensive understanding of the phylogeny and origin, and the genetic diversity of SFTSV, and it could help the studies of other newly discovered tick-borne bunyavirus as reference data and research ideas.

In the complex process of homeostasis of phytohormones cytokinins (CKs), O-glucosylation catalyzed by specific O-glucosyltransferases represents one of important mechanisms of their reversible inactivation. The CK O-glucosyltransferases belong to a highly divergent and polyphyletic multigene superfamily of glycosyltransferases, of which subfamily 1 containing UDP-glycosyltransferases (UGTs) is the largest in the plant kingdom. It contains recently discovered O and P subfamilies present in higher plant species but not in Arabidopsis thaliana. The cis-zeatin O-glucosyltransferase (cisZOG) genes belong to the O subfamily encoding a stereo-specific O-glucosylation of cis-zeatin-type CKs. We studied different homologous genes, their domains and motifs, and performed a phylogenetic reconstruction to elucidate the plant evolution of the cisZOG gene. We found that the cisZOG homologs do not form a clear separate clade, indicating that diversification of the cisZOG gene took place after the diversification of the main angiosperm families, probably within genera or closely related groups. We confirmed that the gene(s) from group O is(are) not present in A. thaliana and is(are) also missing in the family Brassicaceae. However, cisZOG or its metabolites are found among Brassicaceae clade, indicating that remaining genes from other groups (UGT73—group D and UGT85—group G) are able, at least in part, to substitute the function of group O lost during evolution. This study is the first detailed evolutionary evaluation of relationships among different plant ZOGs within angiosperms.

The infraorder Brachyura (true or short-tailed crabs) represents a successful group of marine invertebrates yet with limited genomic resources. Here we report a chromosome-anchored reference genome and transcriptomes of the Chinese mitten crab Eriocheir sinensis, a catadromous crab and invasive species with wide environmental tolerance, strong osmoregulatory capacity and high fertility. We show the expansion of specific gene families in the crab, including F-ATPase, which enhances our knowledge on the adaptive plasticity of this successful invasive species. Our analysis of spatio-temporal transcriptomes and the genome of E. sinensis and other decapods shows that brachyurization development is associated with down-regulation of Hox genes at the megalopa stage when tail shortening occurs. A better understanding of the molecular mechanism regulating sexual development is achieved by integrated analysis of multiple omics. These genomic resources significantly expand the gene repertoire of Brachyura, and provide insights into the biology of this group, and Crustacea in general., Brachyurans, or crabs, are of commercial and ecological importance, but limited genomic resources exist. Here the authors present a chromosome-level genome and expression data for the Chinese mitten crab, to shed light on the biology of this group.

The shoot apical meristems (SAMs) of land plants are crucial for plant growth and organ formation. In several angiosperms, the HAIRY MERISTEM (HAM) genes function as key regulators that control meristem development and stem cell homeostasis. To date, the origin and evolutionary history of the HAM family in land plants remains unclear. Potentially shared and divergent functions of HAM family members from angiosperms and non-angiosperms are also not known. In constructing a comprehensive phylogeny of the HAM family, we show that HAM proteins are widely present in land plants and that HAM proteins originated prior to the divergence of bryophytes. The HAM family was duplicated in a common ancestor of angiosperms, leading to two distinct groups: type I and type II. Type-II HAM members are widely present in angiosperms, whereas type-I HAM members were independently lost in different orders of monocots. Furthermore, HAM members from angiosperms and non-angiosperms (including bryophytes, lycophytes, ferns and gymnosperms) are able to replace the role of the type-II HAM genes in Arabidopsis, maintaining established SAMs and promoting the initiation of new stem cell niches. Our results uncover the conserved functions of HAM family members and reveal the conserved regulatory mechanisms underlying HAM expression patterning in meristems, providing insight into the evolution of key stem cell regulators in land plants.

Chloroplast genome sequences have been used to understand evolutionary events and to infer efficiently phylogenetic relationships. Callitropsis funebris (Cupressaceae) is an endemic species in China. Its phylogenetic position is controversial due to morphological characters similar to those of Cupressus, Callitropsis, and Chamaecyparis. This study used next‐generation sequencing technology to sequence the complete chloroplast genome of Ca. funebris and then constructed the phylogenetic relationship between Ca. funebris and its related species based on a variety of data sets and methods. Simple sequence repeats (SSRs) and adaptive evolution analysis were also conducted. Our results showed that the monophyletic branch consisting of Ca. funebris and Cupressus tonkinensis is a sister to Cupressus, while Callitropsis is not monophyletic; Ca. nootkatensis and Ca. vietnamensis are nested in turn at the base of the monophyletic group Hesperocyparis. The statistical results of SSRs supported the closest relationship between Ca. funebris and Cupressus. By performing adaptive evolution analysis under the phylogenetic background of Cupressales, the Branch model detected three genes and the Site model detected 10 genes under positive selection; and the Branch‐Site model uncovered that rpoA has experienced positive selection in the Ca. funebries branch. Molecular analysis from the chloroplast genome highly supported that Ca. funebris is at the base of Cupressus. Of note, SSR features were found to be able to shed some light on phylogenetic relationships. In short, this chloroplast genomic study has provided new insights into the phylogeny of Ca. funebris and revealed multiple chloroplast genes possibly undergoing adaptive evolution., The complete chloroplast genome sequence of Callitropsis funebris was obtained for the first time. Phylogenetic tree construction based on various methods and simple sequence repeat characteristic analysis showed that Ca. funebri s was located at the base of Cupressus. In the process of Cupressales evolution, positive selection branches were detected in 3 genes, positive selection sites were detected in 10 genes, and positive selection sites were detected in rpoA gene in Ca. funebris.

Photosynthetic organisms have developed diverse antennas composed of chromophorylated proteins to increase photon capture. Cryptophyte algae acquired their photosynthetic organelles (plastids) from a red alga by secondary endosymbiosis. Cryptophytes lost the primary red algal antenna, the red algal phycobilisome, replacing it with a unique antenna composed of αβ protomers, where the β subunit originates from the red algal phycobilisome. The origin of the cryptophyte antenna, particularly the unique α subunit, is unknown. Here we show that the cryptophyte antenna evolved from a complex between a red algal scaffolding protein and phycoerythrin β. Published cryo-EM maps for two red algal phycobilisomes contain clusters of unmodelled density homologous to the cryptophyte-αβ protomer. We modelled these densities, identifying a new family of scaffolding proteins related to red algal phycobilisome linker proteins that possess multiple copies of a cryptophyte-α-like domain. These domains bind to, and stabilise, a conserved hydrophobic surface on phycoerythrin β, which is the same binding site for its primary partner in the red algal phycobilisome, phycoerythrin α. We propose that after endosymbiosis these scaffolding proteins outcompeted the primary binding partner of phycoerythrin β, resulting in the demise of the red algal phycobilisome and emergence of the cryptophyte antenna., Cryptophytes acquired plastids from red algae but replaced the light-harvesting phycobilisome with a unique cryptophyte antenna. Here via analysis of phycobilisome cryo-EM structures, Rathbone et al. propose that the α subunit of the cryptophyte antenna originated from phycobilisome linker proteins

Small RNAs play a major role in the post-transcriptional regulation of gene expression in eukaryotes. Despite the evolutionary importance of streptophyte algae, knowledge on small RNAs in this group of green algae is almost non-existent. We used genome and transcriptome data of 34 algal and plant species, and performed genome-wide analyses of small RNA (miRNA & siRNA) biosynthetic and degradation pathways. The results suggest that Viridiplantae started to evolve plant-like miRNA biogenesis and degradation after the divergence of the Mesostigmatophyceae in the streptophyte algae. We identified two major evolutionary transitions in small RNA metabolism in streptophyte algae; during the first transition, the origin of DCL-New, DCL1, AGO1/5/10 and AGO4/6/9 in the last common ancestor of Klebsormidiophyceae and all other streptophytes could be linked to abiotic stress responses and evolution of multicellularity in streptophytes. During the second transition, the evolution of DCL 2,3,4, and AGO 2,3,7 as well as DRB1 in the last common ancestor of Zygnematophyceae and embryophytes, suggests their possible contribution to pathogen defense and antibacterial immunity. Overall, the origin and diversification of DICER and AGO along with several other small RNA pathway-related genes among streptophyte algae suggested progressive adaptations of streptophyte algae during evolution to a subaerial environment., Wang, Liang et al. conduct a genome-wide investigation into the origin of small RNA pathway-related genes in Viridiplantae. Their findings suggest that streptophyte algae progressively adapted to a subaerial environment through generation of these pathways.

We reconstructed the 2016-2017 Zika virus epidemic in Puerto Rico by using complete genomes to uncover the epidemic's origin, spread, and evolutionary dynamics. Our study revealed that the epidemic was propelled by multiple introductions that spread across the island, intricate evolutionary patterns, and ≈10 months of cryptic transmission.

The DELAY OF GERMINATION1 (DOG1) family genes (DFGs) in Arabidopsis thaliana are involved in seed dormancy, reserve accumulation, and desiccation tolerance. Decoding the molecular evolution of DFGs is key to understanding how these seed programs evolved. This article demonstrates that DFGs have diverged in the four lineages DOG1, DOG1-LIKE4 (DOGL4), DOGL5 and DOGL6, whereas DOGL1, DOGL2 and DOGL3 arose separately within the DOG1 lineage. The systematic DFG nomenclature proposed in this article addresses the current issues of inconsistent DFG annotation and highlights DFG genomic synteny in angiosperms. DFG pseudogenes, or collapsed coding sequences, hidden in the genomes of early-diverging angiosperms are documented here. They suggest ancient birth and loss of DFGs over the course of angiosperm evolution. The proposed models suggest that the origin of DFG diversification dates back to the most recent common ancestor of living angiosperms. The presence of a single form of DFG in nonflowering plants is discussed. Phylogenetic analysis of gymnosperm, lycophyte, and liverwort DFGs and similar genes found in mosses and algae suggests that DFGs diverged from the TGACG motif-binding transcription factor genes before the divergence of the bryophyte lineage.

Recent pangenome studies have revealed a large fraction of the gene content within a species exhibits presence–absence variation (PAV). However, coding regions alone provide an incomplete assessment of functional genomic sequence variation at the species level. Little to no attention has been paid to noncoding regulatory regions in pangenome studies, though these sequences directly modulate gene expression and phenotype. To uncover regulatory genetic variation, we generated chromosome-scale genome assemblies for thirty Arabidopsis thaliana accessions from multiple distinct habitats and characterized species level variation in Conserved Noncoding Sequences (CNS). Our analyses uncovered not only PAV and positional variation (PosV) but that diversity in CNS is nonrandom, with variants shared across different accessions. Using evolutionary analyses and chromatin accessibility data, we provide further evidence supporting roles for conserved and variable CNS in gene regulation. Additionally, our data suggests that transposable elements contribute to CNS variation. Characterizing species-level diversity in all functional genomic sequences may later uncover previously unknown mechanistic links between genotype and phenotype.

Proprotein convertase subtilisin/kexin type 9 (PCSK9) plays a central role in cholesterol homeostasis in humans as a major regulator of LDLR levels. PCSK9 is an intriguing protease in that it does not act by proteolysis but by preventing LDLR recirculation from endosomes to the plasma membrane. This, and the inexistence of any other proteolytic substrate but itself could suggest that PCSK9 is an exquisite example of evolutionary fine-tuning. However, the gene has been lost in several mammalian species, and null alleles are present (albeit at low frequencies) in some human populations without apparently deleterious health effects, raising the possibility that the PCSK9 may have become dispensable in the mammalian lineage. To address this issue, we systematically recovered, assembled, corrected, annotated and analysed publicly available PCSK9 sequences for 420 eutherian species to determine the distribution, frequencies, mechanisms and timing of PCSK9 pseudogenization events, as well as the evolutionary pressures underlying the preservation or loss of the gene. We found a dramatic difference in the patterns of PCSK9 retention and loss between Euarchontoglires—where there is strong pressure for gene preservation—and Laurasiatheria, where multiple independent events have led to PCSK9 loss in most species. These results suggest that there is a fundamental difference in the regulation of cholesterol metabolism between Euarchontoglires and Laurasiatheria, which in turn has important implications for the use of Laurasiatheria species (e.g. pigs) as animal models of human cholesterol-related diseases.

Citrus tristeza virus (CTV) is an important citrus pathogen causing considerable economic loss to citrus production. Knowledge on genetic evolutionary of the CTV population in China remains limited. In this study, 1439 samples were collected from nine citrus-producing areas of China. The coat protein (CP) genes of CTV were amplified by RT-PCR, and sequenced to analyze the genetic evolution. Analysis of the base composition showed an AU preference pattern, with the GC content was lower than AU content. Nine CTV populations were clustered into one clade in neighbor-joining (NJ) tree, indicative of a close phylogenetic relationship among the populations in China. Analysis of molecular variation (AMOVA) revealed that 77.72% genetic variations of CTV populations were observed among populations, with an FST value of 0.223. The values of dN/dS and neutrality test of CP gene were ranged from 0.016 to 0.082 and -1.377 to 1.456, respectively, the results suggesting that all of nine CTV populations were relatively constantly maintained under purifying selection. Our study demonstrated the genetic characteristics and molecular evolution relationship of CTV populations in China, and provided a theoretical basis for scientific control of CTV.

The metallo-β-lactamase fold is an ancient protein structure present in numerous enzyme families responsible for diverse biological processes. The crystal structure of the hyperthermostable crenarchaeal enzyme Igni18 from Ignicoccus hospitalis was solved at 2.3 Å and could resemble a possible first archetype of a multifunctional metallo-β-lactamase. Ancestral enzymes at the evolutionary origin are believed to be promiscuous all-rounders. Consistently, Igni18´s activity can be cofactor-dependently directed from β-lactamase to lactonase, lipase, phosphodiesterase, phosphotriesterase or phospholipase. Its core-domain is highly conserved within metallo-β-lactamases from Bacteria, Archaea and Eukarya and gives insights into evolution and function of enzymes from this superfamily. Structural alignments with diverse metallo-β-lactamase-fold-containing enzymes allowed the identification of Protein Variable Regions accounting for modulation of activity, specificity and oligomerization patterns. Docking of different substrates within the active sites revealed the basis for the crucial cofactor dependency of this enzyme superfamily., The structure of Igni18, ancient enzyme from the Crenarchaean species Ignicoccus hospitalis, is solved and characterized by Pérez-García, Chow and colleagues. This structure provides insight as to the evolution of metallo-beta-lactamases and their functions.

Viroids belong to a very interesting class of molecules attracting researchers in phytopathology and molecular evolution. Here we review recent literature data concerning the genetics of Potato spindle tuber viroid (PSTVd) and the mechanisms related to its pathological effect on the host plants. PSTVd can be transmitted vertically through microspores and macrospores, but not with pollen from another infected plant. The 359 nucleotidelong genomic RNA of PSTVd is highly structured and its 3D-conformation is responsible for interaction with host cellular factors to mediate replication, transport between tissues during systemic infection and the severity of pathological symptoms. RNA replication is prone to errors and infected plants contain a population of mutated forms of the PSTVd genome. Interestingly, at 7 DAI, only 25 % of the newly synthesized RNAs were identical to the master copy, but this proportion increased to up to 70 % at 14 DAI and remained the same afterwards. PSTVd infection induces the immune response in host plants. There are PSTVd strains with a severe, a moderate or a mild pathological effect. Interestingly, viroid replication itself does not necessarily induce strong morphological or physiological symptoms. In the case of PSTVd, disease symptoms may occur due to RNA-interference, which decreases the expression levels of some important cellular regulatory factors, such as, for example, potato StTCP23 from the gibberellic acid pathway with a role in tuber morphogenesis or tomato FRIGIDA-like protein 3 with an early flowering phenotype. This association between the small segments of viroid genomic RNAs complementary to the untranslated regions of cellular mRNAs and disease symptoms provides a way for new resistant cultivars to be developed by genetic editing. To conclude, viroids provide a unique model to reveal the fundamental features of living systems, which appeared early in evolution and still remain undiscovered.Вироиды, небольшие кольцевые молекулы РНК, которые вызывают патогенез у растений, остаются одним из самых необычных биологических объектов, привлекающих внимание не только фитопатологов, но и специалистов в области молекулярной эволюции. В статье приведен обзор последних литературных данных о генетике вироида веретеновидности клубней картофеля (ВВКК) и генетических механизмах формирования патологических состояний у растений-хозяев. ВВКК способен передаваться вертикально (через генеративные клетки зараженного растения), но, в отличие от некоторых других вироидов, не передается через пыльцу от зараженного растения. Большой интерес у исследователей вызывает структура геномной РНК вироида размером 359 нуклеотидов: хорошо известно, что особенности 3D конформации определяют основные параметры взаимодействия с клеточными факторами на стадии репликации, транспорта между различными тканями в процессе системной инфекции, а также степень выраженности симптомов заболевания. При репликации геномной РНК вироидов часто происходят ошибки, приводящие к появлению гетерогенной популяции молекул РНК в тканях зараженного растения. Примечательно, что через 7 дней после инокуляции только 25 % молекул геномной РНК ВВКК соответствовали исходной матрице, использованной для инокуляции, однако эта доля увеличилась до 70 % через 14 дней и далее оставалась на том же уровне. По-видимому, при сохранении у мутантных вариантов геномной РНК способности к репликации вироид обладает высоким потенциалом к отбору эффективных инфекционных форм. ВВКК вызывает у пораженных растений развитие иммунного ответа, механизмы индукции которого недостаточно изучены. Известны сильно- и слабопатогенные штаммы ВВКК, вызывающие разные проявления болезни, фенотипические проявления от которых у пораженных растений в значительной мере различны. Сама по себе репликация вироида не обязательно приводит к выраженным фенотипическим проявлениям, в случае ВВКК они могут быть связаны с участками гомологии между геномной РНК и мРНК некоторых регуляторных генов, например транскрипционного фактора StTCP23 картофеля, участвующего в регуляторном контуре гиббереллиновой кислоты и в контроле морфогенеза клубня. Другой пример – индукция РНК-интерференции против мРНК гена FRIGIDA-like protein 3 у томата, что приводит к раннему цветению. В связи с этим обсуждаются потенциальные способы борьбы с вироидом, основанные на удалении из генома растений таких участков гомологии, расположенных в нетранслируемых областях мРНК и не выполняющих каких-либо функций. В целом вироиды представляют собой уникальную модель для исследования основ организации живых систем, многие из которых возникли на ранних этапах эволюции и остаются до сих пор не выявленными.

The structure and gene sequence of the fish mitochondrial genome are generally considered to be conservative. However, two types of gene arrangements are found in the mitochondrial genome of Anguilliformes. In this paper, we report a complete mitogenome of Muraenesox cinereus (Anguilliformes: Muraenesocidae) with rearrangement phenomenon. The total length of the M. cinereus mitogenome was 17,673 bp, and it contained 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNA genes, and two identical control regions (CRs). The mitochondrial genome of M. cinereus was obviously rearranged compared with the mitochondria of typical vertebrates. The genes ND6 and the conjoint trnE were translocated to the location between trnT and trnP, and one of the duplicated CR was translocated to the upstream of the ND6. The tandem duplication and random loss is most suitable for explaining this mitochondrial gene rearrangement. The Anguilliformes phylogenetic tree constructed based on the whole mitochondrial genome well supports Congridae non-monophyly. These results provide a basis for the future Anguilliformes mitochondrial gene arrangement characteristics and further phylogenetic research.