Selection and gene duplication: a view from the genome

Genome-wide identification of the EIN3/EIL gene family in Ginkgo biloba and functional study of a GbEIL in the ethylene signaling pathway

ETHYLENE-INSENSITIVE3 (EIN3) or EIN3-Like (EIL) proteins, play critical roles in integrating ethylene signaling and physiological regulation in plants by modulating the expression of various downstream genes, such as ethylene-response factors (ERFs). However, little is known about the characteristics of EIN3/EILs in the gymnosperm Ginkgo biloba. In the present study, a genome-wide comparative analysis of Ginkgo EIN3/EIL gene family was performed with those from an array of species, including bryophytes (Physcomitrella patens), gymnosperms (Cycas panzhihuaensis), and angiosperms (Arabidopsis thaliana, Gossypium raimondii, Gossypium hirsutum, Oryza sativa, and Brachypodium distachyon). Within the constructed phylogenetic tree for the 53 EIN3/EILs identified, 5 GbEILs from G. biloba, 2 PpEILs from P. patens, and 3 CpEILs from C. panzhihuaensis were assigned to one cluster, suggesting that their derivation occurred after the split of their ancestors and angiosperms. Although considerable divergence accumulated in amino acid sequences along with the evolutionary process, the specific EIN3_DNA-binding domains were evolutionarily conserved among the 53 EIN3/EILs. Collinearity analysis indicated that whole-genome or segmental duplication and subsequent purifying selection might have prompted the generation and evolution of EIN3/EIL multigene families. Based on the expression patterns of five GbEILs at the four developmental stages of Ginkgo ovules, one GbEIL gene (Gb_03292) was further investigated for its role in mediating ethylene signaling. The functional activity of Gb_03292 was closely related to ethylene signaling, as it complemented the triple response via ectopic expression in ein3eil1 double mutant Arabidopsis. Additionally, GbEIL likely modulates the expression of a Ginkgo ERF (Gb_15517) by directly binding to its promoter. These results demonstrated that the GbEIL gene could have participated in mediating ethylene signal transduction during ovule development in G. biloba. The present study also provides insights into the conservation of ethylene signaling across the gymnosperm G. biloba and angiosperm species.

Evolutionary analysis of species-specific duplications in flatworm genomes

Platyhelminthes, also known as flatworms, is a phylum of bilaterian invertebrates infamous for their parasitic representatives. The classes Cestoda, Monogenea, and Trematoda comprise parasitic helminths inhabiting multiple hosts, including fishes, humans, and livestock, and are responsible for considerable economic damage and burden on human health. As in other animals, the genomes of flatworms have a wide variety of paralogs, genes related via duplication, whose origins could be mapped throughout the evolution of the phylum. Through in-silico analysis, we studied inparalogs, i.e., species-specific duplications, focusing on their biological functions, expression changes, and evolutionary rate. These genes are thought to be key players in the adaptation process of species to each particular niche. Our results showed that genes related with specific functional terms, such as response to stress, transferase activity, oxidoreductase activity, and peptidases, are overrepresented among inparalogs. This trend is conserved among species from different classes, including free-living species. Available expression data from Schistosoma mansoni, a parasite from the trematode class, demonstrated high conservation of expression patterns between inparalogs, but with notable exceptions, which also display evidence of rapid evolution. We discuss how natural selection may operate to maintain these genes and the particular duplication models that fit better to the observations. Our work supports the critical role of gene duplication in the evolution of flatworms, representing the first study of inparalogs evolution at the genome-wide level in this group.

Genomic and Transcriptomic Insights into the Evolution of C4 Photosynthesis in Grasses

C4 photosynthesis has independently evolved over 62 times within 19 angiosperm families. The recurrent evolution of C4 photosynthesis appears to contradict the complex anatomical and biochemical modifications required for the transition from C3 to C4 photosynthesis. In this study, we conducted an integrated analysis of genomics and transcriptomics to elucidate the molecular underpinnings of convergent C4 evolution in the grass family. Our genome-wide exploration of C4-related gene families suggests that the expansion of these gene families may have played an important role in facilitating C4 evolution in the grass family. A phylogenomic synteny network analysis uncovered the emergence of C4 genes in various C4 grass lineages from a common ancestral gene pool. Moreover, through a comparison between non-C4 and C4  PEPCs, we pinpointed 14 amino acid sites exhibiting parallel adaptations. These adaptations, occurring post the BEP-PACMAD divergence, shed light on why all C4 origins in grasses are confined to the PACMAD clade. Furthermore, our study revealed that the ancestor of Chloridoideae grasses possessed a more favorable molecular preadaptation for C4 functions compared to the ancestor of Panicoideae grasses. This molecular preadaptation potentially explains why C4 photosynthesis evolved earlier in Chloridoideae than in Panicoideae and why the C3-to-C4 transition occurred once in Chloridoideae but multiple times in Panicoideae. Additionally, we found that C4 genes share similar cis-elements across independent C4 lineages. Notably, NAD-ME subtype grasses may have retained the ancestral regulatory machinery of the C4  NADP-ME gene, while NADP-ME subtype grasses might have undergone unique cis-element modifications.

Vole genomics links determinate and indeterminate growth of teeth

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Different species of voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars that have similar size and shape, providing alternative models for studying roots. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation. Bulk transcriptomics analyses of embryonic molar development in bank voles also demonstrated conserved patterns of dental gene expression compared to mice, with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars. Our results support ongoing evolution of dental genes across Glires, revealing the complex evolutionary background of convergent evolution for ever-growing molars.

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

INTRODUCTION

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

METHODS

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

RESULTS

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

DISCUSSION

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

Epigenomic divergence correlates with sequence polymorphism in Arabidopsis paralogs

Processes affecting rates of sequence polymorphism are fundamental to the evolution of gene duplicates. The relationship between gene activity and sequence polymorphism can influence the likelihood that functionally redundant gene copies are co-maintained in stable evolutionary equilibria vs other outcomes such as neofunctionalization. Here, we investigate genic variation in epigenome-associated polymorphism rates in Arabidopsis thaliana and consider whether these affect the evolution of gene duplicates. We compared the frequency of sequence polymorphism and patterns of genetic differentiation between genes classified by exon methylation patterns: unmethylated (unM), gene-body methylated (gbM), and transposon-like methylated (teM) states, which reflect divergence in gene expression. We found that the frequency of polymorphism was higher in teM (transcriptionally repressed, tissue-specific) genes and lower in gbM (active, constitutively expressed) genes. Comparisons of gene duplicates were largely consistent with genome-wide patterns - gene copies that exhibit teM accumulate more variation, evolve faster, and are in chromatin states associated with reduced DNA repair. This relationship between expression, the epigenome, and polymorphism may lead to the breakdown of equilibrium states that would otherwise maintain genetic redundancies. Epigenome-mediated polymorphism rate variation may facilitate the evolution of novel gene functions in duplicate paralogs maintained over evolutionary time.

Genome-wide analysis of annexin gene family in Schrenkiella parvula and Eutrema salsugineum suggests their roles in salt stress response

