Molecular evolution of Phytocyanin gene and analysis of expression at different coloring periods in apple (Malus domestica).

BACKGROUND: PC (phytocyanin) is a class of copper-containing electron transfer proteins closely related to plant photosynthesis, abiotic stress responses growth and development in plants, and regulation of the expression of some flavonoids and phenylpropanoids, etc., however, compared with other plants, the PC gene family has not been systematically characterized in apple. RESULTS: A total of 59 MdPC gene members unevenly distributed across 12 chromosomes were identified at the genome-wide level. The proteins of the MdPC family were classified into four subfamilies based on differences in copper binding sites and glycosylation sites: Apple Early nodulin-like proteins (MdENODLs), Apple Uclacyanin-like proteins (MdUCLs), Apple Stellacyanin-like proteins (MdSCLs), and Apple Plantacyanin-like proteins (MdPLCLs). Some MdPC members with similar gene structures and conserved motifs belong to the same group or subfamily. The internal collinearity analysis revealed 14 collinearity gene pairs among members of the apple MdPC gene. Interspecific collinearity analysis showed that apple had 31 and 35 homologous gene pairs with strawberry and grape, respectively. Selection pressure analysis indicated that the MdPC gene was under purifying selection. Prediction of protein interactions showed that MdPC family members interacted strongly with the Nad3 protein. GO annotation results indicated that the MdPC gene also regulated the biosynthesis of phenylpropanoids. Chip data analysis showed that (MdSCL3, MdSCL7 and MdENODL27) were highly expressed in mature fruits and peels. Many cis-regulatory elements related to light response, phytohormones, abiotic stresses and flavonoid biosynthetic genes regulation were identified 2000 bp upstream of the promoter of the MdPC gene, and qRT-PCR results showed that gene members in Group IV (MdSCL1/3, MdENODL27) were up-regulated at all five stages of apple coloring, but the highest expression was observed at the DAF13 (day after fruit bag removal) stage. The gene members in Group II (MdUCL9, MdPLCL3) showed down-regulated or lower expression in the first four stages of apple coloring but up-regulated and highest expression in the DAF 21 stage. CONCLUSION: Herein, one objective of these findings is to provide valuable information for understanding the structure, molecular evolution, and expression pattern of the MdPC gene, another major objective in this study was designed to lay the groundwork for further research on the molecular mechanism of PC gene regulation of apple fruit coloration.

Dynamics of karyotype evolution.

In the evolution of species, the karyotype changes with a timescale of tens to hundreds of thousand years. In the development of cancer, the karyotype often is modified in cancerous cells over the lifetime of an individual. Characterizing these changes and understanding the mechanisms leading to them has been of interest in a broad range of disciplines including evolution, cytogenetics, and cancer genetics. A central issue relates to the relative roles of random vs deterministic mechanisms in shaping the changes. Although it is possible that all changes result from random events followed by selection, many results point to other non-random factors that play a role in karyotype evolution. In cancer, chromosomal instability leads to characteristic changes in the karyotype, in which different individuals with a specific type of cancer display similar changes in karyotype structure over time. Statistical analyses of chromosome lengths in different species indicate that the length distribution of chromosomes is not consistent with models in which the lengths of chromosomes are random or evolve solely by simple random processes. A better understanding of the mechanisms underlying karyotype evolution should enable the development of quantitative theoretical models that combine the random and deterministic processes that can be compared to experimental determinations of the karyotype in diverse settings.

Similar enzymatic functions in distinct bioluminescence systems: evolutionary recruitment of sulfotransferases in ostracod light organs.

Genes from ancient families are sometimes involved in the convergent evolutionary origins of similar traits, even across vast phylogenetic distances. Sulfotransferases are an ancient family of enzymes that transfer sulfate from a donor to a wide variety of substrates, including probable roles in some bioluminescence systems. Here, we demonstrate multiple sulfotransferases, highly expressed in light organs of the bioluminescent ostracod Vargula tsujii, transfer sulfate in vitro to the luciferin substrate, vargulin. We find luciferin sulfotransferases (LSTs) of ostracods are not orthologous to known LSTs of fireflies or sea pansies; animals with distinct and convergently evolved bioluminescence systems compared to ostracods. Therefore, distantly related sulfotransferases were independently recruited at least three times, leading to parallel evolution of luciferin metabolism in three highly diverged organisms. Reuse of homologous genes is surprising in these bioluminescence systems because the other components, including luciferins and luciferases, are completely distinct. Whether convergently evolved traits incorporate ancient genes with similar functions or instead use distinct, often newer, genes may be constrained by how many genetic solutions exist for a particular function. When fewer solutions exist, as in genetic sulfation of small molecules, evolution may be more constrained to use the same genes time and again.

