Phylogenetic Evolution Research Articles

Orphan genes are not a distinct biological entity.

The genome sequencing revolution has revealed that all species possess a large number of unique genes critical for trait variation, adaptation, and evolutionary innovation. One widely used approach to identify such genes consists of detecting protein-coding sequences with no homology in other genomes, termed orphan genes. These genes have been extensively studied, under the assumption that they represent valid proxies for species-specific genes. Here, we critically evaluate taxonomic, phylogenetic, and sequence evolution evidence showing that orphan genes belong to a range of evolutionary ages and thus cannot be assigned to a single lineage. Furthermore, we show that the processes generating orphan genes are substantially more diverse than generally thought and include horizontal gene transfer, transposable element domestication, and overprinting. Thus, orphan genes represent a heterogeneous collection of genes rather than a single biological entity, making them unsuitable as a subject for meaningful investigation of gene evolution and phenotypic innovation.

Genome-wide analysis of the Amorphophallus konjac AkCSLA gene family and its functional characterization in drought tolerance of transgenic arabidopsis

BackgroundAmorphophallus konjac (A. konjac), a perennial tuberous plant, is widely cultivated for its high konjac glucomannan (KGM) content, a heteropolysaccharide with diverse applications. The cellulose synthase-like (CSL) gene family is known to be a group of processive glycan synthases involved in the synthesis of cell-wall polysaccharides and plays an important role in the biological process of KGM. However, in A. konjac the classification, structure, and function of the AkCSLA superfamily have been studied very little.ResultsBioinformatics methods were used to identify the 11 AkCSLA genes from the whole genome of Amorphophallus konjac and to systematically analyze their characteristics, phylogenetic evolution, promoter cis-elements, expression patterns, and subcellular locations. Phylogenetic analysis revealed that the AkCSLA gene family can be divided into three subfamilies (Groups I- III), which have close relationships with Arabidopsis. The promoters of most AkCSLA family members contain MBS elements and ABA response elements. Analysis of expression patterns in different tissues showed that most AkCSLAs are highly expressed in the corms. Notably, PEG6000 induced down-regulation of the expression of most AkCSLAs, including AkCSLA11. Subcellular localization results showed that AkCSLA11 was localized to the plasma membrane, Golgi apparatus and endoplasmic reticulum. Transgenic Arabidopsis experiments demonstrated that overexpression of AkCSLA11 reduced the plant’s drought tolerance. This overexpression also inhibited the expression of drought response genes and altered the sugar components of the cell wall. These findings provide new insights into the response mechanisms of A. konjac to drought stress and may offer potential genetic resources for improving crop drought resistance.ConclusionIn conclusion, the study reveals that the AkCSLA11 gene from A. konjac negatively impacts drought tolerance when overexpressed in Arabidopsis. This discovery provides valuable insights into the mechanisms of plant response to drought stress and may guide future research on crop improvement for enhanced resilience.

Biological Characteristics of the Cytochrome P 450 Family and the Mechanism of Terpinolene Metabolism in Hyalomma asiaticum (Acari: Ixodidae)

The CYP450 enzyme is a superfamily enzyme ubiquitously found in nearly all organisms, playing a vital role in the metabolism of both endogenous and exogenous compounds, and in biosynthesis. Unfortunately, an understanding of its classification, functions, expression characteristics, and other biological traits in Hyalomma asiaticum, a vector for Crimean–Congo Hemorrhagic Fever, as well as of the genes implicated in its natural product metabolism, is lacking. Towards this end, this study has identified 120 H. asiaticum CYP450 genes via transcriptome data in the face of a joint genome threat from terpinolene. The proteins these genes encode are of higher molecular weight, devoid of a signal peptide, and composed of unstable hydrophobic proteins principally containing 1–3 variable transmembrane regions. Phylogenetic evolution classifies these H. asiaticum CYP450 genes into four subfamilies. These genes all encompass complete CYP450 conserved domains, and five specific conserved motifs, albeit with different expression levels. GO and KEGG annotation findings suggest a widespread distribution of these CYP450 genes in many physiological systems, predominantly facilitating lipid metabolism, terpenoid compound metabolism, and polyketone compound metabolism, as well as cofactor and vitamin metabolism at a cellular level. Molecular docking results reveal a hydrophobic interaction between the ARG-103, ARG-104, LEU-106, PHE-109, and ILE-119 amino acid residues in CYP3A8, which is primarily expressed in the fat body, and terpinolene, with a notably up-regulated expression, with affinity = −5.6 kcal/mol. The conservation of these five key amino acid residues varies across 12 tick species, implying differences in terpinolene metabolism efficacy among various tick species. This study thereby fills an existing knowledge gap regarding the biological characteristics of H. asiaticum CYP450 genes and paves the way for further research into the functions of these particular genes.

