Gene duplication and the uniqueness of vertebrate genomes circa 1970-1999

Mu opioids and their receptors: evolution of a concept

Opiates are among the oldest medications available to manage a number of medical problems. Although pain is the current focus, early use initially focused upon the treatment of dysentery. Opium contains high concentrations of both morphine and codeine, along with thebaine, which is used in the synthesis of a number of semisynthetic opioid analgesics. Thus, it is not surprising that new agents were initially based upon the morphine scaffold. The concept of multiple opioid receptors was first suggested almost 50 years ago (Martin, 1967), opening the possibility of new classes of drugs, but the morphine-like agents have remained the mainstay in the medical management of pain. Termed mu, our understanding of these morphine-like agents and their receptors has undergone an evolution in thinking over the past 35 years. Early pharmacological studies identified three major classes of receptors, helped by the discovery of endogenous opioid peptides and receptor subtypes-primarily through the synthesis of novel agents. These chemical biologic approaches were then eclipsed by the molecular biology revolution, which now reveals a complexity of the morphine-like agents and their receptors that had not been previously appreciated.

Hox genes regulation in vertebrates

Hox genes encode transcription factors defining cellular identities along the major and secondary body axes. Their coordinated expression in both space and time is critical for embryonic patterning. Accordingly, Hox genes transcription is tightly controlled at multiple levels, and involves an intricate combination of local and long-range cis-regulatory elements. Recent studies revealed that in addition to transcription factors, dynamic patterns of histone marks and higher-order chromatin structure are important determinants of Hox gene regulation. Furthermore, the emerging picture suggests an involvement of various species of non-coding RNA in targeting activating and repressive complexes to Hox clusters. I review these recent developments and discuss their relevance to the control of Hox gene expression in vivo, as well as to our understanding of transcriptional regulatory mechanisms.

Adult teleost heart expresses two distinct troponin C paralogs: cardiac TnC and a novel and teleost-specific ssTnC in a chamber- and temperature-dependent manner

The teleost-specific whole genome duplication created multiple copies of genes allowing for subfunctionalization of isoforms. In this study, we show that the teleost cardiac Ca²⁺-binding troponin C (TnC) is the product of two distinct genes: cardiac TnC (cTnC, TnnC1a) and a fish-specific slow skeletal TnC (ssTnC, TnnC1b). The ssTnC gene is novel to teleosts as mammals have a single gene commonly referred as cTnC but which is also expressed in slow skeletal muscle. In teleosts, the data strongly indicate that these are two TnC genes are different paralogs. Because we determined that ssTnC exists across many teleosts but not in basal ray-finned fish (e.g., bichir), we propose that these paralogs are the result of an ancestral tandem gene duplication persisting only in teleosts. Quantification of mRNA levels was used to demonstrate distinct expression localization patterns of the paralogs within the chambers of the heart. In the adult zebrafish acclimated at 28°C, ssTnC mRNA levels are twofold greater than cTnC mRNA levels in the atrium, whereas cTnC mRNA was almost exclusively expressed in the ventricle. Meanwhile, rainbow trout acclimated at 5°C showed cTnC mRNA levels in both chambers significantly greater than ssTnC. Distinct responses to temperature acclimation were also quantified in both adult zebrafish and rainbow trout, with mRNA in both chambers shifting to express higher levels of cTnC in 18°C acclimated zebrafish and 5°C acclimated trout. Possible subfunctionalization of TnC isoforms may provide insight into how teleosts achieve physiological versatility in chamber-specific contractile properties.

Sequencing of the sea lamprey (Petromyzon marinus) genome provides insights into vertebrate evolution

Lampreys are representatives of an ancient vertebrate lineage that diverged from our own ∼500 million years ago. By virtue of this deeply shared ancestry, the sea lamprey (P. marinus) genome is uniquely poised to provide insight into the ancestry of vertebrate genomes and the underlying principles of vertebrate biology. Here, we present the first lamprey whole-genome sequence and assembly. We note challenges faced owing to its high content of repetitive elements and GC bases, as well as the absence of broad-scale sequence information from closely related species. Analyses of the assembly indicate that two whole-genome duplications likely occurred before the divergence of ancestral lamprey and gnathostome lineages. Moreover, the results help define key evolutionary events within vertebrate lineages, including the origin of myelin-associated proteins and the development of appendages. The lamprey genome provides an important resource for reconstructing vertebrate origins and the evolutionary events that have shaped the genomes of extant organisms.

Watching the grin fade: tracing the effects of polyploidy on different evolutionary time scales

Polyploidy, or whole-genome duplication (WGD), is a recurrent mutation both in cell lineages and over evolutionary time. By globally changing the relationship between gene copy number and other cellular entities, it can induce dramatic changes at the cellular and phenotypic level. Perhaps surprisingly, then, the insights that these events can bring to understanding other cellular features are not as well appreciated as they could be. In this review, we draw on examples of polyploidy from animals, plants and yeast to explore how investigations of polyploid cells have improved our understanding of the cell cycle, biological network complexity, metabolic phenotypes and tumor biology. We argue that the study of polyploidy across organisms, cell types, and time scales serves not only as a window into basic cell biology, but also as a basis for a predictive biology with applications ranging from crop improvement to treating cancer.

Gene duplications and the early evolution of neural crest development

Neural crest cells are an important cell type present in all vertebrates, and elaboration of the neural crest is thought to have been a key factor in their evolutionary success. Genomic comparisons suggest there were two major genome duplications in early vertebrate evolution, raising the possibility that evolution of neural crest was facilitated by gene duplications. Here, we review the process of early neural crest formation and its underlying gene regulatory network (GRN) as well as the evolution of important neural crest derivatives. In this context, we assess the likelihood that gene and genome duplications capacitated neural crest evolution, particularly in light of novel data arising from invertebrate chordates.

Meiotic maps of sockeye salmon derived from massively parallel DNA sequencing

Background

Meiotic maps are a key tool for comparative genomics and association mapping studies. Next-generation sequencing and genotyping by sequencing are speeding the processes of SNP discovery and the development of new genetic tools, including meiotic maps for numerous species. Currently there are limited genetic resources for sockeye salmon, Oncorhynchus nerka. We develop the first dense meiotic map for sockeye salmon using a combination of novel SNPs found in restriction site associated DNA (RAD tags) and SNPs available from existing expressed sequence tag (EST) based assays.

Results

We discovered and genotyped putative SNPs in 3,430 RAD tags. We removed paralogous sequence variants leaving 1,672 SNPs; these were combined with 53 EST-based SNP genotypes for linkage mapping. The map contained 29 male and female linkage groups, consistent with the haploid chromosome number expected for sockeye salmon. The female map contains 1,057 loci spanning 4,896 cM, and the male map contains 1,118 loci spanning 4,220 cM. Regions of conservation with rainbow trout and synteny between the RAD based rainbow trout map and the sockeye salmon map were established.

