Genome-wide identification and expression analysis of the Dof gene family reveals their involvement in hormone response and abiotic stresses in sunflower (Helianthus annuus L.)

DNA binding with one finger (Dof), plant-specific zinc finger transcription factors, can participate in various physiological and biochemical processes during the life of plants. As one of the most important oil crops in the world, sunflower (Helianthus annuus L.) has significant economic and ornamental value. However, a systematic analysis of H. annuus Dof (HaDof) members and their functions has not been extensively conducted. In this study, we identified 50 HaDof genes that are unevenly distributed on 17 chromosomes of sunflower. We present a comprehensive overview of the HaDof genes, including their chromosome locations, phylogenetic analysis, and expression profile characterization. Phylogenetic analysis classified the 366 Dof members identified from 11 species into four groups (further subdivided into nine subfamilies). Segmental duplications are predominantly contributed to the expansion of sunflower Dof genes, and all segmental duplicate gene pairs are under purifying selection due to strong evolutionary constraints. Furthermore, we observed differential expression patterns for HaDof genes in normal tissues as well as under hormone treatment or abiotic stress conditions by analyzing RNA-seq data from previous studies and RT-qPCR data in our current study. The expression of HaDof04 and HaDof43 were not detected in any samples, which implied that they may be gradually undergoing pseudogenization process. Some HaDof genes, such as HaDof25 and HaDof30, showed responsiveness to exogenous plant hormones, such as kinetin, brassinosteroid, auxin or strigolactone, while others like HaDof15 and HaDof35 may participate in abiotic stress resistance of sunflower seedling. Our study represents the initial step towards understanding the phylogeny and expression characterization of sunflower Dof family genes, which may provide valuable reference information for functional studies on hormone response, abiotic stress resistance, and molecular breeding in sunflower and other species.

Phylogenetic analyses reveal insights into interdomain horizontal gene transfer of microbial lipases

Microbial lipases play a pivotal role in a wide range of biotechnological processes and in the human skin microbiome. However, their evolution remains poorly understood. Accessing the evolutionary process of lipases could contribute to future applications in health and biotechnology. We investigated genetic events associated with the evolutionary trajectory of the microbial family LIP lipases. Using phylogenetic analysis, we identified two distinct horizontal gene transfer (HGT) events from Bacteria to Fungi. Further analysis of human cutaneous mycobiome members such as the lipophilic Malassezia yeasts and CUG-Ser-1 clade (including Candida sp. and other microorganisms associated with cutaneous mycobiota) revealed recent evolutionary processes, with multiple gene duplication events. The Lid region of fungal lipases, crucial for substrate interaction, exhibits varying degrees of conservation among different groups. Our findings suggest the adaptability of the fungal LIP family in various genetic and metabolic contexts and its potential role in niche exploration.

Control of leaf development in the water fern Ceratopteris richardii by the auxin efflux transporter CrPINMa in the CRISPR/Cas9 analysis

BACKGROUND

PIN-FORMED genes (PINs) are crucial in plant development as they determine the directionality of auxin flow. They are present in almost all land plants and even in green algae. However, their role in fern development has not yet been determined. This study aims to investigate the function of CrPINMa in the quasi-model water fern Ceratopteris richardii.

RESULTS

CrPINMa possessed a long central hydrophilic loop and characteristic motifs within it, which indicated that it belonged to the canonical rather than the non-canonical PINs. CrPINMa was positioned in the lineage leading to Arabidopsis PIN6 but not that to its PIN1, and it had undergone numerous gene duplications. CRISPR/Cas9 genome editing had been performed in ferns for the first time, producing diverse mutations including local frameshifts for CrPINMa. Plants possessing disrupted CrPINMa exhibited retarded leaf emergence and reduced leaf size though they could survive and reproduce at the same time. CrPINMa transcripts were distributed in the shoot apical meristem, leaf primordia and their vasculature. Finally, CrPINMa proteins were localized to the plasma membrane rather than other cell parts.

CONCLUSIONS

CRISPR/Cas9 genome editing is feasible in ferns, and that PINs can play a role in fern leaf development.

Origin, evolution, and diversification of inositol 1,4,5-trisphosphate 3-kinases in plants and animals

BACKGROUND

In Eukaryotes, inositol polyphosphates (InsPs) represent a large family of secondary messengers and play crucial roes in various cellular processes. InsPs are synthesized through a series of pohophorylation reactions catalyzed by various InsP kinases in a sequential manner. Inositol 1,4,5-trisphosphate 3-kinase (IP3 3-kinase/IP3K), one member of InsP kinase, plays important regulation roles in InsPs metabolism by specifically phosphorylating inositol 1,4,5-trisphosphate (IP3) to inositol 1,3,4,5-tetrakisphosphate (IP4) in animal cells. IP3Ks were widespread in fungi, plants and animals. However, its evolutionary history and patterns have not been examined systematically.

RESULTS

A total of 104 and 31 IP3K orthologues were identified across 57 plant genomes and 13 animal genomes, respectively. Phylogenetic analyses indicate that IP3K originated in the common ancestor before the divergence of fungi, plants and animals. In most plants and animals, IP3K maintained low-copy numbers suggesting functional conservation during plant and animal evolution. In Brassicaceae and vertebrate, IP3K underwent one and two duplication events, respectively, resulting in multiple gene copies. Whole-genome duplication (WGD) was the main mechanism for IP3K duplications, and the IP3K duplicates have experienced functional divergence. Finally, a hypothetical evolutionary model for the IP3K proteins is proposed based on phylogenetic theory.

CONCLUSION

Our study reveals the evolutionary history of IP3K proteins and guides the future functions of animal, plant, and fungal IP3K proteins.

Frustration can Limit the Adaptation of Promiscuous Enzymes Through Gene Duplication and Specialisation

Virtually all enzymes catalyse more than one reaction, a phenomenon known as enzyme promiscuity. It is unclear whether promiscuous enzymes are more often generalists that catalyse multiple reactions at similar rates or specialists that catalyse one reaction much more efficiently than other reactions. In addition, the factors that shape whether an enzyme evolves to be a generalist or a specialist are poorly understood. To address these questions, we follow a three-pronged approach. First, we examine the distribution of promiscuity in empirical enzymes reported in the BRENDA database. We find that the promiscuity distribution of empirical enzymes is bimodal. In other words, a large fraction of promiscuous enzymes are either generalists or specialists, with few intermediates. Second, we demonstrate that enzyme biophysics is not sufficient to explain this bimodal distribution. Third, we devise a constraint-based model of promiscuous enzymes undergoing duplication and facing selection pressures favouring subfunctionalization. The model posits the existence of constraints between the catalytic efficiencies of an enzyme for different reactions and is inspired by empirical case studies. The promiscuity distribution predicted by our constraint-based model is consistent with the empirical bimodal distribution. Our results suggest that subfunctionalization is possible and beneficial only in certain enzymes. Furthermore, the model predicts that conflicting constraints and selection pressures can cause promiscuous enzymes to enter a 'frustrated' state, in which competing interactions limit the specialisation of enzymes. We find that frustration can be both a driver and an inhibitor of enzyme evolution by duplication and subfunctionalization. In addition, our model predicts that frustration becomes more likely as enzymes catalyse more reactions, implying that natural selection may prefer catalytically simple enzymes. In sum, our results suggest that frustration may play an important role in enzyme evolution.

Uridine diphosphate glucosyltransferases are involved in spinosad resistance in western flower thrips Frankliniella occidentalis (Pergande)

Uridine diphosphate glucosyltransferases (UGTs) play crucial roles in the insect detoxification system and are associated with pesticide resistance. Our previous transcriptomic analysis of spinosad-susceptible (Ivf03) and resistant (NIL-R) Frankliniella occidentalis revealed numerous upregulated UGT genes in the NIL-R strain, suggesting their potential contribution to spinosad resistance. To investigate this hypothesis, here we conducted UGT activity assays and spinosad induction experiments, employing RNA interference (RNAi) techniques for gene function validation. We found significantly elevated UGT activity in the NIL-R strain compared to Ivf03, with 5-nitrouracil showing a substantial synergistic effect on the resistant strain. Eighteen UGT genes were identified in F. occidentalis, with gene expansion and duplication observed within families UGT466, 467, and 468. Ten out of the eighteen UGTs exhibited higher expression levels in NIL-R, specifically FoUGT466B1, FoUGT468A3, and FoUGT468A4 consistently being upregulated across nymphs, males, and females. RNAi-based functional validation targeting these three UGT genes led to increased susceptibility to spinosad in a life stage-, sex-, and dose-dependent manner. These results indicate that UGTs are indeed involved in spinosad resistance in F. occidentalis, and the effects are dependent on life stage, sex, and dose. Therefore, sustainable control for F. occidentalis resistance should always consider these differential responses.

Genome-wide analysis of cellulose synthase gene superfamily in Tectona grandis L.f