Genome-Wide Identification of COMT Gene Family in Maize and its Function in Response to Light.

Maize is a major crop, feed, and industrial material. Caffeic acid-O-methyltransferase (COMT) is a methylase closely associated with lignin biosynthesis and plant growth and resistance. In this study, we identified the COMT gene (ZmCOMT) family in maize and further analyzed its phylogenetic evolution, subcellular localization, and its function in response to light. Thirty-one ZmCOMT genes were identified in the maize genome, which were distributed across eight chromosomes and mainly clustered on chromosome 4. Most ZmCOMT proteins were predicted to localize in the cytoplasm. Ten different conserved motifs were present in most ZmCOMT proteins, and motif1, motif6, and motif7 were highly conserved and present in all ZmCOMT proteins. The photoresponsivity elements were conserved among all members, and ZmCOMT22 and ZmCOMT10 genes responsive to light. This result suggests a potential function for these two genes in lignin biosynthesis which a previous study had linked to light regulation. Jasmonic acid responsive and abscisic acid cis-acting elements were present in the promoter regions of family members, thus the family may be regulated by hormone signaling pathways of maize. In summary, ZmCOMT genes are ancient, and the highly conserved motifs may be significant in survival and evolution of maize. Furthermore, light may influence lignin biosynthesis and photosynthesis through ZmCOMT genes. This research provided theoretical basis for lignin biosynthesis of maize and the potential value of ZmCOMT22 and ZmCOMT10 genes to enhance plant photosynthesis for facing global warming.

Genome-Wide Identification and Expression Analysis of SNAP Gene Family in Wheat

Background/Objectives: The SNAP gene family is a class of proteins containing a SNAP domain, which plays a crucial role in the growth and development of plants. Methods: Bioinformatics methods were used to systematically analyze the gene structure, phylogenetic evolution, chromosomal distribution, physicochemical properties, conserved motifs, and cis-acting elements of the TaSNAP family members. Results: The TaSNAP family comprises members that encode proteins ranging between 120 and 276 amino acids, with isoelectric points spanning from 4.87 to 7.92. Phylogenetic analysis elucidated the categorization of the eight TaSNAP into three distinct subfamilies, wherein members of the same subfamily display marked similarities in their gene structures. Chromosomal mapping revealed the distribution of TaSNAP family members across chromosomes 2A, 2B, 2D, 7A, 7B, and 7D. Utilizing the Plant CARE tool, we identified ten elements linked to plant hormones and four associated with stress responses. Expression analysis via qRT-PCR was performed to assess the levels of the eight TaSNAP genes in various tissues and under diverse abiotic stress conditions. The results indicated heightened expression of most genes in roots compared to spikes. Notably, under ABA stress, the majority of genes exhibited upregulation, whereas certain genes were downregulated under PEG stress, implying a substantial role for SNAP protein in wheat growth and development. Conclusions: This study conducted a comprehensive bioinformatics analysis of each member of the wheat SNAP family, laying a crucial foundation for future functional investigations.

Genome-wide identification of the TRAF gene family in humpback grouper (Cromileptes altivelis) and analysis of their expression in response to Vibrio harveyi challenge