Annexins (Anns) play an important role in plant development, growth and responses to various stresses. Although Ann genes have been characterized in some plants, their role in adaptation mechanisms and tolerance to environmental stresses have not been studied in extremophile plants. In this study, Ann genes in Schrenkiella parvula and Eutrema salsugineum were identified using a genome-wide method and phylogenetic relationships, subcellular distribution, gene structures, conserved residues and motifs and also promoter prediction have been studied through bioinformatics analysis. We identified ten and eight encoding putative Ann genes in S. parvula and E. salsugineum genome respectively, which were divided into six subfamilies according to phylogenetic relationships. By observing conservation in gene structures and protein motifs we found that the majority of Ann members in two extremophile plants are similar. Furthermore, promoter analysis revealed a greater number of GATA, Dof, bHLH and NAC transcription factor binding sites, as well as ABRE, ABRE3a, ABRE4, MYB and Myc cis-acting elements in compare to Arabidopsis thaliana. To gain additional insight into the putative roles of candidate Ann genes, the expression of SpAnn1, SpAnn2 and SpAnn6 in S. parvula was studied in response to salt stress, which indicated that their expression level in shoot increased. Similarly, salt stress induced expression of EsAnn1, 5 and 7, in roots and EsAnn1, 2 and 5 in leaves of E. salsugineum. Our comparative analysis implies that both halophytes have different regulatory mechanisms compared to A. thaliana and suggest SpAnn2 gene play important roles in mediating salt stress.

Architectural groups of a subtelomeric gene family evolve along distinct paths in Candida albicans

Subtelomeres are dynamic genomic regions shaped by elevated rates of recombination, mutation, and gene birth/death. These processes contribute to formation of lineage-specific gene family expansions that commonly occupy subtelomeres across eukaryotes. Investigating the evolution of subtelomeric gene families is complicated by the presence of repetitive DNA and high sequence similarity among gene family members that prevents accurate assembly from whole genome sequences. Here, we investigated the evolution of the telomere-associated (TLO) gene family in Candida albicans using 189 complete coding sequences retrieved from 23 genetically diverse strains across the species. Tlo genes conformed to the 3 major architectural groups (α/β/γ) previously defined in the genome reference strain but significantly differed in the degree of within-group diversity. One group, Tloβ, was always found at the same chromosome arm with strong sequence similarity among all strains. In contrast, diverse Tloα sequences have proliferated among chromosome arms. Tloγ genes formed 7 primary clades that included each of the previously identified Tloγ genes from the genome reference strain with 3 Tloγ genes always found on the same chromosome arm among strains. Architectural groups displayed regions of high conservation that resolved newly identified functional motifs, providing insight into potential regulatory mechanisms that distinguish groups. Thus, by resolving intraspecies subtelomeric gene variation, it is possible to identify previously unknown gene family complexity that may underpin adaptive functional variation.

When LUCA met gnomAD: genetic constraints on universal genes in humans

LUCA, the last universal common ancestor, is the hypothetical most recent common ancestor of the three domains of life which share the universal genes (UG). It seems interesting to evaluate whether the UG phylogeny has had an impact on current Human gene constraints. A list of human homologs of UG was retrieved from the eggNOG database. We analyzed this LUCA gene (LG) group, and a random sample of 500 genes from the gnomAD database (RG group). Gene constraint metrics were retrieved from gnomAD and associations with Mendelian diseases and modes of inheritance were retrieved from OMIM. The LG group consisted of 277 genes and the RG group, 492 (8 genes were in LG group). 38.6% of the genes in the LG group and 25.2% of the genes in the RG group were associated with a Mendelian disease (p < 0.0001). The mode of inheritance was more often autosomal recessive (69.0 vs. 50.5%), and less often autosomal dominant (19.0 vs. 31.3%), or mixed (6.0 vs. 12.1%) for those associated with the LG group (p = 0.048). The LG group was significantly more constrained for missense variants (MOEUF, 0.919 vs. 0.997, p < 0.0001) and was borderline significantly more constrained for loss-of-function variants (LOEUF, 0.872 vs. 0.947, p = 0.051). These results suggest that the UG in humans differs from the rest of the genome in terms of constraints and associated Mendelian diseases. It suggests that phylogenic data can explain some of the characteristics of human genes and could help in interpreting variants.

Characterization of the Calmodulin/Calmodulin-like Protein (CAM/CML) Family in Ginkgo biloba, and the Influence of an Ectopically Expressed GbCML Gene (Gb_30819) on Seedling and Fruit Development of Transgenic Arabidopsis

Calmodulins (CAMs) and calmodulin-like proteins (CMLs) can participate in the regulation of various physiological processes via sensing and decoding Ca²⁺ signals. To reveal the characteristics of the CAM/CML family in Ginkgo biloba, a comprehensive analysis was performed at the genome-wide level. A total of 26 CAMs/CMLs, consisting of 5 GbCAMs and 21 GbCMLs, was identified on 11 out of 12 chromosomes in G. biloba. They displayed a certain degree of multiplicity in their sequences, albeit with conserved EF hands. Collinearity analysis suggested that tandem rather than segmental or whole-genome duplications were likely to play roles in the evolution of the Ginkgo CAM/CML family. Furthermore, GbCAMs/GbCMLs were grouped into higher, lower, and moderate expression in magnitude. The cis-acting regulatory elements involved in phytohormone-responsiveness within GbCAM/GbCML promotors may explain their varied expression profiles. The ectopic expression of a GbCML gene (Gb_30819) in transgenic Arabidopsis led to phenotypes with significantly shortened root length and seedling height, and decreased yields of both pods and seeds. Moreover, an electrophoresis mobility shift assay demonstrated the Ca²⁺-binding activity of Gb_30819 in vitro. Altogether, these results contribute to insights into the characteristics of the evolution and expression of GbCAMs/GbCMLs, as well as evidence for Ca²⁺-CAM/CML pathways functioning within the ancient gymnosperm G. biloba.

Does effective population size affect rates of molecular evolution: Mitochondrial data for host/parasite species pairs in bees suggests not

Adaptive evolutionary theory argues that organisms with larger effective population size (N e) should have higher rates of adaptive evolution and therefore greater capacity to win evolutionary arm races. However, in some certain cases, species with much smaller N e may be able to survive besides their opponents for an extensive evolutionary time. Neutral theory predicts that accelerated rates of molecular evolution in organisms with exceedingly small N e are due to the effects of genetic drift and fixation of slightly deleterious mutations. We test this prediction in two obligate social parasite species and their respective host species from the bee tribe Allodapini. The parasites (genus Inquilina) have been locked into tight coevolutionary arm races with their exclusive hosts (genus Exoneura) for ~15 million years, even though Inquilina exhibit N e that are an order of magnitude smaller than their host. In this study, we compared rates of molecular evolution between host and parasite using nonsynonymous to synonymous substitution rate ratios (dN/dS) of eleven mitochondrial protein-coding genes sequenced from transcriptomes. Tests of selection on mitochondrial genes indicated no significant differences between host and parasite dN/dS, with evidence for purifying selection acting on all mitochondrial genes of host and parasite species. Several potential factors which could weaken the inverse relationship between N e and rate of molecular evolution are discussed.

Congruence between ocean-dispersal modelling and phylogeography explains recent evolutionary history of Cycas species with buoyant seeds