Conclusions

Using RAD sequencing and EST-based SNP assays we successfully generated the first high density linkage map for sockeye salmon.

Molecular phylogeny of OVOL genes illustrates a conserved C2H2 zinc finger domain coupled by hypervariable unstructured regions

OVO-like proteins (OVOL) are members of the zinc finger protein family and serve as transcription factors to regulate gene expression in various differentiation processes. Recent studies have shown that OVOL genes are involved in epithelial development and differentiation in a wide variety of organisms; yet there is a lack of comprehensive studies that describe OVOL proteins from an evolutionary perspective. Using comparative genomic analysis, we traced three different OVOL genes (OVOL1-3) in vertebrates. One gene, OVOL3, was duplicated during a whole-genome-duplication event in fish, but only the copy (OVOL3b) was retained. From early-branching metazoa to humans, we found that a core domain, comprising a tetrad of C2H2 zinc fingers, is conserved. By domain comparison of the OVOL proteins, we found that they evolved in different metazoan lineages by attaching intrinsically-disordered (ID) segments of N/C-terminal extensions of 100 to 1000 amino acids to this conserved core. These ID regions originated independently across different animal lineages giving rise to different types of OVOL genes over the course of metazoan evolution. We illustrated the molecular evolution of metazoan OVOL genes over a period of 700 million years (MY). This study both extends our current understanding of the structure/function relationship of metazoan OVOL genes, and assembles a good platform for further characterization of OVOL genes from diverged organisms.

Impact of foreign DNA integration on tumor biology and on evolution via epigenetic alterations

The insertion of foreign DNA into mammalian genomes can alter their methylation and transcription patterns at remote sites from the locus of foreign DNA integration. The mechanisms leading to these fundamental changes and their frequencies are unknown. Sites and extent of changes in the recipient cells might depend on the location of foreign DNA integration. In the second part of this review, it will be hypothesized that the insertion event itself, for example, of tumor viral DNA via its epigenetic genome-wide consequences, plays an important role in oncogenesis. During evolution, the impact of ancient retrotransposon or retroviral genomes and the ensuing epigenetic alterations in the recipient genomes might have generated cells with completely different transcriptional profiles. Due to the continued presence of the transgenomes these alterations were genetically stable and were selected for or against by the environmental conditions prevalent at the time. These evolutionary effects are very different from those postulated for insertional mutagenesis, added genetic information or regulatory elements placed into the vicinity of cellular functions.

Genome size evolution: sizing mammalian genomes

The study of genome size (GS) and its variation is so fascinating to the scientific community because it constitutes the link between the present-day analytical and molecular studies of the genome and the old trunk of the holistic and synthetic view of the genome. The GS of several taxa vary over a broad range and do not correlate with the complexity of the organisms (the C-value paradox). However, the biology of transposable elements has let us reach a satisfactory view of the molecular mechanisms that give rise to GS variation and novelties, providing a less perplexing view of the significance of the GS (C-enigma). The knowledge of the composition and structure of a genome is a pre-requisite for trying to understand the evolution of the main genome signature: its size. The radiation of mammals provides an approximately 180-million-year test case for theories of how GS evolves. It has been found from data-mining GS databases that GS is a useful cyto-taxonomical instrument at the level of orders/superorders, providing genomic signatures characterizing Monotremata, Marsupialia, Afrotheria, Xenarthra, Laurasiatheria, and Euarchontoglires. A hypothetical ancestral mammalian-like GS of 2.9-3.7 pg has been suggested. This value appears compatible with the average values calculated for the high systematic levels of the extant Monotremata (∼2.97 pg) and Marsupialia (∼4.07 pg), suggesting invasion of mobile DNA elements concurrently with the separation of the older clades of Afrotheria (∼5.5 pg) and Xenarthra (∼4.5 pg) with larger GS, leaving the Euarchontoglires (∼3.4 pg) and Laurasiatheria (∼2.8 pg) genomes with fewer transposable elements. However, the paucity of GS data (546 mammalian species sized from 5,488 living species) for species, genera, and families calls for caution. Considering that mammalian species may be vanished even before they are known, GS data are sorely needed to phenotype the effects brought about by their variation and to validate any hypotheses on GS evolution in mammals.

Evolution of the osteoblast: skeletogenesis in gar and zebrafish

BACKGROUND

Although the vertebrate skeleton arose in the sea 500 million years ago, our understanding of the molecular fingerprints of chondrocytes and osteoblasts may be biased because it is informed mainly by research on land animals. In fact, the molecular fingerprint of teleost osteoblasts differs in key ways from that of tetrapods, but we do not know the origin of these novel gene functions. They either arose as neofunctionalization events after the teleost genome duplication (TGD), or they represent preserved ancestral functions that pre-date the TGD. Here, we provide evolutionary perspective to the molecular fingerprints of skeletal cells and assess the role of genome duplication in generating novel gene functions. We compared the molecular fingerprints of skeletogenic cells in two ray-finned fish: zebrafish (Danio rerio)--a teleost--and the spotted gar (Lepisosteus oculatus)--a "living fossil" representative of a lineage that diverged from the teleost lineage prior to the TGD (i.e., the teleost sister group). We analyzed developing embryos for expression of the structural collagen genes col1a2, col2a1, col10a1, and col11a2 in well-formed cartilage and bone, and studied expression of skeletal regulators, including the transcription factor genes sox9 and runx2, during mesenchymal condensation.

RESULTS

Results provided no evidence for the evolution of novel functions among gene duplicates in zebrafish compared to the gar outgroup, but our findings shed light on the evolution of the osteoblast. Zebrafish and gar chondrocytes both expressed col10a1 as they matured, but both species' osteoblasts also expressed col10a1, which tetrapod osteoblasts do not express. This novel finding, along with sox9 and col2a1 expression in developing osteoblasts of both zebrafish and gar, demonstrates that osteoblasts of both a teleost and a basally diverging ray-fin fish express components of the supposed chondrocyte molecular fingerprint.

CONCLUSIONS

Our surprising finding that the "chondrogenic" transcription factor sox9 is expressed in developing osteoblasts of both zebrafish and gar can help explain the expression of chondrocyte genes in osteoblasts of ray-finned fish. More broadly, our data suggest that the molecular fingerprint of the osteoblast, which largely is constrained among land animals, was not fixed during early vertebrate evolution.

Cytochrome p450 1 genes in birds: evolutionary relationships and transcription profiles in chicken and Japanese quail embryos

Background

Cytochrome P450 1 (CYP1) genes are biomarkers for aryl hydrocarbon receptor (AHR) agonists and may be involved in some of their toxic effects. CYP1s other than the CYP1As are poorly studied in birds. Here we characterize avian CYP1B and CYP1C genes and the expression of the identified CYP1 genes and AHR1, comparing basal and induced levels in chicken and quail embryos.

Methodology/Principal Findings

