The probability of preservation of a newly arisen gene duplicate

A Rigorous Framework to Classify the Postduplication Fate of Paralogous Genes

Gene duplication has a central role in evolution; still, little is known on the fates of the duplicated copies, their relative frequency, and on how environmental conditions affect them. Moreover, the lack of rigorous definitions concerning the fate of duplicated genes hinders the development of a global vision of this process. In this paper, we present a new framework aiming at characterizing and formally differentiating the fate of duplicated genes. Our framework has been tested via simulations, where the evolution of populations has been simulated using aevol, an in silico experimental evolution platform. Our results show several patterns that confirm some of the conclusions from previous studies, while also exhibiting new tendencies; this may open up new avenues to better understand the role of duplications as a driver of evolution.

Evolutionarily conserved 12-oxophytodienoate reductase trans-lncRNA pair affects disease resistance in tea (Camellia sinensis) via the jasmonic acid signaling pathway

Long non-coding RNAs (lncRNAs) have gathered significant attention due to their pivotal role in plant growth, development, and biotic and abiotic stress resistance. Despite this, there is still little understanding regarding the functions of lncRNA in these domains in the tea plant (Camellia sinensis), mainly attributable to the insufficiencies in gene manipulation techniques for tea plants. In this study, we designed a novel strategy to identify evolutionarily conserved trans-lncRNA (ECT-lncRNA) pairs in plants. We used highly consistent base sequences in the exon-overlapping region between trans-lncRNAs and their target gene transcripts. Based on this method, we successfully screened 24 ECT-lncRNA pairs from at least two or more plant species. In tea, as observed in model plants such as Arabidopsis, alfalfa, potatoes, and rice, there exists a trans-lncRNA capable of forming an ECT-lncRNA pair with transcripts of the 12-oxophytodienoate reductase (OPR) family, denoted as the OPRL/OPR pair. Considering evolutionary perspectives, the OPRL gene cluster in each species likely originates from a replication event of the OPR gene cluster. Gene manipulation and gene expression analysis revealed that CsOPRL influences disease resistance by regulating CsOPR expression in tea plants. Furthermore, the knockout of StOPRL1 in Solanum tuberosum led to aberrant growth characteristics and strong resistance to fungal infection. This study provides insights into a strategy for the screening and functional verification of ECT-lncRNA pairs.

Analysing a Group of Homologous BAHD Enzymes Provides Insights into the Evolutionary Transition of Rosmarinic Acid Synthases from Hydroxycinnamoyl-CoA:Shikimate/Quinate Hydroxycinnamoyl Transferases

The interplay of various enzymes and compounds gives rise to the intricate secondary metabolic networks observed today. However, the current understanding of their formation and expansion remains limited. BAHD acyltransferases play important roles in the biosynthesis of numerous significant secondary metabolites. In plants, they are widely distributed and exhibit a diverse range of activities. Among them, rosmarinic acid synthase (RAS) and hydroxycinnamoyl-CoA:shikimate/quinate hydroxycinnamoyl transferase (HCT) have gained significant recognition and have been extensively investigated as prominent members of the BAHD acyltransferase family. Here, we conducted a comprehensive study on a unique group of RAS homologous enzymes in Mentha longifolia that display both catalytic activities and molecular features similar to HCT and Lamiaceae RAS. Subsequent phylogenetic and comparative genome analyses revealed their derivation from expansion events within the HCT gene family, indicating their potential as collateral branches along the evolutionary trajectory, leading to Lamiaceae RAS while still retaining certain ancestral vestiges. This discovery provides more detailed insights into the evolution from HCT to RAS. Our collective findings indicate that gene duplication is the driving force behind the observed evolutionary pattern in plant-specialized enzymes, which probably originated from ancestral enzyme promiscuity and were subsequently shaped by principles of biological adaptation.

Models for the retention of duplicate genes and their biological underpinnings

Gene content in genomes changes through several different processes, with gene duplication being an important contributor to such changes. Gene duplication occurs over a range of scales from individual genes to whole genomes, and the dynamics of this process can be context dependent. Still, there are rules by which genes are retained or lost from genomes after duplication, and probabilistic modeling has enabled characterization of these rules, including their context-dependence. Here, we describe the biology and corresponding mathematical models that are used to understand duplicate gene retention and its contribution to the set of biochemical functions encoded in a genome.

The genome of the ant Tetramorium bicarinatum reveals a tandem organization of venom peptides genes allowing the prediction of their regulatory and evolutionary profiles

BACKGROUND

Venoms have evolved independently over a hundred times in the animal kingdom to deter predators and/or subdue prey. Venoms are cocktails of various secreted toxins, whose origin and diversification provide an appealing system for evolutionary researchers. Previous studies of the ant venom of Tetramorium bicarinatum revealed several Myrmicitoxin (MYRTX) peptides that gathered into seven precursor families suggesting different evolutionary origins. Analysis of the T. bicarinatum genome enabling further genomic approaches was necessary to understand the processes underlying the evolution of these myrmicitoxins.

RESULTS

Here, we sequenced the genome of Tetramorium bicarinatum and reported the organisation of 44 venom peptide genes (vpg). Of the eleven chromosomes that make up the genome of T. bicarinatum, four carry the vpg which are organized in tandem repeats. This organisation together with the ML evolutionary analysis of vpg sequences, is consistent with evolution by local duplication of ancestral genes for each precursor family. The structure of the vpg into two or three exons is conserved after duplication events while the promoter regions are the least conserved parts of the vpg even for genes with highly identical sequences. This suggests that enhancer sequences were not involved in duplication events, but were recruited from surrounding regions. Expression level analysis revealed that most vpg are highly expressed in venom glands, although one gene or group of genes is much more highly expressed in each family. Finally, the examination of the genomic data revealed that several genes encoding transcription factors (TFs) are highly expressed in the venom glands. The search for binding sites (BS) of these TFs in the vpg promoters revealed hot spots of GATA sites in several vpg families.

CONCLUSION

In this pioneering investigation on ant venom genes, we provide a high-quality assembly genome and the annotation of venom peptide genes that we think can fosters further genomic research to understand the evolutionary history of ant venom biochemistry.

Consequences of polyploidy and divergence as revealed by cytogenetic mapping of tandem repeats in African clawed frogs (Xenopus, Pipidae)

Repetitive elements have been identified in several amphibian genomes using whole genome sequencing, but few studies have used cytogenetic mapping to visualize these elements in this vertebrate group. Here, we used fluorescence in situ hybridization and genomic data to map the U1 and U2 small nuclear RNAs and histone H3 in six species of African clawed frog (genus Xenopus), including, from subgenus Silurana, the diploid Xenopus tropicalis and its close allotetraploid relative X. calcaratus and, from subgenus Xenopus, the allotetraploid species X. pygmaeus, X. allofraseri, X. laevis, and X. muelleri. Results allowed us to qualitatively evaluate the relative roles of polyploidization and divergence in the evolution of repetitive elements because our focal species include allotetraploid species derived from two independent polyploidization events - one that is relatively young that gave rise to X. calcaratus and another that is older that gave rise to the other (older) allotetraploids. Our results demonstrated conserved loci number and position of signals in the species from subgenus Silurana; allotetraploid X. calcaratus has twice as many signals as diploid X. tropicalis. However, the content of repeats varied among the other allotetraploid species. We detected almost same number of signals in X. muelleri as in X. calcaratus and same number of signals in X. pygmaeus, X. allofraseri, X. laevis as in the diploid X. tropicalis. Overall, these results are consistent with the proposal that allopolyploidization duplicated these tandem repeats and that variation in their copy number was accumulated over time through reduction and expansion in a subset of the older allopolyploids.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.

Population genomic evidence of selection on structural variants in a natural hybrid zone