Ocean currents play a significant role in driving the long-distance dispersal (LDD), spatial distribution and phylogeographic patterns of many organisms. Integrating phylogeographic analyses and mechanistic ocean current modelling can provide novel insights into the evolutionary history of terrestrial littoral species but has been rarely applied in this context. We focused on a group of Cycas that have buoyant seeds and occupy coastal habitats. By integrating evidence from mechanistic simulations and whole plastomic data, we examined the role of ocean circulation in shaping the phylogeography of these Cycas species. Plastomes of the studied Cycas species showed extreme conservatism, following a post-Pleistocene divergence. Phylogenies revealed three subclades, corresponding to the Pacific Ocean, Sunda Shelf and Indian Ocean. The ocean modelling results indicate that hotspots of seed stranding coincide well with the contemporary distribution of the Cycas species and that drifting trajectories from the three subclades are largely confined to separate regions. These findings suggest that ocean current systems, by driving long-distance dispersal, have shaped the distribution and phylogeography for Cycas with buoyant seeds. This study highlights how the combination of genomic data and ocean drift modelling can help explain phylogeographic patterns and diversity in terrestrial littoral ecosystems.

The expansion and diversity of the CYP75 gene family in Vitaceae

The CYP75 gene family plays an important role in flavonoid biosynthesis in plants. Little is known about the evolution of the gene family within the grape family. Here, we extracted the CYP75 genes from transcriptome data of 15 grape species and 36 representative genomes from other plants to explore the evolutionary history of the CYP75 gene family in Vitaceae. The structure of the CYP75 protein sequences is highly conserved with the variation mainly occurring in the N terminal and the middle region. The evolutionary analyses suggested classifying the CYP75 gene family into three groups in Vitaceae, namely Vitaceae A1, Vitaceae A2 and Vitaceae B. The Vitaceae A1 and A2 belong to the CYP75A subfamily and the Vitaceae B belongs to the CYP75B subfamily. Within the Vitaceae A1, most Vitaceae taxa present only one copy of the CYP75A protein sequence except for Vitis vinifera with a high number of sequences, which might have originated through recent gene duplications after its split from the other species. Vitaceae A2 contain only CYP75A sequences from Vitaceae sister to one from Camellia sinensis, probably representing a relict lineage. The CYP75B proteins were found to be dominated in Vitaceae and other angiosperms. Our results provide important insights into understanding the evolutionary history of the CYP75 gene family in Vitaceae and other angiosperms.

Scale-Free Biology: Integrating Evolutionary and Developmental Thinking

When the history of life on earth is viewed as a history of cell division, all of life becomes a single cell lineage. The growth and differentiation of this lineage in reciprocal interaction with its environment can be viewed as a developmental process; hence the evolution of life on earth can also be seen as the development of life on earth. Here, in reviewing this field, some potentially fruitful research directions suggested by this change in perspective are highlighted. Variation and selection become, for example, bidirectional information flows between scales, while the notions of "cooperation" and "competition" become scale relative. The language of communication, inference, and information processing becomes more useful than the language of causation to describe the interactions of both homogeneous and heterogeneous living systems at any scale. Emerging scale-free theoretical frameworks such as predictive coding and active inference provide conceptual tools for reconceptualizing biology as the study of a unified, multiscale dynamical system.

Genome Organization and Adaptive Potential of Archetypal Organophosphate Degrading Sphingobium fuliginis ATCC 27551

Sphingobium fuliginis ATCC 27551, previously classified as Flavobacterium sp. ATCC 27551, degrades neurotoxic organophosphate insecticides and nerve agents through the activity of a membrane-associated organophosphate hydrolase. This study was designed to determine the complete genome sequence of S. fuliginis ATCC 27551 to unravel its degradative potential and adaptability to harsh environments. The 5,414,624 bp genome with a GC content of 64.4% is distributed between two chromosomes and four plasmids and encodes 5,557 proteins. Of the four plasmids, designated as pSF1, pSF2, pSF3, and pSF4, only two (pSF1 and pSF2) are self-transmissible and contained the complete genetic repertoire for a T4SS. The other two plasmids (pSF3 and pSF4) are mobilizable and both showed the presence of an oriT and relaxase-encoding sequences. The sequence of plasmid pSF3 coincided with the previously determined sequence of pPDL2 and included an opd gene encoding organophosphate hydrolase as a part of the mobile element. About 15,455 orthologous clusters were identified from among the cumulatively annotated genes of 49 Sphingobium species. Phylogenetic analysis done using the core genome consisting of 802 orthologous clusters revealed a close relationship between S. fuliginis ATCC 27551 and bacteria capable of degradation of polyaromatic hydrocarbon compounds. Genes coding for transposases, efflux pumps conferring resistance to heavy metals, and TonR-type outer membrane receptors are selectively enriched in the genome of S. fuliginis ATCC 27551 and appear to contribute to the adaptive potential of the organism to challenging and harsh environments.

Comprehensive analyses of the annexin (ANN) gene family in Brassica rapa, Brassica oleracea and Brassica napus reveals their roles in stress response

Annexins (ANN) are a multigene, evolutionarily conserved family of calcium-dependent and phospholipid-binding proteins that play important roles in plant development and stress resistance. However, a systematic comprehensive analysis of ANN genes of Brassicaceae species (Brassica rapa, Brassica oleracea, and Brassica napus) has not yet been reported. In this study, we identified 13, 12, and 26 ANN genes in B. rapa, B. oleracea, and B. napus, respectively. About half of these genes were clustered on various chromosomes. Molecular evolutionary analysis showed that the ANN genes were highly conserved in Brassicaceae species. Transcriptome analysis showed that different group ANN members exhibited varied expression patterns in different tissues and under different (abiotic stress and hormones) treatments. Meanwhile, same group members from Arabidopsis thaliana, B. rapa, B. oleracea, and B. napus demonstrated conserved expression patterns in different tissues. The weighted gene coexpression network analysis (WGCNA) showed that BnaANN genes were induced by methyl jasmonate (MeJA) treatment and played important roles in jasmonate (JA) signaling and multiple stress response in B. napus.

The evolution of metabolism: How to test evolutionary hypotheses at the genomic level

The origin of primordial metabolism and its expansion to form the metabolic networks extant today represent excellent systems to study the impact of natural selection and the potential adaptive role of novel compounds. Here we present the current hypotheses made on the origin of life and ancestral metabolism and present the theories and mechanisms by which the large chemical diversity of plants might have emerged along evolution. In particular, we provide a survey of statistical methods that can be used to detect signatures of selection at the gene and population level, and discuss potential and limits of these methods for investigating patterns of molecular adaptation in plant metabolism.

Greater genetic and regulatory plasticity of retained duplicates in Epichloë endophytic fungi

Gene duplicates can act as a source of genetic material from which new functions arise. Most duplicated genes revert to single copy genes and only a small proportion are retained. However, it remains unclear why some duplicate genes persist in the genome for an extended time. We investigate this question by analysing retained gene duplicates in the fungal genus Epichloë, ascomycete fungi that form close endophytic symbioses with their host grasses. Retained duplicates within this genus have two independent origins, but both long pre-date the origin and diversification of the genus Epichloë. We find that loss of retained duplicates within the genus is frequent and often associated with speciation. Retained duplicates have faster evolutionary rates (Ka) and show relaxed selection (Ka/Ks) compared to single copy genes. Both features are time-dependent. Through comparison of conspecific strains, we find greater evolutionary rates in coding regions and sequence divergence in regulatory regions of retained duplicates than single copy genes, with this pattern more pronounced for strains adapted to different grass host species. Consistent with this sequence divergence in regulatory regions, transcriptome analyses show greater expression variation of retained duplicates than single copy genes. This suggest that cis-regulatory changes make important contributions to the expression patterns of retained duplicates. Coupled with supporting observations from the model yeast Saccharomyces cerevisiae, these data suggest that genetic robustness and regulatory plasticity are common drivers behind the retention of duplicated genes in fungi.

