Evolutionary Dynamics Research Articles

Genome-wide analysis of XTH gene family in cucumber (Cucumis sativus) against different insecticides to enhance defense mechanism

This study investigates how plants protect themselves, specifically looking at a group of genes called Xyloglucan Endotransglucosylase/Hydrolase (XTH) in cucumber plants. A detailed investigation of CsXTH genes, uncover their evolutionary background, structural motifs, regulatory networks, and functional responsibilities. The investigation utilized various bioinformatics tools and databases to highlight the dynamic interplay between CsXTH genes and abiotic stress responses. The investigation begins with a phylogenetic analysis that reveals the cladding of CsXTH genes into four categories. Further unraveling the conserved motif and domain analyses expose the coordination of xyloglucan modification, highlighting essential motifs and domains shared among these genes. Analyzing cis-regulatory elements helps us understand the complex networks that control the defense mechanisms of cucumber plants. Among these elements, CGTCA-motif and TGACG-motif stand out as key players in interpreting environmental signals and coordinating stress responses. The microRNA Csa-miR169 is becoming a key player in regulating the expression of CsXTH genes, and controlling the biological responses of cucumber plants. Chromosome mapping, Ka/Ks calculations, and protein-protein interaction analysis offer insights into the structural and evolutionary dynamics of CsXTH genes. Additionally, the study explores sub-cellular localization, uncovering the diverse cellular compartments where CsXTH genes are active. Finally, through chrysophanol and physcion, we unveil the adaptable functions of specific CsXTH genes in response to distinct treatments, underscoring their role in orchestrating cucumber's adaptive strategies.

Evolution and regulation of animal sex chromosomes.

Animal sex chromosomes typically carry the upstream sex-determining gene that triggers testis or ovary development and, in some species, are regulated by global dosage compensation in response to functional decay of the Y chromosome. Despite the importance of these pathways, they exhibit striking differences across species, raising fundamental questions regarding the mechanisms underlying their evolutionary turnover. Recent studies of non-model organisms, including insects, reptiles and teleosts, have yielded a broad view of the diversity of sex chromosomes that challenges established theories. Moreover, continued studies in model organisms with recently developed technologies have characterized the dynamics of sex determination and dosage compensation in three-dimensional nuclear space and at single-cell resolution. Here, we synthesize recent insights into sex chromosomes from a variety of species to review their evolutionary dynamics with respect to the canonical model, as well as their diverse mechanisms of regulation.

A new era of mutation rate analyses: Concepts and methods.

The mutation rate is a pivotal biological characteristic, intricately governed by natural selection and historically garnering considerable attention. Recent advances in high-throughput sequencing and analytical methodologies have profoundly transformed our understanding in this domain, ushering in an unprecedented era of mutation rate research. This paper aims to provide a comprehensive overview of the key concepts and methodologies frequently employed in the study of mutation rates. It examines various types of mutations, explores the evolutionary dynamics and associated theories, and synthesizes both classical and contemporary hypotheses. Furthermore, this review comprehensively explores recent advances in understanding germline and somatic mutations in animals and offers an overview of experimental methodologies, mutational patterns, molecular mechanisms, and driving forces influencing variations in mutation rates across species and tissues. Finally, it proposes several potential research directions and pressing questions for future investigations.

Genome-Wide Identification and Characterization of CCT Gene Family from Microalgae to Legumes.