This study aimed to explore Cellulose synthase gene superfamily of teak, and its evolutionary relationship with homologous genes of other woody species. The incidence of evolutionary events like gene duplication and gene loss, influence of the selection pressure, and consequent adaptive functional divergence of the duplicated TgCes gene were assessed alongside it's role in wood coloration. This study identified 39 full-length non-redundant proteins belonging to CesA and Csl gene families. TgCesA and TgCsl proteins with Cellulose synthase domain repeats indicated tandem gene duplication and probable genetic variability, enabling local adaptation. Further, multi-domain protein (MYB-like DNA-binding domain and CesA domain) with maximum introns was also identified indicating gene fusion and formation of complex protein with novel functions. Phylogenetic analysis grouped the genes into seven subfamilies (CesA, CslA, CslC, CslD, CslE, CslG, and CslM) with each undergoing gene duplication and loss along their evolutionary history. Post-species gene duplications and probable neofunctionalization were identified in TgCesA and TgCsl gene families. Each subfamily was found to be under strong purifying selection with a few or no sites under positive selection. Functional divergence analysis further revealed site-specific selective constraints in CesA and Csl genes of the teak Cellulose synthase gene family. Furthermore, protein-protein interaction network analysis identified co-expression of Cellulose synthase gene with flavonoid 3',5'-hydroxylase (F3'5'H, CYP75A), involved in the biosynthesis of xylem anthocyanin compounds, probably responsible for wood coloration. This study thus offers a foundation for future research in wood formation and wood property traits specific to teak and its provenances.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03927-6.

O- and N-Methyltransferases in benzylisoquinoline alkaloid producing plants

BACKGROUND

Secondary metabolites such as benzylisoquinoline alkaloids (BIA) have attracted considerable attention because of their pharmacological properties and potential therapeutic applications. Methyltransferases (MTs) can add methyl groups to alkaloid molecules, altering their physicochemical properties and bioactivity, stability, solubility, and recognition by other cellular components. Five types of O-methyltransferases and two types of N-methyltransferases are involved in BIA biosynthesis.

OBJECTIVE

Since MTs may be the source for the discovery and development of novel biomedical, agricultural, and industrial compounds, we performed extensive molecular and phylogenetic analyses of O- and N-methyltransferases in BIA-producing plants.

METHODS

MTs involved in BIA biosynthesis were isolated from transcriptomes of Berberis koreana and Caulophyllum robustum. We also mined the methyltransferases of Coptis japonica, Papaver somniferum, and Nelumbo nucifera from the National Center for Biotechnology Information protein database. Then, we analyzed the functional motifs and phylogenetic analysis.

RESULT

We mined 42 O-methyltransferases and 8 N-methyltransferases from the five BIA-producing plants. Functional motifs for S-adenosyl-L-methionine-dependent methyltransferases were retained in most methyltransferases, except for the three O-methyltransferases from N. nucifera. Phylogenetic analysis revealed that the methyltransferases were grouped into four clades, I, II, III and IV. The clustering patterns in the phylogenetic analysis suggested a monophyletic origin of methyltransferases and gene duplication within species. The coexistence of different O-methyltransferases in the deep branch subclade might support some cases of substrate promiscuity.

CONCLUSIONS

Methyltransferases may be a source for the discovery and development of novel biomedical, agricultural, and industrial compounds. Our results contribute to further understanding of their structure and reaction mechanisms, which will require future functional studies.

Higher evolutionary dynamics of gene copy number for Drosophila glue genes located near short repeat sequences

BACKGROUND

During evolution, genes can experience duplications, losses, inversions and gene conversions. Why certain genes are more dynamic than others is poorly understood. Here we examine how several Sgs genes encoding glue proteins, which make up a bioadhesive that sticks the animal during metamorphosis, have evolved in Drosophila species.

RESULTS

We examined high-quality genome assemblies of 24 Drosophila species to study the evolutionary dynamics of four glue genes that are present in D. melanogaster and are part of the same gene family - Sgs1, Sgs3, Sgs7 and Sgs8 - across approximately 30 millions of years. We annotated a total of 102 Sgs genes and grouped them into 4 subfamilies. We present here a new nomenclature for these Sgs genes based on protein sequence conservation, genomic location and presence/absence of internal repeats. Two types of glue genes were uncovered. The first category (Sgs1, Sgs3x, Sgs3e) showed a few gene losses but no duplication, no local inversion and no gene conversion. The second group (Sgs3b, Sgs7, Sgs8) exhibited multiple events of gene losses, gene duplications, local inversions and gene conversions. Our data suggest that the presence of short "new glue" genes near the genes of the latter group may have accelerated their dynamics.

CONCLUSIONS

Our comparative analysis suggests that the evolutionary dynamics of glue genes is influenced by genomic context. Our molecular, phylogenetic and comparative analysis of the four glue genes Sgs1, Sgs3, Sgs7 and Sgs8 provides the foundation for investigating the role of the various glue genes during Drosophila life.

Lectin diversity and their positive roles in WSSV replication through regulation of calreticulin expression and inhibiting ALFs expression

In biological evolution, gene duplication (GD) generates new genes to facilitate new functions. C-type lectins (CTLs) in crayfish have been extended by GD to expand their family members. In this study, four CTL genes generated by GD were identified from Procambarus clarkii (PcLec1-4). Among these four genes, PcLec1 can also generate new isoforms with different numbers of tandem repeats through DNA slip mispairing. PcLec1-4 was widely expressed in multiple tissues. The expression levels of PcLec1-4 were upregulated in the intestine of P. clarkii upon white spot syndrome virus (WSSV) challenge at multiple time points. Further analysis indicated that GATA transcription factor regulated PcLec1-4 expression. RNA interference and recombinant PcLec1-4 protein injection experiments suggested that PcLec1-4 promoted the expression of calreticulin (PcCRT) and negatively regulated the expression of antimicrobial peptides, thereby promoting WSSV replication. This study contributes to the understanding of the function of CTLs produced by GD during WSSV invasion in crustaceans.

Functional Divergence and Origin of the Vertebrate Praja Family

The Praja family is an E3 ubiquitin ligase, promoting polyubiquitination and subsequent degradation of substrates. It comprises two paralogs, praja1 and praja2. Prior research suggests these paralogs have undergone functional divergence, with examples, such as their distinct roles in neurite outgrowth. However, the specific evolutionary trajectories of each paralog remain largely unexplored preventing mechanistic understanding of functional differences between paralogs. Here, we investigated the phylogeny and divergence of the vertebrate Praja family through molecular evolutionary analysis. Phylogenetic examination of the vertebrate praja revealed that praja1 and praja2 originated from the common ancestor of placentals via gene duplication, with praja1 evolving at twice the rate of praja2 shortly after the duplication. Moreover, a unique evolutionary trajectory for praja1 relative to other vertebrate Praja was indicated, as evidenced by principal component analysis on GC content, codon usage frequency, and amino acid composition. Subsequent motif/domain comparison revealed conserved N terminus and C terminus in praja1 and praja2, together with praja1-specific motifs, including nuclear localization signal and Ala-Gly-Ser repeats. The nuclear localization signal was demonstrated to be functional in human neuroblastoma SH-SY5Y cells using deletion mutant, while praja2 was exclusively expressed in the nucleus. These discoveries contribute to a more comprehensive understanding of the Praja family's phylogeny and suggest a functional divergence between praja1 and praja2. Specifically, the shift of praja1 into the nucleus implies the degradation of novel substrates located in the nucleus as an evolutionary consequence.

Genome-wide identification of the LRR-RLK gene family in peanut and functional characterization of AhLRR-RLK265 in salt and drought stresses

Leucine-rich repeat receptor-like kinases (LRR-RLKs) play important roles in plant developmental regulations and various stress responses. Peanut (Arachis hypogaea L.) is a worldwide important oil crop; however, no systematic identification or analysis of the peanut LRR-RLK gene family has been reported. In present study, 495 LRR-RLK genes in peanut were identified and analyzed. The 495 AhLRR-RLK genes were classed into 14 groups and 10 subgroups together with their Arabidopsis homologs according to phylogenetic analyses, and 491 of 495 AhLRR-RLK genes unequally located on 20 chromosomes. Analyses of gene structure and protein motif organization revealed similarity in exon/intron and motif organization among members of the same subgroup, further supporting the phylogenetic results. Gene duplication events were found in peanut LRR-RLK gene family via syntenic analysis, which were important in LRR-RLK gene family expansion in peanut. We found that the expression of AhLRR-RLK genes was detected in different tissues using RNA-seq data, implying that AhLRR-RLK genes may differ in function. In addition, Arabidopsis plants overexpressing stress-induced AhLRR-RLK265 displayed lower seed germination rates and root lengths compared to wild-type under exogenous ABA treatment. Notably, overexpression of AhLRR-RLK265 enhanced tolerance to salt and drought stresses in transgenic Arabidopsis. Moreover, the AhLRR-RLK265-OE lines were found to have higher activities of superoxide dismutase (SOD), catalase (CAT), and peroxidase (POD) under salt and drought stress treatments. We believe these results may provide valuable information about the function of peanut LRR-RLK genes for further analysis.

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

