Statistical methods for detecting molecular adaptation

Positive selection and functional diversification of transcription factor Cmr1 homologs in Alternaria

Transcription factor Cmr1 (Colletotrichum melanin regulation 1) and its homologs in several plant fungal pathogens are the regulators of the 1,8-dihydroxynaphthalene (DHN)-melanin biosynthesis pathway and have evolved functional diversification in morphology and pathogenicity. The fungal genus Alternaria comprises the group of "black fungi" that are rich in DHN-melanin in the primary cell wall and septa of the conidia. Some Alternaria species cause many economically important plant diseases worldwide. However, the evolution and function of Cmr1 homologs in Alternaria remain poorly understood. Here, we identified a total of forty-two Cmr1 homologs from forty-two Alternaria spp. and all contained one additional diverse fungal specific transcription factor motif. Phylogenetic analysis indicated the division of these homologs into five major clades and three branches. Dated phylogeny showed the A and D clades diverged latest and earliest, respectively. Molecular evolutionary analyses revealed that three amino acid sites of Cmr1 homologs in Alternaria were the targets of positive selection. Asmr1, the homolog of Cmr1 in the potato early blight pathogen, Alternaria solani was amplified and displayed the sequence conservation at the amino acid level in different A. solani isolates. Asmr1 was further confirmed to have the transcriptional activation activity and was upregulated during the early stage of potato infection. Deletion of asmr1 led to the decreased melanin content and pathogenicity, deformed conidial morphology, and responses to cell wall and fungicide stresses in A. solani. These results suggest positive selection and functional divergence have played a role in the evolution of Cmr1 homologs in Alternaria. KEY POINTS: • Cmr1 homologs were under positive selection in Alternaria species • Asmr1 is a functional transcription factor, involved in spore development, melanin biosynthesis, pathogenicity, and responses to cell wall and fungicide stresses in A. solani • Cmr1 might be used as a potential taxonomic marker of the genus Alternaria.

Human lncRNAs harbor conserved modules embedded in different sequence contexts

We analyzed the structure of human long non-coding RNA (lncRNAs) genes to investigate whether the non-coding transcriptome is organized in modular domains, as is the case for protein-coding genes. To this aim, we compared all known human lncRNA exons and identified 340 pairs of exons with high sequence and/or secondary structure similarity but embedded in a dissimilar sequence context. We grouped these pairs in 106 clusters based on their reciprocal similarities. These shared modules are highly conserved between humans and the four great ape species, display evidence of purifying selection and likely arose as a result of recent segmental duplications. Our analysis contributes to the understanding of the mechanisms driving the evolution of the non-coding genome and suggests additional strategies towards deciphering the functional complexity of this class of molecules.

Evolutionary analysis of species-specific duplications in flatworm genomes

Platyhelminthes, also known as flatworms, is a phylum of bilaterian invertebrates infamous for their parasitic representatives. The classes Cestoda, Monogenea, and Trematoda comprise parasitic helminths inhabiting multiple hosts, including fishes, humans, and livestock, and are responsible for considerable economic damage and burden on human health. As in other animals, the genomes of flatworms have a wide variety of paralogs, genes related via duplication, whose origins could be mapped throughout the evolution of the phylum. Through in-silico analysis, we studied inparalogs, i.e., species-specific duplications, focusing on their biological functions, expression changes, and evolutionary rate. These genes are thought to be key players in the adaptation process of species to each particular niche. Our results showed that genes related with specific functional terms, such as response to stress, transferase activity, oxidoreductase activity, and peptidases, are overrepresented among inparalogs. This trend is conserved among species from different classes, including free-living species. Available expression data from Schistosoma mansoni, a parasite from the trematode class, demonstrated high conservation of expression patterns between inparalogs, but with notable exceptions, which also display evidence of rapid evolution. We discuss how natural selection may operate to maintain these genes and the particular duplication models that fit better to the observations. Our work supports the critical role of gene duplication in the evolution of flatworms, representing the first study of inparalogs evolution at the genome-wide level in this group.

Signals of positive selection in genomes of palearctic Myotis-bats coexisting with a fungal pathogen

Disease can act as a driving force in shaping genetic makeup across populations, even species, if the impacts influence a particularly sensitive part of their life cycles. White-nose disease is caused by a fungal pathogen infecting bats during hibernation. The mycosis has caused massive population declines of susceptible species in North America, particularly in the genus Myotis. However, Myotis bats appear to tolerate infection in Eurasia, where the fungal pathogen has co-evolved with its bat hosts for an extended period of time. Therefore, with susceptible and tolerant populations, the fungal disease provides a unique opportunity to tease apart factors contributing to tolerance at a genomic level to and gain an understanding of the evolution of non-harmful in host-parasite interactions. To investigate if the fungal disease has caused adaptation on a genomic level in Eurasian bat species, we adopted both whole-genome sequencing approaches and a literature search to compile a set of 300 genes from which to investigate signals of positive selection in genomes of 11 Eurasian bats at the codon-level. Our results indicate significant positive selection in 38 genes, many of which have a marked role in responses to infection. Our findings suggest that white-nose syndrome may have applied a significant selective pressure on Eurasian Myotis-bats in the past, which can contribute their survival in co-existence with the pathogen. Our findings provide an insight on the selective pressure pathogens afflict on their hosts using methodology that can be adapted to other host-pathogen study systems.

Teosinte Pollen Drive guides maize diversification and domestication by RNAi

Selfish genetic elements contribute to hybrid incompatibility and bias or 'drive' their own transmission1,2. Chromosomal drive typically functions in asymmetric female meiosis, whereas gene drive is normally post-meiotic and typically found in males. Here, using single-molecule and single-pollen genome sequencing, we describe Teosinte Pollen Drive, an instance of gene drive in hybrids between maize (Zea mays ssp. mays) and teosinte mexicana (Z. mays ssp. mexicana) that depends on RNA interference (RNAi). 22-nucleotide small RNAs from a non-coding RNA hairpin in mexicana depend on Dicer-like 2 (Dcl2) and target Teosinte Drive Responder 1 (Tdr1), which encodes a lipase required for pollen viability. Dcl2, Tdr1 and the hairpin are in tight pseudolinkage on chromosome 5, but only when transmitted through the male. Introgression of mexicana into early cultivated maize is thought to have been critical to its geographical dispersal throughout the Americas3, and a tightly linked inversion in mexicana spans a major domestication sweep in modern maize4. A survey of maize traditional varieties and sympatric populations of teosinte mexicana reveals correlated patterns of admixture among unlinked genes required for RNAi on at least four chromosomes that are also subject to gene drive in pollen from synthetic hybrids. Teosinte Pollen Drive probably had a major role in maize domestication and diversification, and offers an explanation for the widespread abundance of 'self' small RNAs in the germ lines of plants and animals.

Diversification of Toll-like receptor 1 in swamp eel (Monopterus albus)

Toll-like receptor 1 (TLR1) is a pattern recognition receptor that plays critical roles in triggering immune activation via detecting bacterial lipoproteins and lipopeptides. In this study, the genetic characteristic of TLR1 was studied for an important aquaculture fish, swamp eel Monopterus albus. The eel has been seriously threatened by infectious diseases. However, a low level of genetic heterogeneity in the fish that has resulted from a demographic bottleneck presents further challenges in breeding for disease resistance. A comparison with the homologue of closely related species M. javanensis revealed that amino acid replacement (nonsynonymous) but not silent (synonymous) differences have accumulated nonrandomly over the coding sequences of the receptors at the early stage of their phylogenetic split. The combined results from comparative analyses of nonsynonymous-to-synonymous polymorphisms showed that the receptor has undergone significant diversification in M. albus driven by adaptive selection likely after the genetic bottleneck. Some of the changes reported here have taken place in the structures mediating heterodimerization with co-receptor TLR2, ligand recognition, and/or formation of active signaling complex with adaptor, which highlighted key structural elements and strategies of TLR1 in arms race against exogenous challenges. The findings of this study will add to the knowledge base of genetic engineering and breeding for disease resistance in the eel.

Genome-wide analysis of the bZIP gene family in Cinnamomum camphora ('Gantong 1') reveals the putative function in anthocyanin biosynthesis