The CCT (CO, COL and TOC1) gene family has been elucidated to be involved in the functional differentiation of the products in various plant species, but their specific mechanisms are poorly understood. In the present investigation, we conducted a genome-wide identification and phylogenetic analysis of CCT genes from microalgae to legumes. A total of 700 non-redundant members of the CCT gene family from 30 species were identified through a homology search. Phylogenetic clustering with Arabidopsis and domain conservation analysis categorized the CCT genes into three families. Multiple sequence alignment showed that the CCT domain contains important amino acid residues, and each CCT protein contains 24 conserved motifs, as demonstrated by the motif analysis. Whole-genome/segment duplication, as well as tandem duplication, are considered to be the driving forces in the evolutionary trajectory of plant species. This comprehensive investigation into the proliferation of the CCT gene family unveils the evolutionary dynamics whereby WGD/segment duplication is the predominant mechanism contributing to the expansion of the CCT genes. Meanwhile, the examination of the gene expression patterns revealed that the expression patterns of CCT genes vary in different tissues and at different developmental stages of plants, with high expression in leaves, which is consistent with the molecular regulation of flowering in photosynthesis by CCT. Based on the protein-protein interaction analysis of CCT genes in model plants, we propose that the CCT gene family synergistically regulates plant development and flowering with light-signaling factors (PHYs and PIFs) and MYB family transcription factors. Understanding the CCT gene family's molecular evolution enables targeted gene manipulation for enhanced plant traits, including optimized flowering and stress resistance.

Deciphering the molecular toolkit: regulatory elements governing shell biomineralization in marine molluscs.

The intricate process of shell biomineralization in marine molluscs is governed by a complex interplay of regulatory elements, encompassing secretomes, transporters, and noncoding RNA. This review delves into recent advancements in understanding these regulatory mechanisms, emphasizing their significance in elucidating the functions and evolutionary dynamics of the molluscan shell biomineralization process. Central to this intricate orchestration are secretomes with diverse functional domains, selectively exported to the extrapallial space, which directly regulate crystal growth and morphology. Transporters are crucial for substrate transportation in the calcification and maintenance of cellular homeostasis. Beyond proteins and transporters, noncoding RNA molecules are integral components influencing shell biomineralization. This review underscores the nonnegligible roles played by these genetic elements at the molecular level. To comprehend the complexity of biomineralization in mollusc, we explore the origin and evolutionary history of regulatory elements, primarily secretomes. While some elements have recently evolved, others are ancient genes that have been co-opted into the biomineralization toolkit. These elements undergo structural and functional evolution through rapidly evolving repetitive low-complexity domains and domain gain/loss/rearrangements, ultimately shaping a distinctive set of secretomes characterized by both conserved features and evolutionary innovations. This comprehensive review enhances our understanding of molluscan biomineralization at the molecular and genetic levels.

Rising above the noise: the influence of population dynamics on the evolution of acoustic signaling

Acoustic signaling is employed by many sexually reproducing species to select for mates and enhance fitness. However, signaling in dense populations can create an auditory background, or chorus, which may interfere with a signal receiver’s phonotactic selectivity, or the ability to distinguish individual signals. Feedback between the strength of an individual’s signal, phonotactic selectivity, and population size, may interact in complex ways to impact the evolution of signaling within a population, potentially leading to the emergence of silence. Here we formulate a general model that captures the dynamic feedback between individual acoustic signalers, phonotactic selectivity, and the population-level chorus to explore the eco-evolutionary dynamics of an acoustic trait within a population. We find that population dynamics have a significant influence on the evolutionary dynamics of the signaling trait, and that very sharp transitions separate conspicuous from silent populations. Our framework also reveals that increased phonotactic selectivity promotes the stability of signaling populations, and that transitions from signaling to silence are prone to hysteresis. We suggest that understanding the relationship between factors influencing population size, such as environmental productivity, as well as factors influencing phonotactic selectivity, such as anthropogenic noise, are central to understanding the complex mosaic of acoustically signaling and silent populations.

Evaluating Stakeholders’ Decisions in A Blockchain-Based Recycling Construction Waste Project: A Hybrid Evolutionary Game and System Dynamics Approach

