Evolution of Conserved Non-Coding Sequences Within the Vertebrate Hox Clusters Through the Two-Round Whole Genome Duplications Revealed by Phylogenetic Footprinting Analysis

As a result of two-round whole genome duplications, four or more paralogous Hox clusters exist in vertebrate genomes. The paralogous genes in the Hox clusters show similar expression patterns, implying shared regulatory mechanisms for expression of these genes. Previous studies partly revealed the expression mechanisms of Hox genes. However, cis-regulatory elements that control these paralogous gene expression are still poorly understood. Toward solving this problem, the authors searched conserved non-coding sequences (CNSs), which are candidates of cis-regulatory elements. When comparing orthologous Hox clusters of 19 vertebrate species, 208 intergenic conserved regions were found. The authors then searched for CNSs that were conserved not only between orthologous clusters but also among the four paralogous Hox clusters. The authors found three regions that are conserved among all the four clusters and eight regions that are conserved between intergenic regions of two paralogous Hox clusters. In total, 28 CNSs were identified in the paralogous Hox clusters, and nine of them were newly found in this study. One of these novel regions bears a RARE motif. These CNSs are candidates for gene expression regulatory regions among paralogous Hox clusters. The authors also compared vertebrate CNSs with amphioxus CNSs within the Hox cluster, and found that two CNSs in the HoxA and HoxB clusters retain homology with amphioxus CNSs through the two-round whole genome duplications.     

Comparative genomic analysis reveals evolutionary characteristics and patterns of microRNA clusters in vertebrates

MicroRNAs (miRNAs) are a class of small non-coding RNAs that can play important regulatory roles in many important biological processes. Although clustering patterns of miRNA clusters have been uncovered in animals, the origin and evolution of miRNA clusters in vertebrates are still poorly understood. Here, we performed comparative genomic analyses to construct 51 sets of orthologous miRNA clusters (SOMCs) across seven test vertebrate species, a collection of miRNA clusters from two or more species that are likely to have evolved from a common ancestral miRNA cluster, and used these to systematically examine the evolutionary characteristics and patterns of miRNA clusters in vertebrates. We found that miRNA clusters are continuously generated, and most of them tend to be conserved and maintained in vertebrate genomes, although some adaptive gains and losses of miRNA cluster have occurred during evolution. Furthermore, miRNA clusters appeared relatively early in the evolutionary history might suffer from more complicated adaptive gain-and-loss than those young miRNA clusters. Detailed analysis showed that genomic duplication events of ancestral miRNAs or miRNA clusters are likely to be major driving force and apparently contribute to origin and evolution of miRNA clusters. Comparison of conserved with lineage-specific miRNA clusters revealed that the contribution of duplication events for the formation of miRNA cluster appears to be more important for conserved miRNA clusters than lineage-specific. Our study provides novel insights for further exploring the origins and evolution of miRNA clusters in vertebrates at a genome scale.     

Recognition of unknown conserved alternatively spliced exons

The split structure of most mammalian protein-coding genes allows for the potential to produce multiple different mRNA and protein isoforms from a single gene locus through the process of alternative splicing (AS). We propose a computational approach called UNCOVER based on a pair hidden Markov model to discover conserved coding exonic sequences subject to AS that have so far gone undetected. Applying UNCOVER to orthologous introns of known human and mouse genes predicts skipped exons or retained introns present in both species, while discriminating them from conserved noncoding sequences. The accuracy of the model is evaluated on a curated set of genes with known conserved AS events. The prediction of skipped exons in the approximately 1% of the human genome represented by the ENCODE regions leads to more than 50 new exon candidates. Five novel predicted AS exons were validated by RT-PCR and sequencing analysis of 15 introns with strong UNCOVER predictions and lacking EST evidence. These results imply that a considerable number of conserved exonic sequences and associated isoforms are still completely missing from the current annotation of known genes. UNCOVER also identifies a small number of candidates for conserved intron retention.     

Going nuclear: gene family evolution and vertebrate phylogeny reconciled

Gene duplications have been common throughout vertebrate evolution, introducing paralogy and so complicating phylogenetic inference from nuclear genes. Reconciled trees are one method capable of dealing with paralogy, using the relationship between a gene phylogeny and the phylogeny of the organisms containing those genes to identify gene duplication events. This allows us to infer phylogenies from gene families containing both orthologous and paralogous copies. Vertebrate phylogeny is well understood from morphological and palaeontological data, but studies using mitochondrial sequence data have failed to reproduce this classical view. Reconciled tree analysis of a database of 118 vertebrate gene families supports a largely classical vertebrate phylogeny.     

Expression of three <Emphasis Type="Italic">spalt </Emphasis>(<Emphasis Type="Italic">sal</Emphasis>) gene homologues in zebrafish embryos