Basic leucine zipper (bZIP) transcription factors (TFs) regulate plant development, growth, and secondary metabolism. The formation of red bark of new ornamental cultivar 'Gantong 1' is regulated mainly by anthocyanin anabolism. However, it is unclear whether and which bZIP TFs are involved in this process. We identified 89 genes encoding CcbZIP TFs in Cinnamomum camphora genome that could be divided into 14 subfamilies with similar gene structures and conserved motifs. CcbZIP38 and CcbZIP57 were highly conserved compared to HY5 in Arabidopsis thaliana and they were highly expressed in the bark and leaves of 'Gantong 1' at different stages. The target gene enrichment analysis showed that indicating indirect involvement of CcbZIP38 and CcbZIP57 in the regulation of anthocyanin synthesis. Our study contributes to understanding the molecular mechanism of anthocyanin synthesis regulation by CcbZIP TFs and provides a theoretical basis for genetic improvement of ornamental traits in C. camphora.

Evolutionary Dynamics of Accelerated Antiviral Resistance Development in Hypermutator Herpesvirus

Antiviral therapy is constantly challenged by the emergence of resistant pathogens. At the same time, experimental approaches to understand and predict resistance are limited by long periods required for evolutionary processes. Here, we present a herpes simplex virus 1 mutant with impaired proofreading capacity and consequently elevated mutation rates. Comparing this hypermutator to parental wild type virus, we study the evolution of antiviral drug resistance in vitro. We model resistance development and elucidate underlying genetic changes against three antiviral substances. Our analyzes reveal no principle difference in the evolutionary behavior of both viruses, adaptive processes are overall similar, however significantly accelerated for the hypermutator. We conclude that hypermutator viruses are useful for modeling adaptation to antiviral therapy. They offer the benefit of expedited adaptation without introducing apparent bias and can therefore serve as an accelerator to predict natural evolution.

The First Three Mitochondrial Genomes for the Characterization of the Genus Egeirotrioza (Hemiptera: Triozidae) and Phylogenetic Implications

(1) Background: Mitochondrial genomes are important markers for the study of phylogenetics and systematics. Triozidae includes some primary pests of Populus euphratica. The phylogenetic relationships of this group remain controversial due to the lack of molecular data. (2) Methods: Mitochondria of Egeirotrioza Boselli were sequenced and assembled. We analyzed the sequence length, nucleotide composition, and evolutionary rate of Triozidae, combined with the 13 published mitochondrial genomes. (3) Results: The evolutionary rate of protein-coding genes was as follows: ATP8 > ND6 > ND5 > ND2 > ND4 > ND4L > ND1 > ND3 > APT6 > CYTB > COX3 > COX2 > COX1. We reconstructed the phylogenetic relationships of Triozidae based on 16 triozid mitochondrial genomes (thirteen ingroups and three outgroups) using the maximum likelihood (ML) and Bayesian inference (BI) approaches. The phylogenetic analysis of the 16 Triozidae mitochondrial genomes showed that Egeirotrioza was closely related to Leptynoptera. (4) Conclusions: We have identified 13 PCGs, 22 tRNAs, 2 rRNAs, and 1 control region (CR) of all newly sequenced mitochondrial genomes, which were the mitochondrial gene type in animals. The results of this study provide valuable genomic information for the study of psyllid species.

Methylglyoxal detoxifying gene families in tomato: Genome-wide identification, evolution, functional prediction, and transcript profiling

Methylglyoxal (MG) is a highly cytotoxic molecule produced in all biological systems, which could be converted into non-toxic D-lactate by an evolutionarily conserved glyoxalase pathway. Glutathione-dependent glyoxalase I (GLYI) and glyoxalase II (GLYII) are responsible for the detoxification of MG into D-lactate in sequential reactions, while DJ-1 domain containing glyoxalase III (GLYIII) catalyzes the same reaction in a single step without glutathione dependency. Afterwards, D-lactate dehydrogenase (D-LDH) converts D-lactate into pyruvate, a metabolically usable intermediate. In the study, a comprehensive genome-wide investigation has been performed in one of the important vegetable plants, tomato to identify 13 putative GLYI, 4 GLYII, 3 GLYIII (DJ-1), and 4 D-LDH genes. Expression pattern analysis using microarray data confirmed their ubiquitous presence in different tissues and developmental stages. Moreover, stress treatment of tomato seedlings and subsequent qRT-PCR demonstrated upregulation of SlGLYI-2, SlGLYI-3, SlGLYI-6A, SlGLYII-1A, SlGLYII-3B, SlDJ-1A, SlDLDH-1 and SlDLDH-4 in response to different abiotic stresses, whereas SlGLYI-6B, SlGLYII-1B, SlGLYII-3A, SlDJ-1D and SlDLDH-2 were downregulated. Expression data also revealed SlGLYII-1B, SlGLYI-1A, SlGLYI-2, SlDJ-1D, and SlDLDH-4 were upregulated in response to various pathogenic infections, indicating the role of MG detoxifying enzymes in both plant defence and stress modulation. The functional characterization of each of these members could lay the foundation for the development of stress and disease-resistant plants promoting sustainable agriculture and production.

Evolutionary and phylogenetic insights from the mitochondrial genomic analysis of Diceraeus melacanthus and D. furcatus (Hemiptera: Pentatomidae)

The mitochondrial genomes of D. melacanthus and D. furcatus were sequenced and used to investigate the phylogenetic relationships with 54 species of Pentatomidae. Their mitogenomes are 17,197 and 15,444 bp-long, respectively, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22/21 transfer RNA genes, with conserved gene arrangement. Leu, Lys, and Ser were the most common amino acids in their PCGs. PCGs evolutionary analysis indicated their mitogenomes are under purifying selection, and the most conserved genes are from the cytochrome complex, reinforcing their suitability as markers for molecular taxonomy. We identified 490 mtSSRs in 56 Pentatomidae species, with large variation and a positive correlation between mtSSR number and genome size. Three mtSSRs were identified in each Diceraeus species. Only the mtSSR in the nad6 (D. melacanthus) and nad4 (D. furcatus) appear to have application as molecular markers for species characterization. Phylogenetic analysis confirmed the monophyly of Pentatomidae. However, our analysis challenged the monophyly of Pentatominae and Podopinae. We also detected unexpected relationships among some tribes and genera, highlighting the complexity of the internal taxonomic structure of Pentatomidae. Both Diceraeus species were grouped in the same clade with the remaining Carpocorini analyzed.

Genetic Variability Among U.S.-Sentinel Cotton Plot Cotton Leafroll Dwarf Virus and Globally Available Reference Isolates Based on ORF0 Diversity

The aphid-transmitted polerovirus, cotton leafroll dwarf virus (CLRDV), first characterized from symptomatic cotton plants in South America, has been identified in commercial cotton plantings in the United States. Here, the CLRDV intraspecific diversity was investigated by comparative sequence analysis of the most divergent CLRDV coding region, ORF0/P0. Bayesian analysis of ORF0 sequences for U.S. and reference populations resolved three well-supported sister clades comprising one U.S. and two South American lineages. Principal component analysis (PCA) identified seven statistically supported intraspecific populations. The Bayesian phylogeny and PCA dendrogram-inferred relationships were congruent. Population analysis of ORF0 sequences indicated most lineages have evolved under negative selection, albeit certain sites/isolates evolved under positive selection. Both U.S. and South American isolates exhibited extensive ORF0 diversity. At least two U.S. invasion foci were associated with their founder populations in Alabama-Georgia and eastern Texas. The Alabama-Georgia founder is implicated as the source of recent widespread expansion and establishment of secondary disease foci throughout the southeastern-central United States. Based on the geographically restricted distribution, spread of another extant Texas population appeared impeded by a population bottleneck. Extant CLRDV isolates represent several putative introductions potentially associated with catastrophic weather events dispersing viruliferous cotton aphids of unknown origin(s).

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.

Extrinsic and intrinsic drivers of natural killer cell clonality

Clonal expansion of antigen-specific lymphocytes is the fundamental mechanism enabling potent adaptive immune responses and the generation of immune memory. Accompanied by pronounced epigenetic remodeling, the massive proliferation of individual cells generates a critical mass of effectors for the control of acute infections, as well as a pool of memory cells protecting against future pathogen encounters. Classically associated with the adaptive immune system, recent work has demonstrated that innate immune memory to human cytomegalovirus (CMV) infection is stably maintained as large clonal expansions of natural killer (NK) cells, raising questions on the mechanisms for clonal selection and expansion in the absence of re-arranged antigen receptors. Here, we discuss clonal NK cell memory in the context of the mechanisms underlying clonal competition of adaptive lymphocytes and propose alternative selection mechanisms that might decide on the clonal success of their innate counterparts. We propose that the integration of external cues with cell-intrinsic sources of heterogeneity, such as variegated receptor expression, transcriptional states, and somatic variants, compose a bottleneck for clonal selection, contributing to the large size of memory NK cell clones.