Genomic analysis and phylogenetic characterization of Himalayan snow trout, Schizothorax esocinus based on mitochondrial protein-coding genes.

BACKGROUND: Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes. METHODS: Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis. RESULTS: The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%. CONCLUSION: The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species.

Genome-wide identification, molecular evolution and gene expression of P450 gene family in Cyrtotrachelus buqueti.

BACKGROUND: Insect Cytochrome P450 monooxygenase (CYPs or P450s) plays an important role in detoxifying insecticides, causing insect populations to develop resistance. However, the molecular functions of P450 gene family in Cyrtotrachelus buqueti genome are still lacking. RESULTS: In this study, 71 CbuP450 genes have been identified. The amino acids length of CbuP450 proteins was between 183 aa ~ 1041 aa. They are proteins with transmembrane domains. The main component of their secondary structure is α-helix and random coils. Phylogenetic analysis showed that C. buqueti and Rhynchophorus ferrugineus were the most closely related. This gene family has 29 high-frequency codons, which tend to use A/T bases and A/T ending codons. Gene expression analysis showed that CbuP450_23 in the female adult may play an important role on high temperature resistance, and CbuP450_17 in the larval may play an important role on low temperature tolerance. CbuP450_10, CbuP450_17, CbuP450_23, CbuP450_10, CbuP450_16, CbuP450_20, CbuP450_23 and CbuP450_ 29 may be related to the regulation of bamboo fiber degradation genes in C. buqueti. Protein interaction analysis indicates that most CbuP450 proteins are mainly divided into three aspects: encoding the biosynthesis of ecdysteroids, participating in the decomposition of synthetic insecticides, metabolizing insect hormones, and participating in the detoxification of compounds. CONCLUSIONS: We systematically analyzed the gene and protein characteristics, gene expression, and protein interactions of CbuP450 gene family, revealing the key genes involved in the stress response of CbuP450 gene family in the resistance of C. buqueti to high or low temperature stress, and identified the key CbuP450 proteins involved in important life activity metabolism. These results provided a reference for further research on the function of P450 gene family in C. buqueti.

Characterizing genetic variation on the Z chromosome in Schistosoma japonicum reveals host-parasite co-evolution.

BACKGROUND: Schistosomiasis is a neglected tropical disease that afflicts millions of people worldwide; it is caused by Schistosoma, the only dioecious flukes with ZW systems. Schistosoma japonicum is endemic to Asia; the Z chromosome of S. japonicum comprises one-quarter of the entire genome. Detection of positive selection using resequencing data to understand adaptive evolution has been applied to a variety of pathogens, including S. japonicum. However, the contribution of the Z chromosome to evolution and adaptation is often neglected. METHODS: We obtained 1,077,526 high-quality SNPs on the Z chromosome in 72 S. japonicum using re-sequencing data publicly. To examine the faster Z effect, we compared the sequence divergence of S. japonicum with two closely related species, Schistosoma haematobium and S. mansoni. Genetic diversity was compared between the Z chromosome and autosomes in S. japonicum by calculating the nucleotide diversity (π) and Dxy values. Population structure was also assessed based on PCA and structure analysis. Besides, we employed multiple methods including Tajima's D, FST, iHS, XP-EHH, and CMS to detect positive selection signals on the Z chromosome. Further RNAi knockdown experiments were performed to investigate the potential biological functions of the candidate genes. RESULTS: Our study found that the Z chromosome of S. japonicum showed faster evolution and more pronounced genetic divergence than autosomes, although the effect may be smaller than the variation among genes. Compared with autosomes, the Z chromosome in S. japonicum had a more pronounced genetic divergence of sub-populations. Notably, we identified a set of candidate genes associated with host-parasite co-evolution. In particular, LCAT exhibited significant selection signals within the Taiwan population. Further RNA interference experiments suggested that LCAT is necessary for S. japonicum survival and propagation in the definitive host. In addition, we identified several genes related to the specificity of the intermediate host in the C-M population, including Rab6 and VCP, which are involved in adaptive immune evasion to the host. CONCLUSIONS: Our study provides valuable insights into the adaptive evolution of the Z chromosome in S. japonicum and further advances our understanding of the co-evolution of this medically important parasite and its hosts.

