The Many Nuanced Evolutionary Consequences of Duplicated Genes

A chromosome-scale genome of Rhus chinensis Mill. provides new insights into plant-insect interaction and gallotannins biosynthesis

Rhus chinensis Mill., an economically valuable Anacardiaceae species, is parasitized by the galling aphid Schlechtendalia chinensis, resulting in the formation of the Chinese gallnut (CG). Here, we report a chromosomal-level genome assembly of R. chinensis, with a total size of 389.40 Mb and scaffold N50 of 23.02 Mb. Comparative genomic and transcriptome analysis revealed that the enhanced structure of CG and nutritional metabolism contribute to improving the adaptability of R. chinensis to S. chinensis by supporting CG and galling aphid growth. CG was observed to be abundant in hydrolysable tannins (HT), particularly gallotannin and its isomers. Tandem repeat clusters of dehydroquinate dehydratase/shikimate dehydrogenase (DQD/SDH) and serine carboxypeptidase-like (SCPL) and their homologs involved in HT production were determined as specific to HT-rich species. The functional differentiation of DQD/SDH tandem duplicate genes and the significant contraction in the phenylalanine ammonia-lyase (PAL) gene family contributed to the accumulation of gallic acid and HT while minimizing the production of shikimic acid, flavonoids, and condensed tannins in CG. Furthermore, we identified one UDP glucosyltransferase (UGT84A), three carboxylesterase (CXE), and six SCPL genes from conserved tandem repeat clusters that are involved in gallotannin biosynthesis and hydrolysis in CG. We then constructed a regulatory network of these genes based on co-expression and transcription factor motif analysis. Our findings provide a genomic resource for the exploration of the underlying mechanisms of plant-galling insect interaction and highlight the importance of the functional divergence of tandem duplicate genes in the accumulation of secondary metabolites.

Mutational biases favor complexity increases in protein interaction networks after gene duplication

Biological systems can gain complexity over time. While some of these transitions are likely driven by natural selection, the extent to which they occur without providing an adaptive benefit is unknown. At the molecular level, one example is heteromeric complexes replacing homomeric ones following gene duplication. Here, we build a biophysical model and simulate the evolution of homodimers and heterodimers following gene duplication using distributions of mutational effects inferred from available protein structures. We keep the specific activity of each dimer identical, so their concentrations drift neutrally without new functions. We show that for more than 60% of tested dimer structures, the relative concentration of the heteromer increases over time due to mutational biases that favor the heterodimer. However, allowing mutational effects on synthesis rates and differences in the specific activity of homo- and heterodimers can limit or reverse the observed bias toward heterodimers. Our results show that the accumulation of more complex protein quaternary structures is likely under neutral evolution, and that natural selection would be needed to reverse this tendency.

Paralog transcriptional differentiation in the D. melanogaster-specific gene family Sdic across populations and spermatogenesis stages

How recently originated gene copies become stable genomic components remains uncertain as high sequence similarity of young duplicates precludes their functional characterization. The tandem multigene family Sdic is specific to Drosophila melanogaster and has been annotated across multiple reference-quality genome assemblies. Here we show the existence of a positive correlation between Sdic copy number and total expression, plus vast intrastrain differences in mRNA abundance among paralogs, using RNA-sequencing from testis of four strains with variable paralog composition. Single cell and nucleus RNA-sequencing data expose paralog expression differentiation in meiotic cell types within testis from third instar larva and adults. Additional RNA-sequencing across synthetic strains only differing in their Y chromosomes reveal a tissue-dependent trans-regulatory effect on Sdic: upregulation in testis and downregulation in male accessory gland. By leveraging paralog-specific expression information from tissue- and cell-specific data, our results elucidate the intraspecific functional diversification of a recently expanded tandem gene family.

Protein diversification through post-translational modifications, alternative splicing, and gene duplication

Proteins provide the basis for cellular function. Having multiple versions of the same protein within a single organism provides a way of regulating its activity or developing novel functions. Post-translational modifications of proteins, by means of adding/removing chemical groups to amino acids, allow for a well-regulated and controlled way of generating functionally distinct protein species. Alternative splicing is another method with which organisms possibly generate new isoforms. Additionally, gene duplication events throughout evolution generate multiple paralogs of the same genes, resulting in multiple versions of the same protein within an organism. In this review, we discuss recent advancements in the study of these three methods of protein diversification and provide illustrative examples of how they affect protein structure and function.