Three homologues of the Drosophilaregion-specific homeotic gene spalt (sal) have been isolated in zebrafish, sall1a, sall1b and sall3. Phylogenetic analysis of these genes against known salDNA sequences showed zebrafish sall1aand sall1b to be orthologous to other vertebrate sal-1 genes and zebrafish sall3to be orthologous to other vertebrate sal-3 genes, except Xenopus sall3. Phylogenetic reconstruction suggests that zebrafish sall1a and sall1bresulted from a gene duplication event occurring prior to the divergence of the ray-finned and lobe-finned fish lineages. Analysis of the expression pattern of the zebrafish sal genes shows that sall1a and sall3 share expression domains with both orthologous and non-orthologous vertebrate sal genes. Both are expressed in various regions of the CNS, including in primary motor neurons. Outside of the CNS, sall1a expression is observed in the otic vesicle (ear), heart and in a discrete region of the pronephric ducts. These analyses indicate that orthologies between zebrafish sal genes and other vertebrate sal genes do not imply equivalence of expression pattern and, therefore, that biological functions are not entirely conserved. However we suggest that, like other vertebrate sal genes, zebrafish sal genes have a role in neural development. Also, expression of zebrafish sall1a in the otic vesicle, heart sac and the pronephric ducts of zebrafish embryos is possibly consistent with some of the abnormalities seen in Sall1-deficient mice and in Townes-Brocks Syndrome, a human disorder which is caused by mutations in the human spalt gene SALL1.     

The duplicated B-class MADS-box genes display dualistic characters in orchid floral organ identity and growth

Orchidaceae are an excellent model to examine perianth development because of their sophisticated floral architecture. In this study, we identified 24 APETALA3 (AP3)-like and 13 PISTILLA (PI)-like genes from 11 species of orchids and characterized them into four AP3- and two PI-duplicated homologs. The first duplication event in AP3 homologs occurring in the early evolutionary history of the Orchidaceae gave rise to AP3A and AP3B clades. Further duplication events resulted in four subclades, namely AP3A1, AP3A2, AP3B1 and AP3B2, during the evolution of Orchidaceae. The AP3 paralogous genes were expressed throughout inflorescence and floral bud development. From the in situ hybridization results, we noticed that the transition timings from ubiquitous to constrained expression in floral organs for both clades are different. The transition point of expression of the AP3A clade (clades 3 and 4) was at the late floral organ primordia stage. In contrast, that for the AP3B clade (clades 1 and 2) was not observed until the late inflorescence and floral bud stages. In addition, the AP3 orthologous genes revealed diverse expression patterns in various species of orchids, whereas the PI homologs were uniformly expressed in all floral whorls. AP3A2 orthologs play a noticeable role in lip formation because of their exclusive expression in the lip. Further evidence comes from the ectopic expression of AP3A2 detected in the lip-like petals extending from the lip in four sets of peloric mutants. Finally, a Homeotic Orchid Tepal (HOT) model is proposed, in which dualistic characters of duplicated B-class MADS-box genes are involved in orchid perianth development and growth.     

蒺藜苜蓿花器官特异基因的表达分析

蒺藜苜蓿(Medicago truncatula G.)花器官特异表达基因是参与其花器官形成与发育的重要基因。筛选蒺藜苜蓿的花器官特异表达基因, 寻找这类基因在其他模式植物中的直系同源基因, 并将其表达模式在不同植物间进行比较, 有利于深入的理解这类基因在蒺藜苜蓿花器官发育中的功能。根据蒺藜苜蓿表达谱, 并以其PISTILLATA(PI)基因为模板, 文章筛选了97个蒺藜苜蓿花器官特异表达基因(Ratio≥10, 且Z≥7.9)。通过同源比对, 确定了这类基因在拟南芥(Arabidopsis thaliana L.)、大豆(Glycine max L.)、百脉根(Lotus japonicus L.)和水稻(Oryza sativa L.)中的直系同源基因。对这类基因在5种植物中的表达量、表达部位和功能进行比较, 发现进化关系较近的植物, 直系同源基因的表达变异较小, 而进化关系较远的植物, 直系同源基因的表达变异较大。进一步对表达分化较大的直系同源基因进行启动子分析, 发现不同植物中直系同源基因表达模式的变化与启动子中调控元件的特性有关。     

蒺藜苜蓿花器官特异基因的表达分析

