Origins and functional diversification of salinity‐responsive Na+, K+ ATPase α1 paralogs in salmonids

The Salmoniform whole‐genome duplication is hypothesized to have facilitated the evolution of anadromy, but little is known about the contribution of paralogs from this event to the physiological performance traits required for anadromy, such as salinity tolerance. Here, we determined when two candidate, salinity‐responsive paralogs of the Na⁺, K⁺ ATPase α subunit (α1a and α1b) evolved and studied their evolutionary trajectories and tissue‐specific expression patterns. We found that these paralogs arose during a small‐scale duplication event prior to the Salmoniform, but after the teleost, whole‐genome duplication. The ‘freshwater paralog’ (α1a) is primarily expressed in the gills of Salmoniformes and an unduplicated freshwater sister species (Esox lucius) and experienced positive selection in the freshwater ancestor of Salmoniformes and Esociformes. Contrary to our predictions, the ‘saltwater paralog’ (α1b), which is more widely expressed than α1a, did not experience positive selection during the evolution of anadromy in the Coregoninae and Salmonine. To determine whether parallel mutations in Na⁺, K⁺ ATPase α1 may contribute to salinity tolerance in other fishes, we studied independently evolved salinity‐responsive Na⁺, K⁺ ATPase α1 paralogs in Anabas testudineus and Oreochromis mossambicus. We found that a quarter of the mutations occurring between salmonid α1a and α1b in functionally important sites also evolved in parallel in at least one of these species. Together, these data argue that paralogs contributing to salinity tolerance evolved prior to the Salmoniform whole‐genome duplication and that strong selection and/or functional constraints have led to parallel evolution in salinity‐responsive Na⁺, K⁺ ATPase α1 paralogs in fishes.     

Evolutionary Analysis of the Two-Component Systems in <Emphasis Type="Italic">Pseudomonas aeruginosa</Emphasis> PAO1

Gene organization and functional motif analyses of the 123 two-component system (2CS) genes in Pseudomonas aeruginosa PAO1 were carried out. In addition, NJ and ML trees for the sensor kinases and the response regulators were constructed, and the distances measured and comparatively analyzed. It was apparent that more than half of the sensor-regulator gene pairs, especially the 2CSs with OmpR-like regulators, are derivatives of a common ancestor and have most likely co-evolved through gene pair duplication. Several of the 2CS pairs, especially those with NarL-like regulators, however, appeared to be relatively divergent. This is supportive of the recruitment model, in which a sensor gene and regulator gene with different phylogenetic history are assembled to form a 2CS. Correlation of the classification of sensor kinases and response regulators provides further support for these models. Upon comparison of the phylogenetic trees comprised of sensors and regulators, we have identified six congruent clades, which represent the group of the most recently duplicated 2CS gene pairs. Analyses of the congruent 2CS pairs of each of the clades revealed that certain paralogous 2CS pairs may carry a redundant function even after a gene duplication event. Nevertheless, comparative analysis of the putative promoter regions of the paralogs suggested that functional redundancy could be prevented by a differential control. Both codon usage and G+C content of these 2CS genes were found to be comparable with those of the P. aeruginosa genome, suggesting that they are not newly acquired genes.Reviewing Editor: Dr. Martin Kreitman     

Comparative genomic analysis of CIPK gene family in <Emphasis Type="Italic">Arabidopsis</Emphasis> and <Emphasis Type="Italic">Populus</Emphasis>

Calcium serves as a second messenger in various signal transduction pathways in plants. CBL-interacting protein kinases (CIPKs), which have a variety of functions, are involved in calcium signal transduction. Previous, the studies on CIPK family members focused on Arabidopsis and rice. Here, we present a comparative genomic analysis of the CIPK gene family in Arabidopsis and poplar, a model tree species. Twenty-seven potential CIPKs were identified from poplar using genome-wide analysis. Like the CIPK gene family from Arabidopsis, CIPK genes from poplar were also divided into intron-free and intron-harboring groups. In the intron-harboring group, the intron distribution of CIPKs is rather conserved during the genome evolutionary process. Many homologous gene pairs were found in the CIPK gene family, indicating duplication events might contribute to the amplification of this gene family. The phylogenetic comparison of CIPKs in combination with intron distribution analysis revealed that CIPK genes from both Arabidopsis and poplar might have an ancient origin, which formed earlier than the separation of these two eudicot species. Our genomic and bioinformatic analysis will provide an important foundation for further functional dissection of the CBL-CIPK signaling network in poplars. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Functional divergence of gene duplicates through ectopic recombination

