Fundamental patterns of structural evolution revealed by chromosome-length genomes of cactophilic Drosophila.

A thorough understanding of adaptation and speciation requires model organisms with both a history of ecological and phenotypic study as well as a complete set of genomic resources. In particular, high-quality genome assemblies of ecological model organisms are needed to assess the evolution of genome structure and its role in adaptation and speciation. Here, we generate new genomes of cactophilic Drosophila, a crucial model clade for understanding speciation and ecological adaptation in xeric environments. We generated chromosome-level genome assemblies and complete annotations for seven populations across D. mojavensis, D. arizonae, and D. navojoa. We use this data first to establish the most robust phylogeny for this clade to date, and to assess patterns of molecular evolution across the phylogeny, showing concordance with a priori hypotheses regarding adaptive genes in this system. We then show that structural evolution occurs at constant rate across the phylogeny, varies by chromosome, and is correlated with molecular evolution. These results advance the understanding of the D. mojavensis clade by demonstrating core evolutionary genetic patterns and integrating those patterns to generate new gene-level hypotheses regarding adaptation. Our data are presented in a new public database (cactusflybase.arizona.edu), providing one of the most in-depth resources for the analysis of inter- and intraspecific evolutionary genomic data. Furthermore, we anticipate that the patterns of structural evolution identified here will serve as a baseline for future comparative studies to identify the factors that influence the evolution of genome structure across taxa.

Immunophenotypic and molecular changes during progression of papillary urothelial carcinoma.

PURPOSE: Urothelial carcinoma has various molecular subtypes, each with different tumor characteristics. Although it is known that molecular changes occur during tumor progression, little is known about the specifics of these changes. In this study, we performed transcriptional analysis to understand the molecular changes during tumor progression. MATERIALS AND METHODS: Formalin-fixed, paraffin-embedded tumor tissues were obtained from 12 patients with muscle-invasive bladder cancer (MIBC). The invasive and non-invasive papillary areas were identified in papillary urothelial carcinoma specimens. Immunohistochemistry (IHC) and mRNA sequencing were performed for each tumor area. RESULTS: Patients with CK5/6-negative and CK20-positive non-invasive papillary areas were selected and classified into the IHC switch subgroup (CK5/6-positive and CK20-negative in the invasive area) and the IHC unchanged subgroup (CK5/6-negative and CK20-positive in the invasive area) according to the IHC results of the invasive area. We identified differences in the mRNA expression between the non-invasive papillary and invasive areas of the papillary MIBC tissue samples. In both the non-invasive papillary and invasive areas, the IHC switch subgroup showed basal subtype gene expression, while the IHC unchanged subgroup demonstrated luminal subtype gene expression. CONCLUSIONS: The non-invasive papillary area showed a gene expression pattern similar to that of the invasive area. Therefore, even if the non-invasive papillary area exhibits a luminal phenotype on IHC, it can have a basal subtype gene expression depending on the invasive area.

Origin and evolution of microvilli.

BACKGROUND INFORMATION: Microvilli are finger-like, straight, and stable cellular protrusions that are filled with F-actin and present a stereotypical length. They are present in a broad range of cell types across the animal tree of life and mediate several fundamental functions, including nutrient absorption, photosensation, and mechanosensation. Therefore, understanding the origin and evolution of microvilli is key to reconstructing the evolution of animal cellular form and function. Here, we review the current state of knowledge on microvilli evolution and perform a bioinformatic survey of the conservation of genes encoding microvillar proteins in animals and their unicellular relatives. RESULTS: We first present a detailed description of mammalian microvilli based on two well-studied examples, the brush border microvilli of enterocytes and the stereocilia of hair cells. We also survey the broader diversity of microvilli and discuss similarities and differences between microvilli and filopodia. Based on our bioinformatic survey coupled with carefully reconstructed molecular phylogenies, we reconstitute the order of evolutionary appearance of microvillar proteins. We document the stepwise evolutionary assembly of the "molecular microvillar toolkit" with notable bursts of innovation at two key nodes: the last common filozoan ancestor (correlated with the evolution of microvilli distinct from filopodia) and the last common choanozoan ancestor (correlated with the emergence of inter-microvillar adhesions). CONCLUSION AND SIGNIFICANCE: We conclude with a scenario for the evolution of microvilli from filopodia-like ancestral structures in unicellular precursors of animals.