Shaping gene expression and its evolution by chromatin architecture and enhancer activity.

Transcriptional regulation plays a pivotal role in orchestrating the intricate genetic programs governing embryonic development. The expression of developmental genes relies on the combined activity of several cis-regulatory elements (CREs), such as enhancers and silencers, which can be located at long linear distances from the genes that they regulate and that interact with them through establishment of chromatin loops. Mutations affecting their activity or interaction with their target genes can lead to developmental disorders and are thought to have importantly contributed to the evolution of the animal body plan. The income of next-generation-sequencing approaches has allowed identifying over a million of sequences with putative regulatory potential in the human genome. Characterizing their function and establishing gene-CREs maps is essential to decode the logic governing developmental gene expression and is one of the major challenges of the post-genomic era. Chromatin 3D organization plays an essential role in determining how CREs specifically contact their target genes while avoiding deleterious off-target interactions. Our understanding of these aspects has greatly advanced with the income of chromatin conformation capture techniques and fluorescence microscopy approaches to visualize the organization of DNA elements in the nucleus. Here we will summarize relevant aspects of how the interplay between CRE activity and chromatin 3D organization regulates developmental gene expression and how it relates to pathological conditions and the evolution of animal body plan.

Specificities and commonalities of the Planctomycetes plasmidome.

Plasmids, despite their critical role in antibiotic resistance and modern biotechnology, are understood in only a few bacterial groups in terms of their natural ecological dynamics. The bacterial phylum Planctomycetes, known for its unique molecular and cellular biology, has a largely unexplored plasmidome. This study offers a thorough exploration of the diversity of natural plasmids within Planctomycetes, which could serve as a foundation for developing various genetic research tools for this phylum. Planctomycetes plasmids encode a broad range of biological functions and appear to have coevolved significantly with their host chromosomes, sharing many homologues. Recent transfer events of insertion sequences between cohabiting chromosomes and plasmids were also observed. Interestingly, 64% of plasmid genes are distantly related to either chromosomally encoded genes or have homologues in plasmids from other bacterial groups. The planctomycetal plasmidome is composed of 36% exclusive proteins. Most planctomycetal plasmids encode a replication initiation protein from the Replication Protein A family near a putative iteron-containing replication origin, as well as active type I partition systems. The identification of one conjugative and three mobilizable plasmids suggests the occurrence of horizontal gene transfer via conjugation within this phylum. This comprehensive description enhances our understanding of the plasmidome of Planctomycetes and its potential implications in antibiotic resistance and biotechnology.

Genetic diversity and evolution of G12P[6] DS-1-like and G12P[9] AU-1-like Rotavirus strains in Brazil.

In the early 2000s, the global emergence of rotavirus (RVA) G12P[8] genotype was noted, while G12P[6] and G12P[9] combinations remained rare in humans. This study aimed to characterize and phylogenetically analyze three Brazilian G12P[9] and four G12P[6] RVA strains from 2011 to 2020, through RT-PCR and sequencing, in order to enhance our understanding of the genetic relationship between human and animal-origin RVA strains. G12P[6] strains displayed a DS-1-like backbone, showing a distinct genetic clustering. G12P[6] IAL-R52/2020, IAL-R95/2020 and IAL-R465/2019 strains clustered with 2019 Northeastern G12P[6] Brazilian strains and a 2018 Benin strain, whereas IAL-R86/2011 strain grouped with 2010 Northern G12P[6] Brazilian strains and G2P[4] strains from the United States and Belgium. These findings suggest an African genetic ancestry and reassortments with co-circulating American strains sharing the same DS-1-like constellation. No recent zoonotic reassortment was observed, and the DS-1-like constellation detected in Brazilian G12P[6] strains does not seem to be genetically linked to globally reported intergenogroup G1/G3/G9/G8P[8] DS-1-like human strains. G12P[9] strains exhibited an AU-1-like backbone with two different genotype-lineage constellations: IAL-R566/2011 and IAL-R1151/2012 belonged to a VP3/M3.V Lineage, and IAL-R870/2013 to a VP3/M3.II Lineage, suggesting two co-circulating strains in Brazil. This genetic diversity is not observed elsewhere, and the VP3/M3.II Lineage in G12P[9] strains seems to be exclusive to Brazil, indicating its evolution within the country. All three G12P[9] AU-1-like strains were closely relate to G12P[9] strains from Paraguay (2006-2007) and Brazil (2010). Phylogenetic analysis also highlighted that all South American G12P[9] AU-1-like strains had a common origin and supports the hypothesis of their importation from Asia, with no recent introduction from globally circulating G12P[9] strains or reassortments with local G12 strains P[8] or P[6]. Notably, certain genes in the Brazilian G12P[9] AU-1-like strains share ancestry with feline/canine RVAs (VP3/M3.II, NSP4/E3.IV and NSP2/N3.II), whereas NSP1/A3.VI likely originated from artiodactyls, suggesting a history of zoonotic transmission with human strains. This genomic data adds understanding to the molecular epidemiology of G12P[6] and G12P[9] RVA strains in Brazil, offering insights into their genetic diversity and evolution.

