Gene conversion and the evolution of protocadherin gene cluster diversity

DNA damage drives antigen diversification through mosaic Variant Surface Glycoprotein (VSG) formation in Trypanosoma brucei

Antigenic variation, using large genomic repertoires of antigen-encoding genes, allows pathogens to evade host antibody. Many pathogens, including the African trypanosome Trypanosoma brucei, extend their antigenic repertoire through genomic diversification. While evidence suggests that T. brucei depends on the generation of new variant surface glycoprotein (VSG) genes to maintain a chronic infection, a lack of experimentally tractable tools for studying this process has obscured its underlying mechanisms. Here, we present a highly sensitive targeted sequencing approach for measuring VSG diversification. Using this method, we demonstrate that a Cas9-induced DNA double-strand break within the VSG coding sequence can induce VSG recombination with patterns identical to those observed during infection. These newly generated VSGs are antigenically distinct from parental clones and thus capable of facilitating immune evasion. Together, these results provide insight into the mechanisms of VSG diversification and an experimental framework for studying the evolution of antigen repertoires in pathogenic microbes.

EvolClustDB: Exploring Eukaryotic Gene Clusters with Evolutionarily Conserved Genomic Neighbourhoods

Conservation of gene neighbourhood over evolutionary distances is generally indicative of shared regulation or functional association among genes. This concept has been broadly exploited in prokaryotes but its use on eukaryotic genomes has been limited to specific functional classes, such as biosynthetic gene clusters. We here used an evolutionary-based gene cluster discovery algorithm (EvolClust) to pre-compute evolutionarily conserved gene neighbourhoods, which can be searched, browsed and downloaded in EvolClustDB. We inferred ∼35,000 cluster families in 882 different species in genome comparisons of five taxonomically broad clades: Fungi, Plants, Metazoans, Insects and Protists. EvolClustDB allows browsing through the cluster families, as well as searching by protein, species, identifier or sequence. Visualization allows inspecting gene order per species in a phylogenetic context, so that relevant evolutionary events such as gain, loss or transfer, can be inferred. EvolClustDB is freely available, without registration, at http://evolclustdb.org/.

The HUSH complex controls brain architecture and protocadherin fidelity

The HUSH (human silencing hub) complex contains the H3K9me3 binding protein M-phase phosphoprotein 8 (MPP8) and recruits the histone methyltransferase SETDB1 as well as Microrchidia CW-type zinc finger protein 2 (MORC2). Functional and mechanistic studies of the HUSH complex have hitherto been centered around SETDB1 while the in vivo functions of MPP8 and MORC2 remain elusive. Here, we show that genetic inactivation of Mphosph8 or Morc2a in the nervous system of mice leads to increased brain size, altered brain architecture, and behavioral changes. Mechanistically, in both mouse brains and human cerebral organoids, MPP8 and MORC2 suppress the repetitive-like protocadherin gene cluster in an H3K9me3-dependent manner. Our data identify MPP8 and MORC2, previously linked to silencing of repetitive elements via the HUSH complex, as key epigenetic regulators of protocadherin expression in the nervous system and thereby brain development and neuronal individuality in mice and humans.

Conversion between duplicated genes generated by polyploidization contributes to the divergence of poplar and willow

BACKGROUND

Gene conversion has an important effect on duplicate genes produced by polyploidization. Poplar (Populus trichocarpa) and willow (Salix brachista) are leading models and excellent green plants in the Salicaceae. Although much attention has been paid to the evolution of duplicated genes in poplar and willow, the role of conversion between duplicates generated from polyploidization remains poorly understood.

RESULTS

Here, through genomic synteny analyses, we identified duplicate genes generated by the Salicaceae common tetraploidization (SCT) in the poplar and willow genomes. We estimated that at least 0.58% and 0.25% of poplar and willow duplicates were affected by whole-gene conversion after the poplar-willow divergence, with more (5.73% and 2.66%) affected by partial-gene conversion. Moreover, we found that the converted duplicated genes were unevenly distributed on each chromosome in the two genomes, and the well-preserved homoeologous chromosome regions may facilitate the conversion of duplicates. Notably, we found that conversion maintained the similarity of duplicates, likely contributing to the conservation of certain sequences, but is essentially accelerated the rate of evolution and increased species divergence. In addition, we found that converted duplicates tended to have more similar expression patterns than nonconverted duplicates. We found that genes associated with multigene families were preferentially converted. We also found that the genes encoding conserved structural domains associated with specific traits exhibited a high frequency of conversion.

CONCLUSIONS

Extensive conversion between duplicate genes generated from the SCT contributes to the diversification of the family Salicaceae and has had long-lasting effects on those genes with important biological functions.

Mouse models for the study of clustered protocadherins

Since their first description, the clustered protocadherins (cPcdhs) have sparked interest for their potential to generate diverse cell-surface recognition cues and their widespread expression in the nervous system. Through the use of mouse models, we have learned a great deal about the functions served by cPcdhs, and how their molecular diversity is regulated. cPcdhs are essential contributors to a host of processes during neural circuit formation, including neuronal survival, dendritic and axonal branching, self-avoidance and targeting, and synapse formation. Their expression is controlled by the interplay of epigenetic marks with proximal and distal elements involving high order DNA looping, regulating transcription factor binding. Here, we will review various mouse models targeting the cPcdh locus and how they have been instructive in uncovering the regulation and function of the cPcdhs.

Illegitimate Recombination between Duplicated Genes Generated from Recursive Polyploidizations Accelerated the Divergence of the Genus Arachis

The peanut (Arachis hypogaea L.) is the leading oil and food crop among the legume family. Extensive duplicate gene pairs generated from recursive polyploidizations with high sequence similarity could result from gene conversion, caused by illegitimate DNA recombination. Here, through synteny-based comparisons of two diploid and three tetraploid peanut genomes, we identified the duplicated genes generated from legume common tetraploidy (LCT) and peanut recent allo-tetraploidy (PRT) within genomes. In each peanut genome (or subgenomes), we inferred that 6.8–13.1% of LCT-related and 11.3–16.5% of PRT-related duplicates were affected by gene conversion, in which the LCT-related duplicates were the most affected by partial gene conversion, whereas the PRT-related duplicates were the most affected by whole gene conversion. Notably, we observed the conversion between duplicates as the long-lasting contribution of polyploidizations accelerated the divergence of different Arachis genomes. Moreover, we found that the converted duplicates are unevenly distributed across the chromosomes and are more often near the ends of the chromosomes in each genome. We also confirmed that well-preserved homoeologous chromosome regions may facilitate duplicates’ conversion. In addition, we found that these biological functions contain a higher number of preferentially converted genes, such as catalytic activity-related genes. We identified specific domains that are involved in converted genes, implying that conversions are associated with important traits of peanut growth and development.

Evolutionary Perspective and Expression Analysis of Intronless Genes Highlight the Conservation of Their Regulatory Role