Structural variants (SVs) can promote speciation by directly causing reproductive isolation or by suppressing recombination across large genomic regions. Whereas examples of each mechanism have been documented, systematic tests of the role of SVs in speciation are lacking. Here, we take advantage of long-read (Oxford nanopore) whole-genome sequencing and a hybrid zone between two Lycaeides butterfly taxa (L. melissa and Jackson Hole Lycaeides) to comprehensively evaluate genome-wide patterns of introgression for SVs and relate these patterns to hypotheses about speciation. We found >100,000 SVs segregating within or between the two hybridizing species. SVs and SNPs exhibited similar levels of genetic differentiation between species, with the exception of inversions, which were more differentiated. We detected credible variation in patterns of introgression among SV loci in the hybrid zone, with 562 of 1419 ancestry-informative SVs exhibiting genomic clines that deviated from null expectations based on genome-average ancestry. Overall, hybrids exhibited a directional shift towards Jackson Hole Lycaeides ancestry at SV loci, consistent with the hypothesis that these loci experienced more selection on average than SNP loci. Surprisingly, we found that deletions, rather than inversions, showed the highest skew towards excess ancestry from Jackson Hole Lycaeides. Excess Jackson Hole Lycaeides ancestry in hybrids was also especially pronounced for Z-linked SVs and inversions containing many genes. In conclusion, our results show that SVs are ubiquitous and suggest that SVs in general, but especially deletions, might disproportionately affect hybrid fitness and thus contribute to reproductive isolation.

Divergent expression of paralogous genes by modification of shared enhancer activity through a promoter-proximal silencer

The duplication of genes and their associated cis-regulatory elements, or enhancers, is a key contributor to genome evolution and biological complexity. Moreover, many paralogs, particularly tandem duplicates, are fixed for long periods of time under the control of shared enhancers. However, in most cases, the mechanism by which gene expression and function diverge following duplication is not known. Here, we dissect the regulation and function of the paralogous nubbin/pdm2 genes during wing development in Drosophila melanogaster. We show that these paralogs play a redundant role in the wing and that their expression relies on a single shared wing enhancer. However, the two genes differ in their ability to respond to this enhancer, with nub responding in all wing progenitor cells and pdm2 only in a small subset. This divergence is a result of a pdm2-specific silencer element at the pdm2 promoter that receives repressive input from the transcription factor Rotund. Repression through this silencer also depends on nub, allowing pdm2 to fully respond to the wing enhancer when nub expression is perturbed and functional compensation to occur. Thus, expression divergence downstream of a shared enhancer arises as a consequence of silencing the promoter of one paralog.

The logic of virus evolution

Viruses are obligate intracellular parasites. Despite their dependence on host cells, viruses are evolutionarily autonomous, with their own genomes and evolutionary trajectories locked in arms races with the hosts. Here, we discuss a simple functional logic to explain virus macroevolution that appears to define the course of virus evolution. A small core of virus hallmark genes that are responsible for genome replication apparently descended from primordial replicators, whereas most virus genes, starting with those encoding capsid proteins, were subsequently acquired from hosts. The oldest of these acquisitions antedate the last universal cellular ancestor (LUCA). Host gene capture followed two major routes: convergent recruitment of genes with functions that directly benefit virus reproduction and exaptation when host proteins are repurposed for unique virus functions. These forms of host protein recruitment by viruses result in different levels of similarity between virus and host homologs, with the exapted ones often changing beyond easy recognition.

A subfunctionalisation model of gene family evolution predicts balanced tree shapes

We consider a subfunctionalisation model of gene family evolution. A family of n genes that perform z functions is represented by an n×z binary matrix Y_t where a 1 in the ijth position indicates that gene i can perform function j. Y_t evolves according to a continuous time Markov chain (CTMC) that represents the processes of gene duplication, coding region loss and regulatory region loss with the restriction that each function is protected by selection, meaning that each column in the matrix must contain at least one 1. We generate gene trees based on the CTMC {Y_t,t⩾0}. We analyse the long-run behaviour of the model and specify the conditions where we expect gene trees to continue to grow and where we expect them to have a stable number of genes. We show that different choices of rate parameters for the processes of duplication and loss lead to different tree shapes as measured by two common tree-shape statistics: the β-statistic for measuring tree balance and the γ-statistic for assessing diversification rate. We use an extension of β that is estimated from sets of trees. This extension is less biased compared to using the average β value for individual trees. When the rate of gene duplication dominates the rate of gene loss, the process is unstable and the distribution of tree shapes is close to following the uniform ranked tree shape (URT) distribution. However, when the process is stable, gene trees are predicted to have positive values of β indicating balanced trees and negative values of γ indicating that diversification occurs more towards the root of the tree. The results of our analyses suggest that comparing the tree-shape statistics of empirical gene-trees to the predictions presented here will provide a test of the subfunctionalisation model.

Evolutionary Genomics Reveals Multiple Functions of Arylalkylamine N-Acetyltransferase in Fish

As an important hormone, melatonin participates in endocrine regulation of diverse functions in vertebrates. Its biosynthesis is catalyzed by four cascaded enzymes, among them, arylalkylamine N-acetyltransferase (AANAT) is the most critical one. Although only single aanat gene has been identified in most groups of vertebrates, researchers including us have determined that fish have the most diverse of aanat genes (aanat1a, aanat1b, and aanat2), playing various potential roles such as seasonal migration, amphibious aerial vision, and cave or deep-sea adaptation. With the rapid development of genome and transcriptome sequencing, more and more putative sequences of fish aanat genes are going to be available. Related phylogeny and functional investigations will enrich our understanding of AANAT functions in various fish species.

The Role of Interspecific Hybridisation in Adaptation and Speciation: Insights From Studies in Senecio

Hybridisation is well documented in many species, especially plants. Although hybrid populations might be short-lived and do not evolve into new lineages, hybridisaiton could lead to evolutionary novelty, promoting adaptation and speciation. The genus Senecio (Asteraceae) has been actively used to unravel the role of hybridisation in adaptation and speciation. In this article, we first briefly describe the process of hybridisation and the state of hybridisation research over the years. We then discuss various roles of hybridisation in plant adaptation and speciation illustrated with examples from different Senecio species, but also mention other groups of organisms whenever necessary. In particular, we focus on the genomic and transcriptomic consequences of hybridisation, as well as the ecological and physiological aspects from the hybrids' point of view. Overall, this article aims to showcase the roles of hybridisation in speciation and adaptation, and the research potential of Senecio, which is part of the ecologically and economically important family, Asteraceae.

Genome-wide identification of soybean Shaker K+ channel gene family and functional characterization of GmAKT1 in transgenic Arabidopsis thaliana under salt and drought stress

Potassium is a major cationic nutrient involved in numerous physiological processes in plants. The uptake of K+ is mediated by K+ channels and transporters, and the Shaker K+ channel gene family plays an essential role in K+ uptake and stress resistance in plants. However, little is known regarding this family in soybean. In this study, 14 members of the Shaker K+ channel gene family were identified in soybean and were classified into five groups. Protein domain analysis revealed that Shaker K+ channel gene members have an ion transport domain (ion trans), a cyclic nucleotide-binding domain, ankyrin repeat domains, and a dimerization domain in the potassium ion channel. Quantitative real-time polymerase chain reaction analysis indicated that the expression of eight genes (notably GmAKT1) in soybean leaves and roots was significantly increased in response to salt and drought stress. Furthermore, the overexpression of GmAKT1 in Arabidopsis enhanced root length, K+ concentration, and fresh/dry weight ratio compared with wild-type plants subjected to salt and drought stress; this suggests that GmAKT1 improves the tolerance of soybean to abiotic stress. Our results provide important insight into the characterization of Shaker K+ channel gene family members in soybean and highlight the function of GmAKT1 in soybean plants under salt and drought stress.

Detecting Selection on Segregating Gene Duplicates in a Population