Whole-genome Duplications and the Long-term Evolution of Gene Regulatory Networks in Angiosperms

Angiosperms have a complex history of whole-genome duplications (WGDs), with varying numbers and ages of WGD events across clades. These WGDs have greatly affected the composition of plant genomes due to the biased retention of genes belonging to certain functional categories following their duplication. In particular, regulatory genes and genes encoding proteins that act in multiprotein complexes have been retained in excess following WGD. Here, we inferred protein-protein interaction (PPI) networks and gene regulatory networks (GRNs) for seven well-characterized angiosperm species and explored the impact of both WGD and small-scale duplications (SSDs) in network topology by analyzing changes in frequency of network motifs. We found that PPI networks are enriched in WGD-derived genes associated with dosage-sensitive intricate systems, and strong selection pressures constrain the divergence of WGD-derived genes at the sequence and PPI levels. WGD-derived genes in network motifs are mostly associated with dosage-sensitive processes, such as regulation of transcription and cell cycle, translation, photosynthesis, and carbon metabolism, whereas SSD-derived genes in motifs are associated with response to biotic and abiotic stress. Recent polyploids have higher motif frequencies than ancient polyploids, whereas WGD-derived network motifs tend to be disrupted on the longer term. Our findings demonstrate that both WGD and SSD have contributed to the evolution of angiosperm GRNs, but in different ways, with WGD events likely having a more significant impact on the short-term evolution of polyploids.

Dosage balance acts as a time-dependent selective barrier to subfunctionalization

BACKGROUND

Gene duplication is an important process for genome expansion, sometimes allowing for new gene functions to develop. Duplicate genes can be retained through multiple processes, either for intermediate periods of time through processes such as dosage balance, or over extended periods of time through processes such as subfunctionalization and neofunctionalization.

RESULTS

Here, we built upon an existing subfunctionalization Markov model by incorporating dosage balance to describe the interplay between subfunctionalization and dosage balance to explore selective pressures on duplicate copies. Our model incorporates dosage balance using a biophysical framework that penalizes the fitness of genetic states with stoichiometrically imbalanced proteins. These imbalanced states cause increased concentrations of exposed hydrophobic surface areas, which cause deleterious mis-interactions. We draw comparison between our Subfunctionalization + Dosage-Balance Model (Sub + Dos) and the previous Subfunctionalization-Only (Sub-Only) Model. This comparison includes how the retention probabilities change over time, dependent upon the effective population size and the selective cost associated with spurious interaction of dosage-imbalanced partners. We show comparison between Sub-Only and Sub + Dos models for both whole-genome duplication and small-scale duplication events.

CONCLUSION

These comparisons show that following whole-genome duplication, dosage balance serves as a time-dependent selective barrier to the subfunctionalization process, by causing an overall delay but ultimately leading to a larger portion of the genome retained through subfunctionalization. This higher percentage of the genome that is ultimately retained is caused by the alternative competing process, nonfunctionalization, being selectively blocked to a greater extent. In small-scale duplication, the reverse pattern is seen, where dosage balance drives faster rates of subfunctionalization, but ultimately leads to a smaller portion of the genome retained as duplicates. This faster rate of subfunctionalization is because the dosage balance of interacting gene products is negatively affected immediately after duplication and the loss of a duplicate restores the stoichiometric balance. Our findings provide support that the subfunctionalization of genes that are susceptible to dosage balance effects, such as proteins involved in complexes, is not a purely neutral process. With stronger selection against stoichiometrically imbalanced gene partners, the rates of subfunctionalization and nonfunctionalization slow; however, this ultimately leads to a greater proportion of subfunctionalized gene pairs.

Targeted, programmable, and precise tandem duplication in the mammalian genome