Role of paralogs in the sex-bias transcriptional and metabolic regulation of the brain-placental axis in mice

INTRODUCTION

Duplicated genes or paralogs play important roles in the adaptive function of eukaryotic genomes. Animal studies have shown evidence for the functional role of paralogs in pregnancy, but our knowledge about the role of paralogs in the fetoplacental regulation remains limited. In particular, if fetoplacental metabolic regulation is modulated by differential expression of paralogs remains unexamined.

METHODS

In this study, gene expression profiles of day-15 placenta and fetal brain were compared to identify families or groups of paralogous genes expressed in the placenta and brain of male versus female fetuses in mice. A Bayesian modeling was applied to infer directional relationship of transcriptional variation of the paralogs relative to the phylogenetic variation of the genes in each family. Gas chromatography-mass spectrometry (GC-MS) was used to perform untargeted metabolomics analysis of day-15 placenta and fetal brain of both sexes.

RESULTS

We identified paralog groups that were expressed in a sex and/or tissue biased manner between the placenta and fetal brain. Bayesian modeling showed evidence for directional relationship between expression and phylogeny of specific paralogs. These relationships were sex specific. GC-MS analysis identified metabolites that were expressed in a sex-bias manner between the placenta and fetal brain. By performing integrative analysis of the metabolomics and gene expression data, we showed that specific groups of metabolites and paralogous genes were expressed in a coordinated manner between the placenta and fetal brain.

DISCUSSION

The findings of this study collectively suggest that paralogs play an influential role in the regulation of the brain-placental axis in mice.

Expectations of duplicate gene retention under the gene duplicability hypothesis

BACKGROUND

Gene duplication is an important process in evolution. What causes some genes to be retained after duplication and others to be lost is a process not well understood. The most prevalent theory is the gene duplicability hypothesis, that something about the function and number of interacting partners (number of subunits of protein complex, etc.), determines whether copies have more opportunity to be retained for long evolutionary periods. Some genes are also more susceptible to dosage balance effects following WGD events, making them more likely to be retained for longer periods of time. One would expect these processes that affect the retention of duplicate copies to affect the conditional probability ratio after consecutive whole genome duplication events. The probability that a gene will be retained after a second whole genome duplication event (WGD2), given that it was retained after the first whole genome duplication event (WGD1) versus the probability a gene will be retained after WGD2, given it was lost after WGD1 defines the probability ratio that is calculated.

RESULTS

Since duplicate gene retention is a time heterogeneous process, the time between the events (t1) and the time since the most recent event (t2) are relevant factors in calculating the expectation for observation in any genome. Here, we use a survival analysis framework to predict the probability ratio for genomes with different values of t1 and t2 under the gene duplicability hypothesis, that some genes are more susceptible to selectable functional shifts, some more susceptible to dosage compensation, and others only drifting. We also predict the probability ratio with different values of t1 and t2 under the mutational opportunity hypothesis, that probability of retention for certain genes changes in subsequent events depending upon how they were previously retained. These models are nested such that the mutational opportunity model encompasses the gene duplicability model with shifting duplicability over time. Here we present a formalization of the gene duplicability and mutational opportunity hypotheses to characterize evolutionary dynamics and explanatory power in a recently developed statistical framework.

CONCLUSIONS

This work presents expectations of the gene duplicability and mutational opportunity hypotheses over time under different sets of assumptions. This expectation will enable formal testing of processes leading to duplicate gene retention.

The structured coalescent in the context of gene copy number variation

The Structured Coalescent was introduced to describe the coalescent process in spatially subdivided populations with migration. Here, we re-interpret migration routes of individuals in the original model as "migration routes" of single genes in tandemly arranged gene arrays. A gene copy may change its position within the array via unequal recombination. Hence, in a coalescent framework, two copies sampled from two chromosomes may coalesce only if they are at exactly homologous positions. Otherwise, one or multiple recombination events have to occur before they can coalesce, thereby increasing mean coalescence time and expected genetic diversity among the copies in a gene array. We explicitly calculate the transition probabilities on these routes backward in time. We simulate the structured coalescent with migration and coalescence rates informed by the unequal recombination process of gene copies. With this novel interpretation of population structure models we determine coalescence times and expected genetic diversity in samples of orthologous and paralogous copies from a gene family. As a case study, we discuss the site frequency spectrum of a small gene family in the two scenarios of high and of no gene copy number variation among individuals. These examples underline the significance of our model, since standard test-statistics may lead to misinterpretations when analyzing sequence data of multi-copy genes due to their different expected genetic diversity.

Unusual Evolution of Cephalopod Tryptophan Indole-Lyases

Tryptophan indole-lyase (TIL), a pyridoxal-5-phosphate-dependent enzyme, catalyzes the hydrolysis of L-tryptophan (L-Trp) to indole and ammonium pyruvate. TIL is widely distributed among bacteria and bacterial TILs consist of a D2-symmetric homotetramer. On the other hand, TIL genes are also present in several metazoans. Cephalopods have two TILs, TILα and TILβ, which are believed to be derived from a gene duplication that occurred before octopus and squid diverged. However, both TILα and TILβ individually contain disruptive amino acid substitutions for TIL activity, and neither was active when expressed alone. When TILα and TILβ were coexpressed, however, they formed a heterotetramer that exhibited low TIL activity. The loss of TIL activity of the heterotetramer following site-directed mutagenesis strongly suggests that the active heterotetramer contains the TILα/TILβ heterodimer. Metazoan TILs generally have lower kcat values for L-Trp than those of bacterial TILs, but such low TIL activity may be rather suitable for metazoan physiology, where L-Trp is in high demand. Therefore, reduced activity may have been a less likely target for purifying selection in the evolution of cephalopod TILs. Meanwhile, the unusual evolution of cephalopod TILs may indicate the difficulty of post-gene duplication evolution of enzymes with catalytic sites contributed by multiple subunits, such as TIL.

Advances in genome sequencing reveal changes in gene content that contribute to arthropod macroevolution

Current sequencing technology allows for the relatively affordable generation of highly contiguous genomes. Technological advances have made it possible for researchers to investigate the consequences of diverse sorts of genomic variants, such as gene gain and loss. With the extraordinary number of high-quality genomes now available, we take stock of how these genomic variants impact phenotypic evolution. We take care to point out that the identification of genomic variants of interest is only the first step in understanding their impact. Painstaking lab or fieldwork is still required to establish causal relationships between genomic variants and phenotypic evolution. We focus mostly on arthropod research, as this phylum has an impressive degree of phenotypic diversity and is also the subject of much evolutionary genetics research. This article is intended to both highlight recent advances in the field and also to be a primer for learning about evolutionary genetics and genomics.

IS26 mediated blaCTX-M-65 amplification in Escherichia coli increase the antibiotic resistance to cephalosporin in vivo

OBJECTIVES

To characterize two Escherichia coli strains isolated from a patient pre- and post-treatment, using β-lactams and β-lactam/β-lactamase inhibitor combinations (BLBLIs).

METHODS

A combination of antibiotic susceptibility testing (AST) with whole genome sequencing using Illumina and Oxford Nanopore platforms. Long-read sequencing and reverse transcription-quantitative PCR were performed to determine the copy numbers and expression levels of antibiotic resistance genes (ARGs), respectively. Effect on fitness costs were assessed by growth rate determination.

RESULTS

The strain obtained from the patient after the antibiotic treatment (XH989) exhibited higher resistance to cefepime, BLBLIs and quinolones compared with the pre-treatment strain (XH987). Sequencing revealed IS26-mediated duplications of a IS26-fosA3-blaCTX-M-65 plasmid-embedded element in strain XH989. Long-read sequencing (7.4 G data volume) indicated a variation in copy numbers of blaCTX-M-65 within one single culture of strain XH989. Increased copy numbers of the IS26-fosA3-blaCTX-M-65 element were correlated with higher CTX-M-65 expression level and did not impose fitness costs, while facilitating faster growth under high antibiotic concentrations.

CONCLUSION

Our study is an example from the clinic how BLBLIs and β-lactams exposure in vivo possibly promoted the amplification of an IS26-multiple drug resistance (MDR) region. The observation of a copy number variation seen with the blaCTX-M-65 gene in the plasmid of the post-treatment strain expands our knowledge of insertion sequence dynamics and evolution during treatment.

Genome distance and phylogenetic inference accommodating gene duplication, loss and new gene input