Gene duplication is a fundamental process that has the potential to drive phenotypic differences between populations and species. While evolutionarily neutral changes have the potential to affect phenotypes, detecting selection acting on gene duplicates can uncover cases of adaptive diversification. Existing methods to detect selection on duplicates work mostly inter-specifically and are based upon selection on coding sequence changes, here we present a method to detect selection directly on a copy number variant segregating in a population. The method relies upon expected relationships between allele (new duplication) age and frequency in the population dependent upon the effective population size. Using both a haploid and a diploid population with a Moran Model under several population sizes, the neutral baseline for copy number variants is established. The ability of the method to reject neutrality for duplicates with known age (measured in pairwise dS value) and frequency in the population is established through mathematical analysis and through simulations. Power is particularly good in the diploid case and with larger effective population sizes, as expected. With extension of this method to larger population sizes, this is a tool to analyze selection on copy number variants in any natural or experimentally evolving population. We have made an R package available at https://github.com/peterbchi/CNVSelectR/ which implements the method introduced here.

How Important Are Structural Variants for Speciation?

Understanding the genetic basis of reproductive isolation is a central issue in the study of speciation. Structural variants (SVs); that is, structural changes in DNA, including inversions, translocations, insertions, deletions, and duplications, are common in a broad range of organisms and have been hypothesized to play a central role in speciation. Recent advances in molecular and statistical methods have identified structural variants, especially inversions, underlying ecologically important traits; thus, suggesting these mutations contribute to adaptation. However, the contribution of structural variants to reproductive isolation between species—and the underlying mechanism by which structural variants most often contribute to speciation—remain unclear. Here, we review (i) different mechanisms by which structural variants can generate or maintain reproductive isolation; (ii) patterns expected with these different mechanisms; and (iii) relevant empirical examples of each. We also summarize the available sequencing and bioinformatic methods to detect structural variants. Lastly, we suggest empirical approaches and new research directions to help obtain a more complete assessment of the role of structural variants in speciation.

Two polygalacturonase-inhibiting proteins (VrPGIP) of Vigna radiata confer resistance to bruchids (Callosobruchus spp.)

Bruchids (Callosobruchus spp.) are destructive storage pests of mung beans (Vigna radiata). Bruchids infest mature seeds during storage and in the field causing heavy losses. Bruchid resistance in mung bean has been characterized as a dominant trait controlled by a single gene. Several independent mapping studies showed that the Br locus on chromosome 5 was a key quantitative trait loci (QTL) involved in bruchid resistance. Two polygalacturonase-inhibitor protein (PGIP) family genes, VrPGIP1 and VrPGIP2, located in the Br locus may be the primary genes responsible for bruchid resistance in mung bean but no experimental proof is available. We isolated the VrPGIP1 and VrPGIP2 genes from bruchid resistant mung bean cultivar V2802 and purified the proteins by prokaryotic expression. Both VrPGIP1 and VrPGIP2 had polygalacturonase inhibitor activity and both of the PGIP proteins conferred resistance to bruchids in an artificial seed test system. VrPGIPs can inhibit the enzyme activity of polygalacturonase present in males, females and fourth instar larvae of C. maculatus. These results demonstrated that VrPGIP1 and VrPGIP2 play a critical role in bruchid resistance probably through inhibiting polygalacturonase activity.

The new chimeric chiron genes evolved essential roles in zebrafish embryonic development by regulating NAD+ levels

The origination of new genes is important for generating genetic novelties for adaptive evolution and biological diversity. However, their potential roles in embryonic development, evolutionary processes into ancient networks, and contributions to adaptive evolution remain poorly investigated. Here, we identified a novel chimeric gene family, the chiron family, and explored its genetic basis and functional evolution underlying the adaptive evolution of Danioninae fishes. The ancestral chiron gene originated through retroposition of nampt in Danioninae 48-54 million years ago (Mya) and expanded into five duplicates (chiron1-5) in zebrafish 1-4 Mya. The chiron genes (chirons) likely originated in embryonic development and gradually extended their expression in the testis. Functional experiments showed that chirons were essential for zebrafish embryo development. By integrating into the NAD⁺ synthesis pathway, chirons could directly catalyze the NAD⁺ rate-limiting reaction and probably impact two energy metabolism genes (nmnat1 and naprt) to be under positive selection in Danioninae fishes. Together, these results mainly demonstrated that the origin of new chimeric chiron genes may be involved in adaptive evolution by integrating and impacting the NAD⁺ biosynthetic pathway. This coevolution may contribute to the physiological adaptation of Danioninae fishes to widespread and varied biomes in Southeast Asian.

Two novel venom proteins underlie divergent parasitic strategies between a generalist and a specialist parasite

Parasitoids are ubiquitous in natural ecosystems. Parasitic strategies are highly diverse among parasitoid species, yet their underlying genetic bases are poorly understood. Here, we focus on the divergent adaptation of a specialist and a generalist drosophilid parasitoids. We find that a novel protein (Lar) enables active immune suppression by lysing the host lymph glands, eventually leading to successful parasitism by the generalist. Meanwhile, another novel protein (Warm) contributes to a passive strategy by attaching the laid eggs to the gut and other organs of the host, leading to incomplete encapsulation and helping the specialist escape the host immune response. We find that these diverse parasitic strategies both originated from lateral gene transfer, followed with duplication and specialization, and that they might contribute to the shift in host ranges between parasitoids. Our results increase our understanding of how novel gene functions originate and how they contribute to host adaptation.

Expression Partitioning of Duplicate Genes at Single Cell Resolution in Arabidopsis Roots

Gene duplication is a key evolutionary phenomenon, prevalent in all organisms but particularly so in plants, where whole genome duplication (WGD; polyploidy) is a major force in genome evolution. Much effort has been expended in attempting to understand the evolution of duplicate genes, addressing such questions as why some paralog pairs rapidly return to single copy status whereas, in other pairs, both paralogs are retained and may diverge in expression pattern or function. The effect of a gene - its site of expression and thus the initial locus of its function - occurs at the level of a cell comprising a single cell type at a given state of the cell's development. Using Arabidopsis thaliana single cell transcriptomic data we categorized patterns of expression for 11,470 duplicate gene pairs across 36 cell clusters comprising nine cell types and their developmental states. Among these 11,470 pairs, 10,187 (88.8%) had at least one copy expressed in at least one of the 36 cell clusters. Pairs produced by WGD more often had both paralogs expressed in root cells than did pairs produced by small scale duplications. Three quarters of gene pairs expressed in the 36 cell clusters (7,608/10,187) showed extreme expression bias in at least one cluster, including 352 cases of reciprocal bias, a pattern consistent with expression subfunctionalization. More than twice as many pairs showed reciprocal expression bias between cell states than between cell types or between roots and leaves. A group of 33 gene pairs with reciprocal expression bias showed evidence of concerted divergence of gene networks in stele vs. epidermis. Pairs with both paralogs expressed without bias were less likely to have paralogs with divergent mutant phenotypes; such bias-free pairs showed evidence of preservation by maintenance of dosage balance. Overall, we found considerable evidence of shifts in gene expression following duplication, including in >80% of pairs encoding 7,653 genes expressed ubiquitously in all root cell types and states for which we inferred the polarity of change.

When a foreign gene meets its native counterpart: computational biophysics analysis of two PgiC loci in the grass Festuca ovina