A survey of the Sli gene in wild and cultivated potato

Inbred-hybrid breeding of diploid potatoes necessitates breeding lines that are self-compatible. One way of incorporating self-compatibility into incompatible cultivated potato (Solanum tuberosum) germplasm is to introduce the S-locus inhibitor gene (Sli), which functions as a dominant inhibitor of gametophytic self-incompatibility. To learn more about Sli diversity and function in wild species relatives of cultivated potato, we obtained Sli gene sequences that extended from the 5'UTR to the 3'UTR from 133 individuals from 22 wild species relatives of potato and eight diverse cultivated potato clones. DNA sequence alignment and phylogenetic trees based on genomic and protein sequences show that there are two highly conserved groups of Sli sequences. DNA sequences in one group contain the 533 bp insertion upstream of the start codon identified previously in self-compatible potato. The second group lacks the insertion. Three diploid and four polyploid individuals of wild species collected from geographically disjointed localities contained Sli with the 533 bp insertion. For most of the wild species clones examined, however, Sli did not have the insertion. Phylogenetic analysis indicated that Sli sequences with the insertion, in wild species and in cultivated clones, trace back to a single origin. Some diploid wild potatoes that have Sli with the insertion were self-incompatible and some wild potatoes that lack the insertion were self-compatible. Although there is evidence of positive selection for some codon positions in Sli, there is no evidence of diversifying selection at the gene level. In silico analysis of Sli protein structure did not support the hypothesis that amino acid changes from wild-type (no insertion) to insertion-type account for changes in protein function. Our study demonstrated that genetic factors besides the Sli gene must be important for conditioning a switch in the mating system from self-incompatible to self-compatible in wild potatoes.

Evolutionary analysis reveals the role of a non-catalytic domain of peptidyl arginine deiminase 2 in transcriptional regulation

Peptidyl arginine deiminases (PADIs) catalyze protein citrullination, a post-translational conversion of arginine to citrulline. The most widely expressed member of this family, PADI2, regulates cellular processes that impact several diseases. We hypothesized that we could gain new insights into PADI2 function through a systematic evolutionary and structural analysis. Here, we identify 20 positively selected PADI2 residues, 16 of which are structurally exposed and maintain PADI2 interactions with cognate proteins. Many of these selected residues reside in non-catalytic regions of PADI2. We validate the importance of a prominent loop in the middle domain that encompasses PADI2 L162, a residue under positive selection. This site is essential for interaction with the transcription elongation factor (P-TEFb) and mediates the active transcription of the oncogenes c-MYC, and CCNB1, as well as impacting cellular proliferation. These insights could be key to understanding and addressing the role of the PADI2 c-MYC axis in cancer progression.

An evolutionary timeline of the oxytocin signaling pathway

Oxytocin is a neuropeptide associated with both psychological and somatic processes like parturition and social bonding. Although oxytocin homologs have been identified in many species, the evolutionary timeline of the entire oxytocin signaling gene pathway has yet to be described. Using protein sequence similarity searches, microsynteny, and phylostratigraphy, we assigned the genes supporting the oxytocin pathway to different phylostrata based on when we found they likely arose in evolution. We show that the majority (64%) of genes in the pathway are 'modern'. Most of the modern genes evolved around the emergence of vertebrates or jawed vertebrates (540 - 530 million years ago, 'mya'), including OXTR, OXT and CD38. Of those, 45% were under positive selection at some point during vertebrate evolution. We also found that 18% of the genes in the oxytocin pathway are 'ancient', meaning their emergence dates back to cellular organisms and opisthokonta (3500-1100 mya). The remaining genes (18%) that evolved after ancient and before modern genes were classified as 'medium-aged'. Functional analyses revealed that, in humans, medium-aged oxytocin pathway genes are highly expressed in contractile organs, while modern genes in the oxytocin pathway are primarily expressed in the brain and muscle tissue.

Evolutionary genetics of pulmonary anatomical adaptations in deep-diving cetaceans

BACKGROUND

Cetaceans, having experienced prolonged adaptation to aquatic environments, have undergone evolutionary changes in their respiratory systems. This process of evolution has resulted in the emergence of distinctive phenotypic traits, notably the abundance of elastic fibers and thickened alveolar walls in their lungs, which may facilitate alveolar collapse during diving. This structure helps selective exchange of oxygen and carbon dioxide, while minimizing nitrogen exchange, thereby reducing the risk of DCS. Nevertheless, the scientific inquiry into the mechanisms through which these unique phenotypic characteristics govern the diving behavior of marine mammals, including cetaceans, remains unresolved.

RESULTS

This study entails an evolutionary analysis of 42 genes associated with pulmonary fibrosis across 45 mammalian species. Twenty-one genes in cetaceans exhibited accelerated evolution, featuring specific amino acid substitutions in 14 of them. Primarily linked to the development of the respiratory system and lung morphological construction, these genes play a crucial role. Moreover, among marine mammals, we identified eight genes undergoing positive selection, and the evolutionary rates of three genes significantly correlated with diving depth. Specifically, the SFTPC gene exhibited convergent amino acid substitutions. Through in vitro cellular experiments, we illustrated that convergent amino acid site mutations in SFTPC contribute positively to pulmonary fibrosis in marine mammals, and the presence of this phenotype can induce deep alveolar collapse during diving, thereby reducing the risk of DCS during diving.

CONCLUSIONS

The study unveils pivotal genetic signals in cetaceans and other marine mammals, arising through evolution. These genetic signals may influence lung characteristics in marine mammals and have been linked to a reduced risk of developing DCS. Moreover, the research serves as a valuable reference for delving deeper into human diving physiology.

Properties and Mechanisms of Deletions, Insertions, and Substitutions in the Evolutionary History of SARS-CoV-2

SARS-CoV-2 has accumulated many mutations since its emergence in late 2019. Nucleotide substitutions leading to amino acid replacements constitute the primary material for natural selection. Insertions, deletions, and substitutions appear to be critical for coronavirus's macro- and microevolution. Understanding the molecular mechanisms of mutations in the mutational hotspots (positions, loci with recurrent mutations, and nucleotide context) is important for disentangling roles of mutagenesis and selection. In the SARS-CoV-2 genome, deletions and insertions are frequently associated with repetitive sequences, whereas C>U substitutions are often surrounded by nucleotides resembling the APOBEC mutable motifs. We describe various approaches to mutation spectra analyses, including the context features of RNAs that are likely to be involved in the generation of recurrent mutations. We also discuss the interplay between mutations and natural selection as a complex evolutionary trend. The substantial variability and complexity of pipelines for the reconstruction of mutations and the huge number of genomic sequences are major problems for the analyses of mutations in the SARS-CoV-2 genome. As a solution, we advocate for the development of a centralized database of predicted mutations, which needs to be updated on a regular basis.

Genetic mechanism of body size variation in groupers: Insights from phylotranscriptomics

Animal body size variation is of particular interest in evolutionary biology, but the genetic basis remains largely unknown. Previous studies have shown the presence of two parallel evolutionary genetic clusters within the fish genus Epinephelus with evident divergence in body size, providing an excellent opportunity to investigate the genetic basis of body size variation in vertebrates. Herein, we performed phylotranscriptomic analysis and reconstructed the phylogeny of 13 epinephelids originating from the South China Sea. Two genetic clades with an estimated divergence time of approximately 15.4 million years ago were correlated with large and small body size, respectively. A total of 180 rapidly evolving genes and two positively selected genes were identified between the two groups. Functional enrichment analyses of these candidate genes revealed distinct enrichment categories between the two groups. These pathways and genes may play important roles in body size variation in groupers through complex regulatory networks. Based on our results, we speculate that the ancestors of the two divergent groups of groupers may have adapted to different environments through habitat selection, leading to genetic variations in metabolic patterns, organ development, and lifespan, resulting in body size divergence between the two locally adapted populations. These findings provide important insights into the genetic mechanisms underlying body size variation in groupers and species differentiation.

Increasing intra- and inter-subtype HIV diversity despite declining HIV incidence in Uganda