Reversions mask the contribution of adaptive evolution in microbiomes.

When examining bacterial genomes for evidence of past selection, the results depend heavily on the mutational distance between chosen genomes. Even within a bacterial species, genomes separated by larger mutational distances exhibit stronger evidence of purifying selection as assessed by dN/dS, the normalized ratio of nonsynonymous to synonymous mutations. Here, we show that the classical interpretation of this scale dependence, weak purifying selection, leads to problematic mutation accumulation when applied to available gut microbiome data. We propose an alternative, adaptive reversion model with opposite implications for dynamical intuition and applications of dN/dS. Reversions that occur and sweep within-host populations are nearly guaranteed in microbiomes due to large population sizes, short generation times, and variable environments. Using analytical and simulation approaches, we show that adaptive reversion can explain the dN/dS decay given only dozens of locally fluctuating selective pressures, which is realistic in the context of Bacteroides genomes. The success of the adaptive reversion model argues for interpreting low values of dN/dS obtained from long timescales with caution as they may emerge even when adaptive sweeps are frequent. Our work thus inverts the interpretation of an old observation in bacterial evolution, illustrates the potential of mutational reversions to shape genomic landscapes over time, and highlights the importance of studying bacterial genomic evolution on short timescales.

Molecular evolution and genotype shift of Porcine circoviruses type 2 in Vietnam.

Porcine Circovirus Type 2 (PCV2), a significant pathogen in the global swine industry, causes Porcine Circovirus Associated Diseases (PCVAD), contributing to substantial economic losses. This study investigates the genetic diversity and evolutionary dynamics of PCV2 in Vietnam from 2007 to 2023. We sequenced and analyzed 47 PCV2 genomes isolated from swine farms across Vietnam between 2022 and 2023, revealing predominant circulation of PCV2d (80.85%) followed by PCV2b (19.15%). Phylogenetic analysis identified PCV2 genotypes PCV2a, PCV2b, PCV2d, PCV2g, and PCV2h circulating in Vietnam, with PCV2d emerging as the most prevalent genotype. Comparison with historical data highlighted genotype shifts from PCV2b to PCV2d in 2014. Interestingly, PCV2h genotype was mainly observed between 2008 and 2012 but have not been detected since 2014. Regional analysis indicated varied PCV2 epidemiological patterns between northern and southern Vietnam. Amino acid substitutions within the capsid protein were identified, predominantly in antigenic regions critical for immune recognition. Positive selection analysis identified multiple sites under evolutionary pressure, indicating ongoing adaptation of Vietnamese PCV2 strains. These findings enhance understanding of PCV2 dynamics in Vietnam and underscore the importance of continuous surveillance and adaptive management strategies in controlling PCV2-associated diseases in swine populations.

Evolutionary algorithms simulating molecular evolution: a new field proposal.

The genetic blueprint for the essential functions of life is encoded in DNA, which is translated into proteins-the engines driving most of our metabolic processes. Recent advancements in genome sequencing have unveiled a vast diversity of protein families, but compared with the massive search space of all possible amino acid sequences, the set of known functional families is minimal. One could say nature has a limited protein "vocabulary." A major question for computational biologists, therefore, is whether this vocabulary can be expanded to include useful proteins that went extinct long ago or have never evolved (yet). By merging evolutionary algorithms, machine learning, and bioinformatics, we can develop highly customized "designer proteins." We dub the new subfield of computational evolution, which employs evolutionary algorithms with DNA string representations, biologically accurate molecular evolution, and bioinformatics-informed fitness functions, Evolutionary Algorithms Simulating Molecular Evolution.

Selection in molecular evolution.