Duplicative horizontal gene transfer may bring two previously separated homologous genes together, which may raise questions about the interplay between the gene products. One such gene pair is the “native” PgiC1 and “foreign” PgiC2 in the perennial grass Festuca ovina. Both PgiC1 and PgiC2 encode cytosolic phosphoglucose isomerase, a dimeric enzyme whose proper binding is functionally essential. Here, we use biophysical simulations to explore the inter-monomer binding of the two homodimers and the heterodimer that can be produced by PgiC1 and PgiC2 in F. ovina. Using simulated native-state ensembles, we examine the structural properties and binding tightness of the dimers. In addition, we investigate their ability to withstand dissociation when pulled by a force. Our results suggest that the inter-monomer binding is tighter in the PgiC2 than the PgiC1 homodimer, which could explain the more frequent occurrence of the foreign PgiC2 homodimer in dry habitats. We further find that the PgiC1 and PgiC2 monomers are compatible with heterodimer formation; the computed binding tightness is comparable to that of the PgiC1 homodimer. Enhanced homodimer stability and capability of heterodimer formation with PgiC1 are properties of PgiC2 that may contribute to the retaining of the otherwise redundant PgiC2 gene.

A Screen for Gene Paralogies Delineating Evolutionary Branching Order of Early Metazoa

The evolutionary diversification of animals is one of Earth’s greatest marvels, yet its earliest steps are shrouded in mystery. Animals, the monophyletic clade known as Metazoa, evolved wildly divergent multicellular life strategies featuring ciliated sensory epithelia. In many lineages epithelial sensoria became coupled to increasingly complex nervous systems. Currently, different phylogenetic analyses of single-copy genes support mutually-exclusive possibilities that either Porifera or Ctenophora is sister to all other animals. Resolving this dilemma would advance the ecological and evolutionary understanding of the first animals and the evolution of nervous systems. Here we describe a comparative phylogenetic approach based on gene duplications. We computationally identify and analyze gene families with early metazoan duplications using an approach that mitigates apparent gene loss resulting from the miscalling of paralogs. In the transmembrane channel-like (TMC) family of mechano-transducing channels, we find ancient duplications that define separate clades for Eumetazoa (Placozoa + Cnidaria + Bilateria) vs. Ctenophora, and one duplication that is shared only by Eumetazoa and Porifera. In the Max-like protein X (MLX and MLXIP) family of bHLH-ZIP regulators of metabolism, we find that all major lineages from Eumetazoa and Porifera (sponges) share a duplicated gene pair that is sister to the single-copy gene maintained in Ctenophora. These results suggest a new avenue for deducing deep phylogeny by choosing rather than avoiding ancient gene paralogies.

A Developmental Switch Generating Phenotypic Plasticity Is Part of a Conserved Multi-gene Locus

Switching between alternative complex phenotypes is often regulated by "supergenes," polymorphic clusters of linked genes such as in butterfly mimicry. In contrast, phenotypic plasticity results in alternative complex phenotypes controlled by environmental influences rather than polymorphisms. Here, we show that the developmental switch gene regulating predatory versus non-predatory mouth-form plasticity in the nematode Pristionchus pacificus is part of a multi-gene locus containing two sulfatases and two α-N-acetylglucosaminidases (nag). We provide functional characterization of all four genes, using CRISPR-Cas9-based reverse genetics, and show that nag genes and the previously identified eud-1/sulfatase have opposing influences. Members of the multi-gene locus show non-overlapping neuronal expression and epistatic relationships. The locus architecture is conserved in the entire genus Pristionchus. Interestingly, divergence between paralogs is counteracted by gene conversion, as inferred from phylogenies and genotypes of CRISPR-Cas9-induced mutants. Thus, we found that physical linkage accompanies regulatory linkage between switch genes controlling plasticity in P. pacificus.

Paralog dependency indirectly affects the robustness of human cells

The protective redundancy of paralogous genes partly relies on the fact that they carry their functions independently. However, a significant fraction of paralogous proteins may form functionally dependent pairs, for instance, through heteromerization. As a consequence, one could expect these heteromeric paralogs to be less protective against deleterious mutations. To test this hypothesis, we examined the robustness landscape of gene loss-of-function by CRISPR-Cas9 in more than 450 human cell lines. This landscape shows regions of greater deleteriousness to gene inactivation as a function of key paralog properties. Heteromeric paralogs are more likely to occupy such regions owing to their high expression and large number of protein-protein interaction partners. Further investigation revealed that heteromers may also be under stricter dosage balance, which may also contribute to the higher deleteriousness upon gene inactivation. Finally, we suggest that physical dependency may contribute to the deleteriousness upon loss-of-function as revealed by the correlation between the strength of interactions between paralogs and their higher deleteriousness upon loss of function.

Genes Encoding Teleost Fish Ligands and Associated Receptors Remained in Duplicate More Frequently than the Rest of the Genome

Signaling through ligand/receptor interactions is a widespread mechanism across all living taxa. During evolution, however, there has been a diversification in multigene families and changes in their interaction patterns. Among the events that led to the creation of new genes is the whole-genome duplication, which made possible some major innovations. Teleost fishes descended from a common ancestor which underwent one such whole-genome duplication.

In our study, we investigated the effect of complete genome duplication on the evolution of ligand–receptor pairs in teleosts. We selected ten teleost species and used bioinformatics programs and phylogenetic tools in order to study the evolution of the human ligands and receptors that have orthologous genes in fishes, as well as the rest of the fish genomes.

We established that since the complete duplication of the fish genomes, the conservation in duplicate copy of ligand and receptor genes is higher than expected. However, the ligand/receptor pair partners did not necessarily evolve in the same way, and a lot of situations occurred in which one of the partners returned in singleton copy when the other one was maintained in duplicate. This suggests that changes in interaction partners may have taken place during the evolution of teleosts. Moreover, the fate of the ligands and receptor coding genes is partly congruent with the phylogeny of teleosts. However, some incongruences can be observed. We suggest that these incongruences are correlated to the environment.

Genome-wide identification and expression analysis reveal the potential function of ethylene responsive factor gene family in response to Botrytis cinerea infection and ovule development in grapes (Vitis vinifera L.)

The prevention of Botrytis cinerea infection and the study of grape seedlessness are very important for grape industries. Finding correlated regulatory genes is an important approach towards understanding their molecular mechanisms. Ethylene responsive factor (ERF) gene family play critical roles in defence networks and the growth of plants. To date, no large-scale study of the ERF proteins associated with pathogen defence and ovule development has been performed in grape (Vitis vinifera L.). In the present study, we identified 113 ERF genes (VvERF) and named them based on their chromosome locations. The ERF genes could be divided into 11 groups based on a multiple sequence alignment and a phylogenetic comparison with homologues from Arabidopsis thaliana. Synteny analysis and Ka/Ks ratio calculation suggested that segmental and tandem duplications contributed to the expansion of the ERF gene family. The evolutionary relationships between the VvERF genes were investigated by exon-intron structure characterisation, and an analysis of the cis-acting regulatory elements in their promoters suggested potential regulation after stress or hormone treatments. Expression profiling after infection with the fungus, B. cinerea, indicated that ERF genes function in responses to pathogen attack. In addition, the expression levels of most ERF genes were much higher during ovule development in seedless grapes, suggesting a role in ovule abortion related to seedlessness. Taken together, these results indicate that VvERF proteins are involved in responses to Botrytis cinerea infection and in grape ovule development. This information may help guide strategies to improve grape production.

In silico genome-wide identification and comprehensive characterization of the BES1 gene family in soybean

The BES1 transcription factor family play a central role in brassinosteroid signaling pathway that regulates a wide range of plant growth and developmental processes, as well as resistances to various stresses. However, no comprehensive study of the BES1 gene family in soybean has been reported. In this work, 16 GmBES1-like genes were identified in soybean, which could be divided into two clades based on their phylogenetic relationships, gene structures and motif compositions. We then examined their duplication status and evolutionary models. The result showed that most of the GmBES1-like genes have duplicated counterparts generated from the recent Glycine WGD event, and these genes are originated from 6 distinct ancestors before the Gamma WGT event. We further studied the expression profiles of GmBES1-like genes, and found their spatio-temporal and stressed expression patterns varied tremendously. For example, GmBES1-5 and GmBES1-6 were highly expressed in almost every sample, whereas GmBES1-7 and GmBES1-8 were not expressed. Additionally, interaction network analysis revealed the presence of 3 clusters between GmBES1-like genes and other associated genes, implying that they have both the conserved and divergent functions. Lastly, we analyzed the genetic diversity of GmBES1-like genes in 302 resequenced wild, landrace and improved soybean accessions. It showed that most of these genes are well conserved, and they are not changed during domestication and improvement. These results provide insights into the characterization of GmBES1 family and lay the foundation for further functional study of such genes.