The structure of eukaryotic genes is generally a combination of exons interrupted by intragenic non-coding DNA regions (introns) removed by RNA splicing to generate the mature mRNA. A fraction of genes, however, comprise a single coding exon with introns in their untranslated regions or are intronless genes (IGs), lacking introns entirely. The latter code for essential proteins involved in development, growth, and cell proliferation and their expression has been proposed to be highly specialized for neuro-specific functions and linked to cancer, neuropathies, and developmental disorders. The abundant presence of introns in eukaryotic genomes is pivotal for the precise control of gene expression. Notwithstanding, IGs exempting splicing events entail a higher transcriptional fidelity, making them even more valuable for regulatory roles. This work aimed to infer the functional role and evolutionary history of IGs centered on the mouse genome. IGs consist of a subgroup of genes with one exon including coding genes, non-coding genes, and pseudogenes, which conform approximately 6% of a total of 21,527 genes. To understand their prevalence, biological relevance, and evolution, we identified and studied 1,116 IG functional proteins validating their differential expression in transcriptomic data of embryonic mouse telencephalon. Our results showed that overall expression levels of IGs are lower than those of MEGs. However, strongly up-regulated IGs include transcription factors (TFs) such as the class 3 of POU (HMG Box), Neurog1, Olig1, and BHLHe22, BHLHe23, among other essential genes including the β-cluster of protocadherins. Most striking was the finding that IG-encoded BHLH TFs fit the criteria to be classified as microproteins. Finally, predicted protein orthologs in other six genomes confirmed high conservation of IGs associated with regulating neural processes and with chromatin organization and epigenetic regulation in Vertebrata. Moreover, this study highlights that IGs are essential modulators of regulatory processes, such as the Wnt signaling pathway and biological processes as pivotal as sensory organ developing at a transcriptional and post-translational level. Overall, our results suggest that IG proteins have specialized, prevalent, and unique biological roles and that functional divergence between IGs and MEGs is likely to be the result of specific evolutionary constraints.

Wiring the Brain by Clustered Protocadherin Neural Codes

There are more than a thousand trillion specific synaptic connections in the human brain and over a million new specific connections are formed every second during the early years of life. The assembly of these staggeringly complex neuronal circuits requires specific cell-surface molecular tags to endow each neuron with a unique identity code to discriminate self from non-self. The clustered protocadherin (Pcdh) genes, which encode a tremendous diversity of cell-surface assemblies, are candidates for neuronal identity tags. We describe the adaptive evolution, genomic structure, and regulation of expression of the clustered Pcdhs. We specifically focus on the emerging 3-D architectural and biophysical mechanisms that generate an enormous number of diverse cell-surface Pcdhs as neural codes in the brain.

Cephalopod Biology: At the Intersection Between Genomic and Organismal Novelties

Cephalopods are resourceful marine predators that have fascinated generations of researchers as well as the public owing to their advanced behavior, complex nervous system, and significance in evolutionary studies. Recent advances in genomics have accelerated the pace of cephalopod research. Many traditional areas focusing on evolution, development, behavior, and neurobiology, primarily on the morphological level, are now transitioning to molecular approaches. This review addresses the recent progress and impact of genomic and other molecular resources on research in cephalopods. We outline several key directions in which significant progress in cephalopod research is expected and discuss its impact on our understanding of the genetic background behind cephalopod biology and beyond.

Writing, Reading, and Translating the Clustered Protocadherin Cell Surface Recognition Code for Neural Circuit Assembly

The ability of neurites of individual neurons to distinguish between themselves and neurites from other neurons and to avoid self (self-avoidance) plays a key role in neural circuit assembly in both invertebrates and vertebrates. Similarly, when individual neurons of the same type project into receptive fields of the brain, they must avoid each other to maximize target coverage (tiling). Counterintuitively, these processes are driven by highly specific homophilic interactions between cell surface proteins that lead to neurite repulsion rather than adhesion. Among these proteins in vertebrates are the clustered protocadherins (Pcdhs), and key to their function is the generation of enormous cell surface structural diversity. Here we review recent advances in understanding how a Pcdh cell surface code is generated by stochastic promoter choice; how this code is amplified and read by homophilic interactions between Pcdh complexes at the surface of neurons; and, finally, how the Pcdh code is translated to cellular function, which mediates self-avoidance and tiling and thus plays a central role in the development of complex neural circuits. Not surprisingly, Pcdh mutations that diminish homophilic interactions lead to wiring defects and abnormal behavior in mice, and sequence variants in the Pcdh gene cluster are associated with autism spectrum disorders in family-based genetic studies in humans.

Germline Structural Variations Are Preferential Sites of DNA Replication Timing Plasticity during Development

The DNA replication timing program is modulated throughout development and is also one of the main factors influencing the distribution of mutation rates across the genome. However, the relationship between the mutagenic influence of replication timing and its developmental plasticity remains unexplored. Here, we studied the distribution of copy number variations (CNVs) and single nucleotide polymorphisms across the zebrafish genome in relation to changes in DNA replication timing during embryonic development in this model vertebrate species. We show that CNV sites exhibit strong replication timing plasticity during development, replicating significantly early during early development but significantly late during more advanced developmental stages. Reciprocally, genomic regions that changed their replication timing during development contained a higher proportion of CNVs than developmentally constant regions. Developmentally plastic CNV sites, in particular those that become delayed in their replication timing, were enriched for the clustered protocadherins, a set of genes important for neuronal development that have undergone extensive genetic and epigenetic diversification during zebrafish evolution. In contrast, single nucleotide polymorphism sites replicated consistently early throughout embryonic development, highlighting a unique aspect of the zebrafish genome. Our results uncover a hitherto unrecognized interface between development and evolution.

In silico Identification and Expression of Protocadherin Gene Family in Octopus vulgaris

Connecting millions of neurons to create a functional neural circuit is a daunting challenge. Vertebrates developed a molecular system at the cell membrane to allow neurons to recognize each other by distinguishing self from non-self through homophilic protocadherin interactions. In mammals, the protocadherin gene family counts about 50 different genes. By hetero-multimerization, protocadherins are capable of generating an impressive number of molecular interfaces. Surprisingly, in the California two-spot octopus, Octopus bimaculoides, an invertebrate belonging to the Phylum Mollusca, over 160 protocadherins (PCDHs) have been identified. Here we briefly discuss the role of PCDHs in neural wiring and conduct a comparative study of the protocadherin gene family in two closely related octopus species, Octopus vulgaris and O. bimaculoides. A first glance at the expression patterns of protocadherins in O. vulgaris is also provided. Finally, we comment on PCDH evolution in the light of invertebrate nervous system plasticity.

Contrasting patterns of coding and flanking region evolution in mammalian keratin associated protein-1 genes

Mammalian genomes contain a number of duplicated genes, and sequence identity between these duplicates can be maintained by purifying selection. However, between-duplicate recombination can also maintain sequence identity between copies, resulting in a pattern known as concerted evolution where within-genome repeats are more similar to each other than to orthologous repeats in related species. Here we investigated the tandemly-repeated keratin-associated protein 1 (KAP1) gene family, KRTAP1, which encodes proteins that are important components of hair and wool in mammals. Comparison of eutherian mammal KRTAP1 gene repeats within and between species shows a strong pattern of concerted evolution. However, in striking contrast to the coding regions of these genes, we find that the flanking regions have a divergent pattern of evolution. This contrast in evolutionary pattern transitions abruptly near the start and stop codons of the KRTAP1 genes. We reveal that this difference in evolutionary patterns is not explained by conventional purifying selection, nor is it likely a consequence of codon adaptation or reverse transcription of KRTAP1-n mRNA. Instead, the evidence suggests that these contrasting patterns result from short-tract gene conversion events that are biased to the KRTAP1 coding region by selection and/or differential sequence divergence. This work demonstrates the power that gene conversion has to finely shape the evolution of repetitive genes, and provides another distinctive pattern of contrasting evolutionary outcomes that results from gene conversion. A greater emphasis on exploring the evolution of multi-gene eukaryotic families will reveal how common different contrasting evolutionary patterns are in gene duplicates.