TRAF (Tumor necrosis factor receptor-associated factor) proteins are key mediators of signal transduction in cell signaling and immune regulation within the toll-like receptor (TLR) and tumor necrosis factor (TNFR) superfamily. Despite the importance of TRAF genes in teleost innate immunity, study on their functions in C. altivelis is limited. This study utilized bioinformatics methods to identify and named eight TRAF genes (CaTRAF2a, CaTRAF2a-like, CaTRAF2b, CaTRAF3, CaTRAF4a, CaTRAF5, CaTRAF6 and CaTRAF7) in C. altivelis. Phylogenetic, syntenic and molecular evolution revealed that all CaTRAF members were evolutionarily conserved in teleost. Domain analysis indicated the presence of a conserved N-terminal RING finger domain in all CaTRAF proteins. Most CaTRAF proteins also featured a MATH domain at the C-terminal, with the exception of CaTRAF7 which contained seven repeat WD40 domains. In addition, qRT-PCR was used to detect the expression patterns of nine different tissues and eight different embryonic development stages of healthy fish, and it was found that there were spatial and tissue specificities among the members. HE staining revealed evident pathological lesions in the tissues post V. harveyi infection. Atrophy and significant bending of the gill lamellae were observed in the gills, while irregular cell shapes, increased fat vacuoles, and enlarged cell volume were noted in the liver. Intestinal tissues displayed thickening of the muscle layer, elongation of intestinal villi, and increased folds. Moreover, the expression of TRAF gene changed significantly after V. harveyi infection. These results would help to clarify the molecular role of CaTRAF gene in the regulation of immune and inflammatory responses in C. altivelis.

Identification, expression, and function analysis of Rbpms2 splicing variants in Japanese flounder gonad

Rbpms2, an RNA-binding protein with multiple splicing (Rbpms), can interact with RNAs to involve oocyte development, thereby influencing female sex differentiation in vertebrates. Here, two splicing variants of the Rbpms2 gene from Japanese flounder (Paralichthys olivaceus) were identified, namely Rbpms2.1 and Rbpms2.2. The two variants exhibited 98.22 % amino acid homology, both featuring an RNA recognition motif (RRM) domain spanning positions 98–170 amino acids. They were relatively conserved throughout phylogenetic evolution. Differently, the C-terminal region of the Rbpms2.1 contains five additional sequential amino acids (–VRDQP–) compared to Rbpms2.2. The real-time qPCR results demonstrated that Rbpms2.1 and Rbpms2.2 had relatively abundant expression in the gonads of adult Japanese flounder, with higher expression levels in the ovary compared to the testis (P < 0.05). In situ hybridization results showed strong positive expression of Rbpms2 mRNA in oocytes at stages I-III during the V stage of ovarian development. In the testis atstage IV, the expression of Rbpms2 mRNA was mainly concentrated on primary spermatocytes. Importantly, Rbpms2 binding sites were found in the 3′UTR, 5′UTR, and ORF regions of the sex-related genes including dmrt1, sox9, amh, foxl2, and wnt4. siRNA interference and overexpression analysis of Rbpms2.1 and Rbpms2.2 in primary cells of the ovary and testis showed that Rbpms2 can repress the expression of male-related genes (dmrt1, sox9, and amh) and significantly promote the expression of female-related genes (foxl2 and wnt4). Our results revealed that Rbpms2 may play a critical role by targeting the sex-related genes in the gonad development of Japanese flounder.

First insight into the genomes of the Pulmonaria officinalis group (Boraginaceae) provided by repeatome analysis and comparative karyotyping