蒺藜苜蓿(Medicago truncatula G)花器官特异表达基因是参与其花器官形成与发育的重要基因。筛选蒺藜苜蓿的花器官特异表达基因,寻找这类基因在其他模式植物中的直系同源基因,并将其表达模式在不同植物间进行比较,有利于深入的理解这类基因在蒺藜苜蓿花器官发育中的功能。根据蒺藜苜蓿表达谱,并以其PISTILLAZA(PI)基因为模板,文章筛选了97个蒺藜苜蓿花器官特异表达基因(Ratio≥10,且Z≥7.9).通过同源比对,确定了这类基因在拟南芥(Arabidopsis thaliana L.)、大豆(Glycinemax L.)、百脉根(Lotusjaponicus L.)和水稻(Oryzasativa L.)中的直系同源基因。对这类基因在5种植物中的表达量、表达部位和功能进行比较,发现进化关系较近的植物,直系同源基因的表达变异较小,而进化关系较远的植物,直系同源基因的表达变异较大。进一步对表达分化较大的直系同源基因进行启动子分析,发现不同植物中直系同源基因表达模式的变化与启动子中调控元件的特性有关。     

Perm-seq: Mapping Protein-DNA Interactions in Segmental Duplication and Highly Repetitive Regions of Genomes with Prior-Enhanced Read Mapping

Segmental duplications and other highly repetitive regions of genomes contribute significantly to cells’ regulatory programs. Advancements in next generation sequencing enabled genome-wide profiling of protein-DNA interactions by chromatin immunoprecipitation followed by high throughput sequencing (ChIP-seq). However, interactions in highly repetitive regions of genomes have proven difficult to map since short reads of 50–100 base pairs (bps) from these regions map to multiple locations in reference genomes. Standard analytical methods discard such multi-mapping reads and the few that can accommodate them are prone to large false positive and negative rates. We developed Perm-seq, a prior-enhanced read allocation method for ChIP-seq experiments, that can allocate multi-mapping reads in highly repetitive regions of the genomes with high accuracy. We comprehensively evaluated Perm-seq, and found that our prior-enhanced approach significantly improves multi-read allocation accuracy over approaches that do not utilize additional data types. The statistical formalism underlying our approach facilitates supervising of multi-read allocation with a variety of data sources including histone ChIP-seq. We applied Perm-seq to 64 ENCODE ChIP-seq datasets from GM12878 and K562 cells and identified many novel protein-DNA interactions in segmental duplication regions. Our analysis reveals that although the protein-DNA interactions sites are evolutionarily less conserved in repetitive regions, they share the overall sequence characteristics of the protein-DNA interactions in non-repetitive regions.     

Adaptive evolution of the uncoupling protein 1 gene contributed to the acquisition of novel nonshivering thermogenesis in ancestral eutherian mammals

Homeotherms possess various physiological mechanisms to maintain their body temperature, thus allowing them to adapt to various environments. Under cold conditions, most eutherian mammals upregulate heat production in brown adipose tissue (BAT), and uncoupling protein (UCP) 1 is an essential factor in BAT thermogenesis. The evolutionary origin of UCP1 was believed to have been a specific event occurring in eutherian lineages. Recently, however, the UCP1 ortholog was found in fishes, which uncovers a more ancient origin of this gene than previously believed. Here we investigate the evolutionary process of UCP1 by comparative genomic approach. We found that UCP1 evolved rapidly by positive Darwinian selection in the common ancestor of eutherians, although this gene arose in the ancestral vertebrate, since the orthologous genes were shared among most of the vertebrate species. Adaptive evolution occurred after the divergence between eutherians and marsupials, which is consistent with the fact that BAT has been found only in eutherians. Our findings indicate that positive Darwinian selection acted on UCP1 contributed to the acquisition of an efficient mechanism for body temperature regulation in primitive eutherians. Phylogenetic reconstruction of UCP1 with two paralogs (UCP2 and UCP3) among vertebrate species revealed that the gene duplication events which produced these three genes occurred in the common ancestor of vertebrates much earlier than the emergence of eutherians. Thus, our data demonstrate that novel gene function can evolve without de novo gene duplication event.     

Census of orthologous genes and self-organizing maps of biologically relevant transcriptional patterns in chickens (Gallus gallus)