We cloned cDNAs of chicken CYP1C1 and quail CYP1B1 and AHR1. CYP1Cs occur in several bird genomes, but we found no CYP1C gene in quail. The CYP1C genomic region is highly conserved among vertebrates. This region also shares some synteny with the CYP1B region, consistent with CYP1B and CYP1C genes deriving from duplication of a common ancestor gene. Real-time RT-PCR analyses revealed similar tissue distribution patterns for CYP1A4, CYP1A5, CYP1B1, and AHR1 mRNA in chicken and quail embryos, with the highest basal expression of the CYP1As in liver, and of CYP1B1 in eye, brain, and heart. Chicken CYP1C1 mRNA levels were appreciable in eye and heart but relatively low in other organs. Basal transcript levels of the CYP1As were higher in quail than in chicken, while CYP1B1 levels were similar in the two species. 3,3′,4,5,5′-Pentachlorobiphenyl induced all CYP1s in chicken; in quail a 1000-fold higher dose induced the CYP1As, but not CYP1B1.

Conclusions/Significance

The apparent absence of CYP1C1 in quail, and weak expression and induction of CYP1C1 in chicken suggest that CYP1Cs have diminishing roles in tetrapods; similar tissue expression suggests that such roles may be met by CYP1B1. Tissue distribution of CYP1B and CYP1C transcripts in birds resembles that previously found in zebrafish, suggesting that these genes serve similar functions in diverse vertebrates. Determining CYP1 catalytic functions in different species should indicate the evolving roles of these duplicated genes in physiological and toxicological processes.

Toll-like receptors in bony fish: from genomics to function

Receptors that recognize conserved pathogen molecules are the first line of cellular innate immunity defense. Toll-like receptors (TLRs) are the best understood of the innate immune receptors that detect infections in mammals. Key features of the fish TLRs and the factors involved in their signaling cascade have high structural similarity to the mammalian TLR system. However, the fish TLRs also exhibit very distinct features and large diversity which is likely derived from their diverse evolutionary history and the distinct environments that they occupy. Six non-mammalian TLRs were identified in fish. TLR14 shares sequence and structural similarity with TLR1 and 2, and the other five (TLR19, 20, 21, 22 and 23) form a cluster of novel TLRs. TLR4 was lost from the genomes of most fishes, and the TLR4 genes found in zebrafish do not recognize the mammalian agonist LPS and are likely paralogous and not orthologous to mammalian TLR4 genes. TLR6 and 10 are also absent from all fish genomes sequenced to date. Of the at least 16 TLR types identified in fish, direct evidence of ligand specificity has only been shown for TLR2, TLR3, TLR5M, TLR5S and TLR22. The common carp TLR2 was shown to recognize the synthetic triacylated lipopeptide Pam(3)CSK(4) and lipopeptides from gram positive bacteria. The membrane-bound TLR5 (TLR5M) signaling in response to flagellin in rainbow trout is amplified through interaction with the soluble form (TLR5S) in a positive loop feedback. In Fugu, TLR3 is localized to the endoplasmic reticulum (ER) and recognizes relatively short dsRNA, while TLR22 has a surveillance function like the human cell-surface TLR3. Genome and gene duplications have been major contributors to the teleost's rich evolutionary history and genomic diversity. Duplicate or multi-copy TLR genes were identified for TLR3 and 7 in common carp, TLR4b, 5, 8 and 20 in zebrafish, TLR8a in rainbow trout and TLR22 in rainbow trout and Atlantic salmon. The main task for current and near-future fish TLRs research is to develop specificity assays to identify the ligands of all fish TLRs, which will advance comparative immunology research and will contribute to our understanding of disease resistance mechanisms in fish and the development of new adjuvants and/or more effective vaccines and therapeutics.

Comparative genomics of duplicate γ-glutamyl transferase genes in teleosts: medaka (Oryzias latipes), stickleback (Gasterosteus aculeatus), green spotted pufferfish (Tetraodon nigroviridis), fugu (Takifugu rubripes), and zebrafish (Danio rerio)

The availability of multiple teleost (bony fish) genomes is providing unprecedented opportunities to understand the diversity and function of gene duplication events using comparative genomics. Here we examine multiple paralogous genes of γ-glutamyl transferase (GGT) in several distantly related teleost species including medaka, stickleback, green spotted pufferfish, fugu, and zebrafish. Through mining genome databases, we have identified multiple GGT orthologs. Duplicate (paralogous) GGT sequences for GGT1 (GGT1 a and b), GGTL1 (GGTL1 a and b), and GGTL3 (GGTL3 a and b) were identified for each species. Phylogenetic analysis suggests that GGTs are ancient proteins conserved across most metazoan phyla and those paralogous GGTs in teleosts likely arose from the serial 3R genome duplication events. A third GGTL1 gene (GGTL1c) was found in green spotted pufferfish; however, this gene is not present in medaka, stickleback, or fugu. Similarly, one or both paralogs of GGTL3 appear to have been lost in green spotted pufferfish, fugu, and zebrafish. Syntenic relationships were highly maintained between duplicated teleost chromosomes, among teleosts and across ray-finned (Actinopterygii) and lobe-finned (Sarcopterygii) species. To assess subfunction partitioning, six medaka GGT genes were cloned and assessed for developmental and tissue-specific expression. On the basis of these data, we propose a modification of the "duplication-degeneration-complementation" model of subfunction partitioning where quantitative differences rather than absolute differences in gene expression are observed between gene paralogs. Our results demonstrate that multiple GGT genes have been retained within teleost genomes. Questions remain, however, regarding the functional roles of multiple GGTs in these species.

Evolution of ABC transporters by gene duplication and their role in human disease

The ATP-binding cassette (ABC) transporter genes represent the largest family of transporters and these genes are abundant in the genome of all vertebrates. Through analysis of the genome sequence databases we have characterized the full complement of ABC genes from several mammals and other vertebrates. Multiple gene duplication and deletion events were identified in ABC genes in different lineages indicating that the process of gene evolution is still ongoing. Gene duplication resulting in either gene birth or gene death plays a major role in the evolution of the vertebrate ABC genes. The understanding of this mechanism is important in the context of human health because these ABC genes are associated with human disease, involving nearly all organ systems of the body. In addition, ABC genes play an important role in the development of drug resistance in cancer cells. Future genetic, functional, and evolutionary studies of ABC transporters will provide important insight into human and animal biology.

Coexistence of diploid, triploid and tetraploid crucian carp (Carassius auratus) in natural waters

Background

Crucian carp (abbreviated CC) belongs to the genus of Carassius within the family of Cyprinidae. It has been one of the most important freshwater species for Chinese aquaculture and is especially abundant in the Dongting water system of Hunan province. CC used to be considered as all diploid forms. However, coexistence of diploid (abbreviated 2nCC), triploid (abbreviated 3nCC) and tetraploid crucian carp (abbreviated 4nCC) population of the Dongting water system was first found by our recently researches.