The EIL transcription factor family in soybean: Genome-wide identification, expression profiling and genetic diversity analysis

The ETHYLENE INSENSITIVE3-LIKE (EIL) transcription factor family plays a critical role in the ethylene signaling pathway, which regulates a broad spectrum of plant growth and developmental processes, as well as defenses to myriad stresses. Although genome-wide analysis of this family has been carried out for several plant species, no comprehensive analysis of the EIL gene family in soybean has been reported so far. Furthermore, there are few studies on the functions of EIL genes in soybean. In this study, we identified 12 soybean (Gm) EIL genes, which we divided into three groups based on their phylogenetic relationships. We then detected their duplication status and found that most of the GmEIL genes have duplicated copies derived from two whole-genome duplication events. These duplicated genes underwent strong negative selection during evolution. We further analyzed the transcript profiles of GmEIL genes using the transcriptome data and found that their spatio-temporal and stress expression patterns varied considerably. For example, GmEIL1-GmEIL5 were found to be strongly expressed in almost every sample, while GmEIL8-GmEIL12 exhibited low expression, or were not expressed at all. Additionally, these genes showed different responses to dehydration, salinity and phosphate starvation. Finally, we surveyed genetic variations of these genes in 302 resequenced wild soybeans, landraces and improved soybean cultivars. Our data showed that most GmEIL genes are well conserved, and are not modified in domesticated or improved cultivars. Together, these findings provide a potentially valuable resource for characterizing the GmEIL gene family and lay the basis for further elucidation of their molecular mechanisms.

Comprehensive Genomic Analysis and Expression Profiling of Diacylglycerol Kinase (DGK) Gene Family in Soybean (Glycine max) under Abiotic Stresses

Diacylglycerol kinase (DGK) is an enzyme that plays a pivotal role in abiotic and biotic stress responses in plants by transforming the diacylglycerol into phosphatidic acid. However, there is no report on the characterization of soybean DGK genes in spite of the availability of the soybean genome sequence. In this study, we performed genome-wide analysis and expression profiling of the DGK gene family in the soybean genome. We identified 12 DGK genes (namely GmDGK1-12) which all contained conserved catalytic domains with protein lengths and molecular weights ranging from 436 to 727 amino acids (aa) and 48.62 to 80.93 kDa, respectively. Phylogenetic analyses grouped GmDGK genes into three clusters-cluster I, cluster II, and cluster III-which had three, four, and five genes, respectively. The qRT-PCR analysis revealed significant GmDGK gene expression levels in both leaves and roots coping with polyethylene glycol (PEG), salt, alkali, and salt/alkali treatments. This work provides the first characterization of the DGK gene family in soybean and suggests their importance in soybean response to abiotic stress. These results can serve as a guide for future studies on the understanding and functional characterization of this gene family.

Magic roundabout is an endothelial-specific ohnolog of ROBO1 which neo-functionalized to an essential new role in angiogenesis

BACKGROUND

Magic roundabout (ROBO4) is an unusual endothelial-specific paralog of the family of neuronally-expressed axon guidance receptors called roundabouts. Endothelial cells (ECs), whose uninterrupted sheet delimits the lumen of all vertebrate blood vessels and which are absent from invertebrate species, are a vertebrate-specific evolutionary novelty.

RESULTS

Herein, the evolutionary mechanism of the duplication, retention and divergence of ROBO4 was investigated for the first time. Phylogenetic analyses carried out suggested that ROBO4 is a fast-evolving paralog of ROBO1 formed at the base of vertebrates. The ancestral expression pattern was neuronal. ROBO4 dramatically shifted its expression and became exceptionally specific to ECs. The data-mining of FANTOM5 and ENCODE reveals that ROBO4's endothelial expression arises from a single transcription start site (TSS), conserved in mouse, controlled by a proximal promoter with a complex architecture suggestive of regulatory neo-functionalization. (An analysis of promoter probabilities suggested the architecture was not due to a chance arrangement of TFBSes). Further evidence for the neo-functionalization of ROBO4 comes from the analysis of its protein interactions, the rates of protein evolution, and of positively selected sites.

CONCLUSIONS

The neo-functionalization model explains why ROBO4 protein acquired new context-specific biological functions in the control of angiogenesis. This endothelial-specific roundabout receptor is an illustrative example of the emergence of an essential vertebrate molecular novelty and an endothelial-specific signaling sub-network through 2R-WGD. The emergence of novel cell types, such as ECs, might be a neglected evolutionary force contributing to the high rate of retention of duplicates post-2R-WGD. Crucially, expression neo-functionalization to evolutionarily novel sites of expression conceptually extends the classical model of neo-functionalization.

Duplicated flavonoid 3'-hydroxylase and flavonoid 3', 5'-hydroxylase genes in barley genome

Background

Anthocyanin compounds playing multiple biological functions can be synthesized in different parts of barley (Hordeum vulgare L.) plant. The diversity of anthocyanin molecules is related with branching the pathway to alternative ways in which dihydroflavonols may be modified either with the help of flavonoid 3'-hydroxylase (F3'H) or flavonoid 3',5'-hydroxylase (F3'5'H)-the cytochrome P450-dependent monooxygenases. The F3'H and F3'5'H gene families are among the least studied anthocyanin biosynthesis structural genes in barley. The aim of this study was to identify and characterise duplicated copies of the F3'H and F3'5'H genes in the barley genome.

Results

Four copies of the F3'5'H gene (on chromosomes 4HL, 6HL, 6HS and 7HS) and two copies of the F3'H gene (on chromosomes 1HL and 6HS) were identified in barley genome. These copies have either one or two introns. Amino acid sequences analysis demonstrated the presence of the flavonoid hydroxylase-featured conserved motifs in all copies of the F3'H and F3'5'H genes with the exception of F3'5'H-3 carrying a loss-of-function mutation in a conservative cytochrome P450 domain. It was shown that the divergence between F3'H and F3'5'H genes occurred 129 million years ago (MYA) before the emergence of monocot and dicot plant species. The F3'H copy approximately occurred 80 MYA; the appearance of F3'5'H copies occurred 8, 36 and 91 MYA. qRT-PCR analysis revealed the tissue-specific activity for some copies of the studied genes. The F3'H-1 gene was transcribed in aleurone layer, lemma and pericarp (with an increased level in the coloured pericarp), whereas the F3'H-2 gene was expressed in stems only. The F3'5'H-1 gene was expressed only in the aleurone layer, and in a coloured aleurone its expression was 30-fold higher. The transcriptional activity of F3'5'H-2 was detected in different tissues with significantly higher level in uncoloured genotype in contrast to coloured ones. The F3'5'H-3 gene expressed neither in stems nor in aleurone layer, lemma and pericarp. The F3'5'H-4 gene copy was weakly expressed in all tissues analysed.

Conclusion

F3'H and F3'5'H-coding genes involved in anthocyanin synthesis in H. vulgare were identified and characterised, from which the copies designated F3'H-1, F3'H-2, F3'5'H-1 and F3'5'H-2 demonstrated tissue-specific expression patterns. Information on these modulators of the anthocyanin biosynthesis pathway can be used in future for manipulation with synthesis of diverse anthocyanin compounds in different parts of barley plant. Finding both the copies with tissue-specific expression and a copy undergoing pseudogenization demonstrated rapid evolutionary events tightly related with functional specialization of the duplicated members of the cytochrome P450-dependent monooxygenases gene families.