To promote efficient construction waste recycling and reuse, a novel waste management approach based on blockchain technology was introduced to the industry. However, adopting blockchain platforms in construction waste recycling and reuse may impact the behavioral strategies of stakeholders and impede the prediction of the specific impacts of stakeholders’ decisions. Accordingly, this study addresses two primary questions: (1) What are the collaborative framework and the behavioral evolution trends of multiple stakeholders within the context of blockchain? (2) How can the behavioral strategies of multiple stakeholders be systematically coordinated to achieve efficient construction waste recycling and reuse driven by blockchain? To answer these questions, a tripartite game model combined with system dynamics was constructed. In this model, we aimed to elucidate the internal organizational framework, analyze the dynamic evolution process, and assess the influence of decisions made by multiple stakeholders at the individual level. It also offers corresponding policy recommendations for efficient construction waste recycling and reuse driven by blockchain at the system level. This study offers three innovations. First, it considers the decision-making of multiple stakeholders as an interdependent and coevolutionary process to overcome the defects of analyzing only one type of participant. Second, in contrast to the static analysis method, it employs a dynamic system approach to deeply analyze the evolving structures of blockchain-based projects. Third, it provides a theoretical framework for the practical implementation of blockchain-driven platforms in managing construction waste recycling and reuse, thus fostering effective policy development and management practices. This framework aims to promote sustainable development in construction waste recycling and reuse projects in China as well as globally.

Genomic evolution of ST228 SCCmec-I MRSA 10 years after a major nosocomial outbreak.

In this study, we investigated the genomic changes in a major methicillin-resistant Staphylococcus aureus (MRSA) clone following a significant outbreak at a hospital. Whole-genome sequencing of MRSA isolates was utilized to explore the genomic evolution of post-outbreak MRSA strains. The epidemicity of the clone declined over time, coinciding with the introduction of multimodal infection control measures. A genome-wide association study (GWAS) identified multiple genes significantly associated with either high or low epidemic success, indicating alterations in mobilome, virulence, and defense mechanisms. Random Forest models pinpointed a gene related to fibrinogen binding as the most influential predictor of epidemicity. The decline of the MRSA clone may be attributed to various factors, including the implementation of new infection control measures, single nucleotide polymorphisms accumulation, and the genetic drift of a given clone. This research underscores the complex dynamics of MRSA clones, emphasizing the multifactorial nature of their evolution. The decline in epidemicity seems linked to alterations in the clone's genetic profile, with a probable shift towards decreased virulence and adaptation to long-term carriage. Understanding the genomic basis for the decline of epidemic clones is crucial to develop effective strategies for their surveillance and management, as well as to gain insights into the evolutionary dynamics of pathogen genomes.

Evolution in fossil time series reconciles observations in micro- and macroevolution.

Extrapolating microevolutionary models does not always provide satisfactory explanations for phenotypic diversification on million-year time scales. For example, short-term evolutionary change is often modeled assuming a fixed adaptive landscape, but macroevolutionary changes are likely to involve changes in the adaptive landscape itself. A better understanding of how the adaptive landscape changes across different time intervals and how these changes cause populations to evolve has the potential to narrow the gap between micro- and macroevolution. Here, we analyze two fossil diatom time series of exceptional quality and resolution covering time intervals of a few hundred thousand years using models that account for different behaviors of the adaptive landscape. We find that one of the lineages evolves on a randomly and continuously changing landscape, whereas the other lineage evolves on a landscape that shows a rapid shift in the position of the adaptive peak of a magnitude that is typically associated with species-level differentiation. This suggests phenotypic evolution beyond generational timescales may be a consequence of both gradual and sudden repositioning of adaptive peaks. Both lineages are showing rapid and erratic evolutionary change and are constantly readapting towards the optimal trait state, observations that align with evolutionary dynamics commonly observed in contemporary populations. The inferred trait evolution over a span of a few hundred thousand years in these two lineages is therefore chimeric in the sense that it combines components of trait evolution typically observed on both short and long timescales.

Parallel HIV-1 evolutionary dynamics in humans and rhesus macaques who develop broadly neutralizing antibodies.