Teleost Fish-Specific Preferential Retention of Pigmentation Gene-Containing Families After Whole Genome Duplications in Vertebrates

Vertebrate pigmentation is a highly diverse trait mainly determined by neural crest cell derivatives. It has been suggested that two rounds (1R/2R) of whole-genome duplications (WGDs) at the basis of vertebrates allowed changes in gene regulation associated with neural crest evolution. Subsequently, the teleost fish lineage experienced other WGDs, including the teleost-specific Ts3R before teleost radiation and the more recent Ss4R at the basis of salmonids. As the teleost lineage harbors the highest number of pigment cell types and pigmentation diversity in vertebrates, WGDs might have contributed to the evolution and diversification of the pigmentation gene repertoire in teleosts. We have compared the impact of the basal vertebrate 1R/2R duplications with that of the teleost-specific Ts3R and salmonid-specific Ss4R WGDs on 181 gene families containing genes involved in pigmentation. We show that pigmentation genes (PGs) have been globally more frequently retained as duplicates than other genes after Ts3R and Ss4R but not after the early 1R/2R. This is also true for non-pigmentary paralogs of PGs, suggesting that the function in pigmentation is not the sole key driver of gene retention after WGDs. On the long-term, specific categories of PGs have been repeatedly preferentially retained after ancient 1R/2R and Ts3R WGDs, possibly linked to the molecular nature of their proteins (e.g., DNA binding transcriptional regulators) and their central position in protein-protein interaction networks. Taken together, our results support a major role of WGDs in the diversification of the pigmentation gene repertoire in the teleost lineage, with a possible link with the diversity of pigment cell lineages observed in these animals compared to other vertebrates.

Expansion of gene clusters, circular orders, and the shortest Hamiltonian path problem

Clusters of paralogous genes such as the famous HOX cluster of developmental transcription factors tend to evolve by stepwise duplication of its members, often involving unequal crossing over. Gene conversion and possibly other mechanisms of concerted evolution further obfuscate the phylogenetic relationships. As a consequence, it is very difficult or even impossible to disentangle the detailed history of gene duplications in gene clusters. In this contribution we show that the expansion of gene clusters by unequal crossing over as proposed by Walter Gehring leads to distinctive patterns of genetic distances, namely a subclass of circular split systems. Furthermore, when the gene cluster was left undisturbed by genome rearrangements, the shortest Hamiltonian paths with respect to genetic distances coincide with the genomic order. This observation can be used to detect ancient genomic rearrangements of gene clusters and to distinguish gene clusters whose evolution was dominated by unequal crossing over within genes from those that expanded through other mechanisms.

Regulation of neural circuit formation by protocadherins

The protocadherins (Pcdhs), which make up the most diverse group within the cadherin superfamily, were first discovered in the early 1990s. Data implicating the Pcdhs, including ~60 proteins encoded by the tandem Pcdha, Pcdhb, and Pcdhg gene clusters and another ~10 non-clustered Pcdhs, in the regulation of neural development have continually accumulated, with a significant expansion of the field over the past decade. Here, we review the many roles played by clustered and non-clustered Pcdhs in multiple steps important for the formation and function of neural circuits, including dendrite arborization, axon outgrowth and targeting, synaptogenesis, and synapse elimination. We further discuss studies implicating mutation or epigenetic dysregulation of Pcdh genes in a variety of human neurodevelopmental and neurological disorders. With recent structural modeling of Pcdh proteins, the prospects for uncovering molecular mechanisms of Pcdh extracellular and intracellular interactions, and their role in normal and disrupted neural circuit formation, are bright.

Role of Ectopic Gene Conversion in the Evolution of a Candida krusei Pleiotropic Drug Resistance Transporter Family

Gene duplications enable the evolution of novel gene function, but strong positive selection is required to preserve advantageous mutations in a population. This is because frequent ectopic gene conversions (EGCs) between highly similar, tandem-duplicated, sequences, can rapidly remove fate-determining mutations by replacing them with the neighboring parent gene sequences. Unfortunately, the high sequence similarities between tandem-duplicated genes severely hamper empirical studies of this important evolutionary process, because deciphering their correct sequences is challenging. In this study, we employed the eukaryotic model organism Saccharomyces cerevisiae to clone and functionally characterize all 30 alleles of an important pair of tandem-duplicated multidrug efflux pump genes, ABC1 and ABC11, from seven strains of the diploid pathogenic yeast Candida krusei Discovery and functional characterization of their closest ancestor, C. krusei ABC12, helped elucidate the evolutionary history of the entire gene family. Our data support the proposal that the pleiotropic drug resistance (PDR) transporters Abc1p and Abc11p have evolved by concerted evolution for ∼134 MY. While >90% of their sequences remained identical, very strong purifying selection protected six short DNA patches encoding just 18 core amino acid (aa) differences in particular trans membrane span (TMS) regions causing two distinct efflux pump functions. A proline-kink change at the bottom of Abc11p TMS3 was possibly fate determining. Our data also enabled the first empirical estimates for key parameters of eukaryotic gene evolution, they provided rare examples of intron loss, and PDR transporter phylogeny confirmed that C. krusei belongs to a novel, yet unnamed, third major Saccharomycotina lineage.

Characterization of a Single Genomic Locus Encoding the Clustered Protocadherin Receptor Diversity in Xenopus tropicalis

Clustered protocadherins (cPcdhs) constitute the largest subgroup of the cadherin superfamily, and in mammals are grouped into clusters of α-, β-, and γ-types. Tens of tandemly arranged paralogous Pcdh genes of the Pcdh clusters generate a substantial diversity of receptor isoforms. cPcdhs are known to have important roles in neuronal development, and genetic alterations of cPcdhs have been found to be associated with several neurological diseases. Here, we present a first characterization of cPcdhs in Xenopus tropicalis. We determined and annotated all cPcdh isoforms, revealing that they are present in a single chromosomal locus. We validated a total of 96 isoforms, which we show are organized in three distinct clusters. The X. tropicalis cPcdh locus is composed of one α- and two distinct γ-Pcdh clusters (pcdh-γ1 and pcdh-γ2). Bioinformatics analyses assisted by genomic BAC clone sequencing showed that the X. tropicalis α- and γ-Pcdhs are conserved at the cluster level, but, unlike mammals, X. tropicalis does not contain a β-Pcdh cluster. In contrast, the number of γ-Pcdh isoforms has expanded, possibly due to lineage-specific gene duplications. Interestingly, the number of X. tropicalis α-Pcdhs is identical between X. tropicalis and mouse. Moreover, we find highly conserved as well as novel promoter elements potentially involved in regulating the cluster-specific expression of cPcdh isoforms. This study provides important information for the understanding of the evolutionary history of cPcdh genes and future mechanistic studies. It provides an annotated X. tropicalis cPcdh genomic map and a first molecular characterization essential for functional and comparative studies.

Functional test of PCDHB11, the most human-specific neuronal surface protein

Background

Brain-expressed proteins that have undergone functional change during human evolution may contribute to human cognitive capacities, and may also leave us vulnerable to specifically human diseases, such as schizophrenia, autism or Alzheimer’s disease. In order to search systematically for those proteins that have changed the most during human evolution and that might contribute to brain function and pathology, all proteins with orthologs in chimpanzee, orangutan and rhesus macaque and annotated as being expressed on the surface of cells in the human central nervous system were ordered by the number of human-specific amino acid differences that are fixed in modern populations.

Results