Tandem duplications are frequent structural variations of the genome and play important roles in genetic disease and cancer. However, interpreting the phenotypic consequences of tandem duplications remains challenging, in part owing to the lack of genetic tools to model such variations. Here, we developed a strategy, tandem duplication via prime editing (TD-PE), to create targeted, programmable, and precise tandem duplication in the mammalian genome. In this strategy, we design a pair of in trans prime editing guide RNAs (pegRNAs) for each targeted tandem duplication, which encode the same edits but prime the single-stranded DNA (ssDNA) extension in opposite directions. The reverse transcriptase (RT) template of each extension is designed homologous to the target region of the other single guide RNA (sgRNA) to promote the reannealing of the edited DNA strands and the duplication of the fragment in between. We showed that TD-PE produced robust and precise in situ tandem duplications of genomic fragments ranging from ∼50 bp to ∼10 kb, with a maximal efficiency up to 28.33%. By fine-tuning the pegRNAs, we achieved simultaneous targeted duplication and fragment insertion. Finally, we successfully produced multiple disease-relevant tandem duplications, showing the general utility of TD-PE in genetic research.

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.

Humanized yeast to model human biology, disease and evolution

For decades, budding yeast, a single-cellular eukaryote, has provided remarkable insights into human biology. Yeast and humans share several thousand genes despite morphological and cellular differences and over a billion years of separate evolution. These genes encode critical cellular processes, the failure of which in humans results in disease. Although recent developments in genome engineering of mammalian cells permit genetic assays in human cell lines, there is still a need to develop biological reagents to study human disease variants in a high-throughput manner. Many protein-coding human genes can successfully substitute for their yeast equivalents and sustain yeast growth, thus opening up doors for developing direct assays of human gene function in a tractable system referred to as 'humanized yeast'. Humanized yeast permits the discovery of new human biology by measuring human protein activity in a simplified organismal context. This Review summarizes recent developments showing how humanized yeast can directly assay human gene function and explore variant effects at scale. Thus, by extending the 'awesome power of yeast genetics' to study human biology, humanizing yeast reinforces the high relevance of evolutionarily distant model organisms to explore human gene evolution, function and disease.

Retention of duplicated genes in evolution

Gene duplication is a prevalent phenomenon across the tree of life. The processes that lead to the retention of duplicated genes are not well understood. Functional genomics approaches in model organisms, such as yeast, provide useful tools to test the mechanisms underlying retention with functional redundancy and divergence of duplicated genes, including fates associated with neofunctionalization, subfunctionalization, back-up compensation, and dosage amplification. Duplicated genes may also be retained as a consequence of structural and functional entanglement. Advances in human gene editing have enabled the interrogation of duplicated genes in the human genome, providing new tools to evaluate the relative contributions of each of these factors to duplicate gene retention and the evolution of genome structure.

Learning the local landscape of protein structures with convolutional neural networks

One fundamental problem of protein biochemistry is to predict protein structure from amino acid sequence. The inverse problem, predicting either entire sequences or individual mutations that are consistent with a given protein structure, has received much less attention even though it has important applications in both protein engineering and evolutionary biology. Here, we ask whether 3D convolutional neural networks (3D CNNs) can learn the local fitness landscape of protein structure to reliably predict either the wild-type amino acid or the consensus in a multiple sequence alignment from the local structural context surrounding site of interest. We find that the network can predict wild type with good accuracy, and that network confidence is a reliable measure of whether a given prediction is likely going to be correct or not. Predictions of consensus are less accurate and are primarily driven by whether or not the consensus matches the wild type. Our work suggests that high-confidence mis-predictions of the wild type may identify sites that are primed for mutation and likely targets for protein engineering.

Whole-genome sequencing, genome mining, metabolic reconstruction and evolution of pentachlorophenol and other xenobiotic degradation pathways in Bacillus tropicus strain AOA-CPS1