Human immunodeficiency virus (HIV)-1 exhibits remarkable genetic diversity. For this reason, an effective HIV-1 vaccine must elicit antibodies that can neutralize many variants of the virus. While broadly neutralizing antibodies (bnAbs) have been isolated from HIV-1 infected individuals, a general understanding of the virus-antibody coevolutionary processes that lead to their development remains incomplete. We performed a quantitative study of HIV-1 evolution in two individuals who developed bnAbs. We observed strong selection early in infection for mutations affecting HIV-1 envelope glycosylation and escape from autologous strain-specific antibodies, followed by weaker selection for bnAb resistance later in infection. To confirm our findings, we analyzed data from rhesus macaques infected with viruses derived from the same two individuals. We inferred remarkably similar fitness effects of HIV-1 mutations in humans and macaques. Moreover, we observed a striking pattern of rapid HIV-1 evolution, consistent in both humans and macaques, that precedes the development of bnAbs. Our work highlights strong parallels between infection in rhesus macaques and humans, and it reveals a quantitative evolutionary signature of bnAb development.

Selected Lark Mitochondrial Genomes Provide Insights into the Evolution of Second Control Region with Tandem Repeats in Alaudidae (Aves, Passeriformes).

The control region (CR) regulates the replication and transcription of the mitochondrial genome (mitogenome). Some avian mitogenomes possess two CRs, and the second control region (CR2) may enhance replication and transcription; however, the CR2 in lark mitogenome appears to be undergoing loss and is accompanied by tandem repeats. Here, we characterized six lark mitogenomes from Alaudala cheleensis, Eremophila alpestris, Alauda razae, and Calandrella cinerea and reconstructed the phylogeny of Passerida. Through further comparative analysis among larks, we traced the evolutionary process of CR2. The mitochondrial gene orders were conserved in all published lark mitogenomes, with Cytb-trnT-CR1-trnP-ND6-trnE-remnant CR2 with tandem repeat-trnF-rrnS. Phylogenetic analysis revealed Alaudidae and Panuridae are sister groups at the base of Sylvioidea, and sporadic losses of CR2 may occur in their common ancestor. CR sequence and phylogeny analysis indicated CR2 tandem repeats were generated within CR2, originating in the ancestor of all larks, rather than inherited from CR1. The secondary structure comparison of tandem repeat units within and between species suggested slipped-strand mispairing and DNA turnover as suitable models for explaining the origin and evolution of these repeats. This study reveals the evolutionary process of the CR2 containing tandem repeat in Alaudidae, providing reference for understanding the evolutionary characteristics and dynamics of tandem repeats.

Red blood cell dynamics during malaria infection challenge the assumptions of mathematical models of infection dynamics.

For decades, mathematical models have been used to understand the course and outcome of malaria infections (i.e., infection dynamics) and the evolutionary dynamics of the parasites that cause them. A key conclusion of these models is that red blood cell (RBC) availability is a fundamental driver of infection dynamics and parasite trait evolution. The extent to which this conclusion holds will in part depend on model assumptions about the host-mediated processes that regulate RBC availability i.e., removal of uninfected RBCs and supply of RBCs. Diverse mathematical functions have been used to describe host-mediated RBC supply and clearance, but it remains unclear whether they adequately capture the dynamics of RBC supply and clearance during infection. Here, we use a unique dataset, comprising time-series measurements of erythrocyte (i.e., mature RBC) and reticulocyte (i.e., newly supplied RBC) densities during Plasmodium chabaudi malaria infection, and a quantitative data-transformation scheme to elucidate whether RBC dynamics conform to common model assumptions. We found that RBC clearance and supply are not well described by mathematical functions commonly used to model these processes. Furthermore, the temporal dynamics of both processes vary with parasite growth rate in a manner again not captured by existing models. Together, these finding suggest that new model formulations are required if we are to explain and ultimately predict the within-host population dynamics and evolution of malaria parasites.

Mean-field interacting multi-type birth–death processes with a view to applications in phylodynamics