Results

We examined the ploidy level and compared biological characteristics in different ploidy CC. In reproductive mode, 2nCC was bisexual generative and 4nCC generated all-female offspring by gynogenesis. However, 3nCC generated progenies in two different ways. 3nCC produced bisexual triploid offspring fertilized with 3nCC spermatozoa, while it produced all-female triploid offspring by gynogenesis when its ova were activated by heterogenous spermatozoa. The complete mitochondrial DNA of three different ploidy fishes was sequenced and analyzed, suggesting no significant differences. Interestingly, microchromosomes were found only in 3nCC, which were concluded to be the result of hybridization. Allogenetic DNA fragments of Sox genes were obtained in 3nCC and 4nCC, which were absent in 2nCC. Phylogenetics analysis based on Sox4 gene indicated 3nCC and 4nCC formed a separate group from 2nCC.

Conclusions

In summary, this is the first report of the co-existence of three types of different ploidy crucian carps in natural waters in China. It was proved that the coexistence of different ploidy CC was reproductively maintained. We further hypothesized that 3nCC and 4nCC were allopolyploids that resulted from hybridization. The different ploidy CC population we obtained in this study possesses great significance for the study of polyploidization and the evolution of vertebrates.

Evolutionary mode and functional divergence of vertebrate NMDA receptor subunit 2 genes

Background

Ionotropic glutamate receptors in the central nervous system play a major role in numerous brain functions including learning and memory in many vertebrate species. NR2 subunits have been regarded as rate-limiting molecules in controlling the optimal N-methyl-D-aspartate (NMDA) receptor's coincidence-detection property and subsequent learning and memory function across multi-species. However, its evolutionary mode among vertebrate species remains unclear.

Results

With extensive analysis of phylogeny, exon structure, protein domain, paralogon and synteny, we demonstrated that two-round genome duplication generated quartet GRIN2 genes and the third-round fish-specific genome duplication generated extra copies of fish GRIN2 genes. In addition, in-depth investigation has enabled the identification of three novel genes, GRIN2C_Gg, GRIN2D-1_Ol and GRIN2D-2_Tr in the chicken, medaka and fugu genome, respectively. Furthermore, we showed functional divergence of NR2 genes mostly occurred at the first-round duplication, amino acid residues located at the N-terminal Lig_chan domain were responsible for type I functional divergence between these GRIN2 subfamilies and purifying selection has been the prominent natural pressure operating on these diversified GRIN2 genes.

Conclusion and Significance

These findings provide intriguing subjects for testing the 2R and 3R hypothesis and we expect it could provide new insights into the underlying evolution mechanisms of cognition in vertebrate.

The evolution of vertebrate tetraspanins: gene loss, retention, and massive positive selection after whole genome duplications

BACKGROUND

The vertebrate tetraspanin family has many features which make it suitable for preserving the imprint of ancient sequence evolution and amenable for phylogenomic analysis. So we believe that an in-depth analysis of the tetraspanin evolution not only provides more complete understanding of tetraspanin biology, but offers new insights into the influence of the two rounds of whole genome duplication (2R-WGD) at the origin of vertebrates.

RESULTS

A detailed phylogeny of vertebrate tetraspanins was constructed by using multiple lines of information, including sequence-based phylogenetics, key structural features, intron configuration and genomic synteny. In particular, a total of 38 modern tetraspanin ortholog lineages in bony vertebrates have been identified and subsequently classified into 17 ancestral lineages existing before 2R-WGD. Based on this phylogeny, we found that the ohnolog retention rate of tetraspanins after 2R-WGD was three times as the average (a rate similar to those of transcription factors and protein kinases). This high rate didn't increase the tetrapanin family size, but changed the family composition, possibly by displacing vertebrate-specific gene lineages with the lineages conserved across deuterostomes. We also found that the period from 2R-WGD to recent time is controlled by gene losses. Meanwhile, positive selection has been detected on 80% of the branches right after 2R-WGDs, which declines significantly on both magnitude and extensity on the following speciation branches. Notably, the loss of mammalian RDS2 is accompanied by strong positive selection on mammalian ROM1, possibly due to gene loss-induced compensatory evolution.

CONCLUSIONS

First, different from transcription factors and kinases, high duplicate retention rate after 2R-WGD didn't increase the tetraspanin family size but just reshaped the family composition. Second, the evolution of tetraspanins right after 2R-WGD had been impacted by a massive wave of gene loss and positive selection on coding sequences. Third, the lingering effect of 2R-WGD on tetraspanin gene loss and positive selection might last for 300-400 million years.

Comparative genomics in teleost species: Knowledge transfer by linking the genomes of model and non-model fish species

Comparative genomics is a powerful tool to transfer knowledge coming from model fish species to non-model fish species of economic or/and evolutionary interest. Such transfer is of importance as functional studies either are difficult to perform with most non-model species. The first comparative map constructed using the human and the chimpanzee genome allowed the identification of putative orthologues. Although comparative mapping in teleosts is still in its infancy, five model teleost genomes from different orders have been fully sequenced to date and the sequencing of several commercially important species are also underway or near completion. The accessibility of these whole genome sequences and rapid developments in genomics of fish species are paving the way towards new and valuable research in comparative genetics and genomics. With the accumulation of information in model species, the genetic and genomic characterization of non-model, but economically, physiologically or evolutionary important species is now feasible. Furthermore, comparison of low coverage gene maps of non-model fish species against fully sequenced fish species will enhance the efficiency of candidate gene identification projected for quantitative trait loci (QTL) scans for traits of special interest.

Noncoding sequences near duplicated genes evolve rapidly

Gene expression divergence and chromosomal rearrangements have been put forward as major contributors to phenotypic differences between closely related species. It has also been established that duplicated genes show enhanced rates of positive selection in their amino acid sequences. If functional divergence is largely due to changes in gene expression, it follows that regulatory sequences in duplicated loci should also evolve rapidly. To investigate this hypothesis, we performed likelihood ratio tests (LRTs) on all noncoding loci within 5 kb of every transcript in the human genome and identified sequences with increased substitution rates in the human lineage since divergence from Old World Monkeys. The fraction of rapidly evolving loci is significantly higher nearby genes that duplicated in the common ancestor of humans and chimps compared with nonduplicated genes. We also conducted a genome-wide scan for nucleotide substitutions predicted to affect transcription factor binding. Rates of binding site divergence are elevated in noncoding sequences of duplicated loci with accelerated substitution rates. Many of the genes associated with these fast-evolving genomic elements belong to functional categories identified in previous studies of positive selection on amino acid sequences. In addition, we find enrichment for accelerated evolution nearby genes involved in establishment and maintenance of pregnancy, processes that differ significantly between humans and monkeys. Our findings support the hypothesis that adaptive evolution of the regulation of duplicated genes has played a significant role in human evolution.