Gene duplication stimulates evolutionary innovation as the resulting paralogs acquire mutations that lead to sub‐ or neofunctionalization. A comprehensive in silico analysis of paralogs in Saccharomyces cerevisiae reveals that duplicates of cell‐surface and subtelomeric genes also undergo ectopic recombination, which leads to new chimaeric alleles. Mimicking such intergenic recombination events in the FLO (flocculation) family of cell‐surface genes shows that chimaeric FLO alleles confer different adhesion phenotypes than the parental genes. Our results indicate that intergenic recombination between paralogs can generate a large set of new alleles, thereby providing the raw material for evolutionary adaptation and innovation.     

Paralogs are revealed by proportion of heterozygotes and deviations in read ratios in genotyping‐by‐sequencing data from natural populations

Whole‐genome duplications have occurred in the recent ancestors of many plants, fish, and amphibians, resulting in a pervasiveness of paralogous loci and the potential for both disomic and tetrasomic inheritance in the same genome. Paralogs can be difficult to reliably genotype and are often excluded from genotyping‐by‐sequencing (GBS) analyses; however, removal requires paralogs to be identified which is difficult without a reference genome. We present a method for identifying paralogs in natural populations by combining two properties of duplicated loci: (i) the expected frequency of heterozygotes exceeds that for singleton loci, and (ii) within heterozygotes, observed read ratios for each allele in GBS data will deviate from the 1:1 expected for singleton (diploid) loci. These deviations are often not apparent within individuals, particularly when sequence coverage is low; but, we postulated that summing allele reads for each locus over all heterozygous individuals in a population would provide sufficient power to detect deviations at those loci. We identified paralogous loci in three species: Chinook salmon (Oncorhynchus tshawytscha) which retains regions with ongoing residual tetrasomy on eight chromosome arms following a recent whole‐genome duplication, mountain barberry (Berberis alpina) which has a large proportion of paralogs that arose through an unknown mechanism, and dusky parrotfish (Scarus niger) which has largely rediploidized following an ancient whole‐genome duplication. Importantly, this approach only requires the genotype and allele‐specific read counts for each individual, information which is readily obtained from most GBS analysis pipelines.     

Global investigation of the co‐evolution of MIRNA genes and microRNA targets during soybean domestication

Although the selection of coding genes during plant domestication has been well studied, the evolution of MIRNA genes (MIRs) and the interaction between microRNAs (miRNAs) and their targets in this process are poorly understood. Here, we present a genome‐wide survey of the selection of MIRs and miRNA targets during soybean domestication and improvement. Our results suggest that, overall, MIRs have higher evolutionary rates than miRNA targets. Nonetheless, they do demonstrate certain similar evolutionary patterns during soybean domestication: MIRs and miRNA targets with high expression and duplication status, and with greater numbers of partners, exhibit lower nucleotide divergence than their counterparts without these characteristics, suggesting that expression level, duplication status, and miRNA–target interaction are essential for evolution of MIRs and miRNA targets. Further investigation revealed that miRNA–target pairs that are subjected to strong purifying selection have greater similarities than those that exhibited genetic diversity. Moreover, mediated by domestication and improvement, the similarities of a large number of miRNA–target pairs in cultivated soybean populations were increased compared to those in wild soybeans, whereas a small number of miRNA–target pairs exhibited decreased similarity, which may be associated with the adoption of particular domestication traits. Taken together, our results shed light on the co‐evolution of MIRs and miRNA targets during soybean domestication.     

The evolution of teleost pigmentation and the fish‐specific genome duplication

