Search

Gene duplication is a major force driving evolutionary innovation. A classic example is generating new animal toxins via duplication of physiological protein-encoding genes and recruitment into venom. While this process drives the innovation of many animal venoms, reverse recruitment of toxins into nonvenomous cells remains unresolved. Using comparative genomics, we find members of the Membrane Attack Complex and Perforin Family (MAC) have been recruited into venom-injecting cells (cnidocytes), in soft and stony corals and sea anemones, suggesting that the ancestral MAC was a cnidocyte expressed toxin. Further investigation into the model sea anemone Nematostella vectensis reveals that three members have undergone Nematostella -specific duplications leading to their reverse recruitment into endomesodermal cells. Furthermore, simultaneous knockdown of all three endomesodermally expressed MACs leads to mis-development, supporting that these paralogs have nonvenomous function. By resolving the evolutionary history and function of MACs in Nematostella , we provide the first proof for reverse recruitment from venom to organismal development. [ABSTRACT FROM AUTHOR]

Horseshoe crabs represent the most ancestral chelicerate lineage characterized by marine ecology and the possession of lateral compound eyes. While considered living fossils, recent studies reported an unusual number of Pax6 genes in the Atlantic horseshoe crab Limulus polyphemus. Pax genes encode ancient metazoan transcription factors, which comprise seven subfamilies. Among these, the members of the Pax6 subfamily confer critical functions in the development of the head, the visual system, and further body plan components. Arthropods are generally characterized by two Pax6 subfamily homologs that were discovered in Drosophila and named eyeless (ey) and twin of eyeless (toy). However, whole genome sequence searches uncovered three homologs of ey and two homologs of toy in L. polyphemus. These numbers are explained by the occurrence of likely three whole genome duplications in the lineage to the last common ancestor of L. polyphemus and the three additional members of the extant horseshoe crab species cluster. Here, we report that all five L. polyphemus Pax6 paralogs are conserved in the approximately 135-million-year-old horseshoe crab species cluster and that they evolve under strong purifying selection. Largely homogenous protein sequence diversification rates of ey and toy paralogs suggest subfunctionalization as the likeliest preservation trajectory. However, our studies further revealed evidence that the horseshoe crab ey1 and ey2 paralogs share a derived splice isoform that encodes a unique five amino acid-long insertion in helix 3 of the homeodomain. This suggests that the exceptional expansion of the horseshoe crab Pax6 gene family repertoire was also associated with regulatory diversification and possibly innovation.

Horseshoe crabs represent the most ancestral chelicerate lineage characterized by marine ecology and the possession of lateral compound eyes. While considered living fossils, recent studies reported an unusual number of Pax6 genes in the Atlantic horseshoe crab Limulus polyphemus. Pax genes encode ancient metazoan transcription factors, which comprise seven subfamilies. Among these, the members of the Pax6 subfamily confer critical functions in the development of the head, the visual system, and further body plan components. Arthropods are generally characterized by two Pax6 subfamily homologs that were discovered in Drosophila and named eyeless (ey) and twin of eyeless (toy). However, whole genome sequence searches uncovered three homologs of ey and two homologs of toy in L. polyphemus. These numbers are explained by the occurrence of likely three whole genome duplications in the lineage to the last common ancestor of L. polyphemus and the three additional members of the extant horseshoe crab species cluster. Here, we report that all five L. polyphemus Pax6 paralogs are conserved in the approximately 135-million-year-old horseshoe crab species cluster and that they evolve under strong purifying selection. Largely homogenous protein sequence diversification rates of ey and toy paralogs suggest subfunctionalization as the likeliest preservation trajectory. However, our studies further revealed evidence that the horseshoe crab ey1 and ey2 paralogs share a derived splice isoform that encodes a unique five amino acid-long insertion in helix 3 of the homeodomain. This suggests that the exceptional expansion of the horseshoe crab Pax6 gene family repertoire was also associated with regulatory diversification and possibly innovation. [ABSTRACT FROM AUTHOR]

Aquaporin (Aqp) 10 is a member of the aquaglyceroporin subfamily of water channels, and human Aqp10 is permeable to solutes such as glycerol, urea, and boric acid. Tetrapods have a single aqp10 gene, whereas ray-finned fishes have paralogs of this gene through tandem duplication, whole-genome duplication, and subsequent deletion. A previous study on Aqps in the Japanese pufferfish Takifugu rubripes showed that one pufferfish paralog, Aqp10.2b, was permeable to water and glycerol, but not to urea and boric acid. To understand the functional differences of Aqp10s between humans and pufferfish from an evolutionary perspective, we analyzed Aqp10s from an amphibian (Xenopus laevis) and a lobe-finned fish (Protopterus annectens) and Aqp10.1 and Aqp10.2 from several ray-finned fishes (Polypterus senegalus , Lepisosteus oculatus , Danio rerio , and Clupea pallasii). The expression of tetrapod and lobe-finned fish Aqp10s and Aqp10.1-derived Aqps in ray-finned fishes in Xenopus oocytes increased the membrane permeabilities to water, glycerol, urea, and boric acid. In contrast, Aqp10.2-derived Aqps in ray-finned fishes increased water and glycerol permeabilities, whereas those of urea and boric acid were much weaker than those of Aqp10.1-derived Aqps. These results indicate that water, glycerol, urea, and boric acid permeabilities are plesiomorphic activities of Aqp10s and that the ray-finned fish-specific Aqp10.2 paralogs have secondarily reduced or lost urea and boric acid permeability. [ABSTRACT FROM AUTHOR]

Much evidence exists suggesting the presence of genetic functional diversification in plants, though literature associated with the role of functional diversification in the evolution of the plant secondary cell wall (SCW) has sparsely been compiled and reviewed in a recent context. This review aims to elucidate, through the examination of gene phylogenies associated with its biosynthesis and maintenance, the role of functional diversification in shaping the critical, dynamic, and characteristic organelle, the secondary cell wall. It will be asserted that gene families resulting from gene duplication and subsequent functional divergence are present and are heavily involved in SCW biosynthesis and maintenance. Furthermore, diversification will be presented as a significant driver behind the evolution of the many functional characteristics of the SCW. The structure and function of the plant cell wall and its constituents will first be explored, followed by a discussion on the phenomenon of gene duplication and the resulting genetic functional divergence that can emerge. Finally, the major constituents of the SCW and their individual relationships with duplication and divergence will be reviewed to the extent of current knowledge on the subject. [ABSTRACT FROM AUTHOR]

Hypoxia-inducible factor 1 (HIF-1) is a crucial transcriptional factor that can restore oxygen balance in the body by regulating multiple vital activities. Two HIF-1α copies were retained in cyprinid fish after experiencing a teleost-specific genome duplication. How the "divergent collaboration" of HIF-1αA and HIF-1αB proceeds in regulating mitophagy and apoptosis under hypoxic stress in cells of cyprinid fish remains unclear. In this study, zebrafish HIF-1αA/B expression plasmids were constructed and transfected into the epithelioma papulosum cyprini cells and were subjected to hypoxic stress. HIF-1αA induced apoptosis through promoting ROS generation and mitochondrial depolarization when cells were subjected to oxygen deficiency. Conversely, HIF-1αB was primarily responsible for mitophagy induction, prompting ATP production to mitigate apoptosis. HIF-1αA did not induce mitophagy in the mitochondria and lysosomes co-localization assay but it was involved in the regulation of different mitophagy pathways. Over-expression of HIF-1αA increased the expression of bnip3, fundc1, Beclin1, and foxo3, suggesting it has a dual role in mitochondrial autophagy and cell death. Each duplicated copy also experienced functional divergence and target shifting in the regulation of complexes in the mitochondrial electron transport chain (ETC). Our findings shed light on the post-subfunctionalization function of HIF-1αA and HIF-1αB in zebrafish to fine-tune regulation of mitophagy and apoptosis following hypoxia exposure. [ABSTRACT FROM AUTHOR]