Evolutionary genomics of socially polymorphic populations of Pogonomyrmex californicus.

BACKGROUND: Social insects vary considerably in their social organization both between and within species. In the California harvester ant, Pogonomyrmex californicus (Buckley 1867), colonies are commonly founded and headed by a single queen (haplometrosis, primary monogyny). However, in some populations in California (USA), unrelated queens cooperate not only during founding (pleometrosis) but also throughout the life of the colony (primary polygyny). The genetic architecture and evolutionary dynamics of this complex social niche polymorphism (haplometrosis vs pleometrosis) have remained unknown. RESULTS: We provide a first analysis of its genomic basis and evolutionary history using population genomics comparing individuals from a haplometrotic population to those from a pleometrotic population. We discovered a recently evolved (< 200 k years), 8-Mb non-recombining region segregating with the observed social niche polymorphism. This region shares several characteristics with supergenes underlying social polymorphisms in other socially polymorphic ant species. However, we also find remarkable differences from previously described social supergenes. Particularly, four additional genomic regions not in linkage with the supergene show signatures of a selective sweep in the pleometrotic population. Within these regions, we find for example genes crucial for epigenetic regulation via histone modification (chameau) and DNA methylation (Dnmt1). CONCLUSIONS: Altogether, our results suggest that social morph in this species is a polygenic trait involving a potential young supergene. Further studies targeting haplo- and pleometrotic individuals from a single population are however required to conclusively resolve whether these genetic differences underlie the alternative social phenotypes or have emerged through genetic drift.

Characterization of the genetic variation and evolutionary divergence of the CLEC18 family.

BACKGROUND: The C-type lectin family 18 (CLEC18) with lipid and glycan binding capabilities is important to metabolic regulation and innate immune responses against viral infection. However, human CLEC18 comprises three paralogous genes with highly similar sequences, making it challenging to distinguish genetic variations, expression patterns, and biological functions of individual CLEC18 paralogs. Additionally, the evolutionary relationship between human CLEC18 and its counterparts in other species remains unclear. METHODS: To identify the sequence variation and evolutionary divergence of human CLEC18 paralogs, we conducted a comprehensive analysis using various resources, including human and non-human primate reference genome assemblies, human pangenome assemblies, and long-read-based whole-genome and -transcriptome sequencing datasets. RESULTS: We uncovered paralogous sequence variants (PSVs) and polymorphic variants (PVs) of human CLEC18 proteins, and identified distinct signatures specific to each CLEC18 paralog. Furthermore, we unveiled a novel segmental duplication for human CLEC18A gene. By comparing CLEC18 across human and non-human primates, our research showed that the CLEC18 paralogy probably occurred in the common ancestor of human and closely related non-human primates, and the lipid-binding CAP/SCP/TAPS domain of CLEC18 is more diverse than its glycan-binding CTLD. Moreover, we found that certain amino acids alterations at variant positions are exclusive to human CLEC18 paralogs. CONCLUSIONS: Our findings offer a comprehensive profiling of the intricate variations and evolutionary characteristics of human CLEC18.

Epistasis facilitates functional evolution in an ancient transcription factor.