Lysine demethylases KDM6A and UTY: The X and Y of histone demethylation

Histone demethylases remove transcriptional repressive marks from histones in the nucleus. KDM6A (also known as UTX) is a lysine demethylase which acts on the trimethylated lysine at position 27 in histone 3. The KDM6A gene is located on the X chromosome but escapes X inactivation even though it is not located in the pseudoautosomal region. There is a homologue of KDM6A on the Y chromosome, known as UTY. UTY was thought to have lost its demethylase activity and to represent a non-functional remnant of the ancestral KDM6A gene. However, results with knockout mice suggest that the gene is expressed and the protein performs some function within the cell. Female mice with homozygous deletion of Kdm6a do not survive, but hemizygous males are viable, attributed to the presence of the Uty gene. KDM6A is mutated in the human condition Kabuki syndrome type 2 (OMIM 300867) and in many cases of cancer. The amino acid sequence of KDM6A has been conserved across animal phyla, although it is only found on the X chromosome in eutherian mammals. In this review, we reanalyse existing data from various sources (protein sequence comparison, evolutionary genetics, transcription factor binding and gene expression analysis) to determine the function, expression and evolution of KDM6A and UTY and show that UTY has a functional role similar to KDM6A in metabolism and development.

The Essentiality Status of Mouse Duplicate Gene Pairs Correlates with Developmental Co-Expression Patterns

During the evolution of multicellular eukaryotes, gene duplication occurs frequently to generate new genes and/or functions. A duplicated gene may have a similar function to its ancestral gene. Therefore, it may be expected that duplicated genes are less likely to be critical for the survival of an organism, since there are multiple copies of the gene rendering each individual copy redundant. In this study, we explored the developmental expression patterns of duplicate gene pairs and the relationship between development co-expression and phenotypes resulting from the knockout of duplicate genes in the mouse. We define genes that generate lethal phenotypes in single gene knockout experiments as essential genes. We found that duplicate gene pairs comprised of two essential genes tend to be expressed at different stages of development, compared to duplicate gene pairs with at least one non-essential member, showing that the timing of developmental expression affects the ability of one paralogue to compensate for the loss of the other. Gene essentiality, developmental expression and gene duplication are thus closely linked.

Evolutionary and functional analysis of the plant-specific NADPH oxidase gene family in Brassica rapa L

NADPH oxidases (NOXs) have been known as respiratory burst oxidase homologues (RBOHs) in plants. To characterize the evolutionary relationships and functions of RBOHs in Brassica rapa, 134 RBOH homologues were identified from 13 plant species, including 14 members (namely BrRBOH01-14) from B. rapa. There presented 47 gene-pairs among 14 BrRBOHs and other RBOHs, consisting of five pairs within B. rapa, and 15 pairs between B. rapa and Arabidopsis thaliana. Together with phylogenetic analysis, the results suggested that whole-genome duplication might have played an important role in BrRBOH gene expansion, and these duplication events occurred after the divergence of the eudicot and the monocot lineages examined. Furthermore, gene expression of RBOHs in both A. thaliana and B. rapa were assayed via qRT-PCR. An RBOH gene, BrRBOH13 in B. rapa, was transformed into wild-type Arabidopsis plants. The transgenic lines with the overexpressed level of BrRBOH13 conferred to be more tolerant to heavy metal lead (0.05 mM) than wild-type plants. Overall, this integrated analysis at genome-wide level has provided some information on the evolutionary relationships among plant-specific NOXs and the coordinated diversification of gene structure and function in B. rapa.

Effector variation at tandem gene arrays in tissue-dwelling coccidia: who needs antigenic variation anyway?

Locus expansion and diversification is pervasive in apicomplexan genomes and is predominantly found in loci encoding secreted proteins that interact with factors outside of the parasite. Key for understanding the impact of each of these loci on the host requires identification and functional characterization of their protein products, but these repetitive loci often are refractory to genome assembly. In this review we focus on Toxoplasma gondii and its nearest relatives to highlight the known impact of duplicated and diversified loci on our understanding of the host-pathogen molecular arms race. We describe current tools used for the identification and characterization of these loci, and review the most recent examples of how gene-expansion driven diversification can lead to novel gene functions.

Identification of jumonjiC domain containing gene family among the Oryza species and their expression analysis in FL478, a salt tolerant rice genotype

The jumonji (JMJ)-C domain containing proteins belong to histone demethylases family with the ability to demethylate the tri-methylated histone residues. They act as chromatin regulators to regulate many physiological functions in plants. The present study deals with the characterization of JMJ-C gene family members in wild as well as cultivated rice species and their expression analysis in salt tolerant rice genotype, FL478. The genome wide study identified 151 members belonging to JMJ-C gene family in 11 different Oryza species. We also studied their structure, genomic location, gene duplication events, phylogenetic relationship, in silico expression analysis and identified cis elements in their promoters. We also found a few JMJ-C gene family members in rice which underwent duplication before the whole genome duplication event of the rice. The qRT-PCR based expression profiling revealed that out of the total 15 rice JMJ-C members, two were highly expressed in the flag leaf stage of FL478 under salt treatment. These two candidate JMJ-C members were also found to render salinity tolerance when over-expressed in yeast cells. Thus, the present study helps in further structural as well as functional characterization of JMJ-C genes under salinity stress in Oryza species.

Characterization of the APETALA2/Ethylene-responsive factor (AP2/ERF) transcription factor family in sunflower

One of the most prominent families of genes in plants is the AP2/ERF which play an important role in regulating plant growth and responses to various stresses. In this research, a genome-wide survey was conducted to recognize the AP2/ERF genes in sunflower (Helianthus annuus L.), and a total of 288 HaAP2/ERF was obtained. Phylogenetic analysis divided them into four sub-families, including 248 ERF, 4 RAV and 35 AP2, and one subgroup of the Soloist family. Localization of chromosome, gene structure, the conserved motif, gene ontology, interaction networks, homology modeling, the modeling of cis-regulatory elements and the analysis of events in the duplication of genes were carried out for HaAP2/ERF genes. Finally, 9AP2/ERF genes were chosen to confirm the gene expression of the selected genes in leaf and root tissues in various abiotic stress conditions by qPCR. The results confirmed that AP2/ERFs genes could effectively resist abiotic stress. Also, proline content was studied under drought, salinity, cold and heat stress. The results indicated that proline was increased under abiotic stress. This research has been done for the first time to determine the HaAP2/ERF family, which prepared valuable data for the evolutionary and practical research regarding AP2/ERF in sunflower.

Functional diversification accompanies gene family expansion of MED2 homologs in Candida albicans