HIV incidence has been declining in Africa with scale-up of HIV interventions. However, there is limited data on HIV evolutionary trends in African populations with waning epidemics. We evaluated changes in HIV viral diversity and genetic divergence in southern Uganda over a twenty-five-year period spanning the introduction and scale-up of HIV prevention and treatment programs using HIV sequence and survey data from the Rakai Community Cohort Study, an open longitudinal population-based HIV surveillance cohort. Gag (p24) and env (gp41) HIV data were generated from persons living with HIV (PLHIV) in 31 inland semi-urban trading and agrarian communities (1994 to 2018) and four hyperendemic Lake Victoria fishing communities (2011 to 2018) under continuous surveillance. HIV subtype was assigned using the Recombination Identification Program with phylogenetic confirmation. Inter-subtype diversity was estimated using the Shannon diversity index and intra-subtype diversity with the nucleotide diversity and pairwise TN93 genetic distance. Genetic divergence was measured using root-to-tip distance and pairwise TN93 genetic distance analyses. Evolutionary dynamics were assessed among demographic and behavioral sub-groups, including by migration status. 9,931 HIV sequences were available from 4,999 PLHIV, including 3,060 and 1,939 persons residing in inland and fishing communities, respectively. In inland communities, subtype A1 viruses proportionately increased from 14.3% in 1995 to 25.9% in 2017 (p<0.001), while those of subtype D declined from 73.2% in 1995 to 28.2% in 2017 (p<0.001). The proportion of viruses classified as recombinants significantly increased by more than four-fold. Inter-subtype HIV diversity has generally increased. While p24 intra-subtype genetic diversity and divergence leveled off after 2014, diversity and divergence of gp41 increased through 2017. Inter- and intra-subtype viral diversity increased across all population sub-groups, including among individuals with no recent migration history or extra-community sexual partners. This study provides insights into population-level HIV evolutionary dynamics in declining African HIV epidemics following the scale-up of HIV prevention and treatment programs. Continued molecular surveillance may provide a better understanding of the dynamics driving population HIV evolution and yield important insights for epidemic control and vaccine development.

Quantitative pathogenicity and host adaptation in a fungal plant pathogen revealed by whole-genome sequencing

Knowledge of genetic determinism and evolutionary dynamics mediating host-pathogen interactions is essential to manage fungal plant diseases. Studies on the genetic architecture of fungal pathogenicity often focus on large-effect effector genes triggering strong, qualitative resistance. It is not clear how this translates to predominately quantitative interactions. Here, we use the Zymoseptoria tritici-wheat model to elucidate the genetic architecture of quantitative pathogenicity and mechanisms mediating host adaptation. With a multi-host genome-wide association study, we identify 19 high-confidence candidate genes associated with quantitative pathogenicity. Analysis of genetic diversity reveals that sequence polymorphism is the main evolutionary process mediating differences in quantitative pathogenicity, a process that is likely facilitated by genetic recombination and transposable element dynamics. Finally, we use functional approaches to confirm the role of an effector-like gene and a methyltransferase in phenotypic variation. This study highlights the complex genetic architecture of quantitative pathogenicity, extensive diversifying selection and plausible mechanisms facilitating pathogen adaptation.

Combining extracellular matrix proteome and phosphoproteome of chickpea and meta-analysis reveal novel proteoforms and evolutionary significance of clade-specific wall-associated events in plant

Extracellular matrix (ECM) plays central roles in cell architecture, innate defense and cell wall integrity (CWI) signaling. During transition to multicellularity, modular domain structures of ECM proteins and proteoforms have evolved due to continuous adaptation across taxonomic clades under different ecological niche. Although this incredible diversity has to some extent been investigated at protein level, extracellular phosphorylation events and molecular evolution of ECM proteoform families remains unexplored. We developed matrisome proteoform atlas in a grain legume, chickpea and performed meta-analyses of 74 plant matrisomes. MS/MS analysis identified 1,424 proteins and 315 phosphoproteins involved in diverse functions. Cross-species ECM protein network identified proteoforms associated with CWI maintenance system. Phylogenetic characterization of eighteen matrix protein families highlighted the role of taxon-specific paralogs and orthologs. Novel information was acquired on gene expansion and loss, co-divergence, sub functionalization and neofunctionalization during evolution. Modular networks of matrix protein families and hub proteins showed higher diversity across taxonomic clades than among organs. Furthermore, protein families differ in nonsynonymous to synonymous substitution rates. Our study pointed towards the matrix proteoform functionality, sequence divergence variation, interactions between wall remodelers and molecular evolution using a phylogenetic framework. This is the first report on comprehensive matrisome proteoform network illustrating presence of CWI signaling proteins in land plants.

Overexpression and nonsynonymous mutations of UDP-glycosyltransferases are potentially associated with pyrethroid resistance in Anopheles funestus

UDP-glycosyltransferases (UGTs) enzymes are pivotal in insecticide resistance by transforming hydrophobic substrates into more hydrophilic forms for efficient cell elimination. This study provides the first comprehensive investigation of Anopheles funestus UGT genes, their evolution, and their association with pyrethroid resistance. We employed a genome-wide association study using pooled sequencing (GWAS-PoolSeq) and transcriptomics on pyrethroid-resistant An. funestus, along with deep-targeted sequencing of UGTs in 80 mosquitoes Africa-wide. UGT310B2 was consistently overexpressed Africa-wide and significant gene-wise Fst differentiation was observed between resistant and susceptible populations: UGT301C2 and UGT302A3 in Malawi, and UGT306C2 in Uganda. Additionally, nonsynonymous mutations in UGT genes were identified. Gene-wise Tajima's D density curves provide insights into population structures within populations across these countries, supporting previous observations. These findings have important implications for current An. funestus control strategies facilitating the prediction of cross-resistance to other UGT-metabolised polar insecticides, thereby guiding more effective and targeted insecticide resistance management efforts.

Comparative plastome analysis and taxonomic classification of snow lotus species (Saussurea, Asteraceae) in Central Asia and Southern Siberia

Four species of Saussurea, namely S. involucrata, S. orgaadayi, S. bogedaensis, and S. dorogostaiskii, are known as the "snow lotus," which are used as traditional medicines in China (Xinjiang), Kyrgyzstan, Mongolia, and Russia (Southern Siberia). These species are threatened globally, because of illegal harvesting and climate change. Furthermore, the taxonomic classification and identification of these threatened species remain unclear owing to limited research. The misidentification of medicinal species can sometimes be harmful to health. Therefore, the phylogenetic and genomic features of these species need to be confirmed. In this study, we sequenced five complete chloroplast genomes and seven nuclear ITS regions of four snow lotus species and other Saussurea species. We further explored their genetic variety, selective pressure at the sequence level, and phylogenetic relationships using the chloroplast genome, nuclear partial DNA sequences, and morphological features. Plastome of the snow lotus species has a conserved structure and gene content similar to most Saussurea species. Two intergenic regions (ndhJ-ndhK and ndhD-psaC) show significantly high diversity among chloroplast regions. Thus, ITS and these markers are suitable for identifying snow lotus species. In addition, we characterized 43 simple sequence repeats that may be useful in future population genetic studies. Analysis of the selection signatures identified three genes (rpoA, ndhB, and ycf2) that underwent positive selection. These genes may play important roles in the adaptation of the snow lotus species to alpine environments. S. dorogostaiskii is close to S. baicalensis and exhibits slightly different adaptation from others. The taxonomic position of the snow lotus species, confirmed by morphological and molecular evidence, is as follows: (i) S. involucrata has been excluded from the Mongolian flora due to misidentification as S. orgaadayi or S. bogedaensis for a long time; (ii) S. dorogostaiskii belongs to section Pycnocephala subgenus Saussurea, whereas other the snow lotus species belong to section Amphilaena subgenus Amphilaena; and (iii) S. krasnoborovii is synonymous of S. dorogostaiskii. This study clarified the speciation and lineage diversification of the snow lotus species in Central Asia and Southern Siberia.

Evolutionary Analysis of Six Gene Families Part of the Reactive Oxygen Species (ROS) Gene Network in Three Brassicaceae Species

Climate change is expected to intensify the occurrence of abiotic stress in plants, such as hypoxia and salt stresses, leading to the production of reactive oxygen species (ROS), which need to be effectively managed by various oxido-reductases encoded by the so-called ROS gene network. Here, we studied six oxido-reductases families in three Brassicaceae species, Arabidopsis thaliana as well as Nasturtium officinale and Eutrema salsugineum, which are adapted to hypoxia and salt stress, respectively. Using available and new genomic data, we performed a phylogenomic analysis and compared RNA-seq data to study genomic and transcriptomic adaptations. This comprehensive approach allowed for the gaining of insights into the impact of the adaptation to saline or hypoxia conditions on genome organization (gene gains and losses) and transcriptional regulation. Notably, the comparison of the N. officinale and E. salsugineum genomes to that of A. thaliana highlighted changes in the distribution of ohnologs and homologs, particularly affecting class III peroxidase genes (CIII Prxs). These changes were specific to each gene, to gene families subjected to duplication events and to each species, suggesting distinct evolutionary responses. The analysis of transcriptomic data has allowed for the identification of genes related to stress responses in A. thaliana, and, conversely, to adaptation in N. officinale and E. salsugineum.