A pentachlorophenol degrading bacterium was isolated from effluent of a wastewater treatment plant in Durban, South Africa, and identified as Bacillus tropicus strain AOA-CPS1 (BtAOA). The isolate degraded 29% of pentachlorophenol (PCP) within 9 days at an initial PCP concentration of 100 mg L^-1 and 62% of PCP when the initial concentration was set at 350 mg L^-1. The whole-genome of BtAOA was sequenced using Pacific Biosciences RS II sequencer with the Single Molecule, Real-Time (SMRT) Link (version 7.0.1.66975) and analysed using the HGAP4-de-novo assembly application. The contigs were annotated at NCBI, RASTtk and PROKKA prokaryotic genome annotation pipelines. The BtAOA genome is comprised of a 5,246,860-bp chromosome and a 58,449-bp plasmid with a GC content of 35.4%. The metabolic reconstruction for BtAOA showed that the organism has been naturally exposed to various chlorophenolic compounds including PCP and other xenobiotics. The chromosome encodes genes for core processes, stress response and PCP catabolic genes. Analogues of PCP catabolic gene (cpsBDCAE, and p450) sequences were identified from the NCBI annotation data, PCR-amplified from the whole genome of BtAOA, cloned into pET15b expression vector, overexpressed in E. coli BL21 (DE3) expression host, purified and characterized. Sequence mining and comparative analysis of the metabolic reconstruction of the BtAOA genome with closely related strains suggests that the operon encoding the first two enzymes in the PCP degradation pathway were acquired from a pre-existing pterin-carbinolamine dehydratase subsystem. The other two enzymes were recruited via horizontal gene transfer (HGT) from the pool of hypothetical proteins with no previous specific function, while the last enzyme was recruited from pre-existing enzymes from the TCA or serine-glyoxalase cycle via HGT events. This study provides a comprehensive understanding of the role of BtAOA in PCP degradation and its potential exploitation for bioremediation of other xenobiotic compounds.

GhN/AINV13 positively regulates cotton stress tolerance by interacting with the 14-3-3 protein

Neutral/alkaline invertases (N/AINVs) are sucrose hydrolases with important roles in plants. In this study, 15, 15, 15, 29, and 30 N/AINVs were identified in the Gossypium species, G. raimondii, G. herbaceum, G. arboreum, G. hirsutum, and G. barbadense, respectively. Along with two previously discovered branches, α and β, a new clade γ was first discovered in our study. Investigation of gene collinearity showed that whole-genome duplication (WGD) and polyploidization were responsible for the expansion of the N/AINV gene family in allopolyploid Gossypium. Moreover, expression patterns revealed that GhN/AINV3/13/17/23/24/28 from the β clade is highly expressed during the period of fiber initiation. The invertase activity of GhN/AINV13 and GhN/AINV23 were confirmed by restoring defects of invertase-deficient yeast mutant SEY2102. Treatments of abiotic stress showed that most GhN/AINVs were induced in response to polyethylene glycol (PEG) or salt stress. A virus-induced gene-silencing (VIGS) experiment and yeast two-hybrid assay demonstrated that GhN/AINV13 may interact with their positive regulators Gh14-3-3 proteins and participate in the fiber initiation or stress tolerance of cotton. Our results provided fundamental information regarding N/AINVs and highlight their potential functions in cotton stress tolerance.

Exploring the complexity of soybean (Glycine max) transcriptional regulation using global gene co-expression networks

MAIN CONCLUSION

We report a soybean gene co-expression network built with data from 1284 RNA-Seq experiments, which was used to identify important regulators, modules and to elucidate the fates of gene duplicates. Soybean (Glycine max (L.) Merr.) is one of the most important crops worldwide, constituting a major source of protein and edible oil. Gene co-expression networks (GCN) have been extensively used to study transcriptional regulation and evolution of genes and genomes. Here, we report a soybean GCN using 1284 publicly available RNA-Seq samples from 15 distinct tissues. We found modules that are differentially regulated in specific tissues, comprising processes such as photosynthesis, gluconeogenesis, lignin metabolism, and response to biotic stress. We identified transcription factors among intramodular hubs, which probably integrate different pathways and shape the transcriptional landscape in different conditions. The top hubs for each module tend to encode proteins with critical roles, such as succinate dehydrogenase and RNA polymerase subunits. Importantly, gene essentiality was strongly correlated with degree centrality and essential hubs were enriched in genes involved in nucleic acids metabolism and regulation of cell replication. Using a guilt-by-association approach, we predicted functions for 93 of 106 hubs without functional description in soybean. Most of the duplicated genes had different transcriptional profiles, supporting their functional divergence, although paralogs originating from whole-genome duplications (WGD) are more often preserved in the same module than those from other mechanisms. Together, our results highlight the importance of GCN analysis in unraveling key functional aspects of the soybean genome, in particular those associated with hub genes and WGD events.

Genomic architecture and evolutionary antagonism drive allelic expression bias in the social supergene of red fire ants