The launch of large-scale chicken expressed sequence tags (EST) projects has placed the chicken in the lead for the number of EST sequences in agriculturally important animals. More than 451,000 chicken ESTs derived from over 158 libraries have been deposited in the NCBI dbEST database as of December 2003. But how many genes these ESTs represent and how they are expressed in different chicken tissues/organs remain undetermined. In the present research, we developed a human gene-based strategy for census of chicken orthologous genes and identification of their expression patterns. Among 34,157 human coding genes used in the study, BLAST analysis revealed that 11,066 genes provisionally matched 248,628 chicken ESTs. Based on the average EST abundance of the orthologous genes, the current public repository of chicken ESTs could represent 20,000 provisional genes. Analysis of gene expression in 14 single tissues/organs showed that approximately 15% of genes were expressed exclusively in single tissue/organ whereas the remaining 85% of genes were co-expressed in two or more tissues/organs. A majority (91.15%) of genes expressed in chicken embryos were also expressed at post-hatch stages, indicating that most genes activated in chicken embryos could serve housekeeping functions. Self-organizing maps (SOM) analysis organized 8807 provisional genes in selected chicken tissues into 98 clusters with each cluster being indicative of common regulatory factors and pathways. A total of 969 provisional orthologous genes were identified as preferentially expressed genes (PEGs) in various chicken tissues/organs (LOD>3.0). No doubt, the present study on gene expression patterns will provide insight into dynamics of metabolic pathways and tissue/organ programming and reprogramming in chickens.     

Genome-wide analysis of WOX gene family in rice, sorghum, maize, Arabidopsis and poplar

WUSCHEL-related homeobox(WOX)genes form a large gene family specifically expressed in plants.They are known to play important roles in regulating the development of plant tissues and organs by determining cell fate.Recent available whole genome sequences allow us to do more comprehensive phylogenetic analysis of the WOX genes in plants.In the present study,we identified 11 and 21 WOXs from sorghum(Sorghum bicolor)and maize(Zea mays),respectively.The 72 WOX genes from rice(Oryza sativa),sorghum,maize,Arabidopsis(Arabidopsis thaliana)and poplar(Populus trichocarpa)were grouped into three well supported clades with nine subgroups according to the amino acid sequences of their homodomains.Their phylogenetic relationship was also supported by the observation of the motifs outside the homodomain.We observed the variation of duplication events among the nine sub-groups between monocots and eudicots,for instance,more gene duplication events of WOXs within subgroup A for monocots,while,less for dicots in this subgroup.Furthermore,we observed the conserved intron/exon structural patterns of WOX genes in rice,sorghum and Arabidopsis.In addition,WUS(Wuschel)-box and EAR(the ERF-associated amphiphilic repression)-like motif were observed to be conserved among several WOX subgroups in these five plants.Comparative analysis of expression patterns of WOX genes in rice and Arabidopsis suggest that the WOX genes play conserved and various roles in plants.This work provides insights into the evolution of the WOX gene family and is useful for future research.     

Relationships among msx gene structure and function in zebrafish and other vertebrates

The zebrafish genome contains at least five msx homeobox genes, msxA, msxB, msxC, msxD, and the newly isolated msxE. Although these genes share structural features common to all Msx genes, phylogenetic analyses of protein sequences indicate that the msx genes from zebrafish are not orthologous to the Msx1 and Msx2 genes of mammals, birds, and amphibians. The zebrafish msxB and msxC are more closely related to each other and to the mouse Msx3. Similarly, although the combinatorial expression of the zebrafish msx genes in the embryonic dorsal neuroectoderm, visceral arches, fins, and sensory organs suggests functional similarities with the Msx genes of other vertebrates, differences in the expression patterns preclude precise assignment of orthological relationships. Distinct duplication events may have given rise to the msx genes of modern fish and other vertebrate lineages whereas many aspects of msx gene functions during embryonic development have been preserved.      

Gene Duplication and Fragment Recombination Drive Functional Diversification of a Superfamily of Cytoplasmic Effectors in Phytophthora sojae

Phytophthora and other oomycetes secrete a large number of putative host cytoplasmic effectors with conserved FLAK motifs following signal peptides, termed crinkling and necrosis inducing proteins (CRN), or Crinkler. Here, we first investigated the evolutionary patterns and mechanisms of CRN effectors in Phytophthora sojae and compared them to two other Phytophthora species. The genes encoding CRN effectors could be divided into 45 orthologous gene groups (OGG), and most OGGs unequally distributed in the three species, in which each underwent large number of gene gains or losses, indicating that the CRN genes expanded after species evolution in Phytophthora and evolved through pathoadaptation. The 134 expanded genes in P. sojae encoded family proteins including 82 functional genes and expressed at higher levels while the other 68 genes encoding orphan proteins were less expressed and contained 50 pseudogenes. Furthermore, we demonstrated that most expanded genes underwent gene duplication or/and fragment recombination. Three different mechanisms that drove gene duplication or recombination were identified. Finally, the expanded CRN effectors exhibited varying pathogenic functions, including induction of programmed cell death (PCD) and suppression of PCD through PAMP-triggered immunity or/and effector-triggered immunity. Overall, these results suggest that gene duplication and fragment recombination may be two mechanisms that drive the expansion and neofunctionalization of the CRN family in P. sojae, which aids in understanding the roles of CRN effectors within each oomycete pathogen.