PCDHB11, a beta-protocadherin homologous to murine cell adhesion proteins, stood out with 12 substitutions and maintained its lead after normalizing for protein size and applying weights for amino acid exchange probabilities. Human PCDHB11 was found to cause homophilic cell adhesion, but at lower levels than shown for other clustered protocadherins. Homophilic adhesion caused by a PCDHB11 with reversion of human-specific changes was as low as for modern human PCDHB11; while neither human nor reverted PCDHB11 adhered to controls, they did adhere to each other. A loss of function in PCDHB11 is unlikely because intra-human variability did not increase relative to the other human beta-protocadherins.

Conclusions

The brain-expressed protein with the highest number of human-specific substitutions is PCDHB11. In spite of its fast evolution and low intra-human variability, cell-based tests on the only proposed function for PCDHB11 did not indicate a functional change.

The octopus genome and the evolution of cephalopod neural and morphological novelties

Coleoid cephalopods (octopus, squid, and cuttlefish) are active, resourceful predators with a rich behavioral repertoire¹. They have the largest nervous systems among the invertebrates² and present other striking morphological innovations including camera-like eyes, prehensile arms, a highly derived early embryogenesis, and the most sophisticated adaptive coloration system among all animals^1,3. To investigate the molecular bases of cephalopod brain and body innovations we sequenced the genome and multiple transcriptomes of the California two-spot octopus, Octopus bimaculoides. We found no evidence for hypothesized whole genome duplications in the octopus lineage^4–6. The core developmental and neuronal gene repertoire of the octopus is broadly similar to that found across invertebrate bilaterians, except for massive expansions in two gene families formerly thought to be uniquely enlarged in vertebrates: the protocadherins, which regulate neuronal development, and the C2H2 superfamily of zinc finger transcription factors. Extensive mRNA editing generates transcript and protein diversity in genes involved in neural excitability, as previously described⁷, as well as in genes participating in a broad range of other cellular functions. We identified hundreds of cephalopod-specific genes, many of which showed elevated expression levels in such specialized structures as the skin, the suckers, and the nervous system. Finally, we found evidence for large-scale genomic rearrangements that are closely associated with transposable element expansions. Our analysis suggests that substantial expansion of a handful of gene families, along with extensive remodeling of genome linkage and repetitive content, played a critical role in the evolution of cephalopod morphological innovations, including their large and complex nervous systems.

Molecular cloning, sequencing and tissue expression of vasotocin and isotocin precursor genes from Ostariophysian catfishes: phylogeny and evolutionary considerations in teleosts

Basic and neutral neurohypophyseal (NH) nonapeptides have evolved from vasotocin (VT) by a gene duplication at the base of the gnathostome lineage. In teleosts, VT and IT are the basic and neutral peptides, respectively. In the present study, VT and IT precursor genes of Heteropneustes fossilis and Clarias batrachus (Siluriformes, Ostariophysi) were cloned and sequenced. The channel catfish Icatalurus punctatus NH precursor sequences were obtained from EST database. The catfish NH sequences were used along with the available Acanthopterygii and other vertebrate NH precursor sequences to draw phylogenetic inference on the evolutionary history of the teleost NH peptides. Synteny analysis of the NH gene loci in various teleost species was done to complement the phylogenetic analysis. In H. fossilis, the NH transcripts were also sequenced from the ovary. The cloned genes and the deduced precursor proteins showed conserved characteristics of the NH nonapeptide precursors. The genes are expressed in brain and ovary (follicular envelope) of H. fossilis with higher transcript abundance in the brain. The addition of the catfish sequences in the phylogenetic analysis revealed that the VT and IT precursors of the species-rich superorders of teleosts have a distinct phylogenetic history with the Acanthopterygii VT and IT precursors sharing a less evolutionary distance and the Ostariophysi VT and IT having a greater evolutionary distance. The genomic location of VT and IT precursors, and synteny analysis of the NH loci lend support to the phylogenetic inference and suggest a footprint of fish- specific whole genome duplication (3R) and subsequent diploidization in the NH loci. The VT and IT precursor genes are most likely lineage-specific paralogs resulting from differential losses of the 3R NH paralogs in the two superorders. The independent yet consistent retention of VT and IT in the two superorders might be directed by a stringent ligand-receptor selectivity.

Distinctive patterns of evolution of the δ-globin gene (HBD) in primates

In most vertebrates, hemoglobin (Hb) is a heterotetramer composed of two dissimilar globin chains, which change during development according to the patterns of expression of α- and β-globin family members. In placental mammals, the β-globin cluster includes three early-expressed genes, ε(HBE)-γ(HBG)-ψβ(HBBP1), and the late expressed genes, δ (HBD) and β (HBB). While HBB encodes the major adult β-globin chain, HBD is weakly expressed or totally silent. Paradoxically, in human populations HBD shows high levels of conservation typical of genes under strong evolutionary constraints, possibly due to a regulatory role in the fetal-to-adult switch unique of Anthropoid primates. In this study, we have performed a comprehensive phylogenetic and comparative analysis of the two adult β-like globin genes in a set of diverse mammalian taxa, focusing on the evolution and functional divergence of HBD in primates. Our analysis revealed that anthropoids are an exception to a general pattern of concerted evolution in placental mammals, showing a high level of sequence conservation at HBD, less frequent and shorter gene conversion events. Moreover, this lineage is unique in the retention of a functional GATA-1 motif, known to be involved in the control of the developmental expression of the β-like globin genes. We further show that not only the mode but also the rate of evolution of the δ-globin gene in higher primates are strictly associated with the fetal/adult β-cluster developmental switch. To gain further insight into the possible functional constraints that have been shaping the evolutionary history of HBD in primates, we calculated dN/dS (ω) ratios under alternative models of gene evolution. Although our results indicate that HBD might have experienced different selective pressures throughout primate evolution, as shown by different ω values between apes and Old World Monkeys + New World Monkeys (0.06 versus 0.43, respectively), these estimates corroborated a constrained evolution for HBD in Anthropoid lineages, which is unlikely to be related to protein function. Collectively, these findings suggest that sequence change at the δ-globin gene has been under strong selective constraints over 65 Myr of primate evolution, likely due to a regulatory role in ontogenic switches of gene expression.

Mammalian keratin associated proteins (KRTAPs) subgenomes: disentangling hair diversity and adaptation to terrestrial and aquatic environments

Background

Adaptation of mammals to terrestrial life was facilitated by the unique vertebrate trait of body hair, which occurs in a range of morphological patterns. Keratin associated proteins (KRTAPs), the major structural hair shaft proteins, are largely responsible for hair variation.

Results

We exhaustively characterized the KRTAP gene family in 22 mammalian genomes, confirming the existence of 30 KRTAP subfamilies evolving at different rates with varying degrees of diversification and homogenization. Within the two major classes of KRTAPs, the high cysteine (HS) subfamily experienced strong concerted evolution, high rates of gene conversion/recombination and high GC content. In contrast, high glycine-tyrosine (HGT) KRTAPs showed evidence of positive selection and low rates of gene conversion/recombination. Species with more hair and of higher complexity tended to have more KRATP genes (gene expansion). The sloth, with long and coarse hair, had the most KRTAP genes (175 with 141 being intact). By contrast, the “hairless” dolphin had 35 KRTAPs and the highest pseudogenization rate (74% relative to the 19% mammalian average). Unique hair-related phenotypes, such as scales (armadillo) and spines (hedgehog), were correlated with changes in KRTAPs. Gene expression variation probably also influences hair diversification patterns, for example human have an identical KRTAP repertoire as apes, but much less hair.

Conclusions