Supergene regions maintain alleles of multiple genes in tight linkage through suppressed recombination. Despite their importance in determining complex phenotypes, our empirical understanding of early supergene evolution is limited. Here we focus on the young ‘social’ supergene of fire ants, a powerful system for disentangling the effects of evolutionary antagonism and suppressed recombination. We hypothesize that gene degeneration and social antagonism shaped the evolution of the fire ant supergene, resulting in distinct patterns of gene expression. We test these ideas by identifying allelic differences between supergene variants, characterizing allelic expression across populations, castes and body parts, and contrasting allelic expression biases with differences in expression between social forms. We find strong signatures of gene degeneration and gene-specific dosage compensation. On this background, a small portion of the genes has the signature of adaptive responses to evolutionary antagonism between social forms.

Red fire ants (Solenopsis invicta) are native to South America, but the species has spread to North America, Australia and New Zealand where it can be an invasive pest. A reason for this species’ invasiveness types of colonies : one with a single egg-laying queen and another with several queens. However, it is not possible to simply add more queens to a colony with one queen. Instead, the number of queens in a colony is controlled genetically, by a chromosome known as the ‘social chromosome’.

Like many other animals, red fire ants are diploid: their cells have two copies of each chromosome, which can carry two different versions of each gene. The social chromosome is no different, and it comes in two variants, SB and Sb. Each ant can therefore have either two SB chromosomes, leading to a colony with a single queen; or one SB chromosome and one Sb chromosome, leading to a colony with multiple queens. Ants with two copies of the Sb variant die when they are young, so the Sb version is inherited in a similar way to how the Y chromosome is passed on in humans. However, the social chromosome in red fire ants appeared less than one million years ago, making it much younger than the human Y chromosome, which is 180 million years old. This makes the social chromosome a good candidate for examining the early evolution of special chromosome variants that are only inherited.

How differences between the SB and the Sb chromosomes are evolving is an open question, however. Perhaps each version of the social chromosome has been optimised through natural selection to one colony type. Another suggestion is that the Sb chromosome has degenerated over time because its genes cannot be ‘reshuffled’ as they would be on normal chromosomes.

Martinez-Ruiz et al. compared genetic variants on the SB and Sb chromosomes, along with their expression in different types of ant colonies. The analysis showed that the Sb variant is in fact breaking down because of the lack of gene shuffling. This loss is compensated by intact copies of the same genes found on the SB variant, which explains why ants with the Sb variant can only survive if they also carry the SB version. Only a handful of genes on the social chromosomes appear to have been optimised by natural selection. Therefore Martinez-Ruiz et al. concluded the differences between the two chromosomes that lead to different colony types are collateral effects of Sb’s inability to reshuffle its genes.

This work reveals how a special chromosome similar to the Y chromosome in humans evolved. It also shows how multiple complex evolutionary forces can shape a species’ genetic makeup and social forms.

Humanization of yeast genes with multiple human orthologs reveals functional divergence between paralogs

Despite over a billion years of evolutionary divergence, several thousand human genes possess clearly identifiable orthologs in yeast, and many have undergone lineage-specific duplications in one or both lineages. These duplicated genes may have been free to diverge in function since their expansion, and it is unclear how or at what rate ancestral functions are retained or partitioned among co-orthologs between species and within gene families. Thus, in order to investigate how ancestral functions are retained or lost post-duplication, we systematically replaced hundreds of essential yeast genes with their human orthologs from gene families that have undergone lineage-specific duplications, including those with single duplications (1 yeast gene to 2 human genes, 1:2) or higher-order expansions (1:>2) in the human lineage. We observe a variable pattern of replaceability across different ortholog classes, with an obvious trend toward differential replaceability inside gene families, and rarely observe replaceability by all members of a family. We quantify the ability of various properties of the orthologs to predict replaceability, showing that in the case of 1:2 orthologs, replaceability is predicted largely by the divergence and tissue-specific expression of the human co-orthologs, i.e., the human proteins that are less diverged from their yeast counterpart and more ubiquitously expressed across human tissues more often replace their single yeast ortholog. These trends were consistent with in silico simulations demonstrating that when only one ortholog can replace its corresponding yeast equivalent, it tends to be the least diverged of the pair. Replaceability of yeast genes having more than 2 human co-orthologs was marked by retention of orthologous interactions in functional or protein networks as well as by more ancestral subcellular localization. Overall, we performed >400 human gene replaceability assays, revealing 50 new human-yeast complementation pairs, thus opening up avenues to further functionally characterize these human genes in a simplified organismal context.