BackgroundThe genus Pulmonaria (Boraginaceae) represents a taxonomically complex group of species in which morphological similarity contrasts with striking karyological variation. The presence of different numbers of chromosomes in the diploid state suggests multiple hybridization/polyploidization events followed by chromosome rearrangements (dysploidy). Unfortunately, the phylogenetic relationships and evolution of the genome, have not yet been elucidated. Our study focused on the P. officinalis group, the most widespread species complex, which includes two morphologically similar species that differ in chromosome number, i.e. P. obscura (2n = 14) and P. officinalis (2n = 16). Ornamental cultivars, morphologically similar to P. officinalis (garden escapes), whose origin is unclear, were also studied. Here, we present a pilot study on genome size and repeatome dynamics of these closely related species in order to gain new information on their genome and chromosome structure.ResultsFlow cytometry confirmed a significant difference in genome size between P. obscura and P. officinalis, corresponding to the number of chromosomes. Genome-wide repeatome analysis performed on genome skimming data showed that retrotransposons were the most abundant repeat type, with a higher proportion of Ty3/Gypsy elements, mainly represented by the Tekay lineage. Comparative analysis revealed no species-specific retrotransposons or striking differences in their copy number between the species. A new set of chromosome-specific cytogenetic markers, represented by satellite DNAs, showed that the chromosome structure in P. officinalis was more variable compared to that of P. obscura. Comparative karyotyping supported the hybrid origin of putative hybrids with 2n = 15 collected from a mixed population of both species and outlined the origin of ornamental garden escapes, presumably derived from the P. officinalis complex.ConclusionsLarge-scale genome size analysis and repeatome characterization of the two morphologically similar species of the P. officinalis group improved our knowledge of the genome dynamics and differences in the karyotype structure. A new set of chromosome-specific cytogenetic landmarks was identified and used to reveal the origin of putative hybrids and ornamental cultivars morphologically similar to P. officinalis.

Prevalence and species identification of trematode metacercariae in Qiqihar, Northeast China.

Fishborne trematode (FBT) is an important group of parasites that are endemic worldwide to a certain extent. However, despite the epidemiological significance, the species and phylogenetic evolution characteristics of FBT metacercariae have not been well studied. In this study, a total of 600 Pseudorasbora parva (P. parva) specimens were collected from Qiqihar, 61.8% (371/600) were found to be infected with trematode metacercariae. A total of three kinds of trematodes metacercariae were obtained, and they were identified as Clonorchis sinensis (C. sinensis), Metorchis orientalis (M. orientalis), and Metorchis taiwanensis (M. taiwanensis) by morphological and phylogenetic analysis with infection rates of 47.7% (286/600), 15.5% (93/600), and 23.7% (142/600), respectively. Meanwhile, a survey of the three trematodes metacercariae showed that the infection rate of C. sinensis metacercariae was the highest in September, up to 66% (66/100), and the lowest in June at 26% (26/100). The infection rate of M. orientalis metacercariae was the highest in October at 26% (26/100) and the lowest in June at 5% (5/100). The infection rate of M. taiwanensis metacercariae was at its peak in November at 36% (36/100) and the lowest in July at 15% (15/100). The co-infection of metacercariae of C. sinensis and M. taiwanensis was the most common and reached a peak in October, and their infection rate was higher in autumn than in summer. The peak of infection intensity of metacercariae for C. sinensis, M. orientalis, and M. taiwanensis were different: C. sinensis was 24/g in September, M. orientalis was 7/g in October, and M. taiwanensis was 10/g in November. From the above results, it was confirmed that three species of trematodes metacercaria played an important role in infection of second intermediate hosts in Qiqihar region. Studying the morphological characteristics and sequencing the ITS2 gene for a phylogenetic tree of them will be useful for future molecular evolution, biology, and ecology of trematode metacercariae.

Herbarium specimens as tools for exploring the evolution of fatty acid-derived natural products in plants.

Plants synthesize natural products via lineage-specific offshoots of their core metabolic pathways, including fatty acid synthesis. Recent studies have shed light on new fatty acid-derived natural products and their biosynthetic pathways in disparate plant species. Inspired by this progress, we set out to develop tools for exploring the evolution of fatty-acid derived products. We sampled multiple species from all major clades of euphyllophytes, including ferns, gymnosperms, and angiosperms (monocots and eudicots), and we show that the compositional profiles (though not necessarily the total amounts) of fatty-acid derived surface waxes from preserved plant specimens are consistent with those obtained from freshly collected tissue in a semi-quantitative and sometimes quantitative manner. We then sampled herbarium specimens representing 57 monocot species to assess the phylogenetic distribution and evolution, of two fatty acid-derived natural products found in that clade: beta-diketones and alkylresorcinols. These chemical data, combined with analyses of 26 monocot genomes, revealed a co-occurrence (though not necessarily a causal relationship) between whole genome duplication and the evolution of diketone synthases from an ancestral alkylresorcinol synthase-like polyketide synthase. Limitations of using herbarium specimen wax profiles as proxies for those of fresh tissue seem likely to include effects from loss of epicuticular wax crystals, effects from preservation techniques, and variation in wax chemical profiles due to genotype or environment. Nevertheless, this work reinforces the widespread utility of herbarium specimens for studying leaf surface waxes (and possibly other chemical classes) and reveals some of the evolutionary history of fatty acid-derived natural products within monocots.