The availability of an ever-increasing diversity of prokaryotic genomes and metagenomes represents a major opportunity to understand and decipher the mechanisms behind the functional diversification of microbial biosynthetic pathways. However, it remains unclear to what extent a pathway producing a specific molecule from a specific precursor can diversify. In this study, we focus on the biosynthesis of ubiquinone (UQ), a crucial coenzyme that is central to the bioenergetics and to the functioning of a wide variety of enzymes in Eukarya and Pseudomonadota (a subgroup of the formerly named Proteobacteria). UQ biosynthesis involves three hydroxylation reactions on contiguous carbon atoms. We and others have previously shown that these reactions are catalyzed by different sets of UQ-hydroxylases that belong either to the iron-dependent Coq7 family or to the more widespread flavin monooxygenase (FMO) family. Here, we combine an experimental approach with comparative genomics and phylogenetics to reveal how UQ-hydroxylases evolved different selectivities within the constrained framework of the UQ pathway. It is shown that the UQ-FMOs diversified via at least three duplication events associated with two cases of neofunctionalization and one case of subfunctionalization, leading to six subfamilies with distinct hydroxylation selectivity. We also demonstrate multiple transfers of the UbiM enzyme and the convergent evolution of UQ-FMOs toward the same function, which resulted in two independent losses of the Coq7 ancestral enzyme. Diversification of this crucial biosynthetic pathway has therefore occurred via a combination of parallel evolution, gene duplications, transfers, and losses. [ABSTRACT FROM AUTHOR]

The diversity of herbivorous insects is attributed to their propensity to specialize on toxic plants. In an evolutionary twist, toxins betray the identity of their bearers when herbivores coopt them as cues for host-plant finding, but the evolutionary mechanisms underlying this phenomenon are poorly understood. We focused on Scaptomyza flava, an herbivorous drosophilid specialized on isothiocyanate (ITC)-producing (Brassicales) plants, and identified Or67b paralogs that were triplicated as mustard-specific herbivory evolved. Using in vivo heterologous systems for the expression of olfactory receptors, we found that S. flava Or67bs, but not the homologs from microbe-feeding relatives, responded selectively to ITCs, each paralog detecting different ITC subsets. Consistent with this, S. flava was attracted to ITCs, as was Drosophila melanogaster expressing S. flava Or67b3 in the homologous Or67b olfactory circuit. ITCs were likely coopted as olfactory attractants through gene duplication and functional specialization (neofunctionalization and subfunctionalization) in S. flava, a recently derived herbivore.

Abstract Background Gene duplication events are critical for the evolution of new gene functions. Aristaless is a major regulator of distinct developmental processes. It is most known for its role during appendage development across animals. However, more recently other distinct biological functions have been described for this gene and its duplicates. Butterflies and moths have two copies of aristaless, aristaless1 (al1) and aristaless2 (al2), as a result of a gene duplication event. Previous work in Heliconius has shown that both copies appear to have novel functions related to wing color patterning. Here we expand our knowledge of the expression profiles associated with both ancestral and novel functions of Al1 across embryogenesis and wing pigmentation. Furthermore, we characterize Al2 expression, providing a comparative framework between gene copies within the same species, allowing us to understand the origin of new functions following gene duplication. Results Our work shows that the expression of both Al1 and Al2 is associated with the ancestral function of sensory appendage (leg, mouth, spines, and eyes) development in embryos. Interestingly, Al1 exhibits higher expression earlier in embryogenesis while the highest levels of Al2 expression are shifted to later stages of embryonic development. Furthermore, Al1 localization appears extranuclear while Al2 co-localizes tightly with nuclei earlier, and then also expands outside the nucleus later in development. Cellular expression of Al1 and Al2 in pupal wings is broadly consistent with patterns observed during embryogenesis. We also describe, for the first time, how Al1 localization appears to correlate with zones of anterior/posterior elongation of the body during embryonic growth, showcasing a possible new function related to Aristaless’ previously described role in appendage extension. Conclusions Overall, our data suggest that while both gene copies play a role in embryogenesis and wing pigmentation, the duplicates have diverged temporally and mechanistically across those functions. Our study helps clarify principles behind sub-functionalization and gene expression evolution associated with developmental functions following gene duplication events.

Abstract 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.

Background: Gene duplication events are critical for the evolution of new gene functions. Aristaless is a major regulator of distinct developmental processes. It is most known for its role during appendage development across animals. However, more recently other distinct biological functions have been described for this gene and its duplicates. Butterflies and moths have two copies of aristaless, aristaless1 (al1) and aristaless2 (al2), as a result of a gene duplication event. Previous work in Heliconius has shown that both copies appear to have novel functions related to wing color patterning. Here we expand our knowledge of the expression profiles associated with both ancestral and novel functions of Al1 across embryogenesis and wing pigmentation. Furthermore, we characterize Al2 expression, providing a comparative framework between gene copies within the same species, allowing us to understand the origin of new functions following gene duplication. Results: Our work shows that the expression of both Al1 and Al2 is associated with the ancestral function of sensory appendage (leg, mouth, spines, and eyes) development in embryos. Interestingly, Al1 exhibits higher expression earlier in embryogenesis while the highest levels of Al2 expression are shifted to later stages of embryonic development. Furthermore, Al1 localization appears extranuclear while Al2 co-localizes tightly with nuclei earlier, and then also expands outside the nucleus later in development. Cellular expression of Al1 and Al2 in pupal wings is broadly consistent with patterns observed during embryogenesis. We also describe, for the first time, how Al1 localization appears to correlate with zones of anterior/posterior elongation of the body during embryonic growth, showcasing a possible new function related to Aristaless' previously described role in appendage extension. Conclusions: Overall, our data suggest that while both gene copies play a role in embryogenesis and wing pigmentation, the duplicates have diverged temporally and mechanistically across those functions. Our study helps clarify principles behind sub-functionalization and gene expression evolution associated with developmental functions following gene duplication events. [ABSTRACT FROM AUTHOR]

The molecular characteristics and tissue disruption of 10 fatty acid–binding protein (fabp) genes in gilthead seabream (Sparus aurata) were investigated, and their expression levels were found in the fish fed diets with different vegetable oil (VO) sources, which may explore the potential function of fabp genes in S. aurata. For this purpose, the open reading frames of fabp genes involved in the transport and ß‐oxidation of fatty acids (FA) were molecularly cloned and characterized. S. aurata was then exposed to a two‐staged feeding trial (the grow‐out period following a wash‐out period) at low water temperatures. In the grow‐out period, the fish were fed diets containing 50% and 100% ratios of various VOs for 60 days, and in the wash‐out period, the fish were fed a diet containing 100% fish oil (FO) for 30 days. It has been determined that (a) S. aurata and vertebrate fabp/FABP genes are orthologues; (b) spatio‐temporal differences in tissue‐specific patterns of fabp genes differ importantly; for instance, the difference between the highest and lowest values reaches 13 × 105‐fold in the fabp10a; and (c) VO‐based diets upregulated fabp transcript levels in the liver and muscle with some exceptions, such as liver fabp11a and muscle fabp7a. Gene expressions of only the hepatic fabp7b and fabp10a genes were diminished at the end of the wash‐out period. In this study, the authors provide further evidence that dietary FAs affect fabp mRNA expressions in S. aurata. This might be useful in the nutritional control of fabp genes to maintain lipid homeostasis in marine fish fed VO‐based diets at low water temperatures. [ABSTRACT FROM AUTHOR]