Rapid sequence and expression divergence suggest selection for novel function in primate-specific KRAB-ZNF genes

Recent segmental duplications (SDs), arising from duplication events that occurred within the past 35-40 My, have provided a major resource for the evolution of proteins with primate-specific functions. KRAB zinc finger (KRAB-ZNF) transcription factor genes are overrepresented among genes contained within these recent human SDs. Here, we examine the structural and functional diversity of the 70 human KRAB-ZNF genes involved in the most recent primate SD events including genes that arose in the hominid lineage. Despite their recent advent, many parent-daughter KRAB-ZNF gene pairs display significant differences in zinc finger structure and sequence, expression, and splicing patterns, each of which could significantly alter the regulatory functions of the paralogous genes. Paralogs that emerged on the lineage to humans and chimpanzees have undergone more evolutionary changes per unit of time than genes already present in the common ancestor of rhesus macaques and great apes. Taken together, these data indicate that a substantial fraction of the recently evolved primate-specific KRAB-ZNF gene duplicates have acquired novel functions that may possibly define novel regulatory pathways and suggest an active ongoing selection for regulatory diversity in primates.

Identification, characterization and genetic mapping of TLR7, TLR8a1 and TLR8a2 genes in rainbow trout (Oncorhynchus mykiss)

Induction of the innate immune pathways is critical for early anti-viral defense but there is limited understanding of how teleost fish recognize viral molecules and activate these pathways. In mammals, Toll-like receptors (TLR) 7 and 8 bind single-stranded RNA of viral origin and are activated by synthetic anti-viral imidazoquinoline compounds. Herein, we identify and describe the rainbow trout (Oncorhynchus mykiss) TLR7 and TLR8 gene orthologs and their mRNA expression. Two TLR7/8 loci were identified from a rainbow trout bacterial artificial chromosome (BAC) library using DNA fingerprinting and genetic linkage analyses. Direct sequencing of two representative BACs revealed intact omTLR7 and omTLR8a1 open reading frames (ORFs) located on chromosome 3 and a second locus on chromosome 22 that contains an omTLR8a2 ORF and a putative TLR7 pseudogene. We used the omTLR8a1/2 nomenclature for the two trout TLR8 genes as phylogenetic analysis revealed that they and all the other teleost TLR8 genes sequenced to date are similar to the zebrafish TLR8a, but are distinct from the zebrafish TLR8b. The duplicated trout loci exhibit conserved synteny with other fish genomes extending beyond the tandem of TLR7/8 genes. The trout TLR7 and 8a1/2 genes are composed of a single large exon similar to all other described TLR7/8 genes. The omTLR7 ORF is predicted to encode a 1049 amino acid (aa) protein with 84% similarity to the Fugu TLR7 and a conserved pattern of predicted leucine-rich repeats (LRR). The omTLR8a1 and omTLR8a2 are predicted to encode 1035- and 1034-aa proteins, respectively, and have 86% similarity to each other. omTLR8a1 is likely the ortholog of the only Atlantic salmon TLR8 gene described to date as they have 95% aa sequence similarity. The tissue expression profiles of omTLR7, omTLR8a1 and omTLR8a2 in healthy trout were highest in spleen tissue followed by anterior and then posterior kidney tissues. Rainbow trout anterior kidney leukocytes produced elevated levels of pro-inflammatory and type I interferon cytokines mRNA in response to stimulation with the human TLR7/8 agonist R848 or the TLR3 agonist poly I:C. Only poly I:C-induced IFN2 transcription was significantly suppressed in the presence of chloroquine, a compound known to block endosomal acidification and inhibit endosomal maturation. The effect of chloroquine on R848-induced cytokine expression was equivocal and so it remains questionable whether rainbow trout recognition of R848 requires endosomal maturation. TLR7 and TLR8a1 expression levels in rainbow trout anterior kidney leukocytes were not affected by poly I:C or R848 treatments, but surprisingly, TLR8a2 expression was moderately down-regulated by R848. The down-regulation of omTLR8a2 may imply that this gene has evolved to a new or altered function in rainbow trout, as often occurs when the two duplicated genes remain active.

Lineage-specific transcription factors and the evolution of gene regulatory networks

Nature is replete with examples of diverse cell types, tissues and body plans, forming very different creatures from genomes with similar gene complements. However, while the genes and the structures of proteins they encode can be highly conserved, the production of those proteins in specific cell types and at specific developmental time points might differ considerably between species. A full understanding of the factors that orchestrate gene expression will be essential to fully understand evolutionary variety. Transcription factor (TF) proteins, which form gene regulatory networks (GRNs) to act in cooperative or competitive partnerships to regulate gene expression, are key components of these unique regulatory programs. Although many TFs are conserved in structure and function, certain classes of TFs display extensive levels of species diversity. In this review, we highlight families of TFs that have expanded through gene duplication events to create species-unique repertoires in different evolutionary lineages. We discuss how the hierarchical structures of GRNs allow for flexible small to large-scale phenotypic changes. We survey evidence that explains how newly evolved TFs may be integrated into an existing GRN and how molecular changes in TFs might impact the GRNs. Finally, we review examples of traits that evolved due to lineage-specific TFs and species differences in GRNs.

An expressed fgf4 retrogene is associated with breed-defining chondrodysplasia in domestic dogs

Retrotransposition of processed mRNAs is a common source of novel sequence acquired during the evolution of genomes. Although the vast majority of retroposed gene copies, or retrogenes, rapidly accumulate debilitating mutations that disrupt the reading frame, a small percentage become new genes that encode functional proteins. By using a multibreed association analysis in the domestic dog, we demonstrate that expression of a recently acquired retrogene encoding fibroblast growth factor 4 (fgf4) is strongly associated with chondrodysplasia, a short-legged phenotype that defines at least 19 dog breeds including dachshund, corgi, and basset hound. These results illustrate the important role of a single evolutionary event in constraining and directing phenotypic diversity in the domestic dog.

Immunoglobulin superfamily is conserved but evolved rapidly and is active in the silkworm, Bombyx mori

Immunoglobulin superfamily (IgSF) proteins are known for their abilities to specifically recognize and adhere to cells. In this paper, we predicted the presence of 133 IgSF proteins in the silkworm (Bombyx mori) genome. Comparison with similar proteins in other model organisms (Caenorhabditis elegans, Drosophila melanogaster, Anopheles gambiae, Apis mellifera and Homo sapiens) indicated that IgSF proteins are conserved but have rapidly evolved from worms to human beings. However, these proteins are well conserved amongst insects. Silkworm microarray-based expression data showed tissue expression of 57 IgSF genes and microbe-induced differential expression of 37 genes. Based on the expression data, we can conclude that the silkworm IgSF is active.

The role of genome and gene regulatory network canalization in the evolution of multi-trait polymorphisms and sympatric speciation

Background