Mitogenomic Characterization and Phylogenetic Placement of African Hind, Cephalopholis taeniops: Shedding Light on the Evolution of Groupers (Serranidae: Epinephelinae)

The global exploration of evolutionary trends in groupers, based on mitogenomes, is currently underway. This research extensively investigates the structure of and variations in Cephalopholis species mitogenomes, along with their phylogenetic relationships, focusing specifically on Cephalopholis taeniops from the Eastern Atlantic Ocean. The generated mitogenome spans 16,572 base pairs and exhibits a gene order analogous to that of the ancestral teleost's, featuring 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an AT-rich control region. The mitogenome of C. taeniops displays an AT bias (54.99%), aligning with related species. The majority of PCGs in the mitogenome initiate with the start codon ATG, with the exceptions being COI (GTG) and atp6 (TTG). The relative synonymous codon usage analysis revealed the maximum abundance of leucine, proline, serine, and threonine. The nonsynonymous/synonymous ratios were <1, which indicates a strong negative selection among all PCGs of the Cephalopholis species. In C. taeniops, the prevalent transfer RNAs display conventional cloverleaf secondary structures, except for tRNA-serine (GCT), which lacks a dihydrouracil (DHU) stem. A comparative examination of conserved domains and sequence blocks across various Cephalopholis species indicates noteworthy variations in length and nucleotide diversity. Maximum likelihood, neighbor-joining, and Bayesian phylogenetic analyses, employing the concatenated PCGs and a combination of PCGs + rRNAs, distinctly separate all Cephalopholis species, including C. taeniops. Overall, these findings deepen our understanding of evolutionary relationships among serranid groupers, emphasizing the significance of structural considerations in mitogenomic analyses.

Evolutionary analysis and biological characterization of a novel alphabaculovirus isolated from Mythimna separata

Baculoviruses are insect-specific pathogens. Novel baculovirus isolates provide new options for the biological control of pests. Therefore, research into the biological characteristics of newly isolated baculoviruses, including accurate classification and nomenclature, is important. In this study, a baculovirus was isolated from Mythimna separata and its complete genome sequence was determined by next-generation sequencing. The double-stranded DNA genome was 153 882 bp in length, encoding 163 open reading frames. The virus was identified as a variant of Mamestra brassicae multiple nucleopolyhedrovirus (MbMNPV) and designated Mamestra brassicae multiple nucleopolyhedrovirus CHN1 (MbMNPV-CHN1) according to ultrastructural analysis, genome comparison and phylogenetic analysis. Phylogenetic inference placed MbMNPV-CHN1 in a clade containing isolates of MacoNPV-A, MacoNPV-B and MbMNPV, which we have designated the Mb-McNPV group. The genomes of isolates in the Mb-McNPV group exhibited a high degree of collinearity with relatively minor differences in the content of annotated open reading frames. The development of codon usage bias in the Mb-McNPV group was affected mainly by natural selection. MbMNPV-CHN1 shows high infectivity against seven species of Lepidoptera. The yield of MbMNPV-CHN1 in the fourth- and fifth-instar M. separata larvae was 6.25×109-1.23×1010 OBs/cadaver. Our data provide insights into the classification, host range and virulence differences among baculoviruses of the Mb-McNPV group, as well as a promising potential new baculoviral insecticide.

A dynamical limit to evolutionary adaptation

Natural selection makes evolutionary adaptation possible even if the overwhelming majority of new mutations are deleterious. However, in rapidly evolving populations where numerous linked mutations occur and segregate simultaneously, clonal interference and genetic hitchhiking can limit the efficiency of selection, allowing deleterious mutations to accumulate over time. This can in principle overwhelm the fitness increases provided by beneficial mutations, leading to an overall fitness decline. Here, we analyze the conditions under which evolution will tend to drive populations to higher versus lower fitness. Our analysis focuses on quantifying the boundary between these two regimes, as a function of parameters such as population size, mutation rates, and selection pressures. This boundary represents a state in which adaptation is precisely balanced by Muller's ratchet, and we show that it can be characterized by rapid molecular evolution without any net fitness change. Finally, we consider the implications of global fitness-mediated epistasis and find that under some circumstances, this can drive populations toward the boundary state, which can thus represent a long-term evolutionary attractor.

Genome-wide in-silico analysis of ethylene biosynthesis gene family in Musa acuminata L. and their response under nutrient stress

Ethylene is a gaseous phytohormone involved in plants' growth and developmental processes, including seed germination, root initiation, fruit ripening, flower and leaf senescence, abscission, and stress responses. Ethylene biosynthesis (EB) gene analysis in response to nitrogen (N) and potassium (K) stress has not yet been conducted in Musa acuminata (banana) roots. The genome mining of banana (Musa acuminata L.) revealed 14 putative 1-aminocyclopropane-1-carboxylate synthase (ACS), 10 1-aminocyclopropane-1-carboxylate oxidase (ACO), and 3 Ethylene overproducer 1 (ETO1) genes. ACS, ACO, and ETO1 proteins possessed amino acid residues ranging from 422-684, 636-2670, and 893-969, respectively, with molecular weight (Mw) ranging from 4.93-7.55 kD, 10.1-8.3 kD and 10.1-10.78 kD. The number of introns present in ACS, ACO, and ETO1 gene sequences ranges from 0-14, 1-6, and 0-6, respectively. The cis-regulatory element analysis revealed the presence of light-responsive, abscisic acid, seed regulation, auxin-responsive, gibberellin element, endosperm-specific, anoxic inducibility, low-temperature responsiveness, salicylic acid responsiveness, meristem-specific and stress-responsive elements. Comprehensive phylogenetic analyses ACS, ACO, and ETO1 genes of Banana with Arabidopsis thaliana revealed several orthologs and paralogs assisting in understanding the putative functions of these genes. The expression profile of Musa acuminata genes in root under normal and low levels of nitrogen and potassium shows that MaACS14 and MaACO6 expressed highly at normal nitrogen supply. MaACS1 expression was significantly upregulated at low potassium levels, whereas, MaACO6 gene expression was significantly downregulated. The functional divergence and site-specific selective pressures on specific gene sequences of banana have been investigated. The bioinformatics-based genome-wide assessment of the family of banana attempted in the present study could be a significant step for deciphering novel ACS, ACO, and ETO1 genes based on genome-wide expression profiling.

Comparative Evolutionary Genomics in Insects

Genome sequencing quality, in terms of both read length and accuracy, is constantly improving. By combining long-read sequencing technologies with various scaffolding techniques, chromosome-level genome assemblies are now achievable at an affordable price for non-model organisms. Insects represent an exciting taxon for studying the genomic underpinnings of evolutionary innovations, due to ancient origins, immense species-richness, and broad phenotypic diversity. Here we summarize some of the most important methods for carrying out a comparative genomics study on insects. We describe available tools and offer concrete tips on all stages of such an endeavor from DNA extraction through genome sequencing, annotation, and several evolutionary analyses. Along the way we describe important insect-specific aspects, such as DNA extraction difficulties or gene families that are particularly difficult to annotate, and offer solutions. We describe results from several examples of comparative genomics analyses on insects to illustrate the fascinating questions that can now be addressed in this new age of genomics research.

The Complete Mitochondrial Genome and Phylogenetic Analysis of the Freshwater Shellfish Novaculina chinensis (Bivalvia: Pharidae)

Razor clams, belonging to the Pharidae and Solenidae families, are ecologically and economically important; however, very little research has been conducted on the Pharidae family. The genus Novaculina is a marine-derived freshwater lineage, and Novaculina chinensis is a rare freshwater species of the Pharidae family. In order to understand the phylogenetic relationships of N. chinensis, we sequenced the mitochondrial genome of the genus Novaculina, which is 16,262 bp in length and consists of 12 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 ribosomal RNA genes (rRNAs). The phylogenetic relationships of 69 Imparidentian mitochondrial genomes (mitogenomes) indicated that N. chineisis is closely related to Sinonovacula constricta of the order Adapedonta. Our study also found that the Ka/Ks ratios of 12 protein-coding genes in the Pharidae family are lower than one, indicating the occurrence of negative purification selection. Morphological observations of the siphons of N. chinensis, Novaculina myanmarensis, and Novaculina gangetica indicate that N. chinensis may be the ancestral clade of the genus Novaculina, which has not been proposed in previous studies. Our study provides useful molecular information on the phylogenetic and evolutionary relationships of Pharidae and also contributes to the conservation and management of the germplasm resources of N. chinensis.