The control of infections by the vertebrate adaptive immune system requires careful modulation to optimize defense and minimize harm to the host. The Fc receptor-like (FCRL) genes encode immunoregulatory molecules homologous to the receptors for the Fc portion of immunoglobulin (FCR). To date, nine different genes (FCRL1–6, FCRLA, FCRLB and FCRLS) have been identified in mammalian organisms. FCRL6 is located at a separate chromosomal position from the FCRL1-5 locus, has conserved synteny in mammals and is situated between the SLAMF8 and DUSP23 genes. Here, we show that this three gene block underwent repeated duplication in Dasypus novemcinctus (nine-banded armadillo) resulting in six FCRL6 copies, of which five appear functional. Among 21 mammalian genomes analyzed, this expansion was unique to D. novemcinctus. Ig-like domains that derive from the five clustered FCRL6 functional gene copies show high structural conservation and sequence identity. However, the presence of multiple non-synonymous amino acid changes that would diversify individual receptor function has led to the hypothesis that FCRL6 endured subfunctionalization during evolution in D. novemcinctus. Interestingly, D. novemcinctus is noteworthy for its natural resistance to the Mycobacterium leprae pathogen that causes leprosy. Because FCRL6 is chiefly expressed by cytotoxic T and NK cells, which are important in cellular defense responses against M. leprae, we speculate that FCRL6 subfunctionalization could be relevant for the adaptation of D. novemcinctus to leprosy. These findings highlight the species-specific diversification of FCRL family members and the genetic complexity underlying evolving multigene families critical for modulating adaptive immune protection. [ABSTRACT FROM AUTHOR]

Keratan sulfate (KS) is a glycosaminoglycan that is enriched in vertebrate cornea, cartilage, and brain. During embryonic development, highly sulfated KS (HSKS) is first detected in the developing notochord and then in otic vesicles; therefore, HSKS has been used as a molecular marker of the notochord. However, its biosynthetic pathways and functional roles in organogenesis are little known. Here, I surveyed developmental expression patterns of genes related to HSKS biosynthesis in Xenopus embryos. Of these genes, the KS chain-synthesizing glycosyltransferase genes, beta-1,3-N-acetylglucosaminyltransferase (b3gnt7) and beta-1,4-galactosyltransferase (b4galt4), are strongly expressed in the notochord and otic vesicles, but also in other tissues. In addition, their notochord expression is gradually restricted to the posterior end at the tailbud stage. In contrast, carbohydrate sulfotransferase (Chst) genes, chst2, chst3, and chst5.1, are expressed in both notochord and otic vesicles, whereas chst1, chst4/5-like, and chst7 are confined to otic vesicles. Because the substrate for Chst1 and Chst3 is galactose, while that for others is N-acetylglucosamine, combinatorial, tissue-specific expression patterns of Chst genes should be responsible for tissue-specific HSKS enrichment in embryos. As expected, loss of function of chst1 led to loss of HSKS in otic vesicles and reduction of their size. Loss of chst3 and chst5.1 resulted in HSKS loss in the notochord. These results reveal that Chst genes are critical for HSKS biosynthesis during organogenesis. Being hygroscopic, HSKS forms “water bags” in embryos to physically maintain organ structures. In terms of evolution, in ascidian embryos, b4galt and chst-like genes are also expressed in the notochord and regulate notochord morphogenesis. Furthermore, I found that a chst-like gene is also strongly expressed in the notochord of amphioxus embryos. These conserved expression patterns of Chst genes in the notochord of chordate embryos suggest that Chst is an ancestral component of the chordate notochord.

Gene duplication is a fundamental evolutionary process which provides opportunities to acquire new gene functions. In the case of the insulin receptors (InRs) in cockroaches and close-related insects, two successive duplications determined the occurrence of three InR genes: InR2 , InR1 and InR3 , the last two forming a sister cluster to InR2. The biological role of each of the gene duplicates and whether they resulted from neofunctionalization or subfunctionalization is still unclear. The analysis of the sequences from different lineages did not detect positive selection as driving the divergence of InR1 and InR3 , discarding neofunctionalization, and suggesting that there is no functional divergence between both gene copies. Using the cockroach Blattella germanica as a model, we have determined that BgInR2 is the gene with the highest expression levels in all the tissues analyzed, both in adult females and males, as well as in nymphs and embryos. BgInR3 is second in expression levels while BgInR1 is expressed at lower levels and only in some tissues. The selective depletion by RNAi of each of the three InRs, analyzed in terms of phenotype and fat body transcriptomic profiles, resulted in essentially redundant effects, with a magnitude approximately proportional to the level of expression of the respective InR. Therefore, the results indicate that the InR duplicates likely experienced a subfunctionalization process, by which the three InRs maintained similar functions but contributing to those functions proportionally to their expression levels. [Display omitted] • InR2 is the most expressed InR in all the analyzed tissues. • No positive selection as driver of the divergence between InR1 and InR3. • Specific depletion produced effects proportional to the respective expression levels. • Subfunctionalization: collaborative function proportional to the expression levels. [ABSTRACT FROM AUTHOR]

Gene duplication provides raw material for functional innovation, but gene duplicability varies considerably. Previous studies have found widespread asymmetrical sequence evolution between paralogs. However, it remains unknown whether the rate of evolution among paralogs affects their propensity of being retained after another round of whole‐genome duplication (WGD). In this study, we investigated gene groups that have experienced two successive WGDs to determine which of two older duplicates with different evolutionary rates was more likely to retain both younger duplicates. To uncouple the measurement of evolutionary rates from any assignment of duplicate or singleton status, we measured the evolutionary rates of singleton genes in out‐lineages but classified these singleton genes according to whether they are retained or not in a crown group of species. We found that genes that retained younger duplicates in the crown group of genomes were more constrained prior to the younger duplication event than those that failed to leave duplicates. In addition, we also found that the retained clades have more genes in out‐lineages. Subsequent analyses showed that genes in the retained clades were expressed more broadly and highly than genes in the singleton clades. We concluded that the set of repeatedly retained genes after two WGDs is biased toward slowly evolving genes in angiosperms, suggesting that the potential of genes for both functional conservation and divergence likely affects their propensity of being retained after WGD in angiosperms. [ABSTRACT FROM AUTHOR]

The duplication and divergence of stress-related genes might help microbes adapt to different environments. However, little is known about this type of subfunctionalization process in microbes. We here examined the evolution and function of PHB polymerases (PhaCs) from the plant growth-promoting, endophytic, and diazotrophic bacterium Azospirillum brasilense Az39. Besides the previously reported chromosomal phaC gene (phaCAz39-Chr) involved in abiotic stress adaptation, we showed that strain Az39 contains an additional copy of the phaC gene (phaCAz39-p4) in an extra-chromosomal replicon. Heterologous expression analysis demonstrated that the phaCAz39-p4 gene codifies for a functional PHB polymerase. Comparative sequence and evolutionary analyses showed an extreme level of divergence between PhaCAz39-p4 and other known PhaCs, including PhaCAz39-Chr, suggesting a divergent evolution of paralogous copies of phaC from strain Az39. Contrary to its paralogous copy, phaCAz39-p4 is specifically induced in planta and is necessary to enhance plant growth, thus empirically supporting the predicted functional divergence. [ABSTRACT FROM AUTHOR]