Different divergence events for three pairs of PEBPs in Gossypium as implied by evolutionary analysis

INTRODUCTION

The phosphatidylethanolamine-binding protein (PEBP) gene family plays a crucial role in seed germination, reproductive transformation, and other important developmental processes in plants, but its distribution in Gossypium genomes or species, evolutionary properties, and the fates of multiple duplicated genes remain unclear.

OBJECTIVES

The primary objectives of this study were to elucidate the distribution and characteristics of PEBP genes in Gossypium, as well as the evolutionary pattern of duplication and deletion, and functional differentiation of PEBPs in plants.

METHODS

Using the PEBP protein sequences in Arabidopsis thaliana as queries, blast alignment was carried out for the identification of PEBP genes in four sequenced cotton species. Using the primers designed according to the PEBP genome sequences, PEBP genes were cloned from 15 representative genomes of Gossypium genus, and the gene structure, CDS sequence, protein sequence and properties were predicted and phylogenetic analysis was performed. Taking PEBP proteins of grape as reference, grouping of orthologous gene, analysis of phylogeny and divergence of PEBPs in nine species were conducted to reconstruct the evolutionary pattern of PEBP genes in plants.

RESULTS

We identified and cloned 160 PEBPs from 15 cotton species, and the phylogenetic analysis showed that the genes could be classified into the following three subfamilies: MFT-like, FT-like and TFL1-like. There were eight single orthologous group (OG) members in each diploid and 16 double OG members in each tetraploid. An analysis of the expression and selective pressure indicated that expression divergence and strong purification selection within the same OG presented in the PEBP gene family.

CONCLUSION

An evolutionary pattern of duplication and deletion of the PEBP family in the evolutionary history of Gossypium was suggested, and three pairs of genes resulted from different whole-genome duplication events.

From "the Worm" to "the Worms" and Back Again: The Evolutionary Developmental Biology of Nematodes

Since the earliest days of research on nematodes, scientists have noted the developmental and morphological variation that exists within and between species. As various cellular and developmental processes were revealed through intense focus on Caenorhabditis elegans, these comparative studies have expanded. Within the genus Caenorhabditis, they include characterization of intraspecific polymorphisms and comparisons of distinct species, all generally amenable to the same laboratory culture methods and supported by robust genomic and experimental tools. The C. elegans paradigm has also motivated studies with more distantly related nematodes and animals. Combined with improved phylogenies, this work has led to important insights about the evolution of nematode development. First, while many aspects of C. elegans development are representative of Caenorhabditis, and of terrestrial nematodes more generally, others vary in ways both obvious and cryptic. Second, the system has revealed several clear examples of developmental flexibility in achieving a particular trait. This includes developmental system drift, in which the developmental control of homologous traits has diverged in different lineages, and cases of convergent evolution. Overall, the wealth of information and experimental techniques developed in C. elegans is being leveraged to make nematodes a powerful system for evolutionary cellular and developmental biology.

Ankyrin domain encoding genes from an ancient horizontal transfer are functionally integrated into Nasonia developmental gene regulatory networks

BACKGROUND

How regulatory networks incorporate additional components and how novel genes are functionally integrated into well-established developmental processes are two important and intertwined questions whose answers have major implications for understanding the evolution of development. We recently discovered a set of lineage-restricted genes with strong and specific expression patterns along the dorsal-ventral (DV) axis of the embryo of the wasp Nasonia that may serve as a powerful system for addressing these questions. We sought to both understand the evolutionary history of these genes and to determine their functions in the Nasonia DV patterning system.

RESULTS

We have found that the novel DV genes are part of a large family of rapidly duplicating and diverging ankyrin domain-encoding genes that originated most likely by horizontal transfer from a prokaryote in a common ancestor of the wasp superfamily Chalcidoidea. We tested the function of those ankyrin-encoding genes expressed along the DV axis and found that they participate in early embryonic DV patterning. We also developed a new wasp model system (Melittobia) and found that some functional integration of ankyrin genes have been preserved for over 90 million years.

CONCLUSIONS

Our results indicate that regulatory networks can incorporate novel genes that then become necessary for stable and repeatable outputs. Even a modest role in developmental networks may be enough to allow novel or duplicate genes to be maintained in the genome and become fully integrated network components.

Extreme Deviations from Expected Evolutionary Rates in Archaeal Protein Families

Origin of new biological functions is a complex phenomenon ranging from single-nucleotide substitutions to the gain of new genes via horizontal gene transfer or duplication. Neofunctionalization and subfunctionalization of proteins is often attributed to the emergence of paralogs that are subject to relaxed purifying selection or positive selection and thus evolve at accelerated rates. Such phenomena potentially could be detected as anomalies in the phylogenies of the respective gene families. We developed a computational pipeline to search for such anomalies in 1,834 orthologous clusters of archaeal genes, focusing on lineage-specific subfamilies that significantly deviate from the expected rate of evolution. Multiple potential cases of neofunctionalization and subfunctionalization were identified, including some ancient, house-keeping gene families, such as ribosomal protein S10, general transcription factor TFIIB and chaperone Hsp20. As expected, many cases of apparent acceleration of evolution are associated with lineage-specific gene duplication. On other occasions, long branches in phylogenetic trees correspond to horizontal gene transfer across long evolutionary distances. Significant deceleration of evolution is less common than acceleration, and the underlying causes are not well understood; functional shifts accompanied by increased constraints could be involved. Many gene families appear to be “highly evolvable,” that is, include both long and short branches. Even in the absence of precise functional predictions, this approach allows one to select targets for experimentation in search of new biology.

Genome-Wide Bioinformatics Analysis of MAPK Gene Family in Kiwifruit (Actinidia Chinensis)

Mitogen activated protein kinase (MAPK) cascades are universal signal transduction modules that play crucial roles in various biotic and abiotic stresses, hormones, cell division, and developmental processes in plants. Mitogen activated protein kinase (MAPK/MPK), being a part of this cascade, performs an important function for further appropriate cellular responses. Although MAPKs have been investigated in several model plants, no systematic analysis has been conducted in kiwifruit (Actinidia chinensis). In the present study, we identified 18 putative MAPKs in the kiwifruit genome. This gene family was analyzed bioinformatically in terms of their chromosome locations, sequence alignment, gene structures, and phylogenetic and conserved motifs. All members possess fully canonical motif structures of MAPK. Phylogenetic analysis indicated that AcMAPKs could be classified into five subfamilies, and these gene motifs in the same group showed high similarity. Gene structure analysis demonstrated that the number of exons in AcMAPK genes ranged from 2 to 29, suggesting large variation among kiwifruit MAPK genes. The expression profiles of these AcMAPK genes were further investigated using quantitative real-time polymerase chain reaction (qRT-PCR), which demonstrated that AcMAPKs were induced or repressed by various biotic and abiotic stresses and hormone treatments, suggesting their potential roles in the biotic and abiotic stress response and various hormone signal transduction pathways in kiwifruit. The results of this study provide valuable insight into the putative physiological and biochemical functions of MAPK genes in kiwifruit.

Evolutionary cell biology traces the rise of the exomer complex in Fungi from an ancient eukaryotic component

Cargo is transported from the trans-Golgi Network to the plasma membrane by adaptor complexes, which are pan-eukaryotic components. However, in yeast, cargo can also be exported by the exomer complex, a heterotetrameric protein complex consisting of two copies of Chs5, and any two members of four paralogous proteins (ChAPs). To understand the larger relevance of exomer, its phylogenetic distribution and function outside of yeast need to be explored. We find that the four ChAP proteins are derived from gene duplications after the divergence of Yarrowia from the remaining Saccharomycotina, with BC8 paralogues (Bch2 and Chs6) being more diverged relative to the BB8 paralogues (Bch1 and Bud7), suggesting neofunctionalization. Outside Ascomycota, a single preduplicate ChAP is present in nearly all Fungi and in diverse eukaryotes, but has been repeatedly lost. Chs5, however, is a fungal specific feature, appearing coincidentally with the loss of AP-4. In contrast, the ChAP protein is a wide-spread, yet uncharacterized, membrane-trafficking component, adding one more piece to the increasingly complex machinery deduced as being present in our ancient eukaryotic ancestor.