Characterization of the Complete Mitochondrial Genome of the Elongate Loach and Its Phylogenetic Implications in Cobitidae

The elongate loach is an endemic fish in China. Previous studies have provided some insights into the mitochondrial genome composition and the phylogenetic relationships of the elongate loach inferred using protein-coding genes (PCGs), yet detailed information about it remains limited. Therefore, in this study we sequenced the complete mitochondrial genome of the elongate loach and analyzed its structural characteristics. The PCGs and mitochondrial genome were used for selective stress analysis and genomic comparative analysis. The complete mitochondrial genome of the elongate loach, together with those of 35 Cyprinidae species, was used to infer the phylogenetic relationships of the Cobitidae family through maximum likelihood (ML) reconstruction. The results showed that the genome sequence has a full length of 16,591 bp, which includes 13 PCGs, 22 transfer RNA genes (tRNA), 2 ribosomal RNA genes (rRNA), and 2 non-coding regions (CR D-loop and light chain sub-chain replication origin OL). Overall, the elongate loach shared the same gene arrangement and composition of the mitochondrial genes with other teleost fishes. The Ka/Ks ratios of all mitochondrial PCGs were less than 1, indicating that all of the PCGs were evolving under purifying selection. Genome comparison analyses showed a significant sequence homology of species of Leptobotia. A significant identity between L. elongata and the other five Leptobotia species was observed in the visualization result, except for L. mantschurica, which lacked the tRNA-Arg gene and had a shorter tRNA-Asp gene. The phylogenetic tree revealed that the Cobitidae species examined here can be grouped into two clades, with the elongate loach forming a sister relationship with L. microphthalma. This study could provide additional inferences for a better understanding of the phylogenetic relationships among Cobitidae species.

Assessing the phylogenetic relationship among varieties of Toona ciliata (Meliaceae) in sympatry with chloroplast genomes

Toona ciliata is an endangered species due to over-cutting and low natural regeneration in China. Its genetic conservation is of an increasing concern. However, several varieties are recognized according to the leaf and flower traits, which complicates genetic conservation of T. ciliata. Here, we sequenced the whole chloroplast genome sequences of three samples for each of four varieties (T. ciliata var. ciliata, T. ciliata var. yunnanensis, T. ciliata var. pubescens, and T. ciliata var. henryi) in sympatry and assessed their phylogenetic relationship at a fine spatial scale. The four varieties had genome sizes ranged from 159,546 to 159,617 bp and had small variations in genome structure. Phylogenomic analysis indicated that the four varieties were genetically well-mixed in branch groups. Genetic diversity from the whole chloroplast genome sequences of 12 samples was low among varieties (average π = 0.0003). Besides, we investigated genetic variation of 58 samples of the four varieties in sympatry using two markers (psaA and trnL-trnF) and showed that genetic differentiation was generally insignificant among varieties (Ф st = 0%-5%). Purifying selection occurred in all protein-coding genes except for the ycf2 gene that was under weak positive selection. Most amino acid sites in all protein-coding genes were under purifying selection except for a few sites that were under positive selection. The chloroplast genome-based phylogeny did not support the morphology-based classification. The overall results implicated that a conservation strategy based on the T. ciliata complex rather than on intraspecific taxon was more appropriate.

Cross-population genetic analysis revealed genetic variation and selection in the Ovar-DRB1 gene of Indian sheep breeds