Neofunctionalization of duplicated gene copies is thought to be an important process underlying the origin of evolutionary novelty and provides an elegant mechanism for the origin of new phenotypic traits. One putative case where a new gene copy has been linked to a novel morphological trait is the origin of the arachnid patella, a taxonomically restricted leg segment. In spiders, the origin of this segment has been linked to the origin of the paralog dachshund-2, suggesting that a new gene facilitated the expression of a new trait. However, various arachnid groups that possess patellae do not have a copy of dachshund-2, disfavoring the direct link between gene origin and trait origin. We investigated the developmental genetic basis for patellar patterning in the harvestman Phalangium opilio, which lacks dachshund-2. Here, we show that the harvestman patella is established by a novel expression domain of the transcription factor extradenticle. Leveraging this definition of patellar identity, we surveyed targeted groups across chelicerate phylogeny to assess when this trait evolved. We show that a patellar homolog is present in Pycnogonida (sea spiders) and various arachnid orders, suggesting a single origin of the patella in the ancestor of Chelicerata. A potential loss of the patella is observed in Ixodida. Our results suggest that the modification of an ancient gene, rather than the neofunctionalization of a new gene copy, underlies the origin of the patella. Broadly, this work underscores the value of comparative data and broad taxonomic sampling when testing hypotheses in evolutionary developmental biology., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence-function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus , we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning approximately 400 million years of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty., Competing Interests: COMPETING INTERESTS The Wisconsin Alumni Research Foundation has filed a patent application on the technologies described herein with J.G.C. and C.T.H. as inventors. Strains are available for non-commercial, academic use under a material transfer agreement. A.R. is a scientific consultant for LifeMine Therapeutics, Inc.

The diversified natriuretic peptide (NP) family, consisting of four CNPs (CNP1-4), ANP, BNP, and VNP, has been identified in the eel. Here, we successfully cloned additional cnp genes from the brain of eel (a basal teleost) and zebrafish (a later branching teleost). The genes were identified as paralogues of cnp4 generated by the third round of whole genome duplication (3R) in the teleost lineage, thereby being named eel cnp4b and zebrafish cnp4-like, respectively. To examine the histological patterns of their expressions, we employed a newly developed in situ hybridization (ISH) chain reaction using short hairpin DNAs, in addition to conventional ISH. Eel cnp4b was expressed in the medulla oblongata, while mRNAs of eel cnp4a (former cnp4) were localized in the preoptic area. In the zebrafish brain, cnp4-like mRNA was undetectable, while the known cnp4 was expressed in both the preoptic area and medulla oblongata. Together with the different mRNA distribution of cnp4a and cnp4b in eel peripheral tissues determined by RT-PCR and ISH, it is suggested that subfunctionalization by duplicated cnp4s in ancestral teleosts has been retained only in basal teleosts. Intriguingly, cnp4b-expressing neurons in the glossopharyngeal-vagal motor complex of the medulla oblongata were co-localized with choline acetyltransferase, suggesting an involvement of Cnp4b in swallowing and respiration functions that are modulated by the vagus. Since teleost Cnp4 is an ortholog of mammalian CNP, the identified localization of teleost Cnp4 will contribute to future studies aimed at deciphering the physiological functions of CNP. [ABSTRACT FROM AUTHOR]

Chromatin remodelers play a fundamental role in the assembly of chromatin, regulation of transcription, and DNA repair. Biochemical and functional characterizations of the CHD family of chromatin remodelers from a variety of model organisms have shown that these remodelers participate in a wide range of activities. However, because the evolutionary history of CHD homologs is unclear, it is difficult to predict which of these activities are broadly conserved and which have evolved more recently in individual eukaryotic lineages. Here, we performed a comprehensive phylogenetic analysis of 8,042 CHD homologs from 1,894 species to create a model for the evolution of this family across eukaryotes with a particular focus on the timing of duplications that gave rise to the diverse copies observed in plants, animals, and fungi. Our analysis confirms that the three major subfamilies of CHD remodelers originated in the eukaryotic last common ancestor, and subsequent losses occurred independently in different lineages. Improved taxon sampling identified several subfamilies of CHD remodelers in plants that were absent or highly divergent in the model plant Arabidopsis thaliana. Whereas the timing of CHD subfamily expansions in vertebrates corresponds to whole genome duplication events, the mechanisms underlying CHD diversification in land plants appear more complicated. Analysis of protein domains reveals that CHD remodeler diversification has been accompanied by distinct transitions in domain architecture, contributing to the functional differences observed between these remodelers. This study demonstrates the importance of proper taxon sampling when studying ancient evolutionary events to prevent misinterpretation of subsequent lineage-specific changes and provides an evolutionary framework for functional and comparative analysis of this critical chromatin remodeler family across eukaryotes. [ABSTRACT FROM AUTHOR]

The control of infections by the vertebrate adaptive immune system requires careful modulation to optimize defense and minimize harm to the host. The Fc receptor-like (FCRL) genes encode immunoregulatory molecules homologous to the receptors for the Fc portion of immunoglobulin (FCR). To date, nine different genes (FCRL1–6, FCRLA, FCRLB and FCRLS) have been identified in mammalian organisms. FCRL6 is located at a separate chromosomal position from the FCRL1-5 locus, has conserved synteny in mammals and is situated between the SLAMF8 and DUSP23 genes. Here, we show that this three gene block underwent repeated duplication in Dasypus novemcinctus (nine-banded armadillo) resulting in six FCRL6 copies, of which five appear functional. Among 21 mammalian genomes analyzed, this expansion was unique to D. novemcinctus. Ig-like domains that derive from the five clustered FCRL6 functional gene copies show high structural conservation and sequence identity. However, the presence of multiple non-synonymous amino acid changes that would diversify individual receptor function has led to the hypothesis that FCRL6 endured subfunctionalization during evolution in D. novemcinctus. Interestingly, D. novemcinctus is noteworthy for its natural resistance to the Mycobacterium leprae pathogen that causes leprosy. Because FCRL6 is chiefly expressed by cytotoxic T and NK cells, which are important in cellular defense responses against M. leprae, we speculate that FCRL6 subfunctionalization could be relevant for the adaptation of D. novemcinctus to leprosy. These findings highlight the species-specific diversification of FCRL family members and the genetic complexity underlying evolving multigene families critical for modulating adaptive immune protection.

From a micro to macro scale of biological organization, macromolecular diversity and biological heterogeneity are fundamental properties of biological systems. Heterogeneity may result from genetic, epigenetic, and non-genetic characteristics (e.g., tissue microenvironment). Macromolecular diversity and biological heterogeneity are tolerated as long as the sustenance and propagation of life are not disrupted. They also provide the raw materials for microevolutionary changes that may help organisms adapt to new selection pressures arising from the environment. Sequence evolution, functional divergence, and positive selection of gene and promoter dosage play a major role in the evolution of life's diversity including complex metabolic networks, which is ultimately reflected in changes in the allele frequency over time. Robustness in evolvable biological systems is conferred by functional redundancy that is often created by macromolecular diversity and biological heterogeneity. The ability to investigate biological macromolecules at an increasingly finer level has uncovered a wealth of information in this regard. Therefore, the dynamics of biological complexity should be taken into consideration in biomedical research. [ABSTRACT FROM AUTHOR]

Chemokines, relatively small secreted proteins, are involved in cell migration and function in various biological events, including immunity, morphogenesis, and disease. Due to their nature, chemokines tend to be a target of hijacking of immunity by virus and therefore show an exceptionally high mutation rate. Xenopus laevis is considered an excellent model to investigate the effect of whole-genome duplication for gene family evolution. Because its allotetraploidization occurred around 17–18 million years ago, ancestral subgenomes L and S were well conserved. Based on the gene model of human and diploid frog Xenopus tropicalis , we identified 52 chemokine genes and 26 chemokine receptors in X. laevis. The retention rate of the gene in the X. laevis L and S subgenomes was 96% (45/47) and 68% (32/47), respectively. We conducted molecular phylogenetic analysis and found clear orthologies in all receptor genes but not in the ligand genes, suggesting rapid divergences of the ligand. dN/dS calculation demonstrated that dN/dS ratio greater than one was observed in the four ligand genes, cxcl8b.1.S , cxcl18.S , ccl21.S , and xcl1.L , but nothing in receptor genes. These results revealed that the whole-genome duplication promotes diversification of chemokine ligands in X. laevis while conserving the genes necessary for homeostasis, suggesting that selective pressure also supports a rapid divergence of the chemokines in amphibians. [ABSTRACT FROM AUTHOR]