Multi-type birth–death processes underlie approaches for inferring evolutionary dynamics from phylogenetic trees across biological scales, ranging from deep-time species macroevolution to rapid viral evolution and somatic cellular proliferation. A limitation of current phylogenetic birth–death models is that they require restrictive linearity assumptions that yield tractable message-passing likelihoods, but that also preclude interactions between individuals. Many fundamental evolutionary processes – such as environmental carrying capacity or frequency-dependent selection – entail interactions, and may strongly influence the dynamics in some systems. Here, we introduce a multi-type birth–death process in mean-field interaction with an ensemble of replicas of the focal process. We prove that, under quite general conditions, the ensemble’s stochastically evolving interaction field converges to a deterministic trajectory in the limit of an infinite ensemble. In this limit, the replicas effectively decouple, and self-consistent interactions appear as nonlinearities in the infinitesimal generator of the focal process. We investigate a special case that is rich enough to model both carrying capacity and frequency-dependent selection while yielding tractable message-passing likelihoods in the context of a phylogenetic birth–death model.

Chloroplast Genome Profiling and Phylogenetic Insights of the “Qixiadaxiangshui” Pear (Pyrus bretschneideri Rehd.1)

The “Qixiadaxiangshui” pear (Pyrus bretschneideri Rehd.1) is a highly valued cultivar known for its crisp texture, abundant juice, and rich aroma. In this study, we reported the first complete chloroplast genome sequence of the “Qixiadaxiangshui” pear, which is 159,885 bp in length with a GC content of 36.58%. The genome exhibits a typical circular quadripartite structure, comprising a large single-copy region (LSC), a small single-copy region (SSC), and a pair of inverted repeat regions (IRs). A total of 131 genes were identified, including 84 protein-coding genes, 8 rRNA genes, and 37 tRNA genes. We also identified 209 simple sequence repeats (SSRs) and several mutation hotspots, such as ndhC-trnM-CAU and trnR-UCU-atpA, which can be applied in molecular identification and phylogenetic studies of Pyrus. Comparative genomic analysis showed high conservation among ten pear cultivars. Phylogenetic analysis indicated that the “Qixiadaxiangshui” pear is closely related to germplasm Dangshansuli, Wonwhang, and Yali, suggesting a recent common ancestor. These findings provided valuable insights into the genetic diversity and evolutionary dynamics of the Pyrus species and contribute to the conservation and breeding of pear germplasm resources.

Comprehensive analysis of the Xya riparia genome uncovers the dominance of DNA transposons, LTR/Gypsy elements, and their evolutionary dynamics

Transposable elements (TEs) are DNA sequences that can move or replicate within a genome, and their study has become increasingly important in understanding genome evolution and function. The Tridactylidae family, including Xya riparia (pygmy mole cricket), harbors a variety of transposable elements (TEs) that have been insufficiently investigated. Further research is required to fully understand their diversity and evolutionary characteristics. Hence, we conducted a comprehensive repeatome analysis of X. riparia species using the chromosome-level assembled genome. The study aimed to comprehensively analyze the abundance, distribution, and age of transposable elements (TEs) in the genome. The results indicated that the genome was 1.67 Gb, with 731.63 Mb of repetitive sequences, comprising 27% of Class II (443.25 Mb), 16% of Class I (268.45 Mb), and 1% of unknown TEs (19.92 Mb). The study found that DNA transposons dominate the genome, accounting for approximately 60% of the total repeat size, with retrotransposons and unknown elements accounting for 37% and 3% of the genome, respectively. The members of the Gypsy superfamily were the most abundant amongst retrotransposons, accounting for 63% of them. The transposable superfamilies (LTR/Gypsy, DNA/nMITE, DNA/hAT, and DNA/Helitron) collectively constituted almost 70% of the total repeat size of all six chromosomes. The study further unveiled a significant linear correlation (Pearson correlation: r = 0.99, p-value = 0.00003) between the size of the chromosomes and the repetitive sequences. The average age of DNA transposon and retrotransposon insertions ranges from 25 My (million years) to 5 My. The satellitome analysis discovered 13 satellite DNA families that comprise about 0.15% of the entire genome. In addition, the transcriptional analysis of TEs found that DNA transposons were more transcriptionally active than retrotransposons. Overall, the study suggests that the genome of X. riparia is complex, characterized by a substantial portion of repetitive elements. These findings not only enhance our understanding of TE evolution within the Tridactylidae family but also provide a foundation for future investigations into the genomic intricacies of related species.