In sheep, MHC variability is studied widely to explore disease association. The aim of the current study was to explore the genetic diversity of Ovar-DRB diversity across sheep breeds of India. Here, Ovar-DRB1 locus was studied across 20 sheep breeds. DRB1 was amplified (301 bp) and sequenced using a PCR-sequence-based typing approach. Results revealed a high degree of heterozygosity across breeds (mean: 73.99%). Overall mean distance for DRB1 was highest in Sangamneri (0.18) and lowest in Madgyal sheep (0.10). There was a higher rate of transition, across breeds. Further, 39 alleles were isolated in different breeds, out of which 10 were new. To allow easy access and use of the immune-polymorphic database, an online database management system was launched (http://www.mhcdbms.in/). Nucleotide content across breeds for the DRB1 region revealed the richness of GC content (59.26%). Wu-Kabat index revealed vast genetic variation across peptide binding sites (PBS) of DRB1. Residues 6, 66, 69, 52, and 81, were polymorphic showing utility for antigen presentation. All breeds were under positive selection for DRB1 locus (dN > dS). Study revealed the importance of DRB locus diversity for beta chain specifically at PBS across sheep breeds of the Indian subcontinent and presented evidence of positive selection for DRB owing to its evolutionary significance.

Identification and characterization of Glycolate oxidase gene family in garden lettuce (Lactuca sativa cv. 'Salinas') and its response under various biotic, abiotic, and developmental stresses

Glycolate oxidase (GLO) is an FMN-containing enzyme localized in peroxisomes and performs in various molecular and biochemical mechanisms. It is a key player in plant glycolate and glyoxylate accumulation pathways. The role of GLO in disease and stress resistance is well-documented in various plant species. Although studies have been conducted regarding the role of GLO genes from spinach on a microbial level, the direct response of GLO genes to various stresses in short-season and leafy plants like lettuce has not been published yet. The genome of Lactuca sativa cultivar 'Salinas' (v8) was used to identify GLO gene members in lettuce by performing various computational analysis. Dual synteny, protein-protein interactions, and targeted miRNA analyses were conducted to understand the function of GLO genes. The identified GLO genes showed further clustering into two groups i.e., glycolate oxidase (GOX) and hydroxyacid oxidase (HAOX). Genes were observed to be distributed unevenly on three chromosomes, and syntenic analysis revealed that segmental duplication was prevalent. Thus, it might be the main reason for GLO gene diversity in lettuce. Almost all LsGLO genes showed syntenic blocks in respective plant genomes under study. Protein-protein interactions of LsGLO genes revealed various functional enrichments, mainly photorespiration, and lactate oxidation, and among biological processes oxidative photosynthetic carbon pathway was highly significant. Results of in-depth analyses disclosed the interaction of GLO genes with other members of the glycolate pathway and the activity of GLO genes in various organs and developmental stages in lettuce. The extensive genome evaluation of GLO gene family in garden lettuce is believed to be a reference for cloning and studying functional analyses of GLO genes and characterizing other members of glycolate/glyoxylate biosynthesis pathway in various plant species.

Teosinte Pollen Drive guides maize diversification and domestication by RNAi

Meiotic drivers subvert Mendelian expectations by manipulating reproductive development to bias their own transmission. Chromosomal drive typically functions in asymmetric female meiosis, while gene drive is normally postmeiotic and typically found in males. Using single molecule and single-pollen genome sequencing, we describe Teosinte Pollen Drive, an instance of gene drive in hybrids between maize (Zea mays ssp. mays) and teosinte mexicana (Zea mays ssp. mexicana), that depends on RNA interference (RNAi). 22nt small RNAs from a non-coding RNA hairpin in mexicana depend on Dicer-Like 2 (Dcl2) and target Teosinte Drive Responder 1 (Tdr1), which encodes a lipase required for pollen viability. Dcl2, Tdr1, and the hairpin are in tight pseudolinkage on chromosome 5, but only when transmitted through the male. Introgression of mexicana into early cultivated maize is thought to have been critical to its geographical dispersal throughout the Americas, and a tightly linked inversion in mexicana spans a major domestication sweep in modern maize. A survey of maize landraces and sympatric populations of teosinte mexicana reveals correlated patterns of admixture among unlinked genes required for RNAi on at least 4 chromosomes that are also subject to gene drive in pollen from synthetic hybrids. Teosinte Pollen Drive likely played a major role in maize domestication and diversification, and offers an explanation for the widespread abundance of "self" small RNAs in the germlines of plants and animals.

Diversity and evolution of the MHC class II DRB gene in the Capra sibirica experienced a demographic fluctuation in China

The major histocompatibility complex (MHC) genes are the most polymorphic genes in vertebrates, and their proteins play a critical role in adaptive immunity for defense against a variety of pathogens. MHC diversity was lost in many species after experiencing a decline in size. To understand the variation and evolution of MHC genes in the Siberian ibex, Capra sibirica, which has undergone a population decline, we analyzed the variation of the second exon of MHC class II DRB genes in samples collected from five geographic localities in Xinjiang, China, that belong to three diverged mitochondrial clades. Consequently, we identified a total of 26 putative functional alleles (PFAs) with 260 bp in length from 43 individuals, and found one (for 27 individuals) to three (for 5 individuals) PFAs per individual, indicating the presence of one or two DRB loci per haploid genome. The Casi-DRB1*16 was the most frequently occurring PFA, Casi-DRB1*22 was found in only seven individuals, 14 PFAs occurred once, 7 PFAs twice, implying high frequency of rare PFAs. Interestingly, more than half (15) of the PFAs were specific to clade I, only two and three PFAs were specific to clades II and III, respectively. So, we assume that the polygamy and sexual segregation nature of this species likely contributed to the allelic diversity of DRB genes. Genetic diversity indices showed that PFAs of clade II were lower in nucleotide, amino acid, and supertype diversity compared to those of the other two clades. The pattern of allele sharing and FST values between the three clades was to some extent in agreement with the pattern observed in mitochondrial DNA divergence. In addition, recombination analyses revealed no evidence for significant signatures of recombination events. Alleles shared by clades III and the other two clades diverged 6 million years ago, and systematic neighbor grids showed Trans-species polymorphism. Together with the PAML and MEME analyses, the results indicated that the DRB gene in C. sibirica evolved under balancing and positive selection. However, by comparison, it can be clearly seen that different populations were under different selective pressures. Our results are valuable in understanding the diversity and evolution of the DRB gene in a mountain living C. sibirica and in making decisions on future long-term protection strategies.

Diagnostic and commensal Staphylococcus pseudintermedius genomes reveal niche adaptation through parallel selection of defense mechanisms

Staphylococcus pseudintermedius is historically understood as a prevalent commensal and pathogen of dogs, though modern clinical diagnostics reveal an expanded host-range that includes humans. It remains unclear whether differentiation across S. pseudintermedius populations is driven primarily by niche-type or host-species. We sequenced 501 diagnostic and commensal isolates from a hospital, veterinary diagnostic laboratory, and within households in the American Midwest, and performed a comparative genomics investigation contrasting human diagnostic, animal diagnostic, human colonizing, pet colonizing, and household-surface S. pseudintermedius isolates. Though indistinguishable by core and accessory gene architecture, diagnostic isolates harbor more encoded and phenotypic resistance, whereas colonizing and surface isolates harbor similar CRISPR defense systems likely reflective of common household phage exposures. Furthermore, household isolates that persist through anti-staphylococcal decolonization report elevated rates of base-changing mutations in - and parallel evolution of - defense genes, as well as reductions in oxacillin and trimethoprim-sulfamethoxazole susceptibility. Together we report parallel niche-specific bolstering of S. pseudintermedius defense mechanisms through gene acquisition or mutation.

PhylteR: Efficient Identification of Outlier Sequences in Phylogenomic Datasets

In phylogenomics, incongruences between gene trees, resulting from both artifactual and biological reasons, can decrease the signal-to-noise ratio and complicate species tree inference. The amount of data handled today in classical phylogenomic analyses precludes manual error detection and removal. However, a simple and efficient way to automate the identification of outliers from a collection of gene trees is still missing. Here, we present PhylteR, a method that allows rapid and accurate detection of outlier sequences in phylogenomic datasets, i.e. species from individual gene trees that do not follow the general trend. PhylteR relies on DISTATIS, an extension of multidimensional scaling to 3 dimensions to compare multiple distance matrices at once. In PhylteR, these distance matrices extracted from individual gene phylogenies represent evolutionary distances between species according to each gene. On simulated datasets, we show that PhylteR identifies outliers with more sensitivity and precision than a comparable existing method. We also show that PhylteR is not sensitive to ILS-induced incongruences, which is a desirable feature. On a biological dataset of 14,463 genes for 53 species previously assembled for Carnivora phylogenomics, we show (i) that PhylteR identifies as outliers sequences that can be considered as such by other means, and (ii) that the removal of these sequences improves the concordance between the gene trees and the species tree. Thanks to the generation of numerous graphical outputs, PhylteR also allows for the rapid and easy visual characterization of the dataset at hand, thus aiding in the precise identification of errors. PhylteR is distributed as an R package on CRAN and as containerized versions (docker and singularity).

Genome-Wide Analysis of bHLH Family Genes and Identification of Members Associated with Cold/Drought-Induced Photoinhibition in Kandelia obovata

Plant basic helix-loop-helix (bHLH) transcription factors play pivotal roles in responding to stress, including cold and drought. However, it remains unclear how bHLH family genes respond to these stresses in Kandelia obovata. In this study, we identified 75 bHLH members in K. obovata, classified into 11 subfamilies and unevenly distributed across its 18 chromosomes. Collineation analysis revealed that segmental duplication primarily drove the expansion of KobHLH genes. The KobHLH promoters were enriched with elements associated with light response. Through RNA-seq, we identified several cold/drought-associated KobHLH genes. This correlated with decreased net photosynthetic rates (Pn) in the leaves of cold/drought-treated plants. Weighted gene co-expression network analysis (WGCNA) confirmed that 11 KobHLH genes were closely linked to photoinhibition in photosystem II (PS II). Among them, four Phytochrome Interacting Factors (PIFs) involved in chlorophyll metabolism were significantly down-regulated. Subcellular localization showed that KobHLH52 and KobHLH30 were located in the nucleus. Overall, we have comprehensively analyzed the KobHLH family and identified several members associated with photoinhibition under cold or drought stress, which may be helpfulfor further cold/drought-tolerance enhancement and molecular breeding through genetic engineering in K. obovata.

Complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens: genome structures, comparative and phylogenetic analysis

The genus Cerastium includes about 200 species that are mostly found in the temperate climates of the Northern Hemisphere. Here we report the complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens. The length of cp genomes ranged from 147,940 to 148,722 bp. Their quadripartite circular structure had the same gene organization and content, containing 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Repeat sequences varied from 16 to 23 per species, with palindromic repeats being the most frequent. The number of identified SSRs ranged from 20 to 23 per species and they were mainly composed of mononucleotide repeats containing A/T units. Based on Ka/Ks ratio values, most genes were subjected to purifying selection. The newly sequenced chloroplast genomes were characterized by a high frequency of RNA editing, including both C to U and U to C conversion. The phylogenetic relationships within the genus Cerastium and family Caryophyllaceae were reconstructed based on the sequences of 71 protein-coding genes. The topology of the phylogenetic tree was consistent with the systematic position of the studied species. All representatives of the genus Cerastium were gathered in a single clade with C. glomeratum sharing the least similarity with the others.

Bioinformatic Analysis Reveals the Role of Translation Elongation Efficiency Optimisation in the Evolution of Ralstonia Genus

Translation efficiency modulates gene expression in prokaryotes. The comparative analysis of translation elongation efficiency characteristics of Ralstonia genus bacteria genomes revealed that these characteristics diverge in accordance with the phylogeny of Ralstonia. The first branch of this genus is a group of bacteria commonly found in moist environments such as soil and water that includes the species R. mannitolilytica, R. insidiosa, and R. pickettii, which are also described as nosocomial infection pathogens. In contrast, the second branch is plant pathogenic bacteria consisting of R. solanacearum, R. pseudosolanacearum, and R. syzygii. We found that the soil Ralstonia have a significantly lower number and energy of potential secondary structures in mRNA and an increased role of codon usage bias in the optimization of highly expressed genes' translation elongation efficiency, not only compared to phytopathogenic Ralstonia but also to Cupriavidus necator, which is closely related to the Ralstonia genus. The observed alterations in translation elongation efficiency of orthologous genes are also reflected in the difference of potentially highly expressed gene' sets' content among Ralstonia branches with different lifestyles. Analysis of translation elongation efficiency characteristics can be considered a promising approach for studying complex mechanisms that determine the evolution and adaptation of bacteria in various environments.

Phylogenomic Analysis of Cytochrome P450 Gene Superfamily and Their Association with Flavonoids Biosynthesis in Peanut (Arachis hypogaea L.)

Cytochrome P450s (CYPs) constitute extensive enzyme superfamilies in the plants, playing pivotal roles in a multitude of biosynthetic and detoxification pathways essential for growth and development, such as the flavonoid biosynthesis pathway. However, CYPs have not yet been systematically studied in the cultivated peanuts (Arachis hypogaea L.), a globally significant cash crop. This study addresses this knowledge deficit through a comprehensive genome-wide analysis, leading to the identification of 589 AhCYP genes in peanuts. Through phylogenetic analysis, all AhCYPs were systematically classified into 9 clans, 43 gene families. The variability in the number of gene family members suggests specialization in biological functions. Intriguingly, both tandem duplication and fragment duplication events have emerged as pivotal drivers in the evolutionary expansion of the AhCYP superfamily. Ka/Ks analysis underscored the substantial influence of strong purifying selection on the evolution of AhCYPs. Furthermore, we selected 21 genes encoding 8 enzymes associated with the flavonoid pathway. The results of quantitative real-time PCR (qRT-PCR) experiments unveiled stage-specific expression patterns during the development of peanut testa, with discernible variations between pink and red testa. Importantly, we identified a direct correlation between gene expression levels and the accumulation of metabolites. These findings offer valuable insights into elucidating the comprehensive functions of AhCYPs and the underlying mechanisms governing the divergent accumulation of flavonoids in testa of different colors.

The Mitochondrial Genome of Linichthys laticeps (Cypriniformes: Cyprinidae): Characterization and Phylogeny

Mitochondrial genomes (mitogenomes) have been widely used in phylogenetic analysis and evolutionary biology. The Labeoninae is the largest subfamily of Cypriniformes and has great economic importance and ecological value. In this study, we sequenced, annotated, and characterized the complete mitogenome of Linichthys laticeps and then constructed the phylogenetic tree with previously published Labeoninae mitogenomes. The mitogenome of L. laticeps was 16,593 bp in length, with an A + T content of 57.1%. The mitogenome contained a standard set of 37 genes and a control region with the same order and orientation of genes as most fish mitogenomes. Each protein-coding gene (PCG) was initiated by an initial ATG codon, excluding COI, that began with a GTG codon. Furthermore, most of the PCGs were terminated by a conventional stop codon (TAA/TAG), while an incomplete termination codon (TA/T) was detected in 7 of the 13 PCGs. Most tRNA genes in L. laticeps were predicted to fold into the typical cloverleaf secondary structures. The Ka/Ks (ω) values for all PCGs were below one. The phylogenetic relationships of 96 Labeoninae mitogenomes indicated that Labeoninae was not a monophyletic group and L. laticeps was closely related to the genera Discogobio and Discocheilus. Overall, our study provided the first complete annotated mitogenome of L. laticeps, which filled a knowledge gap in Labeoninae and extended the understanding of the taxonomy and mitogenomic phylogeny of the subfamily Labeoninae.

Evolution of the viral genomes of SARS-CoV-2 in association with the changes in local condition: a genomic epidemiological study of a suburban city of Japan

Understanding the factors driving the spread and evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) at the local, regional, national, and international levels is important in protecting against future pandemics. By exploring their viral genomes, we attempted to analyse the spread of SARS-CoV-2 and its evolutionary convergence in Kashiwa City, as an example of a representative commuter town in Japan. From September 2020 to January 2023, a total of 47,134 nasopharyngeal swab and saliva specimens were collected from patients in 47 local clinics and hospitals, covering the vast majority of healthcare facilities. All SARS-CoV-2-positive samples were subjected to whole genome sequencing. Based on the analysis of 5,536 identified genomes, all major strains were represented. Unique regional mutations were occasionally identified in each strain. Inspection of these mutations revealed that the overall base substitution rate increased with progressive waves of the pandemic, at an overall rate of 2.56 bases/year. Interestingly, the spread and evolutionary patterns appeared to be distinct between regions and between individual clinics. Further analysis of the synonymous base substitution rate showed that the speed of viral evolution accelerated coincident with the beginning of public vaccination. Comprehensive genomic epidemiological studies, as presented here, should be useful in precisely understanding the pandemic and preparing for possible future pandemics.