A protein's genetic architecture - the set of causal rules by which its sequence produces its functions - also determines its possible evolutionary trajectories. Prior research has proposed that the genetic architecture of proteins is very complex, with pervasive epistatic interactions that constrain evolution and make function difficult to predict from sequence. Most of this work has analyzed only the direct paths between two proteins of interest - excluding the vast majority of possible genotypes and evolutionary trajectories - and has considered only a single protein function, leaving unaddressed the genetic architecture of functional specificity and its impact on the evolution of new functions. Here, we develop a new method based on ordinal logistic regression to directly characterize the global genetic determinants of multiple protein functions from 20-state combinatorial deep mutational scanning (DMS) experiments. We use it to dissect the genetic architecture and evolution of a transcription factor's specificity for DNA, using data from a combinatorial DMS of an ancient steroid hormone receptor's capacity to activate transcription from two biologically relevant DNA elements. We show that the genetic architecture of DNA recognition consists of a dense set of main and pairwise effects that involve virtually every possible amino acid state in the protein-DNA interface, but higher-order epistasis plays only a tiny role. Pairwise interactions enlarge the set of functional sequences and are the primary determinants of specificity for different DNA elements. They also massively expand the number of opportunities for single-residue mutations to switch specificity from one DNA target to another. By bringing variants with different functions close together in sequence space, pairwise epistasis therefore facilitates rather than constrains the evolution of new functions.

Plastid phylogenomics and cytonuclear discordance in Rubioideae, Rubiaceae.

In this study of evolutionary relationships in the subfamily Rubioideae (Rubiaceae), we take advantage of the off-target proportion of reads generated via previous target capture sequencing projects based on nuclear genomic data to build a plastome phylogeny and investigate cytonuclear discordance. The assembly of off-target reads resulted in a comprehensive plastome dataset and robust inference of phylogenetic relationships, where most intratribal and intertribal relationships are resolved with strong support. While the phylogenetic results were mostly in agreement with previous studies based on plastome data, novel relationships in the plastid perspective were also detected. For example, our analyses of plastome data provide strong support for the SCOUT clade and its sister relationship to the remaining members of the subfamily, which differs from previous results based on plastid data but agrees with recent results based on nuclear genomic data. However, several instances of highly supported cytonuclear discordance were identified across the Rubioideae phylogeny. Coalescent simulation analysis indicates that while ILS could, by itself, explain the majority of the discordant relationships, plastome introgression may be the better explanation in some cases. Our study further indicates that plastomes across the Rubioideae are, with few exceptions, highly conserved and mainly conform to the structure, gene content, and gene order present in the majority of the flowering plants.

The promise and pitfalls of synteny in phylogenomics.

Reconstructing the tree of life remains a central goal in biology. Early methods, which relied on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the tree of life, 2 recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in 2 or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. This Essay identifies and discusses specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the tree of life.

Emergent ribozyme behaviors in oxychlorine brines indicate a unique niche for molecular evolution on Mars.

Mars is a particularly attractive candidate among known astronomical objects to potentially host life. Results from space exploration missions have provided insights into Martian geochemistry that indicate oxychlorine species, particularly perchlorate, are ubiquitous features of the Martian geochemical landscape. Perchlorate presents potential obstacles for known forms of life due to its toxicity. However, it can also provide potential benefits, such as producing brines by deliquescence, like those thought to exist on present-day Mars. Here we show perchlorate brines support folding and catalysis of functional RNAs, while inactivating representative protein enzymes. Additionally, we show perchlorate and other oxychlorine species enable ribozyme functions, including homeostasis-like regulatory behavior and ribozyme-catalyzed chlorination of organic molecules. We suggest nucleic acids are uniquely well-suited to hypersaline Martian environments. Furthermore, Martian near- or subsurface oxychlorine brines, and brines found in potential lifeforms, could provide a unique niche for biomolecular evolution.

Two telomere-to-telomere gapless genomes reveal insights into Capsicum evolution and capsaicinoid biosynthesis.