Sexual reproduction has classically been considered as a barrier to the buildup of discrete phenotypic differentiation. This notion has been confirmed by models of sympatric speciation in which a fixed genetic architecture and a linear genotype phenotype mapping were assumed. In this paper we study the influence of a flexible genetic architecture and non-linear genotype phenotype map on differentiation under sexual reproduction.

We use an individual based model in which organisms have a genome containing genes and transcription factor binding sites. Mutations involve single genes or binding sites or stretches of genome. The genome codes for a regulatory network that determines the gene expression pattern and hence the phenotype of the organism, resulting in a non-linear genotype phenotype map. The organisms compete in a multi-niche environment, imposing selection for phenotypic differentiation.

Results

We find as a generic outcome the evolution of discrete clusters of organisms adapted to different niches, despite random mating. Organisms from different clusters are distinct on the genotypic, the network and the phenotypic level. However, the genome and network differences are constrained to a subset of the genome locations, a process we call genotypic canalization. We demonstrate how this canalization leads to an increased robustness to recombination and increasing hybrid fitness. Finally, in case of assortative mating, we explain how this canalization increases the effectiveness of assortativeness.

Conclusion

We conclude that in case of a flexible genetic architecture and a non-linear genotype phenotype mapping, sexual reproduction does not constrain phenotypic differentiation, but instead constrains the genotypic differences underlying it. We hypothesize that, as genotypic canalization enables differentiation despite random mating and increases the effectiveness of assortative mating, sympatric speciation is more likely than is commonly suggested.

Comparative and developmental study of the immune system in Xenopus

Xenopus laevis is the model of choice for evolutionary, comparative, and developmental studies of immunity, and invaluable research tools including MHC-defined clones, inbred strains, cell lines, and monoclonal antibodies are available for these studies. Recent efforts to use Silurana (Xenopus) tropicalis for genetic analyses have led to the sequencing of the whole genome. Ongoing genome mapping and mutagenesis studies will provide a new dimension to the study of immunity. Here we review what is known about the immune system of X. laevis integrated with available genomic information from S. tropicalis. This review provides compelling evidence for the high degree of similarity and evolutionary conservation between Xenopus and mammalian immune systems. We propose to build a powerful and innovative comparative biomedical model based on modern genetic technologies that takes take advantage of X. laevis and S. tropicalis, as well as the whole Xenopus genus. Developmental Dynamics 238:1249-1270, 2009. (c) 2009 Wiley-Liss, Inc.

Two types of antigen receptor systems in vertebrates

Extant jawless vertebrates, represented by lampreys and hagfishes, have innate immune receptors with variable domains structurally resembling those of T/B-cell receptors. However, they appear to lack cardinal elements of adaptive immunity shared by all jawed vertebrates: major histocompatibility complex molecules and T/B-cell receptors. Thus, it was widely believed that adaptive immunity is unique to jawed vertebrates. Recently, this belief was overturned by the discovery of agnathan antigen receptors named variable lymphocyte receptors. These receptors generate diversity in their antigen-binding sites through assembling highly diverse leucine-rich repeat modules. The crystal structures of hagfish variable lymphocyte receptor monomers indicate that they adopt a horseshoe-shaped structure and likely bind antigens through the hypervariable concave surface. Secreted variable lymphocyte receptors form pentamers or tetramers of dimers and bind antigens with high specificity and avidity. The fact that variable lymphocyte receptors are structurally unrelated to T/B-cell receptors indicates that jawed and jawless vertebrates have developed antigen receptors independently.

Zebrafish neurolin-a and -b, orthologs of ALCAM, are involved in retinal ganglion cell differentiation and retinal axon pathfinding

Neurolin-a and Neurolin-b (also called alcam and nlcam, respectively) are zebrafish orthologs of human ALCAM, an adhesion protein of the immunoglobulin superfamily with functions in axon growth and guidance. Within the developing zebrafish retina, onset and progression of Neurolin-a expression parallels the pattern of retinal ganglion cell (RGC) differentiation. By using a morpholino-based knockdown approach, we show that Neurolin-a (but not Neurolin-b) is necessary for a crucial step in RGC differentiation. Without Neurolin-a, a large proportion of RGCs fail to develop, and RGC axons are absent or reduced in number. Subsequently, Neurolin-a is required for RGC survival and for the differentiation of all other retinal neurons. Neurolin-b is expressed later in well-differentiated RGCs and is required for RGC axon pathfinding. Without Neurolin-b, RGC axons grow in highly aberrant routes along the optic tract and/or fail to reach the optic tectum. Thus, the zebrafish Neurolin paralogs are involved in distinct steps of retinotectal development.

The evolution and functional diversification of animal microRNA genes

microRNAs (miRNAs) are an abundant class of approximately 22 nucleotide (nt) regulatory RNAs that are pervasive in higher eukaryotic genomes. In order to fully understand their prominence in genomes, it is necessary to elucidate the molecular mechanisms that can diversify miRNA activities. In this review, we describe some of the many strategies that allow novel miRNA functions to emerge, with particular emphasis on how miRNA genes evolve in animals. These mechanisms include changes in their sequence, processing, or expression pattern; acquisition of miRNA* functionality or antisense processing; and de novo gene birth. The facility and versatility of miRNAs to evolve and change likely underlies how they have become dominant constituents of higher genomes.

Mechanisms for human genomic rearrangements

Genomic rearrangements describe gross DNA changes of the size ranging from a couple of hundred base pairs, the size of an average exon, to megabases (Mb). When greater than 3 to 5 Mb, such changes are usually visible microscopically by chromosome studies. Human diseases that result from genomic rearrangements have been called genomic disorders. Three major mechanisms have been proposed for genomic rearrangements in the human genome. Non-allelic homologous recombination (NAHR) is mostly mediated by low-copy repeats (LCRs) with recombination hotspots, gene conversion and apparent minimal efficient processing segments. NAHR accounts for most of the recurrent rearrangements: those that share a common size, show clustering of breakpoints, and recur in multiple individuals. Non-recurrent rearrangements are of different sizes in each patient, but may share a smallest region of overlap whose change in copy number may result in shared clinical features among different patients. LCRs do not mediate, but may stimulate non-recurrent events. Some rare NAHRs can also be mediated by highly homologous repetitive sequences (for example, Alu, LINE); these NAHRs account for some of the non-recurrent rearrangements. Other non-recurrent rearrangements can be explained by non-homologous end-joining (NHEJ) and the Fork Stalling and Template Switching (FoSTeS) models. These mechanisms occur both in germ cells, where the rearrangements can be associated with genomic disorders, and in somatic cells in which such genomic rearrangements can cause disorders such as cancer. NAHR, NHEJ and FoSTeS probably account for the majority of genomic rearrangements in our genome and the frequency distribution of the three at a given locus may partially reflect the genomic architecture in proximity to that locus. We provide a review of the current understanding of these three models.

Conservation of notochord gene expression across chordates: insights from the Leprecan gene family