Evolution requires selection. Molecular/chemical/preDarwinian evolution is no exception. One molecule must be selected over another for molecular evolution to occur and advance. Evolution, however, has no goal. The laws of physics have no utilitarian desire, intent or proficiency. Laws and constraints are blind to "usefulness." How then were potential multi-step processes anticipated, valued and pursued by inanimate nature? Can orchestration of formal systems be physico-chemically spontaneous? The purely physico-dynamic self-ordering of Chaos Theory and irreversible non-equilibrium thermodynamic "engines of disequilibria conversion" achieve neither orchestration nor formal organization. Natural selection is a passive and after-the-fact-of-life selection. Darwinian selection reduces to the differential survival and reproduction of the fittest already-living organisms. In the case of abiogenesis, selection had to be 1) Active, 2) Pre-Function, and 3) Efficacious. Selection had to take place at the molecular level prior to the existence of non-trivial functional processes. It could not have been passive or secondary. What naturalistic mechanisms might have been at play?

The complete chloroplast genome of Rhododendronambiguum and comparative genomics of related species.

Rhododendron Linnaeus, 1753, the largest genus of woody plants in the Northern Hemisphere, includes some of the most significant species in horticulture. Rhododendronambiguum Hemsl, 1911, a member of subsection Triflora Sleumer 1947, exemplifies typical alpine Rhododendron species. The analysis of the complete chloroplast genome of R.ambiguum offers new insights into the evolution of Rhododendron species and enhances the resolution of phylogenetic relationships. This genome is composed of 207,478 base pairs, including a pair of inverted repeats (IRs) of 47,249 bp each, separated by a large single-copy (LSC) region of 110,367 bp and a small single-copy (SSC) region of 2,613 bp. It contains 110 genes: 77 protein-coding genes, 29 tRNAs, four unique rRNAs (4.5S, 5S, 16S, and 23S), with 16 genes duplicated in the IRs. Comparative analyses reveal substantial diversity in the Rhododendron chloroplast genome structures, identifying a fourth variant pattern. Specifically, four highly divergent regions (trnI-rpoB, ndhE-psaC, rpl32-ndhF, rrn16S-trnI) were noted in the intergenic spacers. Additionally, 76 simple sequence repeats were identified. Positive selection signals were detected in four genes (cemA, rps4, rpl16, and rpl14), evidenced by high Ka/Ks ratios. Phylogenetic reconstruction based on two datasets (shared protein-coding genes and complete chloroplast genomes) suggests that R.ambiguum is closely related to R.concinnum Hemsley, 1889. However, the phylogenetic positions of subsection Triflora Pojarkova, 1952 species remain unresolved, indicating that the use of complete chloroplast genomes for phylogenetic research in Rhododendron requires careful consideration. Overall, our findings provide valuable genetic information that will enhance understanding of the evolution, molecular biology, and genetic improvement of Rhododendron spieces.

Population bottlenecks and sexual recombination shape diatom microevolution.

Diatoms are single-celled organisms that contribute approximately 20% of the global primary production and play a crucial role in biogeochemical cycles and trophic chains. Despite their ecological importance, our knowledge of microevolution is limited. We developed a model using the SLiM evolutionary framework to address this knowledge gap. As a reference, we used the diatom Pseudo-nitzschia multistriata, which has been extensively studied in the Gulf of Naples. Our model recapitulates what we observe in natural populations, with microevolutionary processes that occur annually during a three-stage bloom phase. Interestingly, we found that non-bloom phases allow the population to maintain sex-generated diversity produced during blooms. This finding suggests that non-bloom phases are critical to counteract bloom-related pressures and mitigate genetic divergence at the species level. Moreover, our model showed that despite the consistent genetic differentiation during bloom phases, the population tends to return to pre-bloom states. While our model is limited to neutral dynamics, our study provides valuable insights into diatoms' microevolution, paving the way to explore the ecological implications of the life history dynamics of these organisms.

Robust Sequence Design Space for the Isothermal Exponential Amplification of Short Oligonucleotides.

Advances in isothermal amplification techniques have accelerated development in biosensing applications and the design of complex molecular devices. The exponential amplification reaction technique, or EXPAR, is uniquely positioned to process molecular information from short oligonucleotide strands (≈10 nucleotides length) typically encountered in molecular computing or microRNA detection. Despite its conceptual simplicity (requiring only a template strand and two enzymes), the issue of nonspecific background amplification has hindered broader adoption. In this work, a new system configuration is established at 37 °C to achieve significantly improved performance. Critical sequence motifs responsible for the excellent signal-to-background profile are identified and generalized as a universal adapter design framework. Orthogonal template sequences generated from the framework are implemented for a triplex reaction and successfully evaluated mixtures of multiple-target inputs in a single-step, one-pot format without the need for exogenous agents.