Both pleiotropic connectivity and mutational correlations can restrict the decoupling of traits under divergent selection, but it is unknown which is more important in trait evolution. To address this question, we create a model that permits within-population variation in both pleiotropic connectivity and mutational correlation, and compare their relative importance to trait evolution. Specifically, we developed an individual-based stochastic model where mutations can affect whether a locus affects a trait and the extent of mutational correlations in a population. We find that traits can decouple whether there is evolution in pleiotropic connectivity or mutational correlation, but when both can evolve, then evolution in pleiotropic connectivity is more likely to allow for decoupling to occur. The most common genotype found in this case is characterized by having one locus that maintains connectivity to all traits and another that loses connectivity to the traits under stabilizing selection (subfunctionalization). This genotype is favored because it allows the subfunctionalized locus to accumulate greater effect size alleles, contributing to increasingly divergent trait values in the traits under divergent selection without changing the trait values of the other traits (genetic modularization). These results provide evidence that partial subfunctionalization of pleiotropic loci may be a common mechanism of trait decoupling under regimes of corridor selection. [ABSTRACT FROM AUTHOR]

Chemokines, relatively small secreted proteins, are involved in cell migration and function in various biological events, including immunity, morphogenesis, and disease. Due to their nature, chemokines tend to be a target of hijacking of immunity by virus and therefore show an exceptionally high mutation rate. Xenopus laevis is considered an excellent model to investigate the effect of whole-genome duplication for gene family evolution. Because its allotetraploidization occurred around 17–18 million years ago, ancestral subgenomes L and S were well conserved. Based on the gene model of human and diploid frog Xenopus tropicalis, we identified 52 chemokine genes and 26 chemokine receptors in X. laevis. The retention rate of the gene in the X. laevis L and S subgenomes was 96% (45/47) and 68% (32/47), respectively. We conducted molecular phylogenetic analysis and found clear orthologies in all receptor genes but not in the ligand genes, suggesting rapid divergences of the ligand. dN/dS calculation demonstrated that dN/dS ratio greater than one was observed in the four ligand genes, cxcl8b.1.S, cxcl18.S, ccl21.S, and xcl1.L, but nothing in receptor genes. These results revealed that the whole-genome duplication promotes diversification of chemokine ligands in X. laevis while conserving the genes necessary for homeostasis, suggesting that selective pressure also supports a rapid divergence of the chemokines in amphibians.

The process of gene duplication followed by gene loss or evolution of new functions has been studied extensively, yet the role gene duplication plays in the function and evolution of fertilization proteins is underappreciated. Gene duplication is observed in many fertilization protein families including Izumo, DCST, ZP, and the TFP superfamily. Molecules mediating fertilization are part of larger gene families expressed in a variety of tissues, but gene duplication followed by structural modifications has often facilitated their cooption into a fertilization function. Repeat expansions of functional domains within a gene also provide opportunities for the evolution of novel fertilization protein. ZP proteins with domain repeat expansions are linked to species-specificity in fertilization and TFP proteins that experienced domain duplications were coopted into a novel sperm function. This review outlines the importance of gene duplications and repeat domain expansions in the evolution of fertilization proteins.

BackgroundGene duplication provides raw material for the evolution of functional innovation. We recently developed a phylogenetic method that classifies evolutionary processes driving the retention of duplicate genes by quantifying divergence between their spatial gene expression profiles and that of their single-copy orthologous gene in a closely related sister species.ResultsHere, we apply our classification method to pairs of duplicate genes in eight mammalian genomes, using data from 11 tissues to construct spatial gene expression profiles. We find that young mammalian duplicates are often functionally conserved, and that expression divergence rapidly increases over evolutionary time. Moreover, expression divergence results in increased tissue specificity, with an overrepresentation of expression in male kidney, underrepresentation of expression in female liver, and strong underrepresentation of expression in testis. Thus, duplicate genes acquire a diversity of new tissue-specific functions outside of the testis, possibly contributing to the origin of a multitude of complex phenotypes during mammalian evolution.ConclusionsOur findings reveal that mammalian duplicate genes are initially functionally conserved, and then undergo rapid functional divergence over evolutionary time, acquiring diverse tissue-specific biological roles. These observations are in stark contrast to the much faster expression divergence and acquisition of broad housekeeping roles we previously observed in Drosophila duplicate genes. Due to the smaller effective population sizes of mammals relative to Drosophila, these analyses implicate natural selection in the functional evolution of duplicate genes.

Cotton is an excellent model for studying crop polyploidization and domestication. Chromatin profiling helps to reveal how histone modifications are involved in controlling differential gene expression between A and D subgenomes in allotetraploid cotton. However, the detailed profiling and functional characterization of broad H3K4me3 and H3K27me3 are still understudied in cotton. In this study, we conducted H3K4me3- and H3K27me3-related ChIP-seq followed by comprehensively characterizing their roles in regulating gene transcription in cotton. We found that H3K4me3 and H3K27me3 exhibited active and repressive roles in regulating the expression of genes between A and D subgenomes, respectively. More importantly, H3K4me3 exhibited enrichment level-, position-, and distance-related impacts on expression levels of related genes. Distinct GO term enrichment occurred between A/D-specific and homeologous genes with broad H3K4me3 enrichment in promoters and gene bodies, suggesting that broad H3K4me3-marked genes might have some unique biological functions between A and D subgenome. An anticorrelation between H3K27me3 enrichment and expression levels of homeologous genes was more pronounced in the A subgenome relative to the D subgenome, reflecting distinct enrichment of H3K27me3 in homeologous genes between A and D subgenome. In addition, H3K4me3 and H3K27me3 marks can indirectly influence gene expression through regulatory networks with TF mediation. Thus, our study provides detailed insights into functions of H3K4me3 and H3K27me3 in regulating differential gene expression and subfunctionalization of homeologous genes, therefore serving as a driving force for polyploidization and domestication in cotton. [ABSTRACT FROM AUTHOR]

SUMMARY: Whole genome duplication (WGD) in plants is typically followed by genomic downsizing, where large portions of the new genome are lost. Whether this downsizing is accompanied by increased or decreased evolutionary rates of the remaining genes is poorly known, not least because homeolog pairings are often obscured by chromosomal rearrangement. Here, we use the newly published genome from a sedge, namely Kobresia littledalei, and CRISPR/Cas‐9 editing to investigate how the Rho WGD event 70 million years ago (MYA) affected transcription factor evolutionary rates, fates, and function in rice (Oryza sativa) and sorghum (Sorghum bicolor). We focus on the 30‐member DNA‐binding with one zinc finger (Dof) transcription factor family in both crops due to their agronomic importance. Using the known speciation dates of rice from Kobresia (97 MYA) and sorghum (50 MYA), we find that rates of amino acid substitution in the critical Dof domain region were over twofold higher during the 20‐million‐year period following the WGD than before or afterward. Through comparison of synteny blocks, we report that at least 11% of Dof genes were purged from 70 to 50 MYA, while only 6% have been lost in the most recent 50‐million‐year interval. CRISPR/Cas9 editing revealed widespread fitness‐related defects in flowering and lack of redundancy of paired members, as well as significant differences in expression between gene pairs. Together these findings demonstrate the strength of Dof genes as a model for deep evolutionary study and offer one of the most detailed portraits yet of the Rho WGD impact on a gene lineage. Significance Statement: This study leverages the newly published Kobresia genome to assess evolutionary rates in Dof transcription factors before and after the Rho whole genome duplication (WGD) 70 million years ago (MYA) in monocots and reports dramatic loss of gene members and increased amino acid substitution rates from 70 to 50 MYA, but lower rates in the 27 million years before and 50 million years after the WGD event. CRISPR/Cas9 editing reveals non‐redundancy in remaining duplicate retained pairs of OsDof genes. [ABSTRACT FROM AUTHOR]