Unravelling the novel sex determination genotype with 'ZY' and a distinctive 2.15-2.95 Mb inversion among poplar species through haplotype-resolved genome assembly and comparative genomics analysis.

Populus tomentosa, an indigenous tree species, is widely distributed and cultivated over 1,000,000 km2 in China, contributing significantly to forest production, ecological conservation and urban-rural greening. Although a reference genome is available for P. tomentosa, the intricate interspecific hybrid origins, chromosome structural variations (SVs) and sex determination mechanisms remain confusion and unclear due to its broad and even overlapping geographical distribution, extensive morphological variations and cross infiltration among white poplar species. We conducted a haplotype-resolved de novo assembly of P. tomentosa elite individual GM107, which comprises subgenomes a and b with a total genome size of 714.9 Mb. We then analysed the formation of hybrid species and the phylogenetic evolution and sex differentiation across the entire genus. Phylogenomic analyses suggested that GM107 likely originated from a hybridisation event between P. alba (♀) and P. davidiana (♂) which diverged at approximately 3.8 Mya. A total of 1551 chromosome SVs were identified between the two subgenomes. More noteworthily, a distinctive inversion structure spanning 2.15-2.95 Mb was unveiled among Populus, Tacamahaca, Turaga, Aigeiros poplar species and Salix, highlighting a unique evolutionary feature. Intriguingly, a novel sex genotype of the ZY type, which represents a crossover between XY and ZW systems, was identified and confirmed through both natural and artificial hybrids populations. These novel insights offer significant theoretical value for the study of the species' evolutionary origins and serve as a valuable resource for ecological genetics and forest biotechnology.

Detection of novel coltivirus-related sequences in Haemaphysalis megaspinosa ticks collected from Kanagawa Prefecture, Japan.

Coltiviruses, belonging to the genus Coltivirus within the family Spinareoviridae, are predominantly tick-borne viruses. Some of these species have been implicated in human diseases; however, their diversity, geographical distribution, and evolutionary dynamics remain inadequately. Therefore, this study was undertaken to explore the phylogenetic evolution of coltiviruses and related viruses. Our results revealed the detection of novel coltivirus-related sequences in adult female Haemaphysalis megaspinosa ticks collected from Kanagawa Prefecture, Japan. Molecular phylogenetic analysis revealed a close association between the sequences and the genome sequences of known coltivirus-related viruses, namely Qinghe tick reovirus and Fennes virus. The putative coltivirus-related virus was tentatively designated the Nakatsu tick virus. This study provides insights into the phylogenetic evolution of coltiviruses and related viruses.

Genome-wide identification, phylogeny and expression analysis of Hsf gene family in Verbena bonariensis under low-temperature stress