We hypothesize that differences in KRTAP gene repertoire and gene expression, together with distinct rates of gene conversion/recombination, pseudogenization and positive selection, are likely responsible for micro and macro-phenotypic hair diversification among mammals in response to adaptations to ecological pressures.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-779) contains supplementary material, which is available to authorized users.

Expansion of stochastic expression repertoire by tandem duplication in mouse Protocadherin-α cluster

Tandem duplications are concentrated within the Pcdh cluster throughout vertebrate evolution and as copy number variations (CNVs) in human populations, but the effects of tandem duplication in the Pcdh cluster remain elusive. To investigate the effects of tandem duplication in the Pcdh cluster, here we generated and analyzed a new line of the Pcdh cluster mutant mice. In the mutant allele, a 218-kb region containing the Pcdh-α2 to Pcdh-αc2 variable exons with their promoters was duplicated and the individual duplicated Pcdh isoforms can be disctinguished. The individual duplicated Pcdh-α isoforms showed diverse expression level with stochastic expression manner, even though those have an identical promoter sequence. Interestingly, the 5'-located duplicated Pcdh-αc2, which is constitutively expressed in the wild-type brain, shifted to stochastic expression accompanied by increased DNA methylation. These results demonstrate that tandem duplication in the Pcdh cluster expands the stochastic expression repertoire irrespective of sequence divergence.

Whole-genome duplication in teleost fishes and its evolutionary consequences

Whole-genome duplication (WGD) events have shaped the history of many evolutionary lineages. One such duplication has been implicated in the evolution of teleost fishes, by far the most species-rich vertebrate clade. After initial controversy, there is now solid evidence that such event took place in the common ancestor of all extant teleosts. It is termed teleost-specific (TS) WGD. After WGD, duplicate genes have different fates. The most likely outcome is non-functionalization of one duplicate gene due to the lack of selective constraint on preserving both. Mechanisms that act on preservation of duplicates are subfunctionalization (partitioning of ancestral gene functions on the duplicates), neofunctionalization (assigning a novel function to one of the duplicates) and dosage selection (preserving genes to maintain dosage balance between interconnected components). Since the frequency of these mechanisms is influenced by the genes' properties, there are over-retained classes of genes, such as highly expressed ones and genes involved in neural function. The consequences of the TS-WGD, especially its impact on the massive radiation of teleosts, have been matter of controversial debate. It is evident that gene duplications are crucial for generating complexity and that WGDs provide large amounts of raw material for evolutionary adaptation and innovation. However, it is less clear whether the TS-WGD is directly linked to the evolutionary success of teleosts and their radiation. Recent studies let us conclude that TS-WGD has been important in generating teleost complexity, but that more recent ecological adaptations only marginally related to TS-WGD might have even contributed more to diversification. It is likely, however, that TS-WGD provided teleosts with diversification potential that can become effective much later, such as during phases of environmental change.

[High gene conversion frequency between genes encoding 2-deoxyglucose-6-phosphate phosphatase in 3 Saccharomyces species]

Gene conversions are nonreciprocal sequence exchanges between genes. They are relatively common in Saccharomyces cerevisiae, but few studies have investigated the evolutionary fate of gene conversions or their functional impacts. Here, we analyze the evolution and impact of gene conversions between the two genes encoding 2-deoxyglucose-6-phosphate phosphatase in S. cerevisiae, Saccharomyces paradoxus and Saccharomyces mikatae. Our results demonstrate that the last half of these genes are subject to gene conversions among these three species. The greater similarity and the greater percentage of GC nucleotides in the converted regions, as well as the absence of long regions of adjacent common converted sites, suggest that these gene conversions are frequent and occur independently in all three species. The high frequency of these conversions probably result from the fact that they have little impact on the protein sequences encoded by these genes.

Both positive and negative selection pressures contribute to the polymorphism pattern of the duplicated human CYP21A2 gene

The human steroid 21-hydroxylase gene (CYP21A2) participates in cortisol and aldosterone biosynthesis, and resides together with its paralogous (duplicated) pseudogene in a multiallelic copy number variation (CNV), called RCCX CNV. Concerted evolution caused by non-allelic gene conversion has been described in great ape CYP21 genes, and the same conversion activity is responsible for a serious genetic disorder of CYP21A2, congenital adrenal hyperplasia (CAH). In the current study, 33 CYP21A2 haplotype variants encoding 6 protein variants were determined from a European population. CYP21A2 was shown to be one of the most diverse human genes (HHe=0.949), but the diversity of intron 2 was greater still. Contrary to previous findings, the evolution of intron 2 did not follow concerted evolution, although the remaining part of the gene did. Fixed sites (different fixed alleles of sites in human CYP21 paralogues) significantly accumulated in intron 2, indicating that the excess of fixed sites was connected to the lack of effective non-allelic conversion and concerted evolution. Furthermore, positive selection was presumably focused on intron 2, and possibly associated with the previous genetic features. However, the positive selection detected by several neutrality tests was discerned along the whole gene. In addition, the clear signature of negative selection was observed in the coding sequence. The maintenance of the CYP21 enzyme function is critical, and could lead to negative selection, whereas the presumed gene regulation altering steroid hormone levels via intron 2 might help fast adaptation, which broadly characterizes the genes of human CNVs responding to the environment.

Purifying selection against gene conversions in the folate receptor genes of primates

We characterized the gene conversions between the human folate receptor (FOLR) genes and those of five other primate species. We found 26 gene conversions having an average length of 534 nucleotides. The length of these conversions is correlated with sequence similarity, converted regions have a higher GC-content and the average size of converted regions from a functional donor to another functional donor is significantly smaller than the average size from a functional donor to a pseudogene. Furthermore, the few conversions observed in the FOLR1 and FOLR2 genes did not change any amino acids in their coding regions and did not affect their promoter regions. In contrast, the promoter and coding regions of the FOLR3 gene are frequently converted and these conversions changed many amino acids in marmoset. These results suggest that purifying selection is limiting the functional impact that frequent gene conversions have on functional folate receptor genes.

The genome-wide landscape of copy number variations in the MUSGEN study provides evidence for a founder effect in the isolated Finnish population

Here we characterized the genome-wide architecture of copy number variations (CNVs) in 286 healthy, unrelated Finnish individuals belonging to the MUSGEN study, where molecular background underlying musical aptitude and related traits are studied. By using Illumina HumanOmniExpress-12v.1.0 beadchip, we identified 5493 CNVs that were spread across 467 different cytogenetic regions, spanning a total size of 287.83 Mb (∼9.6% of the human genome). Merging the overlapping CNVs across samples resulted in 999 discrete copy number variable regions (CNVRs), of which ∼6.9% were putatively novel. The average number of CNVs per person was 20, whereas the average size of CNV per locus was 52.39 kb. Large CNVs (>1 Mb) were present in 4% of the samples. The proportion of homozygous deletions in this data set (∼12.4%) seemed to be higher when compared with three other populations. Interestingly, several CNVRs were significantly enriched in this sample set, whereas several others were totally depleted. For example, a CNVR at chr2p22.1 intersecting GALM was more common in this population (P=3.3706 × 10(-44)) than in African and other European populations. The enriched CNVRs, however, showed no significant association with music-related phenotypes. Moreover, the most common CNV locations in world's normal population cohorts (6q14.1, 11q11) were overrepresented in this population. Thus, the genome-wide CNV investigation in this Finnish sample set demonstrated features that are characteristic to isolated populations. Novel CNVRs and the functional implications of CNVs revealed in this study elucidate structural variation present in this population isolate, and may also serve as candidate gene loci for music-related traits.

Protocadherin-17 function in Zebrafish retinal development