With the rapid growth of entire genome data, phylogenomics focuses on analyzing evolutionary histories and relationships of species, i.e., the tree of life. For decades it has been realized that the genome-wide phylogenetic inference can be approached based upon the dynamic pattern of gene content (the presence/absence of gene families), or extended gene content (absence, presence as a single-copy, or duplicates). Those methods, conceptually or technically, invoked the birth-and-death process to model the evolutionary process (gene duplication or gene loss. One common drawback is that the mechanism of new gene input, including de novo origin of new genes and the lateral gene transfer, has not been explicitly considered. In this paper, the author developed a new genome distance approach for genome phylogeny inference under the origin-birth-death stochastic process. The model takes gene duplication, gene loss and new gene input into account simultaneously. Computer simulations found that the two-genome approach is statistically difficult to distinguish between two proliferation parameters, i.e., the rate of gene duplication and the rate of new gene input. Nevertheless, it has also demonstrated the statistical feasibility for using the loss-genome distance to infer the genome phylogeny, which can avoid the large sampling problem. The strategy to study the universal tree of life was discussed and exemplified by an example.

Genomic and transcriptomic survey of bryozoan Hox and ParaHox genes with emphasis on phylactolaemate bryozoans

BACKGROUND

Bryozoans are mostly sessile aquatic colonial invertebrates belonging to the clade Lophotrochozoa, which unites many protostome bilaterian phyla such as molluscs, annelids and brachiopods. While Hox and ParaHox genes have been extensively studied in various lophotrochozoan lineages, investigations on Hox and ParaHox gene complements in bryozoans are scarce.

RESULTS

Herein, we present the most comprehensive survey of Hox and ParaHox gene complements in bryozoans using four genomes and 35 transcriptomes representing all bryozoan clades: Cheilostomata, Ctenostomata, Cyclostomata and Phylactolaemata. Using similarity searches, phylogenetic analyses and detailed manual curation, we have identified five Hox genes in bryozoans (pb, Dfd, Lox5, Lox4 and Post2) and one ParaHox gene (Cdx). Interestingly, we observed lineage-specific duplication of certain Hox and ParaHox genes (Dfd, Lox5 and Cdx) in some bryozoan lineages.

CONCLUSIONS

The bryozoan Hox cluster does not retain the ancestral lophotrochozoan condition but appears relatively simple (includes only five genes) and broken into two genomic regions, characterized by the loss and duplication of serval genes. Importantly, bryozoans share the lack of two Hox genes (Post1 and Scr) with their proposed sister-taxon, Phoronida, which suggests that those genes were missing in the most common ancestor of bryozoans and phoronids.

Cooption of regulatory modules for tektin paralogs during ciliary band formation in a marine annelid larva

Tektins are a highly conserved family of coiled-coil domain containing proteins known to play a role in structure, stability and function of cilia and flagella. Tektin proteins are thought to form filaments which run the length of the axoneme along the inner surface of the A tubule of each microtubule doublet. Phylogenetic analyses suggest that the tektin family arose via duplications from a single tektin gene in a unicellular organism giving rise to four and five tektin genes in bilaterians and in spiralians, respectively. Although tektins are found in most metazoans, little is known about their expression and function outside of a handful of model species. Here we present the first comprehensive study of tektin family gene expression in any animal system, in the spiralian annelid Platynereis dumerilii. This indirect developing species retains a full ancient spiralian complement of five tektin genes. We show that all five tektins are expressed almost exclusively in known ciliary structures following the expression of the motile cilia master regulator foxJ1. The three older bilaterian tektin-1, tektin-2, and tektin-4 genes, show a high degree of spatial and temporal co-regulation, while the spiralian specific tektin-3/5A and tektin-3/5B show a delay in onset of expression in every ciliary structure. In addition, tektin-3/5B transcripts show a restricted subcellular localization to the most apical region near the multiciliary arrays. The exact recapitulation of the sequence of expression and localization of the five tektins at different times during larval development indicates the cooption of a fixed regulatory and cellular program during the formation of each ciliary band and multiciliated cell type in this spiralian.

The ribosomal protein P0A is required for embryo development in rice

BACKGROUND

The P-stalk is a conserved and vital structural element of ribosome. The eukaryotic P-stalk exists as a P0-(P1-P2)2 pentameric complex, in which P0 function as a base structure for incorporating the stalk onto 60S pre-ribosome. Prior studies have suggested that P0 genes are indispensable for survival in yeast and animals. However, the functions of P0 genes in plants remain elusive.

RESULTS

In the present study, we show that rice has three P0 genes predicted to encode highly conserved proteins OsP0A, OsP0B and OsP0C. All of these P0 proteins were localized both in cytoplasm and nucleus, and all interacted with OsP1. Intriguingly, the transcripts of OsP0A presented more than 90% of the total P0 transcripts. Moreover, knockout of OsP0A led to embryo lethality, while single or double knockout of OsP0B and OsP0C did not show any visible defects in rice. The genomic DNA of OsP0A could well complement the lethal phenotypes of osp0a mutant. Finally, sequence and syntenic analyses revealed that OsP0C evolved from OsP0A, and that duplication of genomic fragment harboring OsP0C further gave birth to OsP0B, and both of these duplication events might happen prior to the differentiation of indica and japonica subspecies in rice ancestor.

CONCLUSION

These data suggested that OsP0A functions as the predominant P0 gene, playing an essential role in embryo development in rice. Our findings highlighted the importance of P0 genes in plant development.

Systematic Analysis of Diverse Polynucleotide Kinase Clp1 Family Proteins in Eukaryotes: Three Unique Clp1 Proteins of Trypanosoma brucei

The Clp1 family proteins, consisting of the Clp1 and Nol9/Grc3 groups, have polynucleotide kinase (PNK) activity at the 5' end of RNA strands and are important enzymes in the processing of some precursor RNAs. However, it remains unclear how this enzyme family diversified in the eukaryotes. We performed a large-scale molecular evolutionary analysis of the full-length genomes of 358 eukaryotic species to classify the diverse Clp1 family proteins. The average number of Clp1 family proteins in eukaryotes was 2.3 ± 1.0, and most representative species had both Clp1 and Nol9/Grc3 proteins, suggesting that the Clp1 and Nol9/Grc3 groups were already formed in the eukaryotic ancestor by gene duplication. We also detected an average of 4.1 ± 0.4 Clp1 family proteins in members of the protist phylum Euglenozoa. For example, in Trypanosoma brucei, there are three genes of the Clp1 group and one gene of the Nol9/Grc3 group. In the Clp1 group proteins encoded by these three genes, the C-terminal domains have been replaced by unique characteristics domains, so we designated these proteins Tb-Clp1-t1, Tb-Clp1-t2, and Tb-Clp1-t3. Experimental validation showed that only Tb-Clp1-t2 has PNK activity against RNA strands. As in this example, N-terminal and C-terminal domain replacement also contributed to the diversification of the Clp1 family proteins in other eukaryotic species. Our analysis also revealed that the Clp1 family proteins in humans and plants diversified through isoforms created by alternative splicing.

Evolution of antimicrobial cysteine-rich peptides in plants

KEY MESSAGE

We analyzed the evolutionary pattern of cysteine-rich peptides (CRPs) to infer the relationship between CRP copy number and plant ecotype, and the origin of bi-domains CRPs. Plants produce cysteine-rich peptides (CRPs) that have long-lasting broad-spectrum antimicrobial activity to protect themselves from various groups of pathogens. We analyzed 240 plant genomes, ranging from algae to eudicots, and discovered that CRPs are widely distributed in plants. Our comparative genomics results revealed that CRP genes have been amplified through both whole genome and local tandem duplication. The copy number of these genes varied significantly across lineages and was associated with the plant ecotype. This may be due to their resistance to changing pathogenic environments. The conserved and lineage-specific CRP families contribute to diverse antimicrobial activities. Furthermore, we investigated the unique bi-domain CRPs that result from unequal crossover events. Our findings provide a unique evolutionary perspective on CRPs and insights into their antimicrobial and symbiosis characteristics.

Genome-wide analysis of the GDSL esterase/lipase family genes in Physcomitrium patens and the involvement of GELP31 in spore germination

In plants, the ability to produce hydrophobic substances that would provide protection from dehydration was required for the transition to land. This genome-wide investigation outlines the evolution of GDSL-type esterase/lipase (GELP) proteins in the moss Physcomitrium patens and suggests possible functions of some genes. GELP proteins play roles in the formation of hydrophobic polymers such as cutin and suberin that protect against dehydration and pathogen attack. GELP proteins are also implicated in processes such as pollen development and seed metabolism and germination. The P. patens GELP gene family comprises 48 genes and 14 pseudogenes. Phylogenetic analysis of all P. patens GELP sequences along with vascular plant GELP proteins with reported functions revealed that the P. patens genes clustered within previously identified A, B and C clades. A duplication model predicting the expansion of the GELP gene family within the P. patens lineage was constructed. Expression analysis combined with phylogenetic analysis suggested candidate genes for functions such as defence against pathogens, cutin metabolism, spore development and spore germination. The presence of relatively fewer GELP genes in P. patens may reduce the occurrence of functional redundancy that complicates the characterization of vascular plant GELP genes. Knockout lines of GELP31, which is highly expressed in sporophytes, were constructed. Gelp31 spores contained amorphous oil bodies and germinated late, suggesting (a) role(s) of GELP31 in lipid metabolism in spore development or germination. Future knockout studies of other candidate GELP genes will further elucidate the relationship between expansion of the family and the ability to withstand the harsh land environment.

Genome-wide identification and characterization of glutathione S-transferase gene family in quinoa (Chenopodium quinoa Willd.)

The present investigation was envisaged for large scale in-silico genome wide identification and characterization of glutathione S-transferases (GSTs) in Chenopodium quinoa. In this study, a total of 120 GST genes (CqGSTs) were identified and divided into 11 classes of which tau and phi were highest in numbers. The average protein length of protein was found to be 279.06 with their corresponding average molecular weight of 31,819.4 kDa. The subcellular localization analysis results showed that proteins were centrally localized in the cytoplasm followed by chloroplast, mitochondria and plastids. Structural analysis revealed the presence of 2 -14 exons in CqGST genes. Most of the proteins possessed two exon one intron organization. MEME analysis identified 15 significantly conserved motifs with a width of 6-50 amino acids. Motifs 1, 3, 2, 5, 6, 8, 9 and 13 were found specifically in tau class family; motifs 3, 4, 5, 6, 7 and 9 were found in phi class gene family, while motifs 3, 4, 13 and 14 were found in metaxin class. Multiple sequence alignment revealed highly conserved N-terminus with active site serine (Ser; S) or cysteine (Cys; C) residue for the activation of GSH binding and GST catalytic activity. The gene loci were found to be unevenly distributed across 18 different chromosomes with a maximum of 17 genes located on chromosome number 7. Dominance of alpha helix was followed by coil, extended strand and beta turns. Gene duplication analysis revealed that segmental duplication and purifying type selection were highest in number and found to be main source of expansion of GST gene family. Cis acting regulatory elements analysis showed the presence of 21 different elements involved in stress, hormone and light response and cellular development. The evolutionary relationship of CqGST proteins carried out using maximum likelihood method revealed that all the tau and phi class GSTs were closely associated with those of G. max, O. sativa and A. thaliana. Molecular docking of GST molecules with the fungicide metalaxyl showed that the CqGSTF1 had the lowest binding energy. The comprehensive study of CqGST gene family in quinoa provides groundwork for further functional analysis of CqGST genes in the species at molecular level and has potential applications in plant breeding.

The power of "controllers": Transposon-mediated duplicated genes evolve towards neofunctionalization

Since the discovery of the first transposon by Dr. Barbara McClintock, the prevalence and diversity of transposable elements (TEs) have been gradually recognized. As fundamental genetic components, TEs drive organismal evolution not only by contributing functional sequences (e.g., regulatory elements or "controllers" as phrased by Dr. McClintock) but also by shuffling genomic sequences. In the latter respect, TE-mediated gene duplications have contributed to the origination of new genes and attracted extensive interest. In response to the development of this field, we herein attempt to provide an overview of TE-mediated duplication by focusing on common rules emerging across duplications generated by different TE types. Specifically, despite the huge divergence of transposition machinery across TEs, we identify three common features of various TE-mediated duplication mechanisms, including end bypass, template switching, and recurrent transposition. These three features lead to one common functional outcome, namely, TE-mediated duplicates tend to be subjected to exon shuffling and neofunctionalization. Therefore, the intrinsic properties of the mutational mechanism constrain the evolutionary trajectories of these duplicates. We finally discuss the future of this field including an in-depth characterization of both the duplication mechanisms and functions of TE-mediated duplicates.

Duplicated zebrafish (Danio rerio) inositol phosphatases inpp5ka and inpp5kb diverged in expression pattern and function

One hurdle in the development of zebrafish models of human disease is the presence of multiple zebrafish orthologs resulting from whole genome duplication in teleosts. Mutations in inositol polyphosphate 5-phosphatase K (INPP5K) lead to a syndrome characterized by variable presentation of intellectual disability, brain abnormalities, cataracts, muscle disease, and short stature. INPP5K is a phosphatase acting at position 5 of phosphoinositides to control their homeostasis and is involved in insulin signaling, cytoskeletal regulation, and protein trafficking. Previously, our group and others have replicated the human phenotypes in zebrafish knockdown models by targeting both INPP5K orthologs inpp5ka and inpp5kb. Here, we show that inpp5ka is the more closely related orthologue to human INPP5K. While both inpp5ka and inpp5kb mRNA expression levels follow a similar trend in the developing head, eyes, and tail, inpp5ka is much more abundantly expressed in these tissues than inpp5kb. In situ hybridization revealed a similar trend, also showing unique localization of inpp5kb in the pineal gland and retina indicating different transcriptional regulation. We also found that inpp5kb has lost its catalytic activity against its preferred substrate, PtdIns(4,5)P2. Since most human mutations are missense changes disrupting phosphatase activity, we propose that loss of inpp5ka alone can be targeted to recapitulate the human presentation. In addition, we show that the function of inpp5kb has diverged from inpp5ka and may play a novel role in the zebrafish.

Extensive ITR expansion of the 2022 Mpox virus genome through gene duplication and gene loss

Poxviruses are known to evolve slower than RNA viruses with only 1-2 mutations/genome/year. Rather than single mutations, rearrangements such as gene gain and loss, which have been discussed as a possible driver for host adaption, were described in poxviruses. In 2022 and 2023 the world is being challenged by the largest global outbreak so far of Mpox virus, and the virus seems to have established itself in the human community for an extended period of time. Here, we report five Mpox virus genomes from Germany with extensive gene duplication and loss, leading to the expansion of the ITR regions from 6400 to up to 24,600 bp. We describe duplications of up to 18,200 bp to the opposed genome end, and deletions at the site of insertion of up to 16,900 bp. Deletions and duplications of genes with functions of supposed immune modulation, virulence and host adaption as B19R, B21R, B22R and D10L are described. In summary, we highlight the need for monitoring rearrangements of the Mpox virus genome rather than for monitoring single mutations only.

Phylogenetic and Evolutionary Analysis of Plant Small RNA 2'-O-Methyltransferase (HEN1) Protein Family

HUA ENHANCER 1 (HEN1) is a pivotal mediator in protecting sRNAs from 3'-end uridylation and 3' to 5' exonuclease-mediated degradation in plants. Here, we investigated the pattern of the HEN1 protein family evolutionary history and possible relationships in the plant lineages using protein sequence analyses and conserved motifs composition, functional domain identification, architecture, and phylogenetic tree reconstruction and evolutionary history inference. According to our results, HEN1 protein sequences bear several highly conserved motifs in plant species retained during the evolution from their ancestor. However, several motifs are present only in Gymnosperms and Angiosperms. A similar trend showed for their domain architecture. At the same time, phylogenetic analysis revealed the grouping of the HEN1 proteins in the three main super clads. In addition, the Neighbor-net network analysis result provides some nodes have multiple parents indicating a few conflicting signals in the data, which is not the consequence of sampling error, the effect of the selected model, or the estimation method. By reconciling the protein and species tree, we considered the gene duplications in several given species and found 170 duplication events in the evolution of HEN1 in the plant lineages. According to our analysis, the main HEN1 superclass mostly showed orthologous sequences that illustrate the vertically transmitting of HEN1 to the main lines. However, in both orthologous and paralogs, we predicted insignificant structural deviations. Our analysis implies that small local structural changes that occur continuously during the folds can moderate the changes created in the sequence. According to our results, we proposed a hypothetical model and evolutionary trajectory for the HEN1 protein family in the plant kingdom.

Divergent expression of aristaless1 and aristaless2 during embryonic appendage and pupal wing development in butterflies

BACKGROUND

Gene duplication events are critical for the evolution of new gene functions. Aristaless is a major regulator of distinct developmental processes. It is most known for its role during appendage development across animals. However, more recently other distinct biological functions have been described for this gene and its duplicates. Butterflies and moths have two copies of aristaless, aristaless1 (al1) and aristaless2 (al2), as a result of a gene duplication event. Previous work in Heliconius has shown that both copies appear to have novel functions related to wing color patterning. Here we expand our knowledge of the expression profiles associated with both ancestral and novel functions of Al1 across embryogenesis and wing pigmentation. Furthermore, we characterize Al2 expression, providing a comparative framework between gene copies within the same species, allowing us to understand the origin of new functions following gene duplication.

RESULTS

Our work shows that the expression of both Al1 and Al2 is associated with the ancestral function of sensory appendage (leg, mouth, spines, and eyes) development in embryos. Interestingly, Al1 exhibits higher expression earlier in embryogenesis while the highest levels of Al2 expression are shifted to later stages of embryonic development. Furthermore, Al1 localization appears extranuclear while Al2 co-localizes tightly with nuclei earlier, and then also expands outside the nucleus later in development. Cellular expression of Al1 and Al2 in pupal wings is broadly consistent with patterns observed during embryogenesis. We also describe, for the first time, how Al1 localization appears to correlate with zones of anterior/posterior elongation of the body during embryonic growth, showcasing a possible new function related to Aristaless' previously described role in appendage extension.

CONCLUSIONS

Overall, our data suggest that while both gene copies play a role in embryogenesis and wing pigmentation, the duplicates have diverged temporally and mechanistically across those functions. Our study helps clarify principles behind sub-functionalization and gene expression evolution associated with developmental functions following gene duplication events.

Identification, evolution and expression analyses of the whole genome-wide PEBP gene family in Brassica napus L

BACKGROUND

With the release of genomic data for B.rapa, B.oleracea, and B.napus, research on the genetic and molecular functions of Brassica spp. has entered a new stage. PEBP genes in plants play an important role in the transition to flowering as well as seed development and germination. Molecular evolutionary and functional analyses of the PEBP gene family in B.napus based on molecular biology methods can provide a theoretical basis for subsequent investigations of related regulators.

RESULTS

In this paper, we identified a total of 29 PEBP genes from B.napus that were located on 14 chromosomes and 3 random locations. Most members contained 4 exons and 3 introns; motif 1 and motif 2 were the characteristic motifs of PEBP members. On the basis of intraspecific and interspecific collinearity analyses, it is speculated that fragment replication and genomic replication are the main drivers of for the amplification and evolution of the PEBP gene in the B.napus genome. The results of promoter cis-elements prediction suggest that BnPEBP family genes are inducible promoters, which may directly or indirectly participate in multiple regulatory pathways of plant growth cycle. Furthermore, the tissue-specific expression results show that the expression levels of BnPEBP family genes in different tissues were quite different, but the gene expression organization and patterns of the same subgroup were basically the same. qRT‒PCR revealed certain spatiotemporal patterns in the expression of the PEBP subgroups in roots, stems, leaves, buds, and siliques, was tissue-specific, and related to function.

CONCLUSIONS

A systematic comparative analysis of the B.napus PEBP gene family was carried out at here. The results of gene identification, phylogenetic tree construction, structural analysis, gene duplication analysis, prediction of promoter cis-elements and interacting proteins, and expression analysis provide a reference for exploring the molecular mechanisms of BnPEBP family genes in future research.

Four ciliate-specific expansion events occurred during actin gene family evolution of eukaryotes

Actin gene family is a divergent and ancient eukaryotic cellular cytoskeletal gene family, and participates in many essential cellular processes. Ciliated protists offer us an excellent opportunity to investigate gene family evolution, since their gene families evolved faster in ciliates than in other eukaryotes. Nonetheless, actin gene family is well studied in few model ciliate species but little is known about its evolutionary patterns in ciliates. Here, we analyzed the evolutionary pattern of eukaryotic actin gene family based on genomes/transcriptomes of 36 species covering ten ciliate classes, as well as those of nine non-ciliate eukaryotic species. Results showed: (1) Except for conventional actins and actin-related proteins (Arps) shared by various eukaryotes, at least four ciliate-specific subfamilies occurred during evolution of actin gene family. Expansions of Act2 and ArpC were supposed to have happen in the ciliate common ancestor, while expansions of ActI and ActII may have occurred in the ancestor of Armophorea, Muranotrichea, and Spirotrichea. (2) The number of actin isoforms varied greatly among ciliate species. Environmental adaptability, whole genome duplication (WGD) or segmental duplication events, distinct spatial and temporal patterns of expression might play driving forces for the increasement of isoform numbers. (3) The 'birth and death' model of evolution could explain the evolution of actin gene family in ciliates. And actin genes have been generally under strong negative selection to maintain protein structures and physiological functions. Collectively, we provided meaningful information for understanding the evolution of eukaryotic actin gene family.

Promoter evolution of mammalian gene duplicates

BACKGROUND

Gene duplication is thought to be a central process in evolution to gain new functions. The factors that dictate gene retention following duplication as well paralog gene divergence in sequence, expression and function have been extensively studied. However, relatively little is known about the evolution of promoter regions of gene duplicates and how they influence gene duplicate divergence. Here, we focus on promoters of paralog genes, comparing their similarity in sequence, in the sets of transcription factors (TFs) that bind them, and in their overall promoter architecture.

RESULTS

We observe that promoters of recent duplications display higher sequence similarity between them and that sequence similarity rapidly declines between promoters of more ancient paralogs. In contrast, similarity in cis-regulation, as measured by the set of TFs that bind promoters of both paralogs, does not simply decrease with time from duplication and is instead related to promoter architecture-paralogs with CpG Islands (CGIs) in their promoters share a greater fraction of TFs, while CGI-less paralogs are more divergent in their TF binding set. Focusing on recent duplication events and partitioning them by their duplication mechanism enables us to uncover promoter properties associated with gene retention, as well as to characterize the evolution of promoters of newly born genes: In recent retrotransposition-mediated duplications, we observe asymmetry in cis-regulation of paralog pairs: Retrocopy genes are lowly expressed and their promoters are bound by fewer TFs and are depleted of CGIs, in comparison with the original gene copy. Furthermore, looking at recent segmental duplication regions in primates enable us to compare successful retentions versus loss of duplicates, showing that duplicate retention is associated with fewer TFs and with CGI-less promoter architecture.

CONCLUSIONS

In this work, we profiled promoters of gene duplicates and their inter-paralog divergence. We also studied how their characteristics are associated with duplication time and duplication mechanism, as well as with the fate of these duplicates. These results underline the importance of cis-regulatory mechanisms in shaping the evolution of new genes and their fate following duplication.

The enhanced Valsa canker resistance conferred by MdLecRK-S.4.3 in Pyrus betulifolia can be largely suppressed by PbePUB36

L-type lectin receptor-like kinases (L-LecRKs) act as a sensor of extracellular signals and an initiator for plant immune responses. However, the function of LecRK-S.4 on plant immunity has not been extensively investigated. At present, in the apple (Malus domestica) genome, we identified that MdLecRK-S.4.3, a homologous gene of LecRK-S.4, was differentially expressed during the occursion of Valsa canker. Over-expression of MdLecRK-S.4.3 facilitated the induction of immune response and enhanced the Valsa canker resistance of apple and pear fruit, and 'Duli-G03' (Pyrus betulifolia) suspension cells. On the contrary, the expression of PbePUB36, RLCK XI subfamily member, was significantly repressed in the MdLecRK-S.4.3 overexpressed cell lines. Over-expression of PbePUB36 interfered with the Valsa canker resistance and immune response caused by up-regulation of MdLecRK-S.4.3. Furthermore, MdLecRK-S.4.3 interacted with BAK1 or PbePUB36 in vivo. In conclusion, MdLecRK-S.4.3 activated various immune responses and positively regulate Valsa canker resistance, which could be largely compromised by PbePUB36. MdLecRK-S.4.3 interacted with PbePUB36 and/or MdBAK1 to mediate the immune responses. This finding provides a reference for studying the molecular mechanism of resistance to Valsa canker and resistance breeding.

Slit-Robo expression in the leech nervous system: insights into eyespot evolution

BACKGROUND

Slit and Robo are evolutionarily conserved ligand and receptor proteins, respectively, but the number of slit and robo gene paralogs varies across recent bilaterian genomes. Previous studies indicate that this ligand-receptor complex is involved in axon guidance. Given the lack of data regarding Slit/Robo in the Lophotrochozoa compared to Ecdysozoa and Deuterostomia, the present study aims to identify and characterize the expression of Slit/Robo orthologs in leech development.

RESULTS

We identified one slit (Hau-slit), and two robo genes (Hau-robo1 and Hau-robo2), and characterized their expression spatiotemporally during the development of the glossiphoniid leech Helobdella austinensis. Throughout segmentation and organogenesis, Hau-slit and Hau-robo1 are broadly expressed in complex and roughly complementary patterns in the ventral and dorsal midline, nerve ganglia, foregut, visceral mesoderm and/or endoderm of the crop, rectum and reproductive organs. Before yolk exhaustion, Hau-robo1 is also expressed where the pigmented eye spots will later develop, and Hau-slit is expressed in the area between these future eye spots. In contrast, Hau-robo2 expression is extremely limited, appearing first in the developing pigmented eye spots, and later in the three additional pairs of cryptic eye spots in head region that never develop pigment. Comparing the expression of robo orthologs between H. austinensis and another glossiphoniid leech, Alboglossiphonia lata allows to that robo1 and robo2 operate combinatorially to differentially specify pigmented and cryptic eyespots within the glossiphoniid leeches.

CONCLUSIONS

Our results support a conserved role in neurogenesis, midline formation and eye spot development for Slit/Robo in the Lophotrochozoa, and provide relevant data for evo-devo studies related to nervous system evolution.

Genetic architecture of human brain evolution

Comparative studies of hominids have long sought to identify mutational events that shaped the evolution of the human nervous system. However, functional genetic differences are outnumbered by millions of nearly neutral mutations, and the developmental mechanisms underlying human nervous system specializations are difficult to model and incompletely understood. Candidate-gene studies have attempted to map select human-specific genetic differences to neurodevelopmental functions, but it remains unclear how to contextualize the relative effects of genes that are investigated independently. Considering these limitations, we discuss scalable approaches for probing the functional contributions of human-specific genetic differences. We propose that a systems-level view will enable a more quantitative and integrative understanding of the genetic, molecular and cellular underpinnings of human nervous system evolution.

Extensive sequence duplication in Arabidopsis revealed by pseudo-heterozygosity

BACKGROUND

It is apparent that genomes harbor much structural variation that is largely undetected for technical reasons. Such variation can cause artifacts when short-read sequencing data are mapped to a reference genome. Spurious SNPs may result from mapping of reads to unrecognized duplicated regions. Calling SNP using the raw reads of the 1001 Arabidopsis Genomes Project we identified 3.3 million (44%) heterozygous SNPs. Given that Arabidopsis thaliana (A. thaliana) is highly selfing, and that extensively heterozygous individuals have been removed, we hypothesize that these SNPs reflected cryptic copy number variation.

RESULTS

The heterozygosity we observe consists of particular SNPs being heterozygous across individuals in a manner that strongly suggests it reflects shared segregating duplications rather than random tracts of residual heterozygosity due to occasional outcrossing. Focusing on such pseudo-heterozygosity in annotated genes, we use genome-wide association to map the position of the duplicates. We identify 2500 putatively duplicated genes and validate them using de novo genome assemblies from six lines. Specific examples included an annotated gene and nearby transposon that transpose together. We also demonstrate that cryptic structural variation produces highly inaccurate estimates of DNA methylation polymorphism.

CONCLUSIONS

Our study confirms that most heterozygous SNP calls in A. thaliana are artifacts and suggest that great caution is needed when analyzing SNP data from short-read sequencing. The finding that 10% of annotated genes exhibit copy-number variation, and the realization that neither gene- nor transposon-annotation necessarily tells us what is actually mobile in the genome suggests that future analyses based on independently assembled genomes will be very informative.

The evolution of the green-light-sensitive visual opsin genes (RH2) in teleost fishes

Vertebrates have four visual cone opsin classes that mediate sensitivity from ultraviolet to red wavelengths of light. The rhodopsin-like 2 (RH2) opsin is sensitive to the central mostly green part of the spectrum. While lost in some terrestrial vertebrates (mammals), the RH2 opsin gene has proliferated during the evolution of teleost fishes. Here, we investigated the genomes of 132 extant teleosts and found between zero and eight RH2 gene copies per species. The RH2 gene shows a dynamic evolutionary history with repeated gene duplications, gene losses, and gene conversions affecting entire orders, families, and species. At least four ancestral duplications provided the substrate for today's RH2 diversity, with duplications occurring in the common ancestors of Clupeocephala (twice), Neoteleostei, and likely Acanthopterygii as well. Despite these evolutionary dynamics, we identified conserved RH2 synteny in two main gene clusters; the slc6A13/synpr cluster is highly conserved within Percomorpha and also present across most teleosts, including Otomorpha, Euteleostei and in parts in tarpons (Elopomorpha), and the mutSH5 cluster, which is specific for Otomorpha. When comparing the number of visual opsin genes (SWS1, SWS2, RH2, LWS, and total cone opsins) with habitat depth, we found that deeper-dwelling species had less (or none) long-wavelength-sensitive opsins. Using retinal/eye transcriptomes in a phylogenetic representative dataset of 32 species, we show that if present in the genome, RH2 is expressed in most fishes except for some species within the tarpons, characins, and gobies (and Osteoglossomorpha and some other characin species have lost the gene). Those species instead express a green-shifted long-wavelength-sensitive LWS opsin. Our study applies modern genomic and transcriptomic tools within a comparative framework to elucidate the evolutionary history of the visual sensory system in teleost fishes.

Evolution of proteins involved in the final steps of juvenile hormone synthesis

Juvenile hormone (JH), a sesquiterpenoid produced by the insect corpus allatum gland (CA), is a key regulator of insect metamorphosis, reproduction, caste differentiation, and polyphenism. The first part of JH biosynthesis occurs via the universal eukaryotic mevalonate pathway. The final steps involve epoxidation and methylation. However, the sequence of these steps might not be conserved among all insects and Crustacea. Therefore, we used available genomic and transcriptomic data and identified JH acid methyltransferase (JHAMT), analyzed their genomic duplications in selected model organisms, and reconstructed their phylogeny. We have further reconstructed phylogeny of FAMeT proteins and show that evolution of this protein group is more complicated than originally appreciated. The analysis delineates important milestones in the evolution of several JH biosynthetic enzymes in arthropods, reviews major literature data on the last steps of JH synthesis, and defines questions and some hypotheses worth pursuing experimentally.

Identification of a second 16-hydroxytabersonine-O-methyltransferase suggests an evolutionary relationship between alkaloid and flavonoid metabolisms in Catharanthus roseus

The medicinal plant Catharanthus roseus biosynthesizes many important drugs for human health, including the anticancer monoterpene indole alkaloids (MIAs) vinblastine and vincristine. Over the past decades, the continuous increase in pharmaceutical demand has prompted several research groups to characterize MIA biosynthetic pathways for considering future metabolic engineering processes of supply. In line with previous work suggesting that diversification can potentially occur at various steps along the vindoline branch, we were here interested in investigating the involvement of distinct isoforms of tabersonine-16-O-methyltransferase (16OMT) which plays a pivotal role in the MIA biosynthetic pathway. By combining homology searches based on the previously characterized 16OMT1, phylogenetic analyses, functional assays in yeast, and biochemical and in planta characterizations, we identified a second isoform of 16OMT, referred to as 16OMT2. 16OMT2 appears to be a multifunctional enzyme working on both MIA and flavonoid substrates, suggesting that a constrained evolution of the enzyme for accommodating the MIA substrate has probably occurred to favor the apparition of 16OMT2 from an ancestral specific flavonoid-O-methyltransferase. Since 16OMT1 and 16OMT2 displays a high sequence identity and similar kinetic parameters for 16-hydroxytabersonine, we postulate that 16OMT1 may result from a later 16OMT2 gene duplication accompanied by a continuous neofunctionalization leading to an almost complete loss of flavonoid O-methyltransferase activity. Overall, these results participate in increasing our knowledge on the evolutionary processes that have likely led to enzyme co-optation for MIA synthesis.

Comparative gene retention analysis in barley, wild emmer, and bread wheat pangenome lines reveals factors affecting gene retention following gene duplication

BACKGROUND

Gene duplication is a prevalent phenomenon and a major driving force underlying genome evolution. The process leading to the fixation of gene duplicates following duplication is critical to understand how genome evolves but remains fragmentally understood. Most previous studies on gene retention are based on gene duplicate analyses in single reference genome. No population-based comparative gene retention analysis has been performed to date.

RESULTS

Taking advantage of recently published genomic data in Triticeae, we dissected a divergent homogentisate phytyltransferase (HPT2) lineage caught in the middle stage of gene fixation following duplication. The presence/absence of HPT2 in barley (diploid), wild emmer (tetraploid), and bread wheat (hexaploid) pangenome lines appears to be associated with gene dosage constraint and environmental adaption. Based on these observations, we adopted a phylogeny-based orthology inference approach and performed comparative gene retention analyses across barley, wild emmer, and bread wheat. This led to the identification of 326 HPT2-pattern-like genes at whole genome scale, representing a pool of gene duplicates in the middle stage of gene fixation. Majority of these HPT2-pattern-like genes were identified as small-scale duplicates, such as dispersed, tandem, and proximal duplications. Natural selection analyses showed that HPT2-pattern-like genes have experienced relaxed selection pressure, which is generally accompanied with partial positive selection and transcriptional divergence. Functional enrichment analyses showed that HPT2-pattern-like genes are over-represented with molecular-binding and defense response functions, supporting the potential role of environmental adaption during gene retention. We also observed that gene duplicates from larger gene family are more likely to be lost, implying a gene dosage constraint effect. Further comparative gene retention analysis in barley and bread wheat pangenome lines revealed combined effects of species-specific selection and gene dosage constraint.

CONCLUSIONS

Comparative gene retention analyses at the population level support gene dosage constraint, environmental adaption, and species-specific selection as three factors that may affect gene retention following gene duplication. Our findings shed light on the evolutionary process leading to the retention of newly formed gene duplicates and will greatly improve our understanding on genome evolution via duplication.

The Fate of Duplicated Enzymes in Prokaryotes: The Case of Isomerases

The isomerases are a unique enzymatic class of enzymes that carry out a great diversity of chemical reactions at the intramolecular level. This class comprises about 300 members, most of which are involved in carbohydrate and terpenoid/polyketide metabolism. Along with oxidoreductases and translocases, isomerases are one of the classes with the highest ratio of paralogous enzymes. Due to its relatively small number of members, it is plausible to explore it in greater detail to identify specific cases of gene duplication. Here, we present an analysis at the level of individual isomerases and identify different members that seem to be involved in duplication events in prokaryotes. As was suggested in a previous study, there is no homogeneous distribution of paralogs, but rather they accumulate into a few subcategories, some of which differ between Archaea and Bacteria. As expected, the metabolic processes with more paralogous isomerases have to do with carbohydrate metabolism but also with RNA modification (a particular case involving an rRNA-modifying isomerase is thoroughly discussed and analyzed in detail). Overall, our findings suggest that the most common fate for paralogous enzymes is the retention of the original enzymatic function, either associated with a dosage effect or with differential expression in response to changing environments, followed by subfunctionalization and, to a much lesser degree, neofunctionalization, which is consistent with what has been reported elsewhere.

Decipher the molecular evolution and expression patterns of Cupin family genes in oilseed rape

Cupin proteins are involved in plant growth and development as well as in response to various stresses. Here, a total of 173 Cupin genes were identified in Brassica napus, and their molecular evolution and expression patterns were analyzed. These genes were classified into ten groups. Motif and exon-intron structure indicated a high degree of conservation within each group during evolution. BnaCupins were distributed on 19 chromosomes and their expansion is mainly contributed by whole-genome duplication (WGD) and segmental duplication events. BnaCupins have undergone severe purifying selection during a long evolutionary process. Meanwhile, some positive selection sites were identified. Expression patterns and cis-element analysis indicated that BnaCupins play significant roles in plant growth and stress responses. In addition, the expression levels of some BnCupins were significantly altered when treated with different conditions (cold, salt, drought, IAA, ABA, and 6-BA). Some BnaCupin interacting proteins, such as glycosyl hydrolase5 (GHs5), carbohydrate kinase (CHKs), ATP-dependent 6-phosphofructokinase (ATP-PFK), S-adenosylmethionine synthase (S-MAT), and aldolase class II (ALD II), were identified by the protein-protein interaction network. It will contribute to enriching our knowledge of the Cupin gene family in B. napus and provide a basis for further studies of their functions.

Genome-wide identification, evolutionary and expression analysis of the cyclin-dependent kinase gene family in peanut

BACKGROUND

Cyclin-dependent kinases (CDKs) are a predominant group of serine/threonine protein kinases that have multi-faceted functions in eukaryotes. The plant CDK members have well-known roles in cell cycle progression, transcriptional regulation, DNA repair, abiotic stress and defense responses, making them promising targets for developing stress adaptable high-yielding crops. There is relatively sparse information available on the CDK family genes of cultivated oilseed crop peanut and its diploid progenitors.

RESULTS

We have identified 52 putative cyclin-dependent kinases (CDKs) and CDK-like (CDKLs) genes in Arachis hypogaea (cultivated peanut) and total 26 genes in each diploid parent of cultivated peanut (Arachis duranensis and Arachis ipaensis). Both CDK and CDKL genes were classified into eight groups based on their cyclin binding motifs and their phylogenetic relationship with Arabidopsis counterparts. Genes in the same subgroup displayed similar exon-intron structure and conserved motifs. Further, gene duplication analysis suggested that segmental duplication events played major roles in the expansion and evolution of CDK and CDKL genes in cultivated peanuts. Identification of diverse cis-acting response elements in CDK and CDKL genes promoter indicated their potential fundamental roles in multiple biological processes. Various gene expression patterns of CDKs and CDKLs in different peanut tissues suggested their involvement during growth and development. In addition, qRT-PCR analysis demonstrated that most representing CDK and CDKL gene family members were significantly down-regulated under ABA, PEG and mannitol treatments.

CONCLUSIONS

Genome-wide analysis offers a comprehensive understanding of the classification, evolution, gene structure, and gene expression profiles of CDK and CDKL genes in cultivated peanut and their diploid progenitors. Additionally, it also provides cell cycle regulatory gene resources for further functional characterization to enhance growth, development and abiotic stress tolerance.

NANOGP1, a tandem duplicate of NANOG, exhibits partial functional conservation in human naïve pluripotent stem cells

Gene duplication events can drive evolution by providing genetic material for new gene functions, and they create opportunities for diverse developmental strategies to emerge between species. To study the contribution of duplicated genes to human early development, we examined the evolution and function of NANOGP1, a tandem duplicate of the transcription factor NANOG. We found that NANOGP1 and NANOG have overlapping but distinct expression profiles, with high NANOGP1 expression restricted to early epiblast cells and naïve-state pluripotent stem cells. Sequence analysis and epitope-tagging revealed that NANOGP1 is protein coding with an intact homeobox domain. The duplication that created NANOGP1 occurred earlier in primate evolution than previously thought and has been retained only in great apes, whereas Old World monkeys have disabled the gene in different ways, including homeodomain point mutations. NANOGP1 is a strong inducer of naïve pluripotency; however, unlike NANOG, it is not required to maintain the undifferentiated status of human naïve pluripotent cells. By retaining expression, sequence and partial functional conservation with its ancestral copy, NANOGP1 exemplifies how gene duplication and subfunctionalisation can contribute to transcription factor activity in human pluripotency and development.

The road less travelled? Exploring the nuanced evolutionary consequences of duplicated genes

Duplicated genes have long been appreciated as both substrates and catalysts of evolutionary processes. From even the simplest cell to complex multicellular animals and plants, duplicated genes have made immeasurable contributions to the phenotypic evolution of all life on Earth. Not merely drivers of morphological innovation and speciation events, however, gene duplications sculpt the evolution of genetic architecture in ways we are only just coming to understand now we have the experimental tools to do so. As such, the present article revisits our understanding of the ways in which duplicated genes evolve, examining closely the various fates they can adopt in light of recent work that yields insights from studies of paralogues from across the tree of life that challenge the classical framework.

Pan-cancer surveys indicate cell cycle-related roles of primate-specific genes in tumors and embryonic cerebrum

BACKGROUND

Despite having been extensively studied, it remains largely unclear why humans bear a particularly high risk of cancer. The antagonistic pleiotropy hypothesis predicts that primate-specific genes (PSGs) tend to promote tumorigenesis, while the molecular atavism hypothesis predicts that PSGs involved in tumors may represent recently derived duplicates of unicellular genes. However, these predictions have not been tested.

RESULTS

By taking advantage of pan-cancer genomic data, we find the upregulation of PSGs across 13 cancer types, which is facilitated by copy-number gain and promoter hypomethylation. Meta-analyses indicate that upregulated PSGs (uPSGs) tend to promote tumorigenesis and to play cell cycle-related roles. The cell cycle-related uPSGs predominantly represent derived duplicates of unicellular genes. We prioritize 15 uPSGs and perform an in-depth analysis of one unicellular gene-derived duplicate involved in the cell cycle, DDX11. Genome-wide screening data and knockdown experiments demonstrate that DDX11 is broadly essential across cancer cell lines. Importantly, non-neutral amino acid substitution patterns and increased expression indicate that DDX11 has been under positive selection. Finally, we find that cell cycle-related uPSGs are also preferentially upregulated in the highly proliferative embryonic cerebrum.

CONCLUSIONS

Consistent with the predictions of the atavism and antagonistic pleiotropy hypotheses, primate-specific genes, especially those PSGs derived from cell cycle-related genes that emerged in unicellular ancestors, contribute to the early proliferation of the human cerebrum at the cost of hitchhiking by similarly highly proliferative cancer cells.

Characterization and expression profiling of G protein-coupled receptors (GPCRs) in Spodoptera litura (Lepidoptera: Noctuidae)

Spodoptera litura is a highly destructive omnivorous pest, and they caused serious damage to various crops. G protein-coupled receptors (GPCRs) mediate dozens of physiological processes including reproduction, development, life span and behaviors, but the information of these receptors has been lacking in S. litura. Here, we methodically identified 122 GPCRs in S. litura and made an assay of their expression patterns in different tissues. Comparing the identified GPCRs with homologous genes of other insects, it is obvious that the subfamily A2 (biogenic amine receptors) and the subfamily A3 (neuropeptide and protein hormone receptors) of S. litura have expanded to a certain extent, which may be related to the omnivorous nature and drought environment resistance of S. litura. Besides, the large Methuselah (Mth)/Methuselah-like (Mthl) subfamily of S. litura may be involved in many physiological functions such as longevity and stress response. Apart from duplicate receptors, the loss of parathyroid hormone receptor (PTHR) and the bride of sevenless (Boss) receptor in the lepidopteran insects may imply a new pattern of wing formation and energy metabolism in lepidopteran insects. In addition, the high expression level of GPCRs in different tissues reflects the functional diversity of GPCRs regulating. Systemic identification and initial characterization of GPCRs in S. litura provide a basis for further studies to reveal the functions of these receptors in regulating physiology and behavior.