Interregulation of CDKA/CDK1 and the Plant-Specific Cyclin-Dependent Kinase CDKB in Control of the Chlamydomonas Cell Cycle

The cyclin-dependent kinase CDK1 is essential for mitosis in fungi and animals. Plant genomes contain the CDK1 ortholog CDKA and a plant kingdom-specific relative, CDKB. The green alga Chlamydomonas reinhardtii has a long G1 growth period followed by rapid cycles of DNA replication and cell division. We show that null alleles of CDKA extend the growth period prior to the first division cycle and modestly extend the subsequent division cycles, but do not prevent cell division, indicating at most a minor role for the CDK1 ortholog in mitosis in Chlamydomonas. A null allele of cyclin A has a similar though less extreme phenotype. In contrast, both CDKB and cyclin B are essential for mitosis. CDK kinase activity measurements imply that the predominant in vivo complexes are probably cyclin A-CDKA and cyclin B-CDKB. We propose a negative feedback loop: CDKA activates cyclin B-CDKB. Cyclin B-CDKB in turn promotes mitotic entry and inactivates cyclin A-CDKA. Cyclin A-CDKA and cyclin B-CDKB may redundantly promote DNA replication. We show that the anaphase-promoting complex is required for inactivation of both CDKA and CDKB and is essential for anaphase. These results are consistent with findings in Arabidopsis thaliana and may delineate the core of plant kingdom cell cycle control that, compared with the well-studied yeast and animal systems, exhibits deep conservation in some respects and striking divergence in others.

Papain-like cysteine proteases in Carica papaya: lineage-specific gene duplication and expansion

Background

Papain-like cysteine proteases (PLCPs), a large group of cysteine proteases structurally related to papain, play important roles in plant development, senescence, and defense responses. Papain, the first cysteine protease whose structure was determined by X-ray crystallography, plays a crucial role in protecting papaya from herbivorous insects. Except the four major PLCPs purified and characterized in papaya latex, the rest of the PLCPs in papaya genome are largely unknown.

Results

We identified 33 PLCP genes in papaya genome. Phylogenetic analysis clearly separated plant PLCP genes into nine subfamilies. PLCP genes are not equally distributed among the nine subfamilies and the number of PLCPs in each subfamily does not increase or decrease proportionally among the seven selected plant species. Papaya showed clear lineage-specific gene expansion in the subfamily III. Interestingly, all four major PLCPs purified from papaya latex, including papain, chymopapain, glycyl endopeptidase and caricain, were grouped into the lineage-specific expansion branch in the subfamily III. Mapping PLCP genes on chromosomes of five plant species revealed that lineage-specific expansions of PLCP genes were mostly derived from tandem duplications. We estimated divergence time of papaya PLCP genes of subfamily III. The major duplication events leading to lineage-specific expansion of papaya PLCP genes in subfamily III were estimated at 48 MYA, 34 MYA, and 16 MYA. The gene expression patterns of the papaya PLCP genes in different tissues were assessed by transcriptome sequencing and qRT-PCR. Most of the papaya PLCP genes of subfamily III expressed at high levels in leaf and green fruit tissues.

Conclusions

Tandem duplications played the dominant role in affecting copy number of PLCPs in plants. Significant variations in size of the PLCP subfamilies among species may reflect genetic adaptation of plant species to different environments. The lineage-specific expansion of papaya PLCPs of subfamily III might have been promoted by the continuous reciprocal selective effects of herbivore attack and plant defense.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4394-y) contains supplementary material, which is available to authorized users.

Intrinsic adaptive value and early fate of gene duplication revealed by a bottom-up approach

The population genetic mechanisms governing the preservation of gene duplicates, especially in the critical very initial phase, have remained largely unknown. Here, we demonstrate that gene duplication confers per se a weak selective advantage in scenarios of fitness trade-offs. Through a precise quantitative description of a model system, we show that a second gene copy serves to reduce gene expression inaccuracies derived from pervasive molecular noise and suboptimal gene regulation. We then reveal that such an accuracy in the phenotype yields a selective advantage in the order of 0.1% on average, which would allow the positive selection of gene duplication in populations with moderate/large sizes. This advantage is greater at higher noise levels and intermediate concentrations of the environmental molecule, when fitness trade-offs become more evident. Moreover, we discuss how the genome rearrangement rates greatly condition the eventual fixation of duplicates. Overall, our theoretical results highlight an original adaptive value for cells carrying new-born duplicates, broadly analyze the selective conditions that determine their early fates in different organisms, and reconcile population genetics with evolution by gene duplication.

A possible role of DNA methylation in functional divergence of a fast evolving duplicate gene encoding odorant binding protein 11 in the honeybee

Although gene duplication is seen as the main path to evolution of new functions, molecular mechanisms by which selection favours the gain versus loss of newly duplicated genes and minimizes the fixation of pseudo-genes are not well understood. Here, we investigate in detail a duplicate honeybee gene obp11 belonging to a fast evolving insect gene family encoding odorant binding proteins (OBPs). We report that obp11 is expressed only in female bees in rare antennal sensilla basiconica in contrast to its tandem partner obp10 that is expressed in the brain in both females and males (drones). Unlike all other obp genes in the honeybee, obp11 is methylated suggesting that functional diversification of obp11 and obp10 may have been driven by an epigenetic mechanism. We also show that increased methylation in drones near one donor splice site that correlates with higher abundance of a transcript variant encoding a truncated OBP11 protein is one way of controlling its contrasting expression. Our data suggest that like in mammals and plants, DNA methylation in insects may contribute to functional diversification of proteins produced from duplicated genes, in particular to their subfunctionalization by generating complementary patterns of expression.

Comprehensive analysis of correlations among codon usage bias, gene expression, and substitution rate in Arachis duranensis and Arachis ipaënsis orthologs

The relationship between evolutionary rates and gene expression in model plant orthologs is well documented. However, little is known about the relationships between gene expression and evolutionary trends in Arachis orthologs. We identified 7,435 one-to-one orthologs, including 925 single-copy and 6,510 multiple-copy sequences in Arachis duranensis and Arachis ipaënsis. Codon usage was stronger for shorter polypeptides, which were encoded by codons with higher GC contents. Highly expressed coding sequences had higher codon usage bias, GC content, and expression breadth. Additionally, expression breadth was positively correlated with polypeptide length, but there was no correlation between gene expression and polypeptide length. Inferred selective pressure was also negatively correlated with both gene expression and expression breadth in all one-to-one orthologs, while positively but non-significantly correlated with gene expression in sequences with signatures of positive selection. Gene expression levels and expression breadth were significantly higher for single-copy genes than for multiple-copy genes. Similarly, the gene expression and expression breadth in sequences with signatures of purifying selection were higher than those of sequences with positive selective signatures. These results indicated that gene expression differed between single-copy and multiple-copy genes as well as sequences with signatures of positive and purifying selection.

Genome Mutational and Transcriptional Hotspots Are Traps for Duplicated Genes and Sources of Adaptations