Cadherin cell adhesion molecules play crucial roles in vertebrate development including the development of the retina. Most studies have focused on examining functions of classic cadherins (e.g. N-cadherin) in retinal development. There is little information on the function of protocadherins in the development of the vertebrate visual system. We previously showed that protocadherin-17 mRNA was expressed in developing zebrafish retina during critical stages of the retinal development. To gain insight into protocadherin-17 function in the formation of the retina, we analyzed eye development and differentiation of retinal cells in zebrafish embryos injected with protocadherin-17 specific antisense morpholino oligonucleotides (MOs). Protocadherin-17 knockdown embryos (pcdh17 morphants) had significantly reduced eyes due mainly to decreased cell proliferation. Differentiation of several retinal cell types (e.g. retinal ganglion cells) was also disrupted in the pcdh17 morphants. Phenotypic rescue was achieved by injection of protocadherin-17 mRNA. Injection of a vivo-protocadherin-17 MO into one eye of embryonic zebrafish resulted in similar eye defects. Our results suggest that protocadherin-17 plays an important role in the normal formation of the zebrafish retina.

Stronger purifying selection against gene conversions in a pathogenic Saccharomyces cerevisiae strain

Gene conversions most often have no selective impact, but some are selectively disadvantageous whereas others are selectively advantageous. Although gene conversions have been extensively studied in yeasts, very little is known about their selective impact in pathological yeasts. Here, we used the GENECONV software to compare the characteristics of candidate gene conversions found in a pathogenic strain (YJM789) and a nonpathogenic strain (S288c) of Saccharomyces cerevisiae. Interestingly, the pathogenic strain has fewer gene conversions when compared with the nonpathogenic strain. Of the 123 conversions we identified, 27 were identical or similar between the two strains, 62 were specific to the S288c strain, and 34 were specific to the YJM789 strain. Identical and similar conversions likely represent conversions that are under similar levels of purifying selection in both strains. The lower number of gene conversions in most gene families of the pathogenic strain is likely the result of higher purifying selection in this strain. In contrast, the higher number of conversions found in the YRF1 helicase gene family of the pathogenic strain could represent an example of adaptive gene conversions involved in maintaining its telomeres.

Tracing back the nascence of a new sex-determination pathway to the ancestor of bees and ants

In several Hymenoptera, sexual fate is determined by the allelic composition at the complementary sex-determiner locus, a sex-determination mechanism that can strongly affect population dynamics. To date, the molecular identification of complementary sex determiner has only been achieved in the honeybee, where the complementary sex-determiner gene was reported to have arisen from duplication of the feminizer gene. Strikingly, the complementary sex-determiner gene was also proposed to be unique to the honeybee lineage. Here we identify feminizer and complementary sex-determiner orthologues in bumble bees and ants. We further demonstrate that the duplication of feminizer that produced complementary sex determiner occurred before the divergence of Aculeata species (~120 Myr ago). Finally, we provide evidence that the two genes evolved concertedly through gene conversion, complementary sex-determiner evolution being additionally shaped by mosaic patterns of selection. Thus, the complementary sex-determiner gene likely represents the molecular basis for single locus-complementary sex determination in the Aculeata infra-order, and possibly, in the entire Hymenoptera order.

In several Hymenoptera species - ants, bees and wasps - sexual fate is determined by the allelic composition at the complementary sex - determiner locus. This study identifies the honeybee complementary sex - determiner in bumble bee and ant orthologues, previously thought to be unique to the honeybee lineage.

The clustered protocadherins Pcdhα and Pcdhγ form a heteromeric complex in zebrafish

The clustered protocadherin genes encode a diverse collection of neuronal cell surface receptors. These genes have been proposed to play roles in axon targeting, synaptic development and neuronal survival, although their specific cellular roles remain poorly defined. In zebrafish there are four clustered protocadherin genes, two pcdhα clusters and two pcdhγ clusters, that give rise to over 100 distinct proteins, each with a distinct ectodomain (EC). The zebrafish is an excellent model in which to address the function of protocadherins during neural development, as the embryos are transparent, develop rapidly, and are amenable to experimental manipulation. As a first step to investigating the clustered protocadherins during zebrafish development, we have generated antibodies against the common cytodomains of zebrafish Pcdhγ. We compare the distribution of Pcdhγ with Pcdhα and find a similar pan-neuronal pattern, with strong labeling of neurons within all major regions of the central nervous system. Pcdhα and Pcdhγ are particularly enriched in the developing visual system, with strong labeling found in the synaptic layers of the retina, as well as the optic tectum. Consistent with studies in mouse, we find that Pcdhα and Pcdhγ are present in a complex, as they can be co-immunoprecipitated from zebrafish larval extracts. This interaction is direct and occurs through the ECs of these proteins. Using standard bead aggregation assays, we find no evidence for intrinsic adhesive ability by either Pcdhγ or Pcdhα, suggesting that they do not function as cell adhesion molecules.

Genetic variations and haplotype diversity of the UGT1 gene cluster in the Chinese population

Vertebrates require tremendous molecular diversity to defend against numerous small hydrophobic chemicals. UDP-glucuronosyltransferases (UGTs) are a large family of detoxification enzymes that glucuronidate xenobiotics and endobiotics, facilitating their excretion from the body. The UGT1 gene cluster contains a tandem array of variable first exons, each preceded by a specific promoter, and a common set of downstream constant exons, similar to the genomic organization of the protocadherin (Pcdh), immunoglobulin, and T-cell receptor gene clusters. To assist pharmacogenomics studies in Chinese, we sequenced nine first exons, promoter and intronic regions, and five common exons of the UGT1 gene cluster in a population sample of 253 unrelated Chinese individuals. We identified 101 polymorphisms and found 15 novel SNPs. We then computed allele frequencies for each polymorphism and reconstructed their linkage disequilibrium (LD) map. The UGT1 cluster can be divided into five linkage blocks: Block 9 (UGT1A9), Block 9/7/6 (UGT1A9, UGT1A7, and UGT1A6), Block 5 (UGT1A5), Block 4/3 (UGT1A4 and UGT1A3), and Block 3′ UTR. Furthermore, we inferred haplotypes and selected their tagSNPs. Finally, comparing our data with those of three other populations of the HapMap project revealed ethnic specificity of the UGT1 genetic diversity in Chinese. These findings have important implications for future molecular genetic studies of the UGT1 gene cluster as well as for personalized medical therapies in Chinese.

Molecular codes for neuronal individuality and cell assembly in the brain

The brain contains an enormous, but finite, number of neurons. The ability of this limited number of neurons to produce nearly limitless neural information over a lifetime is typically explained by combinatorial explosion; that is, by the exponential amplification of each neuron's contribution through its incorporation into "cell assemblies" and neural networks. In development, each neuron expresses diverse cellular recognition molecules that permit the formation of the appropriate neural cell assemblies to elicit various brain functions. The mechanism for generating neuronal assemblies and networks must involve molecular codes that give neurons individuality and allow them to recognize one another and join appropriate networks. The extensive molecular diversity of cell-surface proteins on neurons is likely to contribute to their individual identities. The clustered protocadherins (Pcdh) is a large subfamily within the diverse cadherin superfamily. The clustered Pcdh genes are encoded in tandem by three gene clusters, and are present in all known vertebrate genomes. The set of clustered Pcdh genes is expressed in a random and combinatorial manner in each neuron. In addition, cis-tetramers composed of heteromultimeric clustered Pcdh isoforms represent selective binding units for cell-cell interactions. Here I present the mathematical probabilities for neuronal individuality based on the random and combinatorial expression of clustered Pcdh isoforms and their formation of cis-tetramers in each neuron. Notably, clustered Pcdh gene products are known to play crucial roles in correct axonal projections, synaptic formation, and neuronal survival. Their molecular and biological features induce a hypothesis that the diverse clustered Pcdh molecules provide the molecular code by which neuronal individuality and cell assembly permit the combinatorial explosion of networks that supports enormous processing capability and plasticity of the brain.