The notochord is a defining character of the chordates, and the T-box transcription factor Brachyury has been shown to be required for notochord development in all chordates examined. In the ascidian Ciona intestinalis, at least 44 notochord genes have been identified as bona fide transcriptional targets of Brachyury. We examined the embryonic expression of a subset of murine orthologs of Ciona Brachyury target genes in the notochord to assess its conservation throughout chordate evolution. We focused on analyzing the Leprecan gene family, which in mouse is composed of three genes, as opposed to the single-copy Ciona gene. We found that all three mouse Leprecan genes are expressed in the notochord. Additionally, while Leprecan expression in C. intestinalis is confined to the notochord, expression of its mouse orthologs includes dorsal root ganglia, limb buds, branchial arches, and developing kidneys. These results have interesting implications for the evolution and development of chordates.

Evolutionary history of the alpha2,8-sialyltransferase (ST8Sia) gene family: tandem duplications in early deuterostomes explain most of the diversity found in the vertebrate ST8Sia genes

BACKGROUND

The animal sialyltransferases, which catalyze the transfer of sialic acid to the glycan moiety of glycoconjugates, are subdivided into four families: ST3Gal, ST6Gal, ST6GalNAc and ST8Sia, based on acceptor sugar specificity and glycosidic linkage formed. Despite low overall sequence identity between each sialyltransferase family, all sialyltransferases share four conserved peptide motifs (L, S, III and VS) that serve as hallmarks for the identification of the sialyltransferases. Currently, twenty subfamilies have been described in mammals and birds. Examples of the four sialyltransferase families have also been found in invertebrates. Focusing on the ST8Sia family, we investigated the origin of the three groups of alpha2,8-sialyltransferases demonstrated in vertebrates to carry out poly-, oligo- and mono-alpha2,8-sialylation.

RESULTS

We identified in the genome of invertebrate deuterostomes, orthologs to the common ancestor for each of the three vertebrate ST8Sia groups and a set of novel genes named ST8Sia EX, not found in vertebrates. All these ST8Sia sequences share a new conserved family-motif, named "C-term" that is involved in protein folding, via an intramolecular disulfide bridge. Interestingly, sequences from Branchiostoma floridae orthologous to the common ancestor of polysialyltransferases possess a polysialyltransferase domain (PSTD) and those orthologous to the common ancestor of oligosialyltransferases possess a new ST8Sia III-specific motif similar to the PSTD. In osteichthyans, we have identified two new subfamilies. In addition, we describe the expression profile of ST8Sia genes in Danio rerio.

CONCLUSION

Polysialylation appeared early in the deuterostome lineage. The recent release of several deuterostome genome databases and paralogons combined with synteny analysis allowed us to obtain insight into events at the gene level that led to the diversification of the ST8Sia genes, with their corresponding enzymatic activities, in both invertebrates and vertebrates. The initial expansion and subsequent divergence of the ST8Sia genes resulted as a consequence of a series of ancient duplications and translocations in the invertebrate genome long before the emergence of vertebrates. A second subset of ST8sia genes in the vertebrate genome arose from whole genome duplication (WGD) R1 and R2. Subsequent selective ST8Sia gene loss is responsible for the characteristic ST8Sia gene expression pattern observed today in individual species.

Transgenesis procedures in Xenopus

Stable integration of foreign DNA into the frog genome has been the purpose of several studies aimed at generating transgenic animals or producing mutations of endogenous genes. Inserting DNA into a host genome can be achieved in a number of ways. In Xenopus, different strategies have been developed which exhibit specific molecular and technical features. Although several of these technologies were also applied in various model organizms, the attributes of each method have rarely been experimentally compared. Investigators are thus confronted with a difficult choice to discriminate which method would be best suited for their applications. To gain better understanding, a transgenesis workshop was organized by the X-omics consortium. Three procedures were assessed side-by-side, and the results obtained are used to illustrate this review. In addition, a number of reagents and tools have been set up for the purpose of gene expression and functional gene analyses. This not only improves the status of Xenopus as a powerful model for developmental studies, but also renders it suitable for sophisticated genetic approaches. Twenty years after the first reported transgenic Xenopus, we review the state of the art of transgenic research, focusing on the new perspectives in performing genetic studies in this species.

The emergence and loss of gonadotropin-releasing hormone in protostomes: orthology, phylogeny, structure, and function

Gonadotropin-releasing hormone (GNRH) is a neuropeptide critical for reproductive activation and maintenance in vertebrates. The recent elucidation of molluscan GNRH-like sequences led to several important questions regarding the evolution of the GNRH family. For instance, are molluscan and chordate GNRHs true orthologs? Has GNRH been retained in most protostomian lineages? What was the function of the ancestral GNRH? The goal of this review is to provide a critical analysis of GNRH evolution based on data available from the known forms of protostomian GNRH. Judging from the orthology between chordate and protostomian GNRH receptors, conservation of several structural motifs on the GNRH peptide, and exon/intron arrangement conserved between protostomian and chordate GNRH genomic sequences, we conclude that chordate and protostomian GNRHs likely share a common ancestor. Based on our analysis of phylogenetic distribution, we also hypothesize that GNRH may have been lost in the ecdysozoan lineage but preserved in lophotrochozoans. Lastly, we propose that the ancestral function of GNRH is to serve as a general neural regulator, and its considerable specialization in reproduction seen in chordates is a consequence of neofunctionalization following gene duplication.

A novel, functional, and highly divergent sex hormone-binding globulin that may participate in the local control of ovarian functions in salmonids

A cDNA encoding for a novel rainbow trout SHBG was identified and characterized. Phylogenetic analysis showed that this novel SHBG, named SHBGb, was a highly divergent paralog of the classical SHBG (SHBGa) form previously known in vertebrates including zebrafish, seabass, and rainbow trout. Using all available sequences, no SHBGb-like sequence could be identified in any fish species besides Atlantic salmon. Rainbow trout SHBGa and SHBGb share only 26% sequence identity at the amino acid level and exhibit totally distinct tissue distribution, thus demonstrating a functional shift of SHBGb. Indeed, shbga mRNA was predominantly expressed in liver and spleen but could not be detected in the ovary, whereas shbgb had a predominant ovarian expression but could not be detected in liver. Despite its high divergence, rainbow trout SHBGb expressed in COS-7 cells could bind estradiol and testosterone with high affinity and specificity. Both rainbow trout shbgb mRNA and proteins were localized to the granulosa cells of vitellogenic ovarian follicles, whereas SHBGb immunoreactivity was also found in theca cells. Finally, shbgb ovarian mRNA expression exhibited a significant drop between late vitellogenesis and oocyte maturation at a time when ovarian aromatase (cyp19a) gene expression and estradiol circulating levels exhibited a dramatic decrease. Together, these observations show that SHBGb is a functional and highly divergent SHBG paralog probably arising from a salmonid-specific duplication of the shbg gene.