Gene duplication facilitates functional diversification and provides greater phenotypic flexibility to an organism. Expanded gene families arise through repeated gene duplication but the extent of functional divergence that accompanies each paralogous gene is generally unexplored because of the difficulty in isolating the effects of single family members. The telomere-associated (TLO) gene family is a remarkable example of gene family expansion, with 14 members in the more pathogenic Candida albicans relative to two TLO genes in the closely-related species C. dubliniensis. TLO genes encode interchangeable Med2 subunits of the major transcriptional regulatory complex Mediator. To identify biological functions associated with each C. albicans TLO, expression of individual family members was regulated using a Tet-ON system and the strains were assessed across a range of phenotypes involved in growth and virulence traits. All TLOs affected multiple phenotypes and a single phenotype was often affected by multiple TLOs, including simple phenotypes such as cell aggregation and complex phenotypes such as virulence in a Galleria mellonella model of infection. No phenotype was regulated by all TLOs, suggesting neofunctionalization or subfunctionalization of ancestral properties among different family members. Importantly, regulation of three phenotypes could be mapped to individual polymorphic sites among the TLO genes, including an indel correlated with two phenotypes, growth in sucrose and macrophage killing. Different selective pressures have operated on the TLO sequence, with the 5’ conserved Med2 domain experiencing purifying selection and the gene/clade-specific 3’ end undergoing extensive positive selection that may contribute to the impact of individual TLOs on phenotypic variability. Therefore, expansion of the TLO gene family has conferred unique regulatory properties to each paralog such that it influences a range of phenotypes. We posit that the genetic diversity associated with this expansion contributed to C. albicans success as a commensal and opportunistic pathogen.

Gene duplication is a rapid mechanism to generate additional sequences for natural selection to act upon and confer greater organismal fitness. If additional copies of the gene are beneficial, this process may be repeated to produce an expanded gene family containing many copies of related sequences. Following duplication, individual gene family members may retain functions of the ancestral gene or acquire new functions through mutation. How functional diversification accompanies expansion into large gene families remains largely unexplored due to the difficulty in assessing individual genes in the presence of the remaining family members. Here, we addressed this question using an inducible promoter to regulate expression of individual genes of the TLO gene family in the commensal yeast and opportunistic pathogen Candida albicans, which encode components of a major transcriptional regulator. Induced expression of individual TLOs affected a wide range of phenotypes such that significant functional overlap occurred among TLO genes and most phenotypes were affected by more than one TLO. Induced expression of individual TLOs did not produce massive phenotypic effects in most cases, suggesting that functional overlap among TLO genes may buffer new mutations that arise. Specific sequence variants among the TLO genes correlated with certain phenotypes and these sequence variants did not necessarily correlate with sequence similarity across the entire gene. Therefore, individual TLO family members evolved specific functional roles following duplication that likely reflect a combination of inherited function and new mutation.

Chromosomal assignment of centromere-specific histone CENH3 genes in rye (Secale cereale L.) and their phylogeny

Centromeres are essential for correct chromosome segregation during cell division and are determined by the presence of centromere-specific histone 3 (CENH3). Most of the diploid plant species, in which the structure and copy number of CENH3 genes have been determined, have this gene as a singleton; however, some cereal species in the tribe Triticeae have been found to have CENH3 in two variants. In this work, using the set of the wheat-rye addition lines we wanted to establish the chromosomal assignment of the CENH3 genes in the cultivated rye, Secale cereale (Linnaeus, 1753), in order to expand our knowledge about synteny conservation in the most important cereal species and about their chromosome evolution. To this end, we have also analyzed data in available genome sequencing databases. As a result, the αCENH3 and βCENH3 forms have been assigned to rye chromosomes 1R and 6R: specifically, the commonest variants αCENH3v1 and βCENH3v1 to chromosome 1R, and the rare variants, αCENH3v2 and probably βCENH3v2, to chromosome 6R. No other CENH3 variants have been found by analysis of the rye genome sequencing databases. Our chromosomal assignment of CENH3 in rye has been found to be the same as that in barley, suggesting that both main forms of CENH3 appeared in a Triticeae species before the barley and wheatrye lineages split.

Genome-wide analysis of DUF221 domain-containing gene family in Oryza species and identification of its salinity stress-responsive members in rice

DUF221 domain-containing genes (DDP genes) play important roles in developmental biology, hormone signalling transduction, and responses to abiotic stress. Therefore to understand their structural and evolutionary relationship, we did a genome-wide analysis of this important gene family in rice. Further, through comparative genomics, DDP genes from Oryza sativa subsp. (indica), nine different wild species of rice and Arabidopsis were also identified. We also found an expansion of the DDP gene families in rice and Arabidopsis which is due to the segmental duplication events in some of the gene family members. In general, a highly purifying selection was found acting on all the deduced paralogous and orthologous DDP gene pairs. The data from microarray and subsequent qRT-PCR analysis revealed that although several OsDDPs were differentially regulated under salinity stress, yet OsDDP6 was upregulated at all the developmental stages in salt tolerant rice genotype, FL478. Interestingly, OsDDP6 was found to be involved in proline metabolism pathway as indicated by protein network analysis. The diverse gene structures, varied transmembrane topologies and the differential expression patterns implied the functional diversity in DDP genes. Therefore, the comprehensive evolutionary analysis of DDP genes from different Oryza species and Arabidopsis performed in this study will provide the basis for further functional validation studies vis-à-vis DDP genes of rice and other plant species.

Genome-wide analysis of the CCCH zinc finger family identifies tissue specific and stress responsive candidates in chickpea (Cicer arietinum L.)

The CCCH zinc finger is a group of proteins characterised by a typical motif consisting of three cysteine residues and one histidine residue. These proteins have been reported to play important roles in regulation of plant growth, developmental processes and environmental responses. In the present study, genome wide analysis of the CCCH zinc finger gene family was carried out in the available chickpea genome. Various bioinformatics tools were employed to predict 58 CCCH zinc finger genes in chickpea (designated CarC3H1-58), which were analysed for their physio-chemical properties. Phylogenetic analysis classified the proteins into 12 groups in which members of a particular group had similar structural organization. Further, the numbers as well as the types of CCCH motifs present in the CarC3H proteins were compared with those from Arabidopsis and Medicago truncatula. Synteny analysis revealed valuable information regarding the evolution of this gene family. Tandem and segmental duplication events were identified and their Ka/Ks values revealed that the CarC3H gene family in chickpea had undergone purifying selection. Digital, as well as real time qRT-PCR expression analysis was performed which helped in identification of several CarC3H members that expressed preferentially in specific chickpea tissues as well as during abiotic stresses (desiccation, cold, salinity). Moreover, molecular characterization of an important member CarC3H45 was carried out. This study provides comprehensive genomic information about the important CCCH zinc finger gene family in chickpea. The identified tissue specific and abiotic stress specific CCCH genes could be potential candidates for further characterization to delineate their functional roles in development and stress.

Co-accumulation of cis-regulatory and coding mutations during the pseudogenization of the Xenopus laevis homoeologs six6.L and six6.S

Common models for the evolution of duplicated genes after genome duplication are subfunctionalization, neofunctionalization, and pseudogenization. Although the crucial roles of cis-regulatory mutations in subfunctionalization are well-documented, their involvement in pseudogenization and/or neofunctionalization remains unclear. We addressed this issue by investigating the evolution of duplicated homeobox genes, six6.L and six6.S, in the allotetraploid frog Xenopus laevis. Based on a comparative expression analysis, we observed similar eye-specific expression patterns for the two loci and their single ortholog in the ancestral-type diploid species Xenopus tropicalis. However, we detected lower levels of six6.S expression than six6.L expression. The six6.S enhancer sequence was more highly diverged from the orthologous enhancer of X. tropicalis than the six6.L enhancer, and showed weaker activity in a transgenic reporter assay. Based on a phylogenetic analysis of the protein sequences, we observed greater divergence between X. tropicalis Six6 and Six6.S than between X. tropicalis Six6 and Six6.L, and the observed mutations were reminiscent of a microphthalmia mutation in human SIX6. Misexpression experiments showed that six6.S has weaker eye-enlarging activity than six6.L, and targeted disruption of six6.L reduced the eye size more significantly than that of six6.S. These results suggest that enhancer attenuation stimulates the accumulation of hypomorphic coding mutations, or vice versa, in one duplicated gene copy and facilitates pseudogenization. We also underscore the value of the allotetraploid genome of X. laevis as a resource for studying latent pathogenic mutations.