BackgroundThe heat shock transcription factor (Hsf) is a crucial regulator of plant stress resistance, playing a key role in plant stress response, growth, and development regulation.ResultsIn this study, we utilized bioinformatics tools to screen 25 VbHsf members, which were named VbHsf1-VbHsf25. We used bioinformatics methods to analyze the sequence structure, physicochemical properties, conserved motifs, phylogenetic evolution, chromosome localization, promoter cis-acting elements, collinearity, and gene expression of Hsf heat shock transcription factor family members under low-temperature stress. The results revealed that the majority of the Hsf genes contained motif1, motif2, and motif3, signifying that these three motifs were highly conserved in the Hsf protein sequence of Verbena bonariensis. Although there were some variations in motif deletion among the members, the domain remained highly conserved. The theoretical isoelectric point ranged from 4.17 to 9.71, with 21 members being unstable proteins and the remainder being stable proteins. Subcellular localization predictions indicated that all members were located in the nucleus. Phylogenetic analysis of the Hsf gene family in V. bonariensis and Arabidopsis thaliana revealed that the Hsf gene family of V. bonariensis could be categorized into three groups, with group A comprising 17 members and group C having at least two members. Among the 25 Hsf members, there were 1–3 exons located on seven chromosome fragments, which were unevenly distributed. Collinearity analysis demonstrated the presence of seven pairs of homologous genes in the VbHsf gene family. The Ka/Ks ratios were less than one, indicating that the VbHsf gene underwent purification selection pressure. Additionally, nine genes in V. bonariensis were found to have collinearity with A. thaliana. Promoter analysis revealed that the promoters of all VbHsf genes contained various types of cis-acting elements related to hormones and stress. Based on RNA-seq data, qRT-PCR analysis of six highly expressed genes was performed, and it was found that VbHsf5, VbHsf14, VbHsf17, VbHsf18, VbHsf20 and VbHsf21 genes were highly expressed at 12 h of low-temperature treatment, and the expression decreased after 24 h, among which VbHsf14 was up-regulated at 12 h of low-temperature by 70-fold.ConclusionsOur study may help reveal the important roles of Hsf in plant development and show insight for the further molecular breeding of V. bonariensis.

Exploring the phylogeny and depth evolution of cusk eels and their relatives (Ophidiiformes: Ophidioidei)

With 289 known species in 51 genera, the ophidiiform family Ophidiidae together with their relatives from the Carapidae (36 species in eight genera) of the same suborder Ophidioidei dominate the deep sea, but some occur also in shallow water habitats. Despite their high species diversity in the deep sea and wide bathymetric distributions, their phylogenetic relationships and evolution remain unexplored due in part to sampling difficulties. Thanks to the biodiversity exploratory program entitled “Tropical Deep-Sea Benthos” and joint efforts between Taiwan and French teams for sampling from different localities across the Indo-West Pacific over the last two decades, we are able to compile comprehensive datasets for investigations. In this study, 59 samples representing 36 of 59 known ophidioid genera are selected and used to construct a multi-gene dataset to infer the phylogenetic relationships of ophidioid fishes and their relatives. Our results reveal that the Ophidiidae forms a paraphyletic group with respect to the Carapidae. The four main clades of Ophidioidei resolved are the (1) clade comprising species from the subfamily Brotulinae; (2) clade that includes species in the genera Acanthonus and Xyelacyba; (3) clade grouping Hypopleuron caninum with species from the family Carapidae; and (4) clade containing the species in the subfamily Brotulotaenilinae, Neobythitinae (in part), and Ophidiinae. Accordingly, we suggest the following new revisions based on our results and proposed morphological diagnoses. The subfamily Brotulinae should be elevated to the family level. The genera Xyelacyba and probably Tauredophidium (unsampled in this study) should be included in the newly established family Acanthonidae with Acanthonus. The families Carapidae and Ophidiidae are re-defined. Our time-calibrated phylogenetic and ancestral depth reconstructions enable us to clarify the evolutionary history of ophidiiform fishes and infer past patterns of species distributions at different depths. While Ophidiiformes is inferred to have originated in shallow waters around 96.25 million years ago (Mya), the common ancestor to the Ophidioidei is inferred to have invaded the deep sea around 90.22 Mya, the dates coinciding with the global anoxic event of the OAE2. The observed bathymetric distribution patterns in Ophidioidei most likely point to the mesopelagic zone as the center of origin and diversification. This was followed by multiple events of depth transitions or range expansions towards either shallower waters or greater depth zones, which were likely triggered by past climate changes during the Paleogene-Neogene.

Genome-wide identification, phylogenetic, structural and functional evolution of thecore components of ABA signaling in plant species: a focus on rice.