An approach for prioritizing candidate genes from RNA-seq using preclinical cocaine self-administration datasets as a test case

RNA-sequencing (RNA-seq) technology has led to a surge of neuroscience research using animal models to probe the complex molecular mechanisms underlying brain function and behavior, including substance use disorders. However, findings from rodent studies often fail to be translated into clinical treatments. Here, we developed a novel pipeline for narrowing candidate genes from preclinical studies by translational potential and demonstrated its utility in 2 RNA-seq studies of rodent self-administration. This pipeline uses evolutionary conservation and preferential expression of genes across brain tissues to prioritize candidate genes, increasing the translational utility of RNA-seq in model organisms. Initially, we demonstrate the utility of our prioritization pipeline using an uncorrected P-value. However, we found no differentially expressed genes in either dataset after correcting for multiple testing with false discovery rate (FDR < 0.05 or <0.1). This is likely due to low statistical power that is common across rodent behavioral studies, and, therefore, we additionally illustrate the use of our pipeline on a third dataset with differentially expressed genes corrected for multiple testing (FDR < 0.05). We also advocate for improved RNA-seq data collection, statistical testing, and metadata reporting that will bolster the field's ability to identify reliable candidate genes and improve the translational value of bioinformatics in rodent research.

Microevolutionary patterns in ecotypes of the symbiotic cyanobacterium UCYN-A revealed from a Northwest Atlantic coastal time series

UCYN-A is a globally important nitrogen-fixing symbiotic microbe often found in colder regions and coastal areas where nitrogen fixation has been overlooked. We present a 3-year coastal Northwest Atlantic time series of UCYN-A by integrating oceanographic data with weekly nifH and16S rRNA gene sequencing and quantitative PCR assays for UCYN-A ecotypes. High UCYN-A relative abundances dominated by A1 to A4 ecotypes reoccurred annually in the coastal Northwest Atlantic. Although UCYN-A was detected every summer/fall, the ability to observe separate ecotypes may be highly dependent on sampling time given intense interannual and weekly variability of ecotype-specific occurrences. Additionally, much of UCYN-A's rarer diversity was populated by short-lived neutral mutational variants, therefore providing insight into UCYN-A's microevolutionary patterns. For instance, rare ASVs exhibited community composition restructuring annually, while also sharing a common connection to a dominant ASV within each ecotype. Our study provides additional perspectives for interpreting UCYN-A intraspecific diversity and underscores the need for high-resolution datasets when deciphering spatiotemporal ecologies within UCYN-A.

Molecular Evolution and Protein Structure Variation of Dkk Family

Dkks have inhibitory effects on the Wnt signaling pathway, which is involved in the development of skin and its appendages and the regulation of hair growth. The nucleotide sequences were compared and analyzed to further investigate the relationship between the structure and function of the Dkk gene family and vertebrate epidermal hair. The analysis of the molecular evolution of the Dkk family revealed that the evolution rate of the genes changed significantly after speciation, with the Aves and Reptilia branches showing accelerated evolution. Additionally, positive selection was observed at specific sites. The tertiary structure of the protein was also predicted. The analysis of the functional divergence of the Dkk family revealed that the functional divergence coefficient of each gene was greater than 0, with most of the functional divergence sites were located in the Cys-2 domain and a few in the Cys-1 domain. This suggests that the amino acid and functional divergence sites may play a role in regulating the binding of the Dkk family to LRP5/6, and thus affect the inhibition of Wnt signaling, leading to different functions of Dkk1, Dkk2, and Dkk4 in the development of skin hair follicles. In addition, the Dkk families of Aves and Reptilia may have undergone adaptive evolution and functional divergence.

Universality of evolutionary trajectories under arbitrary forms of self-limitation and competition

The assumption of constant population size is central in population genetics. It led to a large body of results that is robust to modeling choices and that has proven successful to understand evolutionary dynamics. In reality, allele frequencies and population size are both determined by the interaction between a population and the environment. Relaxing the constant-population assumption has two big drawbacks. It increases the technical difficulty of the analysis, and it requires specifying a mechanism for the saturation of the population size, possibly making the results contingent on model details. Here we develop a framework that encompasses a great variety of systems with an arbitrary mechanism for population growth limitation. By using techniques based on scale separation for stochastic processes, we are able to calculate analytically properties of evolutionary trajectories, such as the fixation probability. Remarkably, these properties assume a universal form with respect to our framework, which depends on only three parameters related to the intergeneration timescale, the invasion fitness, and the carrying capacity of the strains. In other words, different systems, such as Lotka-Volterra or a chemostat model (contained in our framework), share the same evolutionary outcomes after a proper remapping of their parameters. An important and surprising consequence of our results is that the direction of selection can be inverted, with a population evolving to reach lower values of invasion fitness.