Paralogous vitamin D receptors in teleosts: transition of nuclear receptor function

The availability of multiple teleost (bony fish) genomes is providing unprecedented opportunities to understand the diversity and function of gene duplication events using comparative genomics. Here we describe the cloning and functional characterization of two novel vitamin D receptor (VDR) paralogs from the freshwater teleost medaka (Oryzias latipes). VDR sequences were identified through mining of the medaka genome database in which gene organization and structure was determined. Two distinct VDR genes were identified in the medaka genome and mapped to defined loci. Each VDR sequence exhibits unique intronic organization and dissimilar 5' untranslated regions, suggesting they are not isoforms of the same gene locus. Phylogenetic comparison with additional teleosts and mammalian VDR sequences illustrate that two distinct clusters are formed separating aquatic and terrestrial species. Nested within the teleost cluster are two separate clades for VDRalpha and VDRbeta. The topology of teleost VDR sequences is consistent with the notion of paralogous genes arising from a whole genome duplication event prior to teleost radiation. Functional characterization was conducted through the development of VDR expression vectors including Gal4 chimeras containing the yeast Gal4 DNA binding domain fused to the medaka VDR ligand binding domain and full-length protein. The common VDR ligand 1alpha,25-dihydroxyvitamin D3 [1alpha,25(OH)(2)D(3)] resulted in significant transactivation activity with both the Gal4 and full-length constructs of medaka (m) VDRbeta. Comparatively, transactivation of mVDRalpha with 1alpha,25(OH)(2)D(3) was highly attenuated, suggesting a functional divergence between these two nuclear receptor paralogs. We additionally demonstrate through coactivator studies that mVDRalpha is still functional; however, it exhibits a different sensitivity to 1alpha,25(OH)(2)D(3), compared with VDRbeta. These results suggest that in mVDRalpha and VDRbeta have undergone a functional divergence through a process of sub- and/or neofunctionalization of VDR nuclear receptor gene pairs.

Fbxw15/Fbxo12J is an F-box protein-encoding gene selectively expressed in oocytes of the mouse ovary

In recent years, several factors required for follicular assembly and/or early growth of newly formed primordial follicles have been characterized, but additional factors likely remain to be identified. We have used cDNA arrays to compare gene expression in the neonatal mouse ovary at 48 h (when primordial follicles are being assembled) and at 96 h (when early follicular growth is taking place) after birth to that of ovaries collected <24 h after birth (when follicles have not yet been formed). Segregating genes according to their pattern of expression revealed the presence of one cluster of 24 genes for which expression consistently increased at 48 and 96 h. The top increaser in this cluster encodes a approximately 1.5-kb mRNA containing an open reading frame of 1401 bp that encodes a protein of 466 amino acids. The predicted 52.3-kDa protein is a member of the F-box-only (FBXO) protein family, termed FBXW15 or FBXO12J. It has a cytoplasmic localization that includes the endoplasmic reticulum. Expression of Fbxw15/Fbxp12J mRNA is oocyte-specific; the mRNA is first detected on Gestational Day 18, decreasing thereafter to minimal levels on the day of birth. The prevalence of Fbxw15/Fbxp12J mRNA increases again at 48 and 96 h after birth, coinciding with the time of follicular assembly and the initiation of early follicular growth, respectively. The specific expression of Fbxw15/Fbxp12J in oocytes and its developmental pattern of expression suggest a role for this gene in the regulation of oocyte physiology.

Identification and characterization of two closely related histone H4 arginine 3 methyltransferases in Arabidopsis thaliana

Arginine methylation of histone H3 and H4 plays important roles in transcriptional regulation in eukaryotes such as yeasts, fruitflies, nematode worms, fish and mammals; however, less is known in plants. In the present paper, we report the identification and characterization of two Arabidopsis thaliana protein arginine N-methyltransferases, AtPRMT1a and AtPRMT1b, which exhibit high homology with human PRMT1. Both AtPRMT1a and AtPRMT1b methylated histone H4, H2A, and myelin basic protein in vitro. Site-directed mutagenesis of the third arginine (R3) on the N-terminus of histone H4 to lysine (H4R3N) completely abolished the methylation of histone H4. When fused to GFP (green fluorescent protein), both methyltransferases localized to the cytoplasm as well as to the nucleus. Consistent with their subcellular distribution, GST (glutathione transferase) pull-down assays revealed an interaction between the two methyltransferases, suggesting that both proteins may act together in a functional unit. In addition, we demonstrated that AtFib2 (Arabidopsis thaliana fibrillarin 2), an RNA methyltransferase, is a potential substrate for AtPRMT1a and AtPRMT1b, and, furthermore, uncovered a direct interaction between the protein methyltransferase and the RNA methyltransferase. Taken together, our findings implicate AtPRMT1a and AtPRMT1b as H4-R3 protein arginine N-methyltransferases in Arabidopsis and may be involved in diverse biological processes inside and outside the nucleus.

Phylogenomic analyses of KCNA gene clusters in vertebrates: why do gene clusters stay intact?

Background

Gene clusters are of interest for the understanding of genome evolution since they provide insight in large-scale duplications events as well as patterns of individual gene losses. Vertebrates tend to have multiple copies of gene clusters that typically are only single clusters or are not present at all in genomes of invertebrates. We investigated the genomic architecture and conserved non-coding sequences of vertebrate KCNA gene clusters. KCNA genes encode shaker-related voltage-gated potassium channels and are arranged in two three-gene clusters in tetrapods. Teleost fish are found to possess four clusters. The two tetrapod KNCA clusters are of approximately the same age as the Hox gene clusters that arose through duplications early in vertebrate evolution. For some genes, their conserved retention and arrangement in clusters are thought to be related to regulatory elements in the intergenic regions, which might prevent rearrangements and gene loss. Interestingly, this hypothesis does not appear to apply to the KCNA clusters, as too few conserved putative regulatory elements are retained.

Results

We obtained KCNA coding sequences from basal ray-finned fishes (sturgeon, gar, bowfin) and confirmed that the duplication of these genes is specific to teleosts and therefore consistent with the fish-specific genome duplication (FSGD). Phylogenetic analyses of the genes suggest a basal position of the only intron containing KCNA gene in vertebrates (KCNA7). Sistergroup relationships of KCNA1/2 and KCNA3/6 support that a large-scale duplication gave rise to the two clusters found in the genome of tetrapods. We analyzed the intergenic regions of KCNA clusters in vertebrates and found that there are only a few conserved sequences shared between tetrapods and teleosts or between paralogous clusters. The orthologous teleost clusters, however, show sequence conservation in these regions.

Conclusion

The lack of overall conserved sequences in intergenic regions suggests that there are either other processes than regulatory evolution leading to cluster conservation or that the ancestral regulatory relationships among genes in KCNA clusters have been changed together with their regulatory sites.