A genome-wide analysis had identified 642 ABA core component genes from 20 plant species, which were further categorized into three distinct subfamilies. The gene structures and evolutionary relationships of these genes had been characterized. PP2C_1, PP2C_2, and SnRK2_1 had emerged as key players in mediating the ABA signaling transduction pathway, specifically in rice, in response to abiotic stresses. The plant hormone abscisic acid (ABA) is essential for growth, development, and stress response, relying on its core components, pyrabactin resistance, pyrabactin resistance-like, and the regulatory component of ABA receptor (PYR/PYL/RCAR), 2C protein phosphatase (PP2C), sucrose non-fermenting-1-related protein kinase 2 (SnRK2). However, there's a lack of research on their structural evolution and functional differentiation across plants. Our study analyzed the phylogenetic, gene structure, homology, and duplication evolution of this complex in 20 plant species. We found conserved patterns in copy number and homology across subfamilies. Segmental and tandem duplications drove the evolution of these genes, while whole-genome duplication (WGD) expanded PYR/PYL/RCAR and PP2C subfamilies, enhancing environmental adaptation. In rice and Arabidopsis, the PYR/PYL/RCAR, PP2C, and SnRK2 genes showed distinct tissue-specific expression and responded to various stresses. Notably, PP2C_1 and PP2C_2 interacted with SnRK2_1 and were crucial for ABA signaling in rice. These findings offered new insights into ABA signaling evolution, interactions, and integration in green plants, benefiting future research in agriculture, evolutionary biology, ecology, and environmental science.

PEMPS: a phylogenetic software tool to model the evolution of metabolic pathways

BackgroundMetabolic pathways support the enzyme flux that converts input chemicals into energy and cellular building blocks. With a constant rate of input, steady-state flux is achieved when metabolite concentrations and reaction rates remain constant over time. Individual genes undergo mutation, while selection acts on higher level functions of the pathway, such as steady-state flux where applicable. Modeling the evolution of metabolic pathways through mechanistic sets of ordinary differential equations is a piece of the genotype–phenotype map model for interpreting genetic variation and inter-specific differences. Such models can generate distinct compensatory changes and adaptive changes from directional selection, indicating single nucleotide polymorphisms and fixed differences that could affect phenotype. If used for inference, this would ultimately enable detection of selection on metabolic pathways as well as inference of ancestral states for metabolic pathway function.ResultsA software tool for simulating the evolution of metabolic pathways based upon underlying biochemistry, phylogenetics, and evolutionary considerations is presented. The Python program, Phylogenetic Evolution of Metabolic Pathway Simulator (PEMPS), implements a mutation-selection framework to simulate the evolution of the pathway over a phylogeny by interfacing with COPASI to calculate the steady-state flux of the metabolic network, introducing mutations as alterations in parameter values according to a model, and calculating a fitness score and corresponding probability of fixation based on the change in steady-state flux value(s). Results from simulations are consistent with a priori expectations of fixation probabilities and systematic change in model parameters.ConclusionsThe PEMPS program simulates the evolution of a metabolic pathway with a mutation-selection modeling framework based on criteria like steady-state flux that is designed to work with SBML-formatted kinetic models, and Newick-formatted phylogenetic trees. The Python software is run on the Linux command line and is available at https://github.com/nmccloskey/PEMPS.

Challenge accepted: Evolutionary lineages versus taxonomic classification of North American shrub willows (Salix).

The huge diversity of Salix subgenus Chamaetia/Vetrix clade in North America and the lack of phylogenetic resolution within this clade has presented a difficult but fascinating challenge for taxonomists to resolve. Here we tested the existing taxonomic classification with molecular tools. In this study, 132 samples representing 46 species from 22 described sections of shrub willows from the United States and Canada were analyzed and combined with 67 samples from Eurasia. The ploidy levels of the samples were determined using flow cytometry and nQuire. Sequences were produced using a RAD sequencing approach and subsequently analyzed with ipyrad, then used for phylogenetic reconstructions (RAxML, SplitsTree), dating analyses (BEAST, SNAPPER), and character evolution analyses of 14 selected morphological traits (Mesquite). The RAD sequencing approach allowed the production of a well-resolved phylogeny of shrub willows. The resulting tree showed an exclusively North American (NA) clade in sister position to a Eurasian clade, which included some North American endemics. The NA clade began to diversify in the Miocene. Polyploid species appeared in each observed clade. Character evolution analyses revealed that adaptive traits such as habit and adaxial nectaries evolved multiple times independently. The diversity in shrub willows was shaped by an evolutionary radiation in North America. Most species were monophyletic, but the existing sectional classification could not be supported by molecular data. Nevertheless, monophyletic lineages share several morphological characters, which might be useful in the revision of the taxonomic classification of shrub willows.