Genome-wide analysis of Dof transcription factors reveals functional characteristics during development and response to biotic stresses in pepper

The DNA-binding with one zinc finger proteins (Dofs) are a plant-specific family of transcription factors. The Dofs are involved in a variety of biological processes such as phytohormone production, seed development, and environmental adaptation. Dofs have been previously identified in several plants, but not in pepper. We identified 33 putative Dof genes in pepper (CaDofs). To gain an overview of the CaDofs, we analyzed phylogenetic relationships, protein motifs, and evolutionary history. We divided the 33 CaDofs, containing 25 motifs, into four major groups distributed on eight chromosomes. We discovered an expansion of the CaDofs dated to a recent duplication event. Segmental duplication that occurred before the speciation of the Solanaceae lineages was predominant among the CaDofs. The global gene-expression profiling of the CaDofs by RNA-seq analysis showed distinct temporal and pathogen-specific variation during development and response to biotic stresses (two TMV strains, PepMoV, and Phytophthora capsici), suggesting functional diversity among the CaDofs. These results will provide the useful clues into the responses of Dofs in biotic stresses and promote a better understanding of their multiple function in pepper and other species.

Unraveling the message: insights into comparative genomics of the naked mole-rat

Animals have evolved to survive, and even thrive, in different environments. Genetic adaptations may have indirectly created phenotypes that also resulted in a longer lifespan. One example of this phenomenon is the preternaturally long-lived naked mole-rat. This strictly subterranean rodent tolerates hypoxia, hypercapnia, and soil-based toxins. Naked mole-rats also exhibit pronounced resistance to cancer and an attenuated decline of many physiological characteristics that often decline as mammals age. Elucidating mechanisms that give rise to their unique phenotypes will lead to better understanding of subterranean ecophysiology and biology of aging. Comparative genomics could be a useful tool in this regard. Since the publication of a naked mole-rat genome assembly in 2011, analyses of genomic and transcriptomic data have enabled a clearer understanding of mole-rat evolutionary history and suggested molecular pathways (e.g., NRF2-signaling activation and DNA damage repair mechanisms) that may explain the extraordinarily longevity and unique health traits of this species. However, careful scrutiny and re-analysis suggest that some identified features result from incorrect or imprecise annotation and assembly of the naked mole-rat genome: in addition, some of these conclusions (e.g., genes involved in cancer resistance and hairlessness) are rejected when the analysis includes additional, more closely related species. We describe how the combination of better study design, improved genomic sequencing techniques, and new bioinformatic and data analytical tools will improve comparative genomics and ultimately bridge the gap between traditional model and nonmodel organisms.

Multiple evolution of flavonoid 3',5'-hydroxylase