Illuminating Genetic Diversity and Selection Signatures in Matou Goats through Whole-Genome Sequencing Analysis.

(1) Background: Matou goats, native to Hunan and Hubei provinces in China, are renowned for their exceptional meat and skin quality. However, a comprehensive whole-genome-based exploration of the genetic architecture of this breed is scant in the literature. (2) Methods: To address this substantial gap, we used whole-genome sequences of 20 Matou goats and compared them with published genomic data of 133 goats of different breeds across China. This comprehensive investigation sought to assess genetic diversity, population structure, and the presence of genomic selection signals. (3) Results: The whole genome of Matou goat populations yielded a substantial catalog of over 19 million single nucleotide polymorphisms (SNPs), primarily distributed within intergenic and intron regions. The phylogenetic tree analysis revealed distinct clades corresponding to each goat population within the dataset. Notably, this analysis positioned Matou goats in a closer genetic affinity with Guizhou White goats, compared to other recognized goat breeds. This observation was corroborated by principal component analysis (PCA) and admixture analysis. Remarkably, Matou goats exhibited diminished genetic diversity and a notable degree of inbreeding, signifying a reduced effective population size. Moreover, the study employed five selective sweep detection methods (including PI, CLR, PI-Ratio, Fst, and XP-EHH) to screen top signal genes associated with critical biological functions, encompassing cardiomyocytes, immunity, coat color, and meat quality. (4) Conclusions: In conclusion, this study significantly advances our understanding of the current genetic landscape and evolutionary dynamics of Matou goats. These findings underscore the importance of concerted efforts in resource conservation and genetic enhancement for this invaluable breed.

Ecological and evolutionary dynamics of CRISPR-Cas systems in Clostridium botulinum: Insights from genome mining and comparative analysis

Understanding the prevalence and distribution of CRISPR-Cas systems across different strains can illuminate the ecological and evolutionary dynamics of Clostridium botulinum populations. In this study, we conducted genome mining to characterize the CRISPR-Cas systems of C. botulinum strains. Our analysis involved retrieving complete genome sequences of these strains and assessing the diversity, prevalence, and evolution of their CRISPR-Cas systems. Subsequently, we performed an analysis of homology in spacer sequences from identified CRISPR arrays to investigate and characterize the range of targeted phages and plasmids. Additionally, we investigated the evolutionary trajectory of C. botulinum strains under selective pressures from foreign invasive DNA. Our findings revealed that 306 strains possessed complete CRISPR-Cas structures, comprising 58% of the studied C. botulinum strains. Secondary structure prediction of consensus repeats indicated that subtype II-C, with longer stems compared to subtypes ID and IB, tended to form more stable RNA secondary structures. Moreover, protospacer motif analysis demonstrated that strains with subtype IB CRISPR-Cas systems exhibited 5′-CGG-3′, 5′-CC-3′, and 5′-CAT-3′ motifs in the 3′ flanking regions of protospacers. The diversity observed in CRISPR-Cas systems indicated their classification into subtypes IB, ID, II-C, III-B, and III-D. Furthermore, our results showed that systems with subtype ID and III-D frequently harbored similar spacer patterns. Moreover, analysis of spacer sequences homology with phage and prophage genomes highlighted the specific activities exhibited by subtype IB and III-B against phages and plasmids, providing valuable insights into the functional specialization within these systems.

Differentiation in detoxification gene complements, including neofunctionalization of duplicated cytochrome P450 genes, between lineages of cotton bollworm, Helicoverpa armigera.