Gene duplication generates new genetic material, which has been shown to lead to major innovations in unicellular and multicellular organisms. A whole-genome duplication occurred in the ancestor of Saccharomyces yeast species but 92% of duplicates returned to single-copy genes shortly after duplication. The persisting duplicated genes in Saccharomyces led to the origin of major metabolic innovations, which have been the source of the unique biotechnological capabilities in the Baker’s yeast Saccharomyces cerevisiae. What factors have determined the fate of duplicated genes remains unknown. Here, we report the first demonstration that the local genome mutation and transcription rates determine the fate of duplicates. We show, for the first time, a preferential location of duplicated genes in the mutational and transcriptional hotspots of S. cerevisiae genome. The mechanism of duplication matters, with whole-genome duplicates exhibiting different preservation trends compared to small-scale duplicates. Genome mutational and transcriptional hotspots are rich in duplicates with large repetitive promoter elements. Saccharomyces cerevisiae shows more tolerance to deleterious mutations in duplicates with repetitive promoter elements, which in turn exhibit higher transcriptional plasticity against environmental perturbations. Our data demonstrate that the genome traps duplicates through the accelerated regulatory and functional divergence of their gene copies providing a source of novel adaptations in yeast.

Diurnal Cycling Transcription Factors of Pineapple Revealed by Genome-Wide Annotation and Global Transcriptomic Analysis

Circadian clock provides fitness advantage by coordinating internal metabolic and physiological processes to external cyclic environments. Core clock components exhibit daily rhythmic changes in gene expression, and the majority of them are transcription factors (TFs) and transcription coregulators (TCs). We annotated 1,398 TFs from 67 TF families and 80 TCs from 20 TC families in pineapple, and analyzed their tissue-specific and diurnal expression patterns. Approximately 42% of TFs and 45% of TCs displayed diel rhythmic expression, including 177 TF/TCs cycling only in the nonphotosynthetic leaf tissue, 247 cycling only in the photosynthetic leaf tissue, and 201 cycling in both. We identified 68 TF/TCs whose cycling expression was tightly coupled between the photosynthetic and nonphotosynthetic leaf tissues. These TF/TCs likely coordinate key biological processes in pineapple as we demonstrated that this group is enriched in homologous genes that form the core circadian clock in Arabidopsis and includes a STOP1 homolog. Two lines of evidence support the important role of the STOP1 homolog in regulating CAM photosynthesis in pineapple. First, STOP1 responds to acidic pH and regulates a malate channel in multiple plant species. Second, the cycling expression pattern of the pineapple STOP1 and the diurnal pattern of malate accumulation in pineapple leaf are correlated. We further examined duplicate-gene retention and loss in major known circadian genes and refined their evolutionary relationships between pineapple and other plants. Significant variations in duplicate-gene retention and loss were observed for most clock genes in both monocots and dicots.

Natural variation in stress response gene activity in the allopolyploid Arabidopsis suecica

BACKGROUND

Allopolyploids contain genomes composed of more than two complete sets of chromosomes that originate from at least two species. Allopolyploidy has been suggested as an important evolutionary mechanism that can lead to instant speciation. Arabidopsis suecica is a relatively recent allopolyploid species, suggesting that its natural accessions might be genetically very similar to each other. Nonetheless, subtle phenotypic differences have been described between different geographic accessions of A. suecica grown in a common garden.

RESULTS

To determine the degree of genomic similarity between different populations of A. suecica, we obtained transcriptomic sequence, quantified SNP variation within the gene space, and analyzed gene expression levels genome-wide from leaf material grown in controlled lab conditions. Despite their origin from the same progenitor species, the two accessions of A. suecica used in our study show genomic and transcriptomic variation. We report significant gene expression differences between the accessions, mostly in genes with stress-related functions. Among the differentially expressed genes, there are a surprising number of homoeologs coordinately regulated between sister accessions.

CONCLUSIONS

Many of these homoeologous genes and other differentially expressed genes affect transpiration and stomatal regulation, suggesting that they might be involved in the establishment of the phenotypic differences between the two accessions.

The house spider genome reveals an ancient whole-genome duplication during arachnid evolution

BACKGROUND

The duplication of genes can occur through various mechanisms and is thought to make a major contribution to the evolutionary diversification of organisms. There is increasing evidence for a large-scale duplication of genes in some chelicerate lineages including two rounds of whole genome duplication (WGD) in horseshoe crabs. To investigate this further, we sequenced and analyzed the genome of the common house spider Parasteatoda tepidariorum.

RESULTS

We found pervasive duplication of both coding and non-coding genes in this spider, including two clusters of Hox genes. Analysis of synteny conservation across the P. tepidariorum genome suggests that there has been an ancient WGD in spiders. Comparison with the genomes of other chelicerates, including that of the newly sequenced bark scorpion Centruroides sculpturatus, suggests that this event occurred in the common ancestor of spiders and scorpions, and is probably independent of the WGDs in horseshoe crabs. Furthermore, characterization of the sequence and expression of the Hox paralogs in P. tepidariorum suggests that many have been subject to neo-functionalization and/or sub-functionalization since their duplication.

CONCLUSIONS

Our results reveal that spiders and scorpions are likely the descendants of a polyploid ancestor that lived more than 450 MYA. Given the extensive morphological diversity and ecological adaptations found among these animals, rivaling those of vertebrates, our study of the ancient WGD event in Arachnopulmonata provides a new comparative platform to explore common and divergent evolutionary outcomes of polyploidization events across eukaryotes.

The house spider genome reveals an ancient whole-genome duplication during arachnid evolution

Background

The duplication of genes can occur through various mechanisms and is thought to make a major contribution to the evolutionary diversification of organisms. There is increasing evidence for a large-scale duplication of genes in some chelicerate lineages including two rounds of whole genome duplication (WGD) in horseshoe crabs. To investigate this further, we sequenced and analyzed the genome of the common house spider Parasteatoda tepidariorum.

Results

We found pervasive duplication of both coding and non-coding genes in this spider, including two clusters of Hox genes. Analysis of synteny conservation across the P. tepidariorum genome suggests that there has been an ancient WGD in spiders. Comparison with the genomes of other chelicerates, including that of the newly sequenced bark scorpion Centruroides sculpturatus, suggests that this event occurred in the common ancestor of spiders and scorpions, and is probably independent of the WGDs in horseshoe crabs. Furthermore, characterization of the sequence and expression of the Hox paralogs in P. tepidariorum suggests that many have been subject to neo-functionalization and/or sub-functionalization since their duplication.

Conclusions

Our results reveal that spiders and scorpions are likely the descendants of a polyploid ancestor that lived more than 450 MYA. Given the extensive morphological diversity and ecological adaptations found among these animals, rivaling those of vertebrates, our study of the ancient WGD event in Arachnopulmonata provides a new comparative platform to explore common and divergent evolutionary outcomes of polyploidization events across eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0399-x) contains supplementary material, which is available to authorized users.

Molecular Characterisation of a Novel Isoform of Hepatic Antimicrobial Peptide, Hepcidin (Le-Hepc), from Leiognathus equulus and Analysis of Its Functional Properties In Silico

Hepcidin represents a family of cysteine-rich antimicrobial peptides that are mainly expressed in the liver of living organisms. In this study, we have identified and characterised a novel isoform of hepcidin from the common pony fish, Leiognathus equulus (Le-Hepc). A 261-bp fragment cDNA coding for 86 amino acids was obtained. Homologous analysis showed that Le-Hepc belongs to the hepcidin super family and shares sequence identity with other known fish pre-propeptide hepcidin sequences. The ORF encodes for a 24-amino acid (aa) signal peptide coupled to a 36-aa prodomain followed by a 26-aa mature peptide. The mature peptide region has a calculated molecular weight of 2.73 kDa, a net positive charge of +2 and a theoretical pI of 8.23. Phylogenetic analysis of Le-Hepc showed a strong relationship with other fish hepcidin sequences and clustered into HAMP2 group hepcidins. Secondary structural analysis indicated that Le-Hepc mature peptide contains two antiparallel β-sheets strengthened by four disulphide bonds formed by eight conserved cysteine residues. The physicochemical properties of the peptide and its structural parameters are in agreement with characteristic features of an antimicrobial peptide. This is the first report of an antimicrobial peptide from the common pony fish, L. equulus.