Genome-wide identification and expression analysis of the SPL transcription factor family and its response to abiotic stress in Pisum sativum L

Squamous promoter binding protein-like (SPL) genes encode plant-specific transcription factors (TFs) that play essential roles in modulating plant growth, development, and stress response. Pea (Pisum sativum L.) is a coarse grain crop of great importance in food production, biodiversity conservation and molecular genetic research, providing genetic information and nutritional resources for improving agricultural production and promoting human health. However, only limited researches on the structure and functions of SPL genes exist in pea (PsSPLs). In this study, we identified 22 PsSPLs and conducted a genome-wide analysis of their physical characteristics, chromosome distribution, gene structure, phylogenetic evolution and gene expression patterns. As a result, the PsSPLs were unevenly distributed on the seven chromosomes of pea and harbored the SBP domain, which is composed of approximately 76 amino acid residues. The phylogenetic analysis revealed that the PsSPLs clustered into eight subfamilies and showed high homology with SPL genes in soybean. Further analysis showed the presence of segmental duplications in the PsSPLs. The expression patterns of 22 PsSPLs at different tissues, developmental stages and under various stimulus conditions were evaluated by qRT-PCR method. It was found that the expression patterns of PsSPLs from the same subfamily were similar in different tissues, the transcripts of most PsSPLs reached the maximum peak value at 14 days after anthesis in the pod. Abiotic stresses can cause significantly up-regulated PsSPL19 expression with spatiotemporal specificity, in addition, four plant hormones can cause the up-regulated expression of most PsSPLs including PsSPL19 in a time-dependent manner. Therefore, PsSPL19 could be a key candidate gene for signal transduction during pea growth and development, pod formation, abiotic stress and plant hormone response. Our findings should provide insights for the elucidating of development regulation mechanism and breeding for resistance to abiotic stress pea.

Dispersion of SARS-CoV-2 lineage BA.5.1.25 and its descendants in Peru during two COVID-19 waves in 2022

During the third year of the pandemic in Peru, the persistent transmission of SARS-CoV-2 led to the appearance of more transmissible and immune-evasive Omicron sublineages; in that context, the National Genomic Surveillance of SARS-CoV-2 performed by the Peruvian National Institute of Health detected spike mutations in the circulating Omicron BA.5.1.25 sublineage which was later designated as DJ.1 and increased during the fourth COVID-19 wave, this eventually branched into new sublineages. The introduction, emergence, and timing of the most recent common ancestor (tMRCA) of BA.5.1.25 and its descendants (DJ.1, DJ.1.1, DJ.1.2, and DJ.1.3) were investigated in this paper as well as the time lags between their emergence and identification by the Peruvian National Institute of Health. Our findings show that ongoing genomic surveillance of SARS-CoV-2 is critical for understanding its phylogenetic evolution and the emergence of novel variations.

Physical Seed Dormancy in Legumes: Molecular Advances and Perspectives.

Physical dormancy of seeds is a form of dormancy due to the presence of an impermeable seed coat layer, and it represents a feature for plants to adapt to environmental changes over an extended period of phylogenetic evolution. However, in agricultural practice, physical dormancy is problematic. because it prevents timely and uniform seed germination. Therefore, physical dormancy is an important agronomical trait to target in breeding and domestication, especially for many leguminous crops. Compared to the well-characterized physiological dormancy, research progress on physical dormancy at the molecular level has been limited until recent years, due to the lack of suitable research materials. This review focuses on the structure of seed coat, factors affecting physical dormancy, genes controlling physical dormancy, and plants suitable for studying physical dormancy at the molecular level. Our goal is to provide a plethora of information for further molecular research on physical dormancy.