Teleost fishes have evolved a unique complexity and diversity of pigmentation and colour patterning that is unmatched among vertebrates. Teleost colouration is mediated by five different major types of neural‐crest derived pigment cells, while tetrapods have a smaller repertoire of such chromatophores. The genetic basis of teleost colouration has been mainly uncovered by the cloning of pigmentation genes in mutants of zebrafish Danio rerio and medaka Oryzias latipes. Many of these teleost pigmentation genes were already known as key players in mammalian pigmentation, suggesting partial conservation of the corresponding developmental programme among vertebrates. Strikingly, teleost fishes have additional copies of many pigmentation genes compared with tetrapods, mainly as a result of a whole‐genome duplication that occurred 320–350 million years ago at the base of the teleost lineage, the so‐called fish‐specific genome duplication. Furthermore, teleosts have retained several duplicated pigmentation genes from earlier rounds of genome duplication in the vertebrate lineage, which were lost in other vertebrate groups. It was hypothesized that divergent evolution of such duplicated genes may have played an important role in pigmentation diversity and complexity in teleost fishes, which therefore not only provide important insights into the evolution of the vertebrate pigmentary system but also allow us to study the significance of genome duplications for vertebrate biodiversity.     

Characterization of genes with tissue-specific differential expression patterns in Populus

Like many plants, Populus has an evolutionary history in which several, both recent and more ancient, genome duplication events have occurred and, therefore, constitutes an excellent model system for studying the functional evolution of genes. In the present study, we have focused on the properties of genes with tissue-specific differential expression patterns in poplar. We identified the genes by analyzing digital expression profiles derived by mapping 90,000+ expressed sequence tags (ESTs) from 18 sources to the predicted genes of Populus. Our sequence analysis suggests that tissue-specific differentially expressed genes have less diverged paralogs than average, indicating that gene duplication events is an important event in the pathway leading to this type of expression pattern. The functional analysis showed that genes coding for proteins involved in processes of functional importance for the specific tissue(s) in which they are expressed and genes coding for regulatory or responsive proteins are most common among the differentially expressed genes, demonstrating that the expression differentiation process is under strong selective pressure. Thus, our data supports a model where gene duplication followed by gene specialization or expansion of the regulatory and responsive networks leads to tissue-specific differential expression patterns. We have also searched for clustering of genes with similar expression pattern into gene-expression neighborhoods within the Populus genome. However, we could not detect any major clustering among the analyzed genes with highly specific expression patterns. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

A transcriptome‐wide study on the microRNA‐ and the Argonaute 1‐enriched small RNA‐mediated regulatory networks involved in plant leaf senescence

Leaf senescence is an important physiological process during the plant life cycle. However, systemic studies on the impact of microRNAs (miRNAs) on the expression of senescence‐associated genes (SAGs) are lacking. Besides, whether other Argonaute 1 (AGO1)‐enriched small RNAs (sRNAs) play regulatory roles in leaf senescence remains unclear. In this study, a total of 5,123 and 1,399 AGO1‐enriched sRNAs, excluding miRNAs, were identified in Arabidopsis thaliana and rice (Oryza sativa), respectively. After retrieving SAGs from the Leaf Senescence Database, all of the AGO1‐enriched sRNAs and the miRBase‐registered miRNAs of these two plants were included for target identification. Supported by degradome signatures, 200 regulatory pairs involving 120 AGO1‐enriched sRNAs and 40 SAGs, and 266 regulatory pairs involving 64 miRNAs and 42 SAGs were discovered in Arabidopsis. Moreover, 13 genes predicted to interact with some of the above‐identified target genes at protein level were validated as regulated by 17 AGO1‐enriched sRNAs and ten miRNAs in Arabidopsis. In rice, only one SAG was targeted by three AGO1‐enriched sRNAs, and one SAG was targeted by miR395. However, five AGO1‐enriched sRNAs were conserved between Arabidopsis and rice. Target genes conserved between the two plants were identified for three of the above five sRNAs, pointing to the conserved roles of these regulatory pairs in leaf senescence or other developmental procedures. Novel targets were discovered for three of the five AGO1‐enriched sRNAs in rice, indicating species‐specific functions of these sRNA–target pairs. These results could advance our understanding of the sRNA‐involved molecular processes modulating leaf senescence.     