Phylogenetic signatures reveal multilevel selection and fitness costs in SARS-CoV-2.

Background: Large-scale sequencing of SARS-CoV-2 has enabled the study of viral evolution during the COVID-19 pandemic. Some viral mutations may be advantageous to viral replication within hosts but detrimental to transmission, thus carrying a transient fitness advantage. By affecting the number of descendants, persistence times and growth rates of associated clades, these mutations generate localised imbalance in phylogenies. Quantifying these features in closely-related clades with and without recurring mutations can elucidate the tradeoffs between within-host replication and between-host transmission. Methods: We implemented a novel phylogenetic clustering algorithm ( mlscluster, https://github.com/mrc-ide/mlscluster) to systematically explore time-scaled phylogenies for mutations under transient/multilevel selection. We applied this method to a SARS-CoV-2 time-calibrated phylogeny with >1.2 million sequences from England, and characterised these recurrent mutations that may influence transmission fitness across PANGO-lineages and genomic regions using Poisson regressions and summary statistics. Results: We found no major differences across two epidemic stages (before and after Omicron), PANGO-lineages, and genomic regions. However, spike, nucleocapsid, and ORF3a were proportionally more enriched for transmission fitness polymorphisms (TFP)-homoplasies than other proteins. We provide a catalog of SARS-CoV-2 sites under multilevel selection, which can guide experimental investigations within and beyond the spike protein. Conclusions: This study provides empirical evidence for the existence of important tradeoffs between within-host replication and between-host transmission shaping the fitness landscape of SARS-CoV-2. This method may be used as a fast and scalable means to shortlist large sequence databases for sites under putative multilevel selection which may warrant subsequent confirmatory analyses and experimental confirmation.

Viral mutations can potentially carry a transient advantage, being simultaneously favourable for replication within hosts (e.g. by evading host immune responses) and deleterious to transmission (e.g. by having reduced cell binding). To identify such mutations, called transmission fitness polymorphisms (TFPs), we developed a clustering algorithm entitled mlscluster that computes clade-level statistics based on the number of descendants, persistence times, and growth rates of clades carrying a specific mutation in comparison with their immediate sisters without the mutation, which usually are different than expected in the presence of such TFPs. We then applied it to a representative SARS-CoV-2 time-scaled tree with >1 million whole-genome sequences from England. Our statistical analysis suggested approximately constant levels of transient selection across waves driven by very distinct variants. It also showed that genomic regions of known functional significance such as spike, nucleocapsid, and ORF3a were enriched for TFPs. This is the one of the first studies to characterise SARS-CoV-2 recurrent mutations potentially under multilevel selection, providing empirical evidence for the existence of important tradeoffs in selection between intrahost replication and inter-host transmission. Therefore, it provides target mutations for realistic coalescent-based modelling and laboratory-based investigations of their impacts and mechanisms of interaction with human cells.

Conservation and selective pressures shaping baleen whale olfactory receptor genes supports their use of olfaction in the marine environment.

The relative importance of various sensory modalities can shift in response to evolutionary transitions, resulting in changes to underlying gene families encoding their reception systems. The rapid birth-and-death process underlying the evolution of the large olfactory receptor (OR) gene family has accelerated genomic-level change for the sense of smell in particular. The transition from the land to sea in marine mammals is an attractive model for understanding the influence of habitat shifts on sensory systems, with the retained OR repertoire of baleen whales contrasting with its loss in toothed whales. In this study, we examine to what extent the transition from a terrestrial to a marine environment has influenced the evolution of baleen whale OR repertoires. We developed Gene Mining Pipeline (GMPipe) (https://github.com/AprilJauhal/GMPipe), which can accurately identify large numbers of candidate OR genes. GMPipe identified 707 OR sequences from eight baleen whale species. These repertoires exhibited distinct family count distributions compared to terrestrial mammals, including signs of relative expansion in families OR10, OR11 and OR13. While many receptors have been lost or show signs of random drift in baleen whales, others exhibit signs of evolving under purifying or positive selection. Over 85% of OR genes could be sorted into orthologous groups of sequences containing at least four homologous sequences. Many of these groups, particularly from family OR10, presented signs of relative expansion and purifying selective pressure. Overall, our results suggest that the relatively small size of baleen whale OR repertoires result from specialisation to novel olfactory landscapes, as opposed to random drift.

Phylogeography of horseshoe bat sarbecoviruses in Vietnam and neighbouring countries. Implications for the origins of SARS-CoV and SARS-CoV-2.

Previous studies on horseshoe bats (Rhinolophus spp.) have described many coronaviruses related to SARS-CoV (SARSCoVr) in China and only a few coronaviruses related to SARS-CoV-2 (SARSCoV2r) in Yunnan (southern China), Cambodia, Laos and Thailand. Here, we report the results of several field missions carried out in 2017, 2021 and 2022 across Vietnam during which 1218 horseshoe bats were sampled from 19 locations. Sarbecoviruses were detected in 11% of faecal RNA extracts, with much more positives among Rhinolophus thomasi (46%). We assembled 38 Sarbecovirus genomes, including 32 SARSCoVr, four SARSCoV2r, and two recombinants of SARSCoVr and SARSCoV2r (RecSar), one showing a Spike protein very similar to SARS-CoV-2. We detected a bat co-infected with four coronaviruses, including two sarbecoviruses. Our analyses revealed that Sarbecovirus genomes evolve in Vietnam under strong geographical and host constraints. First, we found evidence for a deep separation between viruses from northern Vietnam and those from central and southern Vietnam. Second, we detected only SARSCoVr in Rhinolophus thomasi, both SARSCoVr and SARSCoV2r in Rhinolophus affinis, and only RecSar in Rhinolophus pusillus captured close to the border with China. Third, the bias in favour of Uracil in synonymous third codon positions of SARSCoVr extracted from R. thomasi showed a negative correlation with latitudes. Our results also provided support for an emergence of SARS-CoV in horseshoe bats from northern Yunnan and emergence of SARS-CoV-2 in horseshoe bats from northern Indochina subtropical forests (southern Yunnan, northern Laos and north-western Vietnam).

Genomics of ecological adaptation in Canary Island Descurainia (Brassicaceae) and comparisons with other Brassicaceae.

Oceanic archipelagos provide striking examples of lineages that have radiated over pronounced ecological gradients. Accompanying this diversification, lineages have evolved adaptations allowing survival in extreme environments. Here, we investigate the genomic basis of ecological adaptation in Canary Island Descurainia (Brassicaceae), an island relative of Arabidopsis. The seven endemic species have diversified in situ along an elevational and ecological gradient, from low-elevation scrub to high-elevation sub-alpine desert. We first generated a reference genome for Descurainia millefolia, phylogenetic analysis of which placed it as sister to D. sophioides. Ninety-six gene families were found to be specific to D. millefolia and a further 1087 and 1469 gene families have expanded or contracted in size, respectively, along the D. millefolia branch. We then employed genome re-sequencing to sample 14 genomes across the seven species of Canary Island Descurainia and an outgroup. Phylogenomic analyses were consistent with previous reconstructions of Canary Island Descurainia in resolving low- and high-elevation clades. Using the branch-site dN/dS method, we detected positive selection for 275 genes on the branch separating the low- and high-elevation species and these positively selected genes (PSGs) were significantly enriched for functions related to reproduction and stress tolerance. Comparing PSGs to those in analyses of adaptation to elevation and/or latitude in other Brassicaceae, we found little evidence of widespread convergence and gene reuse, except for two examples, one of which was a significant overlap between Descurainia and Draba nivalis, a species restricted to high latitudes. The study of Canary Island Descurainia suggests that the transition to high-elevation environments such as that found in the high mountains of the Canary Islands involves selection on genes related to reproduction and stress tolerance but that repeated evolution across different lineages that have evolved into similar habitats is limited, indicating substantially different molecular trajectories to adaptation.

Molecular Evolutionary Analyses of the Fusion Genes in Human Parainfluenza Virus Type 4.

The human parainfluenza virus type 4 (HPIV4) can be classified into two distinct subtypes, 4a and 4b. The full lengths of the fusion gene (F gene) of 48 HPIV4 strains collected during the period of 1966-2022 were analyzed. Based on these gene sequences, the time-scaled evolutionary tree was constructed using Bayesian Markov chain Monte Carlo methods. A phylogenetic tree showed that the first division of the two subtypes occurred around 1823, and the most recent common ancestors of each type, 4a and 4b, existed until about 1940 and 1939, respectively. Although the mean genetic distances of all strains were relatively wide, the distances in each subtype were not wide, indicating that this gene was conserved in each subtype. The evolutionary rates of the genes were relatively low (4.41 × 10-4 substitutions/site/year). Moreover, conformational B-cell epitopes were predicted in the apex of the trimer fusion protein. These results suggest that HPIV4 subtypes diverged 200 years ago and the progenies further diverged and evolved.

How Phenograms and Cladograms Became Molecular Phylogenetic Trees.

Tree diagrams are the prevailing form of visualization in biological classification and phylogenetics. Already during the time of the so-called Systematist Wars from the mid-1960s until the 1980s most journal articles and textbooks published by systematists contained tree diagrams. Although this episode of systematics is well studied by historians and philosophers of biology, most analyses prioritize scientific theories over practices and tend to emphasize conflicting theoretical assumptions. In this article, I offer an alternative perspective by viewing the conflict through the lens of representational practices with a case study on tree diagrams that were used by numerical taxonomists (phenograms) and cladists (cladograms). I argue that the current state of molecular phylogenetics should not be interpreted as the result of a competition of views within systematics. Instead, molecular phylogenetics arose independently of systematics and elements of cladistics and phenetics were integrated into the framework of molecular phylogenetics, facilitated by the compatibility of phenetic and cladistic practices with the quantitative approach of molecular phylogenetics. My study suggests that this episode of scientific change is more complex than common narratives of battles and winners or conflicts and compromises. Today, cladograms are still used and interpreted as specific types of molecular phylogenetic trees. While phenograms and cladograms represented different forms of knowledge during the time of the Systematist Wars, today they are both used to represent evolutionary relationships. This indicates that diagrams are versatile elements of scientific practice that can change their meaning, depending on the context of use within theoretical frameworks.

Engineering regioselectivity of glycosyltransferase for efficient polydatin synthesis.

Resveratrol is a promising functional ingredient applied in food products. However, low bioavailability and poor water solubility, which can be improved by glycosylation, hinder its application. A uridine diphosphate-dependent glycosyltransferase (UGT) from Bacillus subtilis 168 (named UGTBS) presents potential application for resveratrol glycosylation; nonetheless, imprecise regioselectivity renders the synthesis of resveratrol-3-O-ß-D-glucoside (polydatin) difficult. Therefore, molecular evolution was applied to UGTBS. A triple mutant Y14I/I62G/M315W was developed for 3-OH glycosylation of resveratrol and polydatin accounted for 91% of the total product. Kinetic determination and molecular docking indicated that the enhancement of hydrogen bond interaction and altered conformation of the binding pocket increases the enzyme's affinity for the 3-OH group, stabilizing the enzyme-substrate intermediate and promoting polydatin formation. Furthermore, a fed-batch cascade reaction by periodic addition of resveratrol was conducted and nearly 20 mM polydatin was obtained. The mutant Y14I/I62G/M315W can be used for polydatin manufacture.

Genome-wide identification and molecular evolution of elongation family of very long chain fatty acids proteins in Cyrtotrachelus buqueti.

To reveal the molecular function of elongation family of very long chain fatty acids(ELO) protein in Cyrtotrachelus buqueti, we have identified 15 ELO proteins from C.buqueti genome. 15 CbuELO proteins were located on four chromosomes. Their isoelectric points ranged from 9.22 to 9.68, and they were alkaline. These CbuELO proteins were stable and hydrophobic. CbuELO proteins had transmembrane movement, and had multiple phosphorylation sites. The secondary structure of CbuELO proteins was mainly α-helix. A total of 10 conserved motifs were identified in CbuELO protein family. Phylogenetic analysis showed that molecular evolutionary relationships of ELO protein family between C. buqueti and Tribolium castaneum was the closest. Developmental transcriptome analysis indicated that CbuELO10, CbuELO13 and CbuELO02 genes were key enzyme genes that determine the synthesis of very long chain fatty acids in pupae and eggs, CbuELO6 and CbuELO7 were that in the male, and CbuELO8 and CbuELO11 were that in the larva. Transcriptome analysis under different temperature conditions indicated that CbuELO1, CbuELO5, CbuELO12 and CbuELO14 participated in regulating temperature stress responses. Transcriptome analysis at different feeding times showed CbuELO12 gene expression level in all feeding time periods was significant downregulation. The qRT-PCR experiment verified expression level changes of CbuELO gene family under different temperature and feeding time conditions. Protein-protein interaction analysis showed that 9 CbuELO proteins were related to each other, CbuELO1, CbuELO4 and CbuELO12 had more than one interaction relationship. These results lay a theoretical foundation for further studying its molecular function during growth and development of C. buqueti.

How did plants evolve the prenylation of specialized phenolic metabolites by means of UbiA prenyltransferases?

Prenylated phenolics occur in over 4000 species in the plant kingdom, most of which are known as specialized metabolites with high chemical diversity. Many of them have been identified as pharmacologically active compounds from various medicinal plants, in which prenyl residues play a key role in these activities. Prenyltransferases (PTs) responsible for their biosynthesis have been intensively studied in the last two decades. These enzymes are membrane-bound proteins belonging to the UbiA superfamily that occurs from bacteria to humans, and in particular those involved in plant specialized metabolism show strict specificities for both substrates and products. This article reviews the enzymatic features of plant UbiA PTs, including C- and O-prenylation, molecular evolution, and application of UbiA PTs in synthetic biology.

Genome-wide identification and molecular evolution of Dof gene family in Camellia oleifera.

DNA binding with one finger(Dof) gene family is a class of transcription factors which play an important role on plant growth and development. Genome-wide identification results indicated that there were 45 Dof genes(ColDof) in C.oleifera genome. All 45 ColDof proteins were non-transmembrane and non-secretory proteins. Phosphorylation site analysis showed that biological function of ColDof proteins were mainly realized by phosphorylation at serine (Ser) site. The secondary structure of 44 ColDof proteins was dominated by random coil, and only one ColDof protein was dominated by α-helix. ColDof genes' promoter region contained a variety of cis-acting elements, including light responsive regulators, gibberellin responsive regulators, abscisic acid responsive regulators, auxin responsive regulators and drought induction responsive regulators. The SSR sites analysis showed that the proportion of single nucleotide repeats and the frequency of A/T in ColDof genes were the largest. Non-coding RNA analysis showed that 45 ColDof genes contained 232 miRNAs. Transcription factor binding sites of ColDof genes showed that ColDof genes had 5793 ERF binding sites, 4381 Dof binding sites, 2206 MYB binding sites, 3702 BCR-BPC binding sites. ColDof9, ColDof39 and ColDof44 were expected to have the most TFBSs. The collinearity analysis showed that there were 40 colinear locis between ColDof proteins and AtDof proteins. Phylogenetic analysis showed that ColDof gene family was most closely related to that of Camellia sinensis var. sinensis cv.Biyun and Camellia lanceoleosa. Protein-protein interaction analysis showed that ColDof34, ColDof20, ColDof28, ColDof35, ColDof42 and ColDof26 had the most protein interactions. The transcriptome analysis of C. oleifera seeds showed that 21 ColDof genes were involved in the growth and development process of C. oleifera seeds, and were expressed in 221 C. oleifera varieties. The results of qRT-PCR experiments treated with different concentrations NaCl and PEG6000 solutions indicated that ColDof1, ColDof2, ColDof14 and ColDof36 not only had significant molecular mechanisms for salt stress tolerance, but also significant molecular functions for drought stress tolerance in C. oleifera. The results of this study provide a reference for further understanding of the function of ColDof genes in C.oleifera.