Molecular Evolution of the Glutathione S-Transferase Family in the Bemisia tabaci Species Complex

The glutathione S-transferase (GST) family plays an important role in the adaptation of herbivorous insects to new host plants and other environmental constrains. The family codes for enzymes that neutralize reactive oxygen species and phytotoxins through the conjugation of reduced glutathione. Here, we studied the molecular evolution of the GST family in Bemisia tabaci, a complex of >35 sibling species, differing in their geographic and host ranges. We tested if some enzymes evolved different functionality, by comparing their sequences in six species, representing five of the six major genetic clades in the complex. Comparisons of the nonsynonymous to synonymous substitution ratios detected positive selection events in 11 codons of 5 cytosolic GSTs. Ten of them are located in the periphery of the GST dimer, suggesting a putative involvement in interactions with other proteins. Modeling the tertiary structure of orthologous enzymes, identified additional 19 mutations in 9 GSTs, likely affecting the enzymes' functionality. Most of the mutation events were found in the environmentally responsive classes Delta and Sigma, indicating a slightly different delta/sigma tool box in each species. At a broader genomic perspective, our analyses indicated a significant expansion of the Delta GST class in B. tabaci and a general association between the diet breadth of hemipteran species and their total number of GST genes. We raise the possibility that at least some of the identified changes improve the fitness of the B. tabaci species carrying them, leading to their better adaptation to specific environments.

Different origins of paralogues of salmonid TNR1 and TNFR2: Characterisation and expression analysis of four TNF receptor genes in rainbow trout Oncorhynchus mykiss

Mammalian TNFR1 and TNFR2 bind TNFα and TNFβ, and provide key communication signals to a variety of cell types during development and immune responses that are crucial for cell survival, proliferation and apoptosis. In teleost fish TNFβ is absent but TNFα has been expanded by the third whole genome duplication (3R WGD) and again by a 4R WGD in some lineages, leading to the four TNFα paralogues known in salmonids. Two paralogues for each of TNFR1 and TNFR2 have been cloned in rainbow trout in this study and are present in other salmonid genomes. Whilst the TNFR2 paralogues were generated via the 4R salmonid WGD, the TNFR1 paralogues arose from a local en bloc duplication. Functional diversification of TNFR paralogues was evidenced by differential gene expression and modulation, upstream ATGs affecting translation, ATTTA motifs in the 3'-UTR regulating mRNA stability, and post-translational modification by N-glycosylation. Trout TNFR are highly expressed in immune tissues/organs, and other tissues, in a gene- and tissue-specific manner. Furthermore, their expression is differentially modulated by PAMPs and cytokines in a cell type- and stimulant-specific manner. Such findings suggest an important role of the TNF/TNFR axis in the immune response and other physiological processes in fish.

Common Carp mef2 Genes: Evolution and Expression

The MEF2 (myocyte enhancer factor 2) family belongs to the MADS-box superfamily of eukaryotic transcription factors. The vertebrate genes compose four distinct subfamilies designated MEF2A, -B, -C, and -D. There are multiple mef2 genes in the common carp (Cyprinus carpio). So far, the embryonic expression patterns of these genes and the evolution of fish mef2 genes have been barely investigated. In this study, we completed the coding information of C. carpio mef2ca2 and mef2d1 genes via gene cloning and presented two mosaic mef2 sequences as evidence for recombination. We also analyzed the phylogenetic relationship and conserved synteny of mef2 genes and proposed a new evolutionary scenario. In our version, MEF2B and the other three vertebrate subfamilies were generated in parallel from the single last ancestor via two rounds of whole genome duplication events that occurred at the dawn of vertebrates. Moreover, we examined the expression patterns of C. carpio mef2 genes during embryogenesis, by using whole-mount in situ hybridization, and found the notochord to be a new expression site for these genes except for mef2ca1&2. Our results thus provide new insights into the evolution and expression of mef2 genes.