Copyright of Journal of Henan Agricultural Sciences is the property of Editorial Board of Journal of Henan Agricultural Sciences and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Proteinaceous pheromones that diversify through gene duplication can result in shifts in courtship cocktails that may serve as a mechanism for reproductive isolation. The molecular evolution of pheromones has been extensively studied in salamanders, but how these genes and associated novel courtship glands have codiversified has not been evaluated. In this study we used transcriptional analyses to examine the relationship between pheromone diversification and gland type in three divergent lineages of plethodontid salamanders. Our results revealed that plethodontid salamanders express up to eight divergent Sodefrin Precursor-like Factor genes (spf, representing both alpha and beta subfamilies) along with Plethodontid Modulating Factor (pmf) and Plethodontid Receptivity Factor (prf). Expression of pheromone genes is tissue specific with pmf, prf, and some spf genes restricted to the mental gland. In contrast, the caudal gland shows strong expression of the other spf genes. We found evidence for punctuated changes in pheromone cocktail composition related to the loss of metamorphosis, and subsequent extreme reduction of the mental gland, in a paedomorphic lineage. Our study provides insight into how pheromone diversification can be partitioned into unique glands, which may lead to cocktail specificity in behavioral modules during courtship. [ABSTRACT FROM AUTHOR]

Abstract Background Convergent and parallel evolution provide unique insights into the mechanisms of natural selection. Some of the most striking convergent and parallel (collectively recurrent) amino acid substitutions in proteins are adaptive, but there are also many that are selectively neutral. Accordingly, genome-wide assessment has shown that recurrent sequence evolution in orthologs is chiefly explained by nearly neutral evolution. For paralogs, more frequent functional change is expected because additional copies are generally not retained if they do not acquire their own niche. Yet, it is unknown to what extent recurrent sequence differentiation is discernible after independent gene duplications in different eukaryotic taxa. Results We develop a framework that detects patterns of recurrent sequence evolution in duplicated genes. This is used to analyze the genomes of 90 diverse eukaryotes. We find a remarkable number of families with a potentially predictable functional differentiation following gene duplication. In some protein families, more than ten independent duplications show a similar sequence-level differentiation between paralogs. Based on further analysis, the sequence divergence is found to be generally asymmetric. Moreover, about 6% of the recurrent sequence evolution between paralog pairs can be attributed to recurrent differentiation of subcellular localization. Finally, we reveal the specific recurrent patterns for the gene families Hint1/Hint2, Sco1/Sco2 and vma11/vma3. Conclusions The presented methodology provides a means to study the biochemical underpinning of functional differentiation between paralogs. For instance, two abundantly repeated substitutions are identified between independently derived Sco1 and Sco2 paralogs. Such identified substitutions allow direct experimental testing of the biological role of these residues for the repeated functional differentiation. We also uncover a diverse set of families with recurrent sequence evolution and reveal trends in the functional and evolutionary trajectories of this hitherto understudied phenomenon.

Cotton is an excellent model for studying crop polyploidization and domestication. Chromatin profiling helps to reveal how histone modifications are involved in controlling differential gene expression between A and D subgenomes in allotetraploid cotton. However, the detailed profiling and functional characterization of broad H3K4me3 and H3K27me3 are still understudied in cotton. In this study, we conducted H3K4me3- and H3K27me3-related ChIP-seq followed by comprehensively characterizing their roles in regulating gene transcription in cotton. We found that H3K4me3 and H3K27me3 exhibited active and repressive roles in regulating the expression of genes between A and D subgenomes, respectively. More importantly, H3K4me3 exhibited enrichment level-, position-, and distance-related impacts on expression levels of related genes. Distinct GO term enrichment occurred between A/D-specific and homeologous genes with broad H3K4me3 enrichment in promoters and gene bodies, suggesting that broad H3K4me3-marked genes might have some unique biological functions between A and D subgenome. An anticorrelation between H3K27me3 enrichment and expression levels of homeologous genes was more pronounced in the A subgenome relative to the D subgenome, reflecting distinct enrichment of H3K27me3 in homeologous genes between A and D subgenome. In addition, H3K4me3 and H3K27me3 marks can indirectly influence gene expression through regulatory networks with TF mediation. Thus, our study provides detailed insights into functions of H3K4me3 and H3K27me3 in regulating differential gene expression and subfunctionalization of homeologous genes, therefore serving as a driving force for polyploidization and domestication in cotton.

Covalent modifications of histone proteins act as epigenetic regulators of gene expression. We report the distribution of two active histone marks (H3K4me3 and H3K36me3) in 14-day leaves in two lines of Brassica rapa L. by chromatin immunoprecipitation sequencing. Both lines were enriched with H3K4me3 and H3K36me3 marks at the transcription start site, and the transcription level of a gene was associated with the level of H3K4me3 and H3K36me3. H3K4me3- and H3K36me3-marked genes showed low tissue-specific gene expression, and genes with both H3K4me3 and H3K36me3 had a high level of expression and were constitutively expressed. Bivalent active and repressive histone modifications such as H3K4me3 and H3K27me3 marks or antagonistic coexistence of H3K36me3 and H3K27me3 marks were observed in some genes. Expression may be susceptible to changes by abiotic and biotic stresses in genes having both H3K4me3 and H3K27me3 marks. We showed that the presence of H3K36me3 marks was associated with different gene expression levels or tissue specificity between paralogous paired genes, suggesting that H3K36me3 might be involved in subfunctionalization of the subgenomes.

Covalent modifications of histone proteins act as epigenetic regulators of gene expression. We report the distribution of two active histone marks (H3K4me3 and H3K36me3) in 14-day leaves in two lines of Brassica rapa L. by chromatin immunoprecipitation sequencing. Both lines were enriched with H3K4me3 and H3K36me3 marks at the transcription start site, and the transcription level of a gene was associated with the level of H3K4me3 and H3K36me3. H3K4me3- and H3K36me3-marked genes showed low tissue-specific gene expression, and genes with both H3K4me3 and H3K36me3 had a high level of expression and were constitutively expressed. Bivalent active and repressive histone modifications such as H3K4me3 and H3K27me3 marks or antagonistic coexistence of H3K36me3 and H3K27me3 marks were observed in some genes. Expression may be susceptible to changes by abiotic and biotic stresses in genes having both H3K4me3 and H3K27me3 marks. We showed that the presence of H3K36me3 marks was associated with different gene expression levels or tissue specificity between paralogous paired genes, suggesting that H3K36me3 might be involved in subfunctionalization of the subgenomes. [ABSTRACT FROM AUTHOR]

Vertebrates, often considered as the most sophisticatedly organized animals, have a number of unique morphological features that ensured their evolutionary stability and success. Many have argued that the genetic basis for these innovations was provided by the genome duplications occurred at the early stages of vertebrate evolution. One of the important results of genome duplication is the emergence of additional copies of regulatory genes. Having been removed from the constraining pressure of natural selection, these copies acquired an opportunity to modify their structure and functions, influencing the organism's development and morphology. The consequence of ancient genomic duplications is, for example, that about 35% of human genes are represented by at least two homologous copies. The idea of genome duplications occurred at the early stages of vertebrate evolution was first proposed in 1970, but questions about the number, the scale (whole-genome or local) and the evolutionary timing of these duplications are still actively debated. In recent years (2018–2020), due to the rapid development of methods for processing of big data of high-throughput genome sequencing in different lineages of vertebrates and their closest relatives, cephalochordates and tunicates, a number of comparative studies aimed at identifying groups of syntheny in genomes of different evolutionary lines and reconstruction of ancestral vertebrate chromosomes were obtained. As a result, several models describing putative scenarios of genome duplications in vertebrate evolution have been proposed. In parallel, detailed laboratory studies focused on the expression and functional properties of different families of regulatory genes were performed for representatives of several vertebrate groups. In these studies, a plethora of new information on the molecular mechanisms of the embryonic development in hitherto poorly studied representatives of the evolutionarily ancient vertebrate linages (e.g. cyclostomes, cartilaginous fishes and sturgeons) have been collected. In this review, we consider modern concepts of the mechanisms and consequences of genomic duplications in the light of recent experimental data and currently proposed models of vertebrate genome evolution. [ABSTRACT FROM AUTHOR]