Identification of CTCF as a master regulator of the clustered protocadherin genes

The brain is a large and complex network of neurons. Specific neuronal connectivity is thought to be based on the combinatorial expression of the 52 protocadherins (Pcdh) membrane adhesion proteins, whereby each neuron expresses only a specific subset. Pcdh genes are arranged in tandem, in a cluster of three families: Pcdhα, Pcdhβ and Pcdhγ. The expression of each Pcdh gene is regulated by a promoter that has a regulatory conserved sequence element (CSE), common to all 52 genes. The mechanism and factors controlling individual Pcdh gene expression are currently unknown. Here we show that the promoter of each Pcdh gene contains a gene-specific conserved control region, termed specific sequence element (SSE), located adjacent and upstream to the CSE and activates transcription together with the CSE. We purified the complex that specifically binds the SSE-CSE region and identified the CCTC binding-factor (CTCF) as a key molecule that binds and activates Pcdh promoters. Our findings point to CTCF as a factor essential for Pcdh expression and probably governing neuronal connectivity.

Evolution of orthologous tandemly arrayed gene clusters

BACKGROUND

Tandemly Arrayed Gene (TAG) clusters are groups of paralogous genes that are found adjacent on a chromosome. TAGs represent an important repertoire of genes in eukaryotes. In addition to tandem duplication events, TAG clusters are affected during their evolution by other mechanisms, such as inversion and deletion events, that affect the order and orientation of genes. The DILTAG algorithm developed in 1 makes it possible to infer a set of optimal evolutionary histories explaining the evolution of a single TAG cluster, from an ancestral single gene, through tandem duplications (simple or multiple, direct or inverted), deletions and inversion events.

RESULTS

We present a general methodology, which is an extension of DILTAG, for the study of the evolutionary history of a set of orthologous TAG clusters in multiple species. In addition to the speciation events reflected by the phylogenetic tree of the considered species, the evolutionary events that are taken into account are simple or multiple tandem duplications, direct or inverted, simple or multiple deletions, and inversions. We analysed the performance of our algorithm on simulated data sets and we applied it to the protocadherin gene clusters of human, chimpanzee, mouse and rat.

CONCLUSIONS

Our results obtained on simulated data sets showed a good performance in inferring the total number and size distribution of duplication events. A limitation of the algorithm is however in dealing with multiple gene deletions, as the algorithm is highly exponential in this case, and becomes quickly intractable.

Gene conversions in the growth hormone gene family of primates: stronger homogenizing effects in the Hominidae lineage

In humans, the growth hormone/chorionic somatomammotropin gene family is composed of five highly similar genes. We characterized the gene conversions that occurred between the growth hormone genes of 11 primate species. We detected 48 conversions using GENECONV and others were only detected using phylogenetic analyses. Gene conversions were detected in all species analyzed, their average size (±standard deviation) is 197.8±230.4 nucleotides, the size of the conversions is correlated with sequence similarity and converted regions are significantly more GC-rich than non-converted regions. Gene conversions have a stronger homogenizing effect in Hominidae genes than in other primate species. They are also less frequent in conserved gene regions and towards functionally important genes. This suggests that the high degree of sequence similarity observed between the growth hormone genes of primate species is a consequence of frequent gene conversions in gene regions which are under little selective constraints.

Complete genome sequence of Treponema paraluiscuniculi, strain Cuniculi A: the loss of infectivity to humans is associated with genome decay

Treponema paraluiscuniculi is the causative agent of rabbit venereal spirochetosis. It is not infectious to humans, although its genome structure is very closely related to other pathogenic Treponema species including Treponema pallidum subspecies pallidum, the etiological agent of syphilis. In this study, the genome sequence of Treponema paraluiscuniculi, strain Cuniculi A, was determined by a combination of several high-throughput sequencing strategies. Whereas the overall size (1,133,390 bp), arrangement, and gene content of the Cuniculi A genome closely resembled those of the T. pallidum genome, the T. paraluiscuniculi genome contained a markedly higher number of pseudogenes and gene fragments (51). In addition to pseudogenes, 33 divergent genes were also found in the T. paraluiscuniculi genome. A set of 32 (out of 84) affected genes encoded proteins of known or predicted function in the Nichols genome. These proteins included virulence factors, gene regulators and components of DNA repair and recombination. The majority (52 or 61.9%) of the Cuniculi A pseudogenes and divergent genes were of unknown function. Our results indicate that T. paraluiscuniculi has evolved from a T. pallidum-like ancestor and adapted to a specialized host-associated niche (rabbits) during loss of infectivity to humans. The genes that are inactivated or altered in T. paraluiscuniculi are candidates for virulence factors important in the infectivity and pathogenesis of T. pallidum subspecies.

Neutral and non-neutral evolution of duplicated genes with gene conversion

Gene conversion is one of the major mutational mechanisms involved in the DNA sequence evolution of duplicated genes. It contributes to create unique patters of DNA polymorphism within species and divergence between species. A typical pattern is so-called concerted evolution, in which the divergence between duplicates is maintained low for a long time because of frequent exchanges of DNA fragments. In addition, gene conversion affects the DNA evolution of duplicates in various ways especially when selection operates. Here, we review theoretical models to understand the evolution of duplicates in both neutral and non-neutral cases. We also explain how these theories contribute to interpreting real polymorphism and divergence data by using some intriguing examples.

Gene conversion in angiosperm genomes with an emphasis on genes duplicated by polyploidization

Angiosperm genomes differ from those of mammals by extensive and recursive polyploidizations. The resulting gene duplication provides opportunities both for genetic innovation, and for concerted evolution. Though most genes may escape conversion by their homologs, concerted evolution of duplicated genes can last for millions of years or longer after their origin. Indeed, paralogous genes on two rice chromosomes duplicated an estimated 60–70 million years ago have experienced gene conversion in the past 400,000 years. Gene conversion preserves similarity of paralogous genes, but appears to accelerate their divergence from orthologous genes in other species. The mutagenic nature of recombination coupled with the buffering effect provided by gene redundancy, may facilitate the evolution of novel alleles that confer functional innovations while insulating biological fitness of affected plants. A mixed evolutionary model, characterized by a primary birth-and-death process and occasional homoeologous recombination and gene conversion, may best explain the evolution of multigene families.

A Bayesian approach for fast and accurate gene tree reconstruction

Recent sequencing and computing advances have enabled phylogenetic analyses to expand to both entire genomes and large clades, thus requiring more efficient and accurate methods designed specifically for the phylogenomic context. Here, we present SPIMAP, an efficient Bayesian method for reconstructing gene trees in the presence of a known species tree. We observe many improvements in reconstruction accuracy, achieved by modeling multiple aspects of evolution, including gene duplication and loss (DL) rates, speciation times, and correlated substitution rate variation across both species and loci. We have implemented and applied this method on two clades of fully sequenced species, 12 Drosophila and 16 fungal genomes as well as simulated phylogenies and find dramatic improvements in reconstruction accuracy as compared with the most popular existing methods, including those that take the species tree into account. We find that reconstruction inaccuracies of traditional phylogenetic methods overestimate the number of DL events by as much as 2-3-fold, whereas our method achieves significantly higher accuracy. We feel that the results and methods presented here will have many important implications for future investigations of gene evolution.