Duplicate gene evolution and expression in the wake of vertebrate allopolyploidization

Background  

The mechanism by which duplicate genes originate – whether by duplication of a whole genome or of a genomic segment – influences their genetic fates. To study events that trigger duplicate gene persistence after whole genome duplication in vertebrates, we have analyzed molecular evolution and expression of hundreds of persistent duplicate gene pairs in allopolyploid clawed frogs (Xenopus and Silurana). We collected comparative data that allowed us to tease apart the molecular events that occurred soon after duplication from those that occurred later on. We also quantified expression profile divergence of hundreds of paralogs during development and in different tissues.     

Dual‐color oligo‐FISH can reveal chromosomal variations and evolution in Oryza species

Fluorescence in situ hybridization using probes based on oligonucleotides (oligo‐FISH) is a useful tool for chromosome identification and karyotype analysis. Here we developed two oligo‐FISH probes that allow the identification of each of the 12 pairs of chromosomes in rice (Oryza sativa). These two probes comprised 25 717 (green) and 25 215 (red) oligos (45 nucleotides), respectively, and generated 26 distinct FISH signals that can be used as a barcode to uniquely label each of the 12 pairs of rice chromosomes. Standard karyotypes of rice were established using this system on both mitotic and meiotic chromosomes. Moreover, dual‐color oligo‐FISH was used to characterize diverse chromosomal abnormalities. Oligo‐FISH analyses using these probes in various wild Oryza species revealed that chromosomes from the AA, BB or CC genomes generated specific and intense signals similar to those in rice, while chromosomes with the EE genome generated less specific signals and the FF genome gave no signal. Together, the oligo‐FISH probes we established will be a powerful tool for studying chromosome variations and evolution in the genus Oryza.     

Comparative genomics using <Emphasis Type="Italic">Fugu</Emphasis> reveals insights into regulatory subfunctionalization

Background  

A major mechanism for the preservation of gene duplicates in the genome is thought to be mediated via loss or modification of cis-regulatory subfunctions between paralogs following duplication (a process known as regulatory subfunctionalization). Despite a number of gene expression studies that support this mechanism, no comprehensive analysis of regulatory subfunctionalization has been undertaken at the level of the distal cis-regulatory modules involved. We have exploited fish-mammal genomic alignments to identify and compare more than 800 conserved non-coding elements (CNEs) that associate with genes that have undergone fish-specific duplication and retention.     

Triplicate genes for mitochondrial ADP/ATP carriers in the aerobic yeast <Emphasis Type="Italic">Yarrowia lipolytica</Emphasis> are regulated differentially in the absence of oxygen

Yarrowia lipolytica is a strictly aerobic fungus, which differs from the extensively studied model yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe with respect to its physiology, genetics and dimorphic growth habit. We isolated and sequenced cDNA and genomic clones (YlAAC1) from Y. lipolytica that encode a mitochondrial ADP/ATP carrier. The YlAAC1 gene can complement the S. cerevisiae aac2 deletion mutant. Southern hybridization, analysis of Yarrowia clones obtained in the course of the Génolevures project, and further sequencing revealed the existence of two paralogs of the YlAAC1 gene, which were named YlAAC2 and YlAAC3, respectively. Phylogenetic analysis showed that YlAAC1 and YlAAC2 were more closely related to each other than to YlAAC3, and are likely to represent the products of a recent gene duplication. All three Y. lipolytica YlAAC genes group together on the phylogenetic tree, suggesting that YlAAC3 is derived from a more ancient duplication within the Y. lipolytica lineage. A similar branching pattern for the three ScAAC paralogs in the facultative anaerobe S. cerevisiae demonstrates that two rounds of duplication of AAC genes occurred independently at least twice in the evolution of hemiascomycetous yeasts. Surprisingly, in both the aerobic Y. lipolytica and the facultative anaerobe S. cerevisiae, the three paralogs are differentially regulated in the absence of oxygen. Apparently, Y. lipolytica can sense hypoxia and down-regulate target genes in response.