Abstract Background Correct timing of flowering is critical for plants to produce enough viable offspring. In Arabidopsis thaliana (Arabidopsis), flowering time is regulated by an intricate network of molecular signaling pathways. Arabidopsis srr1–1 mutants lacking SENSITIVITY TO RED LIGHT REDUCED 1 (SRR1) expression flower early, particularly under short day (SD) conditions (1). SRR1 ensures that plants do not flower prematurely in such non-inductive conditions by controlling repression of the key florigen FT. Here, we have examined the role of SRR1 in the closely related crop species Brassica napus. Results Arabidopsis SRR1 has five homologs in Brassica napus. They can be divided into two groups, where the A02 and C02 copies show high similarity to AtSRR1 on the protein level. The other group, including the A03, A10 and C09 copies all carry a larger deletion in the amino acid sequence. Three of the homologs are expressed at detectable levels: A02, C02 and C09. Notably, the gene copies show a differential expression pattern between spring and winter type accessions of B. napus. When the three expressed gene copies were introduced into the srr1–1 background, only A02 and C02 were able to complement the srr1–1 early flowering phenotype, while C09 could not. Transcriptional analysis of known SRR1 targets in Bna.SRR1-transformed lines showed that CYCLING DOF FACTOR 1 (CDF1) expression is key for flowering time control via SRR1. Conclusions We observed subfunctionalization of the B. napus SRR1 gene copies, with differential expression between early and late flowering accessions of some Bna.SRR1 copies. This suggests involvement of Bna.SRR1 in regulation of seasonal flowering in B. napus. The C09 gene copy was unable to complement srr1–1 plants, but is highly expressed in B. napus, suggesting specialization of a particular function. Furthermore, the C09 protein carries a deletion which may pinpoint a key region of the SRR1 protein potentially important for its molecular function. This is important evidence of functional domain annotation in the highly conserved but unique SRR1 amino acid sequence.

GRAS transcription factors play important roles in many processes of plant development as well as abiotic and biotic stress responses. However, little is known about this gene family in bread wheat (Triticum aestivum), one of the most important crops worldwide. The completion of a quality draft genome allows genome-wide detection and evolutionary analysis of the GRAS gene family in wheat. In this study, 188 TaGRAS genes were detected and divided into 12 subfamilies based on phylogenetic analyses: DELLA, DLT, HAM, LISCL, SCL3, SCL4/7, SCR, SHR, PAT1, Os19, Os4 and LAS. Tandem and segmental duplications are the main contributors to the expansion of TaGRAS, which may contribute to the adaptation of wheat to various environmental conditions. A high rate of homoeolog retention during hexaploidization was detected, suggesting the nonredundancy and biological importance of TaGRAS homoeologs. Systematic analyses of TaGRAS indicated the conserved expression pattern and function of the same subfamily during evolution. In addition, we detected five genes belonging to the LISCL subfamily induced by both biotic and abiotic stresses and they may be potential targets for further research through gene editing. Using degradome and ChIP-seq data, we identified the targets of miR171 and histone modifications and further analyzed the contribution of epigenetic modification to the subfunctionalization of TaGRAS. This study laid a foundation for further functional elucidation of TaGRAS genes.

Gene duplication is an important driver of genomic diversity that can promote adaptive evolution. However, like most mutations, a newly duplicated gene is often deleterious and removed from the genome by drift or natural selection. The early molecular changes that occur soon after duplication therefore may influence the long‐term survival of gene duplicates, but relatively little empirical data exist on the events near the onset of duplication before mutations have time to accumulate. In this study, we contrast gene expression and DNA methylation levels of duplicate genes in the threespine stickleback, Gasterosteus aculeatus, including recently emerged duplications that segregate as copy number variations (CNVs). We find that younger duplicate genes have higher levels of promoter methylation than older genes, and that gene CNVs have higher promoter methylation than non‐CNVs. These results suggest preferential duplication of highly methylated genes or rapid methylation changes soon after duplication. We also find a negative association between methylation and expression, providing a putative role for methylation in suppressing transcription that compensates for increases in gene copy numbers and promoting paralog retention. We propose that methylation contributes to the longevity of young duplicate genes, extending the window of opportunity for functional divergence via mutation. [ABSTRACT FROM AUTHOR]

Learning about the roles that duplicate genes play in the origins of novel phenotypes requires an understanding of how their functions evolve. A previous method for achieving this goal, CDROM, employs gene expression distances as proxies for functional divergence and then classifies the evolutionary mechanisms retaining duplicate genes from comparisons of these distances in a decision tree framework. However, CDROM does not account for stochastic shifts in gene expression or leverage advances in contemporary statistical learning for performing classification, nor is it capable of predicting the parameters driving duplicate gene evolution. Thus, here we develop CLOUD, a multi-layer neural network built on a model of gene expression evolution that can both classify duplicate gene retention mechanisms and predict their underlying evolutionary parameters. We show that not only is the CLOUD classifier substantially more powerful and accurate than CDROM, but that it also yields accurate parameter predictions, enabling a better understanding of the specific forces driving the evolution and long-term retention of duplicate genes. Further, application of the CLOUD classifier and predictor to empirical data from Drosophila recapitulates many previous findings about gene duplication in this lineage, showing that new functions often emerge rapidly and asymmetrically in younger duplicate gene copies, and that functional divergence is driven by strong natural selection. Hence, CLOUD represents a major advancement in classifying retention mechanisms and predicting evolutionary parameters of duplicate genes, thereby highlighting the utility of incorporating sophisticated statistical learning techniques to address long-standing questions about evolution after gene duplication. [ABSTRACT FROM AUTHOR]

People with NR5A1 mutations experience testicular dysgenesis, ovotestes, or adrenal insufficiency, but we do not completely understand the origin of this phenotypic diversity. NR5A1 is expressed in gonadal soma precursor cells before expression of the sex-determining gene SRY. Many fish have two co-orthologs of NR5A1 that likely partitioned ancestral gene subfunctions between them. To explore ancestral roles of NR5A1, we knocked out nr5a1a and nr5a1b in zebrafish. Single-cell RNA-seq identified nr5a1a-expressing cells that co-expressed genes for steroid biosynthesis and the chemokine receptor Cxcl12a in 1-day postfertilization (dpf) embryos, as does the mammalian adrenal-- gonadal (interrenal-gonadal) primordium. In 2dpf embryos, nr5a1a was expressed stronger in the interrenal-gonadal primordium than in the early hypothalamus but nr5a1b showed the reverse. Adult Leydig cells expressed both ohnologs and granulosa cells expressed nr5a1a stronger than nr5a1b. Mutants for nr5a1a lacked the interrenal, formed incompletely differentiated testes, had no Leydig cells, and grew far larger than normal fish. Mutants for nr5a1b formed a disorganized interrenal and their gonads completely disappeared. All homozygous mutant genotypes lacked secondary sex characteristics, including male breeding tubercles and female sex papillae, and had exceedingly low levels of estradiol, 11-ketotestosterone, and cortisol. RNA-seq showed that at 21dpf, some animals were developing as females and others were not, independent of nr5a1 genotype. By 35dpf, all mutant genotypes greatly under-expressed ovary-biased genes. Because adult nr5a1a mutants form gonads but lack an interrenal and conversely, adult nr5a1b mutants lack a gonad but have an interrenal, the adrenal, and gonadal functions of the ancestral nr5a1 gene partitioned between ohnologs after the teleost genome duplication, likely owing to reciprocal loss of ancestral tissue-specific regulatory elements. Identifying such elements could provide hints to otherwise unexplained cases of Differences in Sex Development. [ABSTRACT FROM AUTHOR]