Ectopic gene conversions in the genome of ten hemiascomycete yeast species

We characterized ectopic gene conversions in the genome of ten hemiascomycete yeast species. Of the ten species, three diverged prior to the whole genome duplication (WGD) event present in the yeast lineage and seven diverged after it. We analyzed gene conversions from three separate datasets: paralogs from the three pre-WGD species, paralogs from the seven post-WGD species, and common ohnologs from the seven post-WGD species. Gene conversions have similar lengths and frequency and occur between sequences having similar degrees of divergence, in paralogs from pre- and post-WGD species. However, the sequences of ohnologs are both more divergent and less frequently converted than those of paralogs. This likely reflects the fact that ohnologs are more often found on different chromosomes and are evolving under stronger selective pressures than paralogs. Our results also show that ectopic gene conversions tend to occur more frequently between closely linked genes. They also suggest that the mechanisms responsible for the loss of introns in S. cerevisiae are probably also involved in the gene 3'-end gene conversion bias observed between the paralogs of this species.

Gene conversion and evolution of gene families: an overview

The importance of gene conversion for the evolution of gene families is reviewed. Four problems concerning gene conversion, i.e., concerted evolution, generation of useful variation, deleterious effects, and relation to neofunctionalization, are discussed by surveying reported examples of evolving gene families. Emphasis is given toward understanding interactive effects of gene conversion and natural selection.

New insights into the evolution of metazoan cadherins

Mining newly sequenced genomes of basal metazoan organisms reveals the evolutionary origin of modern protein families. Specific cell-cell adhesion and intracellular communication are key processes in multicellular animals, and members of the cadherin superfamily are essential players in these processes. Mammalian genomes contain over 100 genes belonging to this superfamily. By a combination of tBLASTn and profile hidden Markov model analyses, we made an exhaustive search for cadherins and compiled the cadherin repertoires in key organisms, including Branchiostoma floridae (amphioxus), the sea anemone Nematostella vectensis, and the placozoan Trichoplax adhaerens. Comparative analyses of multiple protein domains within known and novel cadherins enabled us to reconstruct the complex evolution in metazoa of this large superfamily. Five main cadherin branches are represented in the primitive metazoan Trichoplax: classical (CDH), flamingo (CELSR), dachsous (DCHS), FAT, and FAT-like. Classical cadherins, such as E-cadherin, arose from an Urmetazoan cadherin, which progressively lost N-terminal extracellular cadherin repeats, whereas its cytoplasmic domain, which binds the armadillo proteins p120ctn and β-catenin, remained quite conserved from placozoa to man. The origin of protocadherins predates the Bilateria and is likely rooted in an ancestral FAT cadherin. Several but not all protostomians lost protocadherins. The emergence of chordates coincided with a great expansion of the protocadherin repertoire. The evolution of ancient metazoan cadherins points to their unique and crucial roles in multicellular animal life.

Comparative profiling of microRNA expression between neural stem cells and motor neurons in embryonic spinal cord in rat

Neural stem cells' transplantation has been proposed as a future therapy for spinal cord injury. The challenge is how to make proportionally more neural stem cells differentiate into spinal motor neurons. Recent reports reveal that microRNAs play an important role in regulating stem cell self-renewal and differentiation. The aim of this study was to compare the profiling of microRNA expression between neural stem cells and motor neurons and to find candidate targets that direct differentiation of neural stem cells into motor neurons. We performed a parallel isolation and purification of motor neurons and neural stem cells from the same rat embryonic spinal cord sample. With the high-throughput TaqMan low-density array platform, 44 differentially expressed microRNAs were identified (22 specially expressed microRNAs in motor neurons and neural stem cells, respectively). Using bioinformatic methods, clustering, transcriptional regulation and target genes of differential microRNAs were analyzed. Furthermore, miR-126 specially expressed in cultured motor neurons identified by TaqMan low-density array was significantly elevated in choline acetyltransferase-positive neurons differentiated from the neural stem cells. These findings suggest that specially expressed microRNAs may contribute to the directed differentiation of neural stem cells into motor neurons and are potential targets for therapeutic interventions following spinal cord injury.

Regulation of protocadherin gene expression by multiple neuron-restrictive silencer elements scattered in the gene cluster

The clustered protocadherins are a subfamily of neuronal cell adhesion molecules that play an important role in development of the nervous systems in vertebrates. The clustered protocadherin genes exhibit complex expression patterns in the central nervous system. In this study, we have investigated the molecular mechanism underlying neuronal expression of protocadherin genes using the protocadherin gene cluster in fugu as a model. By in silico prediction, we identified multiple neuron-restrictive silencer elements (NRSEs) scattered in the fugu protocadherin cluster and demonstrated that these elements bind specifically to NRSF/REST in vitro and in vivo. By using a transgenic Xenopus approach, we show that these NRSEs regulate neuronal specificity of protocadherin promoters by suppressing their activity in non-neuronal tissues. We provide evidence that protocadherin genes that do not contain an NRSE in their 5' intergenic region are regulated by NRSEs in the regulatory region of their neighboring genes. We also show that protocadherin clusters in other vertebrates such as elephant shark, zebrafish, coelacanth, lizard, mouse and human, contain different sets of multiple NRSEs. Taken together, our data suggest that the neuronal specificity of protocadherin cluster genes in vertebrates is regulated by the NRSE-NRSF/REST system.

Cloning and comparative analyses of the zebrafish Ugt repertoire reveal its evolutionary diversity

UDP-glucuronosyltransferases (Ugts) are a supergene family of phase II drug-metabolizing enzymes that catalyze the conjugation of numerous hydrophobic small molecules with the UDP-glucuronic acid, converting them into hydrophilic molecules. Here, we report the identification and cloning of the complete zebrafish Ugt gene repertoire. We found that the zebrafish genome contains 45 Ugt genes that can be divided into three families: Ugt1, Ugt2, and Ugt5. Both Ugt1 and Ugt2 have two unlinked clusters: a and b. The Ugt1a, Ugt1b, Ugt2a, and Ugt2b clusters each contain variable and constant regions, similar to that of the protocadherin (Pcdh), immunoglobulin (Ig), and T-cell receptor (Tcr) clusters. Cloning the full-length coding sequences confirmed that each of the variable exons is separately spliced to the set of constant exons within each zebrafish Ugt cluster. Comparative analyses showed that both a and b clusters of the zebrafish Ugt1 and Ugt2 genes have orthologs in other teleosts, suggesting that they may be resulted from the "fish-specific" whole-genome duplication event. The Ugt5 genes are a novel family of Ugt genes that exist in teleosts and amphibians. Their entire open reading frames are encoded by single large exons. The zebrafish Ugt1, Ugt2, and Ugt5 genes can generate additional transcript diversity through alternative splicing. Based on phylogenetic analyses, we propose that the ancestral tetrapod and teleost Ugt1 clusters contained multiple Ugt1 paralogs. After speciation, these ancestral Ugt1 clusters underwent lineage-specific gene loss and duplication. The ancestral vertebrate Ugt2 cluster also underwent lineage-specific duplication. The intronless Ugt5 open reading frames may be derived from retrotransposition followed by gene duplication. They have been expanded dramatically in teleosts and have become the most abundant Ugt family in these lineages. These findings have interesting implications regarding the molecular evolution of genes with diversified variable exons in vertebrates.