Multiple F3'5'H evolution from F3'H has occurred in dicotyledonous plants. Efficient pollinator attraction is probably the driving force behind, as this allowed for the synthesis of delphinidin-based blue anthocyanins. The cytochrome P450-dependent monooxygenases flavonoid 3'-hydroxylase (F3'H) and flavonoid 3',5'-hydroxylase (F3'5'H) hydroxylate the B-ring of flavonoids at the 3'- and 3'- and 5'-position, respectively. Their divergence took place early in plant evolution. While F3'H is ubiquitously present in higher plants, the distribution of F3'5'H is scattered. Here, we report that F3'5'H has repeatedly evolved from F3'H precursors at least four times in dicotyledonous plants: In the Asteraceae, we identified F3'5'Hs specific for the subfamilies Cichorioideae and Asteroideae, and additionally an F3'5'H that seems to be specific for the genus Echinops of the subfamily Carduoideae; moreover, characterisation of a sequence from Billardiera heterophylla (formerly Sollya heterophylla) (Pittosporaceae) showed that the independent evolution of an F3'5'H has occurred at least once also in another family. The evolution of F3'5'H from an F3'H precursor represents a gain of enzymatic function, probably triggered by an amino acid change at one position of substrate recognition site 6. The gain of F3'5'H activity allows for the synthesis of delphinidin-based anthocyanins which usually provide the basis for lilac to blue flower colours. Therefore, the need for an efficient pollinator attraction is probably the driving force behind the multiple F3'5'H evolution.

Genome-Wide Analysis of the AP2/ERF Transcription Factors Family and the Expression Patterns of DREB Genes in Moso Bamboo (Phyllostachys edulis)

The AP2/ERF transcription factor family, one of the largest families unique to plants, performs a significant role in terms of regulation of growth and development, and responses to biotic and abiotic stresses. Moso bamboo (Phyllostachys edulis) is a fast-growing non-timber forest species with the highest ecological, economic and social values of all bamboos in Asia. The draft genome of moso bamboo and the available genomes of other plants provide great opportunities to research global information on the AP2/ERF family in moso bamboo. In total, 116 AP2/ERF transcription factors were identified in moso bamboo. The phylogeny analyses indicated that the 116 AP2/ERF genes could be divided into three subfamilies: AP2, RAV and ERF; and the ERF subfamily genes were divided into 11 groups. The gene structures, exons/introns and conserved motifs of the PeAP2/ERF genes were analyzed. Analysis of the evolutionary patterns and divergence showed the PeAP2/ERF genes underwent a large-scale event around 15 million years ago (MYA) and the division time of AP2/ERF family genes between rice and moso bamboo was 15–23 MYA. We surveyed the putative promoter regions of the PeDREBs and showed that largely stress-related cis-elements existed in these genes. Further analysis of expression patterns of PeDREBs revealed that the most were strongly induced by drought, low-temperature and/or high salinity stresses in roots and, in contrast, most PeDREB genes had negative functions in leaves under the same respective stresses. In this study there were two main interesting points: there were fewer members of the PeDREB subfamily in moso bamboo than in other plants and there were differences in DREB gene expression profiles between leaves and roots triggered in response to abiotic stress. The information produced from this study may be valuable in overcoming challenges in cultivating moso bamboo.

Molecular evolution of NASP and conserved histone H3/H4 transport pathway

BACKGROUND

NASP is an essential protein in mammals that functions in histone transport pathways and maintenance of a soluble reservoir of histones H3/H4. NASP has been studied exclusively in Opisthokonta lineages where some functional diversity has been reported. In humans, growing evidence implicates NASP miss-regulation in the development of a variety of cancers. Although a comprehensive phylogenetic analysis is lacking, NASP-family proteins that possess four TPR motifs are thought to be widely distributed across eukaryotes.

RESULTS

We characterize the molecular evolution of NASP by systematically identifying putative NASP orthologs across diverse eukaryotic lineages ranging from excavata to those of the crown group. We detect extensive silent divergence at the nucleotide level suggesting the presence of strong purifying selection acting at the protein level. We also observe a selection bias for high frequencies of acidic residues which we hypothesize is a consequence of their critical function(s), further indicating the role of functional constraints operating on NASP evolution. Our data indicate that TPR1 and TPR4 constitute the most rapidly evolving functional units of NASP and may account for the functional diversity observed among well characterized family members. We also show that NASP paralogs in ray-finned fish have different genomic environments with clear differences in their GC content and have undergone significant changes at the protein level suggesting functional diversification.

CONCLUSION

We draw four main conclusions from this study. First, wide distribution of NASP throughout eukaryotes suggests that it was likely present in the last eukaryotic common ancestor (LECA) possibly as an important innovation in the transport of H3/H4. Second, strong purifying selection operating at the protein level has influenced the nucleotide composition of NASP genes. Further, we show that selection has acted to maintain a high frequency of functionally relevant acidic amino acids in the region that interrupts TPR2. Third, functional diversity reported among several well characterized NASP family members can be explained in terms of quickly evolving TPR1 and TPR4 motifs. Fourth, NASP fish specific paralogs have significantly diverged at the protein level with NASP2 acquiring a NNR domain.

Genome-wide identification, phylogeny, evolution and expression patterns of AP2/ERF genes and cytokinin response factors in Brassica rapa ssp. pekinensis

The AP2/ERF transcription factor family is one of the largest families involved in growth and development, hormone responses, and biotic or abiotic stress responses in plants. In this study, 281 AP2/ERF transcription factor unigenes were identified in Chinese cabbage. These superfamily members were classified into three families (AP2, ERF, and RAV). The ERF family was subdivided into the DREB subfamily and the ERF subfamily with 13 groups (I- XI) based on sequence similarity. Duplication, evolution and divergence of the AP2/ERF genes in B. rapa and Arabidopsis thaliana were investigated and estimated. Cytokinin response factors (CRFs), as a subclade of the AP2/ERF family, are important transcription factors that define a branch point in the cytokinin two-component signal (TCS) transduction pathway. Up to 21 CRFs with a conserved CRF domain were retrieved and designated as BrCRFs. The amino acid sequences, conserved regions and motifs, phylogenetic relationships, and promoter regions of the 21 BrCRFs were analyzed in detail. The BrCRFs broadly expressed in various tissues and organs. The transcripts of BrCRFs were regulated by factors such as drought, high salinity, and exogenous 6-BA, NAA, and ABA, suggesting their involvement in abiotic stress conditions and regulatory mechanisms of plant hormone homeostasis. These results provide new insight into the divergence, variation, and evolution of AP2/ERF genes at the genome-level in Chinese cabbage.

Following gene duplication, paralog interference constrains transcriptional circuit evolution

Most models of gene duplication assume that the ancestral functions of the preduplication gene are independent and can therefore be neatly partitioned between descendant paralogs. However, many gene products, such as transcriptional regulators, are components within cooperative assemblies; here, we show that a natural consequence of duplication and divergence of such proteins can be competitive interference between the paralogs. Our example is based on the duplication of the essential MADS-box transcriptional regulator Mcm1, which is found in all fungi and regulates a large set of genes. We show that a set of historical amino acid sequence substitutions minimized paralog interference in contemporary species and, in doing so, increased the molecular complexity of this gene regulatory network. We propose that paralog interference is a common constraint on gene duplicate evolution, and its resolution, which can generate additional regulatory complexity, is needed to stabilize duplicated genes in the genome.

Afrobatrachian mitochondrial genomes: genome reorganization, gene rearrangement mechanisms, and evolutionary trends of duplicated and rearranged genes

Background

Mitochondrial genomic (mitogenomic) reorganizations are rarely found in closely-related animals, yet drastic reorganizations have been found in the Ranoides frogs. The phylogenetic relationships of the three major ranoid taxa (Natatanura, Microhylidae, and Afrobatrachia) have been problematic, and mitogenomic information for afrobatrachians has not been available. Several molecular models for mitochondrial (mt) gene rearrangements have been proposed, but observational evidence has been insufficient to evaluate them. Furthermore, evolutionary trends in rearranged mt genes have not been well understood. To gain molecular and phylogenetic insights into these issues, we analyzed the mt genomes of four afrobatrachian species (Breviceps adspersus, Hemisus marmoratus, Hyperolius marmoratus, and Trichobatrachus robustus) and performed molecular phylogenetic analyses. Furthermore we searched for two evolutionary patterns expected in the rearranged mt genes of ranoids.

Results

Extensively reorganized mt genomes having many duplicated and rearranged genes were found in three of the four afrobatrachians analyzed. In fact, Breviceps has the largest known mt genome among vertebrates. Although the kinds of duplicated and rearranged genes differed among these species, a remarkable gene rearrangement pattern of non-tandemly copied genes situated within tandemly-copied regions was commonly found. Furthermore, the existence of concerted evolution was observed between non-neighboring copies of triplicated 12S and 16S ribosomal RNA regions.

Conclusions

Phylogenetic analyses based on mitogenomic data support a close relationship between Afrobatrachia and Microhylidae, with their estimated divergence 100 million years ago consistent with present-day endemism of afrobatrachians on the African continent. The afrobatrachian mt data supported the first tandem and second non-tandem duplication model for mt gene rearrangements and the recombination-based model for concerted evolution of duplicated mt regions. We also showed that specific nucleotide substitution and compositional patterns expected in duplicated and rearranged mt genes did not occur, suggesting no disadvantage in employing these genes for phylogenetic inference.

The influence of assortativity on the robustness and evolvability of gene regulatory networks upon gene birth

Gene regulatory networks (GRNs) represent the interactions between genes and gene products, which drive the gene expression patterns that produce cellular phenotypes. GRNs display a number of characteristics that are beneficial for the development and evolution of organisms. For example, they are often robust to genetic perturbation, such as mutations in regulatory regions or loss of gene function. Simultaneously, GRNs are often evolvable as these genetic perturbations are occasionally exploited to innovate novel regulatory programs. Several topological properties, such as degree distribution, are known to influence the robustness and evolvability of GRNs. Assortativity, which measures the propensity of nodes of similar connectivity to connect to one another, is a separate topological property that has recently been shown to influence the robustness of GRNs to point mutations in cis-regulatory regions. However, it remains to be seen how assortativity may influence the robustness and evolvability of GRNs to other forms of genetic perturbation, such as gene birth via duplication or de novo origination. Here, we employ a computational model of genetic regulation to investigate whether the assortativity of a GRN influences its robustness and evolvability upon gene birth. We find that the robustness of a GRN generally increases with increasing assortativity, while its evolvability generally decreases. However, the rate of change in robustness outpaces that of evolvability, resulting in an increased proportion of assortative GRNs that are simultaneously robust and evolvable. By providing a mechanistic explanation for these observations, this work extends our understanding of how the assortativity of a GRN influences its robustness and evolvability upon gene birth.

Genomic organization and transcription of the medaka and zebrafish cellular retinol-binding protein (rbp) genes

In this study, we examined the evolutionary trajectories and the common ancestor of medaka rbp genes by comparing them to the well-studied rbp/RBP genes from zebrafish and other vertebrates. We describe here gene structure, sequence identity, phylogenetic analysis and conserved gene synteny of medaka rbp genes and their putative proteins as well as the tissue-specific distribution of rbp transcripts in adult medaka and zebrafish. Medaka rbp genes consist of four exons separated by three introns that encode putative polypeptides of 134-138 amino acids, a genomic organization characteristic of rbp genes. Medaka Rbp sequences share highest sequence identity and similarity with their orthologs in vertebrates, and form a distinct clade with them in phylogenetic analysis. Conserved gene synteny was evident among medaka, zebrafish and human rbp/RBP genes, which provides compelling evidence that the medaka rbp1, rbp2a, rbp2b, rbp5, rbp7a and rbp7b genes arose from a common ancestor of vertebrates. Moreover, the duplicated rbp2 and rbp7 genes most likely exist owing to a whole-genome duplication (WGD) event specific to the teleost fish lineage. Selection pressure and the nonparametric relative rate test of the medaka and zebrafish duplicated rbp2 and rbp7 genes suggest that these duplicated genes are subjected to purifying selection and one paralog might have evolved at an accelerated rate compared to its sister duplicate since the WGD. The steady-state levels of medaka and zebrafish rbp1, rbp2a, rbp2b and rbp5 transcripts in various tissues suggest that medaka rbp1, rbp2a and rbp2b genes have retained the regulatory elements of an ancestral RBP1 and RBP2 genes, and the medaka rbp5 gene has acquired new function. Furthermore, the tissue-specific regulations of rbp7a and rbp7b genes have diverged markedly in medaka and zebrafish since the teleost-specific WGD.

Genetic interactions of Arabidopsis thaliana damaged DNA binding protein 1B (DDB1B) with DDB1A, DET1, and COP1

Damaged DNA Binding protein 1 (DDB1)–CULLIN4 E3 ubiquitin ligase complexes have been implicated in diverse biological processes in a range of organisms. Arabidopsis thaliana encodes two homologs of DDB1, DDB1A, and DDB1B. In this study we use a viable partial loss of function allele of DDB1B, ddb1b-2, to examine genetic interactions with DDB1A, DET1 and COP1. Although the ddb1b-2 ddb1a double mutant is lethal, ddb1a ddb1b-2/+ and ddb1b-2 ddb1a/+ heterozygotes exhibit few developmental phenotypes but do exhibit decreased tolerance of ultraviolet light. In addition, germination in ddb1a and ddb1a ddb1b-2/+ was found to be sensitive to salt and mannitol, and both DDB1 single mutants as well as the heterozygotes exhibited heat sensitivity. DE-ETIOLATED1 (DET1) and CONSTITUTIVE PHOTOMORPHOGENIC1 (COP1) are negative regulators of light development which interact with DDB1-CUL4 complexes. Although ddb1a strongly enhances det1 phenotypes in both dark- and light-grown seedlings, ddb1b-2 weakly enhanced the det1 short hypocotyl phenotype in the dark, as well as enhancing anthocyanin levels and suppressing the det1 low chlorophyll phenotype in light-grown seedlings. In adults, ddb1a suppresses det1 early flowering and enhances the det1 dwarf phenotype. A similar trend was observed in ddb1b-2 det1 double mutants, although the effects were smaller in magnitude. In cop1 mutants, ddb1b-2 enhanced the cop1-4 short hypocotyl phenotype in dark and light, enhanced anthocyanin levels in cop1-1 in the light, but had no effect in adults. Thus the requirement for DDB1B varies in the course of development, from COP1-specific effects in hypocotyls to DET1-specific in adults.

Genome-wide comparative analysis of annexin superfamily in plants

Most annexins are calcium-dependent, phospholipid-binding proteins with suggested functions in response to environmental stresses and signaling during plant growth and development. They have previously been identified and characterized in Arabidopsis and rice, and constitute a multigene family in plants. In this study, we performed a comparative analysis of annexin gene families in the sequenced genomes of Viridiplantae ranging from unicellular green algae to multicellular plants, and identified 149 genes. Phylogenetic studies of these deduced annexins classified them into nine different arbitrary groups. The occurrence and distribution of bona fide type II calcium binding sites within the four annexin domains were found to be different in each of these groups. Analysis of chromosomal distribution of annexin genes in rice, Arabidopsis and poplar revealed their localization on various chromosomes with some members also found on duplicated chromosomal segments leading to gene family expansion. Analysis of gene structure suggests sequential or differential loss of introns during the evolution of land plant annexin genes. Intron positions and phases are well conserved in annexin genes from representative genomes ranging from Physcomitrella to higher plants. The occurrence of alternative motifs such as K/R/HGD was found to be overlapping or at the mutated regions of the type II calcium binding sites indicating potential functional divergence in certain plant annexins. This study provides a basis for further functional analysis and characterization of annexin multigene families in the plant lineage.

The origin and evolution of new genes

New genes are a major source of genetic innovation in genomes. However, until recently, understanding how new genes originate and how they evolve was hampered by the lack of appropriate genetic datasets. The advent of the genomic era brought about a revolution in the amount of data available to study new genes. For the first time, decades-old theoretical principles could be tested empirically and novel and unexpected avenues of research opened up. This chapter explores how genomic data can and is being used to study both the origin and evolution of new genes and the surprising discoveries made thus far.

Genomic characterization of the European sea bass Dicentrarchus labrax reveals the presence of a novel uncoupling protein (UCP) gene family member in the teleost fish lineage

Background

Uncoupling proteins (UCP) are evolutionary conserved mitochondrial carriers that control energy metabolism and therefore play important roles in several physiological processes such as thermogenesis, regulation of reactive oxygen species (ROS), growth control, lipid metabolism and regulation of insulin secretion. Despite their importance in various physiological processes, their molecular function remains controversial. The evolution and phylogenetic distribution may assist to identify their general biological function and structure-function relationships. The exact number of uncoupling protein genes in the fish genome and their evolution is unresolved.

Results

Here we report the first characterisation of UCP gene family members in sea bass, Dicentrarchus labrax, and then retrace the evolution of the protein family in vertebrates. Four UCP genes that are shared by five other fish species were identified in sea bass genome. Phylogenetic reconstitution among vertebrate species and synteny analysis revealed that UCP1, UCP2 and UCP3 evolved from duplication events that occurred in the common ancestor of vertebrates, whereas the novel fourth UCP originated specifically in the teleost lineage. Functional divergence analysis among teleost species revealed specific amino acid positions that have been subjected to altered functional constraints after duplications.

Conclusions

This work provides the first unambiguous evidence for the presence of a fourth UCP gene in teleost fish genome and brings new insights into the evolutionary history of the gene family. Our results suggest functional divergence among paralogues which might result from long-term and differential selective pressures, and therefore, provide the indication that UCP genes may have diverse physiological functions in teleost fishes. Further experimental analysis of the critical amino acids identified here may provide valuable information on the physiological functions of UCP genes.

The effect of functional compensation among duplicate genes can constrain their evolutionary divergence

Gene duplicates have the inherent property of initially being functionally redundant. This means that they can compensate for the effect of deleterious variation occurring at one or more sister sites. Here, I present data bearing on evolutionary theory that illustrates the manner in which any functional adaptation in duplicate genes is markedly constrained because of the compensatory utility provided by a sustained genetic redundancy. Specifically, a two-locus epistatic model of paralogous genes was simulated to investigate the degree of purifying selection imposed, and whether this would serve to impede any possible biochemical innovation. Three population sizes were considered to see if, as expected, there was a significant difference in any selection for robustness. Interestingly, physical linkage between tandem duplicates was actually found to increase the probability of any neofunctionalization and the efficacy of selection, contrary to what is expected in the case of singleton genes. The results indicate that an evolutionary trade-off often exists between any functional change under either positive or relaxed selection and the need to compensate for failures due to degenerative mutations, thereby guaranteeing the reliability of protein production.