Flowering is a vulnerable, but crucial phase in building crop yield. Proper timing of this period is therefore decisive in obtaining optimal yields. However, genetic regulation of flowering integrates many different environmental signals and is therefore extremely complex. This complexity increases in polyploid crops which carry two or more chromosome sets, like wheat, potato or rapeseed. Here, I summarize the current state of knowledge about flowering time gene copies in rapeseed (Brassica napus), an important oil crop with a complex polyploid history and a close relationship to Arabidopsis thaliana. The current data show a high demand for more targeted studies on flowering time genes in crops rather than in models, allowing better breeding designs and a deeper understanding of evolutionary principles. Over evolutionary time, some copies of rapeseed flowering time genes changed or lost their original role, resulting in subfunctionalization of the respective homologs. For useful applications in breeding, such patterns of subfunctionalization need to be identified and better understood.

In different lineages of C4 plants, the release of CO2 by decarboxylation of a C4 acid near rubisco is catalyzed by NADP-malic enzyme (ME) or NAD-ME, and the facultative use of phosphoenolpyruvate carboxykinase. The co-option of gene lineages during the evolution of C4-NADP-ME has been thoroughly investigated, whereas that of C4-NAD-ME has received less attention. In this work, we aimed at elucidating the mechanism of recruitment of NAD-ME for its function in the C4 pathway by focusing on the eudicot family Cleomaceae. We identified a duplication of NAD-ME in vascular plants that generated the two paralogs lineages: α- and β-NAD-ME. Both gene lineages were retained across seed plants, and their fixation was likely driven by a degenerative process of sub-functionalization, which resulted in a NAD-ME operating primarily as a heteromer of α- and β-subunits. We found most angiosperm genomes maintain a 1:1 β-NAD-ME/α-NAD-ME (β/α) relative gene dosage, but with some notable exceptions mainly due to additional duplications of β-NAD-ME subunits. For example, a significantly high proportion of species with C4-NAD-ME-type photosynthesis have a non-1:1 ratio of β/α. In the Brassicales, we found C4 species with a 2:1 ratio due to a β-NAD-ME duplication (β1 and β2); this was also observed in the C3Tarenaya hassleriana and Brassica crops. In the independently evolved C4 species, Gynandropsis gynandra and Cleome angustifolia, all three genes were affected by C4 evolution with α- and β1-NAD-ME driven by adaptive selection. In particular, the β1-NAD-MEs possess many differentially substituted amino acids compared with other species and the β2-NAD-MEs of the same species. Five of these amino acids are identically substituted in β1-NAD-ME of G. gynandra and C. angustifolia, two of them were identified as positively selected. Using synteny analysis, we established that β-NAD-ME duplications were derived from ancient polyploidy events and that α-NAD-ME is in a unique syntenic context in both Cleomaceae and Brassicaceae. We discuss our hypotheses for the evolution of NAD-ME and its recruitment for C4 photosynthesis. We propose that gene duplications provided the basis for the recruitment of NAD-ME in C4 Cleomaceae and that all members of the NAD-ME gene family have been adapted to fit the C4-biochemistry. Also, one of the β-NAD-ME gene copies was independently co-opted for its function in the C4 pathway.

PurposeThe zebrafish lens is well suited for studies of physiology and development due to its rapid formation in the embryo and genetic accessibility. Aquaporin 0 (AQP0), a lens-specific membrane protein, is required for lens clarity. Zebrafish have two copies of AQP0 (Aqp0a and b), whereas mammals have a single, multifunctional protein. Here we demonstrate a reliable knockdown/rescue system in zebrafish and use it to provide evidence for subfunctionalization of Aqp0a and b, as well as to show that calcium-mediated regulation of Aqp0a in zebrafish lenses is necessary for transparency.MethodsCoinjection of antisense oligonucleotides and DNA rescue constructs into zebrafish embryos, followed by evaluation of the developing fish for cataracts, was used to analyze the functions of Aqp0a and b. The water permeability and regulation characteristics of each rescue protein were tested in a Xenopus oocyte swelling assay.ResultsBoth copies of AQP0 are necessary for lens clarity in the zebrafish, and neither is sufficient. Water permeability is necessary but also insufficient. Phosphorylation and regulation of Aqp0a are required for its function.ConclusionsIn the zebrafish lens, the two closely related AQP0s have acquired distinct functions that are both necessary for lens development and clarity. Regulation of AQP0 water permeability, a well-studied phenomenon in vitro, may be physiologically relevant in the living lens.

Developing soybean lines with high levels of stearic acid is a primary goal of the soybean industry. Most high-stearic-acid soybeans carry different GmSACPD-C mutated alleles. However, due to the dual role of GmSACPD-C in seeds and nodule development, all derived deleterious GmSACPD-C mutant alleles are of extremely poor agronomic value because of defective nodulation. The soybean stearoyl-acyl carrier protein desaturase (GmSACPD) gene family is composed of five members. Comparative genomics analysis indicated that SACPD genes were duplicated and derived from a common ancestor that is still present in chlorophytic algae. Synteny analysis showed the presence of segment duplications between GmSACPD-A/GmSACPD-B , and GmSACPD-C/GmSACPD-D. GmSACPD-E was not contained in any duplicated segment and may be the result of tandem duplication. We developed a TILLING by Target Capture Sequencing (Tilling-by-Sequencing+) technology, a versatile extension of the conventional TILLING by sequencing, and successfully identified 12, 14, and 18 ethyl methanesulfonate mutants at the GmSACPD-A , GmSACPD-B , and GmSACPD-D genes, respectively. Functional analysis of all identified mutants revealed an unprecedented role of GmSACPD-A, GmSACPD-B, and GmSACPD-D in unsaturated fatty acid biosynthesis without affecting nodule development and structure. This discovery will positively impact the development of high-stearic-acid lines to enhance soybean nutritional value without potential developmental tradeoffs. [ABSTRACT FROM AUTHOR]

SUMMARY: Although SQUAMOSA promoter‐binding‐like (SPL) transcription factors are important regulators of development in rice (Oryza sativa), prior assessments of the SPL family have been limited to single genes. A functional comparison across the full gene family in standardized genetic backgrounds has not been reported previously. Here, we demonstrate that the SPL gene family in rice is enriched due to the most recent whole genome duplication (WGD). Notably, 10 of 19 rice SPL genes (52%) cluster in four units that have persisted for at least 50 million years. We show that SPL gene grouping and retention following WGD is widespread in angiosperms, suggesting the conservatism and importance of this gene arrangement. We used Cas9 editing to generate transformation lines for all 19 SPL genes in a common set of backgrounds, and found that knockouts of 14 SPL genes exhibited defects in plant height, 10 exhibited defects in panicle size, and nine had altered grain lengths. We observed subfunctionalization of genes in the paleoduplicated pairs, but little evidence of neofunctionalization. Expression of OsSPL3 was negatively correlated with that of its closest neighbor in its synteny group, OsSPL4, and its sister paired gene, OsSPL12, in the opposing group. Nucleotide diversity was lower in eight of the nine singleton genes in domesticated rice, relative to wild rice, whereas the reverse was true for the paired genes. Together, these results provide functional information on eight previously unexamined OsSPL family members and suggest that paleoduplicate pair redundancy benefits plant survival and innovation. [ABSTRACT FROM AUTHOR]