Here we investigate the evolutionary dynamics of five enzyme superfamilies (CYPs, GSTs, UGTs, CCEs and ABCs) involved in detoxification in Helicoverpa armigera. The reference assembly for an African isolate of the major lineages, H. a. armigera, has 373 genes in the five superfamilies. Most of its CYPs, GSTs, UGTs and CCEs and a few of its ABCs occur in blocks and most of the clustered genes are in subfamilies specifically implicated in detoxification. Most of the genes have orthologues in the reference genome for the Oceania lineage, H. a. conferta. However, clustered orthologues and subfamilies specifically implicated in detoxification show greater sequence divergence and less constraint on non-synonymous differences between the two assemblies than do other members of the five superfamilies. Two duplicated CYPs, which were found in the H. a. armigera but not H. a. conferta reference genome, were also missing in 16 Chinese populations spanning two different lineages of H. a. armigera. The enzyme produced by one of these duplicates has higher activity against esfenvalerate than a previously described chimeric CYP mutant conferring pyrethroid resistance. Various transposable elements were found in the introns of most detoxification genes, generating diverse gene structures. Extensive resequencing data for the Chinese H. a. armigera and H. a. conferta lineages also revealed complex copy number polymorphisms in 17 CCE001s in a cluster also implicated in pyrethroid metabolism, with substantial haplotype differences between all three lineages. Our results suggest that cotton bollworm has a versatile complement of detoxification genes which are evolving in diverse ways across its range.

On the evolutionary origin of discrete phenotypic plasticity.

Phenotypic plasticity provides an attractive strategy for adapting to various environments, but the evolutionary mechanism of the underlying genetic system is poorly understood. We use a simple gene regulatory network model to explore how a species acquires phenotypic plasticity, particularly focusing on discrete phenotypic plasticity, which has been difficult to explain by quantitative genetic models. Our approach employs a population genetic framework that integrates the developmental process, where each individual undergoes growth to develop its phenotype, which subsequently becomes subject to selection pressures. Our model considers two alternative types of environments, with the gene regulatory network including a sensor gene that turns on and off depending on the type of environment. With this assumption, we demonstrate that the system gradually adapts by acquiring the ability to produce two distinct optimum phenotypes under two types of environments without changing genotype, resulting in phenotypic plasticity. We find that the resulting plasticity is often discrete after a lengthy period of evolution. Our results suggest that gene regulatory networks have a notable capacity to flexibly produce various phenotypes in response to environmental changes. This study also shows that the evolutionary dynamics of phenotype may differ significantly between mechanistic-based developmental models and quantitative genetics models, suggesting the utility of incorporating gene regulatory networks into evolutionary models.

Analysis of Codon Usage Bias in Chloroplast Genomes of Dryas octopetala var. asiatica (Rosaceae).

Dryas octopetala var. asiatica, a dwarf shrub belonging to the Rosaceae family and native to Asia, exhibits notable plasticity in photosynthesis in response to temperature variations. However, the codon usage patterns and factors influencing them in the chloroplast genome of this species have not yet been documented. This study sequenced and assembled the complete genome of D. octopetala var. asiatica. The annotated genes in the chloroplast genome were analyzed for codon composition through multivariate statistical methods including a neutrality plot, a parity rule 2 (PR2) bias plot, and an effective number of codons (ENC) plot using CodonW 1.4.2 software. The results indicated that the mean GC content of 53 CDSs was 38.08%, with the average GC content at the third codon base position being 27.80%, suggesting a preference for A/U(T) at the third codon position in chloroplast genes. Additionally, the chloroplast genes exhibited a weak overall codon usage bias (CUB) based on ENC values and other indicators. Correlation analysis showed a significant negative correlation between ENC value and GC2, an extremely positive correlation with GC3, but no correlation with GC1 content. These findings highlight the importance of the codon composition at the third position in influencing codon usage bias. Furthermore, our analysis indicated that the CUB of the chloroplast genome of D. octopetala var. asiatica was primarily influenced by natural selection and other factors. Finally, this study identified UCA, CCU, GCU, AAU, GAU, and GGU as the optimal codons. These results offer a foundational understanding for genetic modification and evolutionary dynamics of the chloroplast genome of D. octopetala var. asiatica.