Chili pepper (Capsicum) is known for its unique fruit pungency due to the presence of capsaicinoids. The evolutionary history of capsaicinoid biosynthesis and the mechanism of their tissue specificity remain obscure due to the lack of high-quality Capsicum genomes. Here, we report two telomere-to-telomere (T2T) gap-free genomes of C. annuum and its wild nonpungent relative C. rhomboideum to investigate the evolution of fruit pungency in chili peppers. We precisely delineate Capsicum centromeres, which lack high-copy tandem repeats but are extensively invaded by CRM retrotransposons. Through phylogenomic analyses, we estimate the evolutionary timing of capsaicinoid biosynthesis. We reveal disrupted coding and regulatory regions of key biosynthesis genes in nonpungent species. We also find conserved placenta-specific accessible chromatin regions, which likely allow for tissue-specific biosynthetic gene coregulation and capsaicinoid accumulation. These T2T genomic resources will accelerate chili pepper genetic improvement and help to understand Capsicum genome evolution.

Installation of LYRM proteins in early eukaryotes to regulate the metabolic capacity of the emerging mitochondrion.

Core mitochondrial processes such as the electron transport chain, protein translation and the formation of Fe-S clusters (ISC) are of prokaryotic origin and were present in the bacterial ancestor of mitochondria. In animal and fungal models, a family of small Leu-Tyr-Arg motif-containing proteins (LYRMs) uniformly regulates the function of mitochondrial complexes involved in these processes. The action of LYRMs is contingent upon their binding to the acylated form of acyl carrier protein (ACP). This study demonstrates that LYRMs are structurally and evolutionarily related proteins characterized by a core triplet of α-helices. Their widespread distribution across eukaryotes suggests that 12 specialized LYRMs were likely present in the last eukaryotic common ancestor to regulate the assembly and folding of the subunits that are conserved in bacteria but that lack LYRM homologues. The secondary reduction of mitochondria to anoxic environments has rendered the function of LYRMs and their interaction with acylated ACP dispensable. Consequently, these findings strongly suggest that early eukaryotes installed LYRMs in aerobic mitochondria as orchestrated switches, essential for regulating core metabolism and ATP production.

Pervasive relaxed selection in termite genomes.

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

Unveiling the A-to-I mRNA editing machinery and its regulation and evolution in fungi.

A-to-I mRNA editing in animals is mediated by ADARs, but the mechanism underlying sexual stage-specific A-to-I mRNA editing in fungi remains unknown. Here, we show that the eukaryotic tRNA-specific heterodimeric deaminase FgTad2-FgTad3 is responsible for A-to-I mRNA editing in Fusarium graminearum. This editing capacity relies on the interaction between FgTad3 and a sexual stage-specific protein called Ame1. Although Ame1 orthologs are widely distributed in fungi, the interaction originates in Sordariomycetes. We have identified key residues responsible for the FgTad3-Ame1 interaction. The expression and activity of FgTad2-FgTad3 are regulated through alternative promoters, alternative translation initiation, and post-translational modifications. Our study demonstrates that the FgTad2-FgTad3-Ame1 complex can efficiently edit mRNA in yeasts, bacteria, and human cells, with important implications for the development of base editors in therapy and agriculture. Overall, this study uncovers mechanisms, regulation, and evolution of RNA editing in fungi, highlighting the role of protein-protein interactions in modulating deaminase function.

Intraspecific and interspecific variations in the synonymous codon usage in mitochondrial genomes of 8 pleurotus strains.

In this study, we investigated the codon bias of twelve mitochondrial core protein coding genes (PCGs) in eight Pleurotus strains, two of which are from the same species. The results revealed that the codons of all Pleurotus strains had a preference for ending in A/T. Furthermore, the correlation between codon base compositions and codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) indices was also detected, implying the influence of base composition on codon bias. The two P. ostreatus species were found to have differences in various base bias indicators. The average effective number of codons (ENC) of mitochondrial core PCGs of Pleurotus was found to be less than 35, indicating strong codon preference of mitochondrial core PCGs of Pleurotus. The neutrality plot analysis and PR2-Bias plot analysis further suggested that natural selection plays an important role in Pleurotus codon bias. Additionally, six to ten optimal codons (ΔRSCU > 0.08 and RSCU > 1) were identified in eight Pleurotus strains, with UGU and ACU being the most widely used optimal codons in Pleurotus. Finally, based on the combined mitochondrial sequence and RSCU value, the genetic relationship between different Pleurotus strains was deduced, showing large variations between them. This research has improved our understanding of synonymous codon usage characteristics and evolution of this important fungal group.