Evolution of coding and non-coding genes in HOX clusters of a marsupial

ABSTRACT: BACKGROUND: The HOX gene clusters are thought to be highly conserved amongst mammals and othervertebrates, but the long non-coding RNAs have only been studied in detail in human andmouse. The sequencing of the kangaroo genome provides an opportunity to use comparativeanalyses to compare the HOX clusters of a mammal with a distinct body plan to those ofother mammals. RESULTS: Here we report a comparative analysis of HOX gene clusters between an Australian marsupialof the kangaroo family and the eutherians. There was a strikingly high level of conservationof HOX gene sequence and structure and non-protein coding genes including the microRNAsmiRNA-196a, miRNA-196b, miRNA-10a and miRNA-10b and the long non-coding RNAsHOTAIR, HOTAIRM1 and HOXA11AS that play critical roles in regulating gene expressionand controlling development. By microRNA deep sequencing and comparative genomicanalyses, two conserved microRNAs (miR-10a and miR-10b) were identified and one newcandidate microRNA with typical hairpin precursor structure that is expressed in bothfibroblasts and testes was found. The prediction of microRNA target analysis showed thatseveral known microRNA targets, such as mir-10, mir-414 and mir-464, were found in thetammar HOX clusters. In addition, several novel and putative miRNAs were identified thatoriginated from elsewhere in the tammar genome and that target the tammar HOXB andHOXD clusters. CONCLUSIONS: This study confirms that the emergence of known long non-coding RNAs in the HOXclusters clearly predate the marsupial-eutherian divergence 160 Ma ago. It also identified anew potentially functional microRNA as well as conserved miRNAs. These non-codingRNAs may participate in the regulation of HOX genes to influence the body plan of thismarsupial.     

Phylogenetic analyses of Zostera species based on rbcL and matK nucleotide sequences: implications for the origin and diversification of seagrasses in Japanese waters

Seagrasses are composed of four families belonging to angiosperms and they are thought to become adaptive to aquatic life independently. Zosteraceae is one such family and because of the relatively high species diversity around Japan and Korea coast areas, the family might have arisen therefrom. To elucidate the origin and evolution of Zosteraceae which consists of three genera, Phyllospadix, Zostera, and Heterozostera, 2.8 kb nucleotide sequences of rbcL and matK genes in the chloroplast genome were examined for various species, including cosmopolitan Z. marina and endemic Z. caulescens. The phylogenetic analysis reveals the following three features. First, based on the synonymous nucleotide substitution rate of the rice chloroplast genome, we estimated the divergence times between Zosteraceae and its closest relative, Potamogetonaceae, and between different genera, Zostera and Phyllospadix, as approximately 100 million years (myr) and 36 myr, respectively, suggesting that Zosteraceae emerged somewhere in the period from 36 myr ago to 100 myr ago. Second, two subgenera of Zostera, Zostera and Zosterella, exhibit their reciprocal monophyly and appear to have differentiated from each other approximately 33 myr ago. However, the third genus Heterozostera branched off only 5 myr ago from the stem lineage leading to Zosterella and this seems too recent in comparison with the ancient divergence of the two subgenera. Third, we estimated the most recent common ancestor of subgenus Zostera as 6 myr. In Z. marina four haplotypes were found in the sample and have diversified in the past 1.5 myr. One haplotype is shared by both sides of the Japan Archipelago and its closely related haplotypes occur also in eastern Pacific Ocean. Based on these phylogeographic analyses, we propose a provisional age related classification of Zosteraceae to argue the origin and evolution.     

Molecular cloning of retrovirus-like genes present in multiple copies in the Syrian hamster genome.

Endogenous retrovirus-like sequences homologous to intracisternal type-A particle (IAP) genes, which are present in the inbred mouse (Mus musculus) genome, were cloned from a Syrian hamster gene library. A typical hamster IAP gene was 7 kb long and segments homologous to long terminal repeat (IAP) sequences present in Mus musculus IAP genes were located at both ends of the gene. Contrary to the pattern found in the Mus musculus IAP genes, the organization of the cloned hamster IAP genes was not markedly polymorphic and deletion was not observed among these cloned genes. A sequence about 0.8 kb long and located close to the 3' end of the hamster IAP gene was well conserved in both IAP gene families, although they showed less overall homology with one another. The reiteration frequency of the hamster IAP genes was calculated to be 950 copies per haploid genome. Since such IAP genes with the above properties were not found in the genome of the Chinese hamster, whose progenitors diverged from those of the Syrian hamster about 7.5 Myr ago, the integration of a huge number of Syrian hamster IAP genes must have occurred subsequent to such divergence.     

The ancient source of a distinct gene family encoding proteins featuring RING and C(3)H zinc-finger motifs with abundant expression in developing brain and nervous system

Intronless genes can arise by germline retrotransposition of a cDNA originating as mRNA from an intron-containing source gene. Previously, we described several members of a family of intronless mammalian genes encoding a novel class of zinc-finger proteins, including one that shows imprinted expression and one that escapes X-inactivation. We report here the identification and characterization of the Makorin ring finger protein 1 gene (MKRN1), a highly transcribed, intron-containing source for this family of genes. Phylogenetic analyses clearly indicate that the MKRN1 gene is the ancestral founder of this gene family. We have identified MKRN1 orthologs from human, mouse, wallaby, chicken, fruitfly, and nematode, underscoring the age and conservation of this gene. The MKRN gene family encodes putative ribonucleoproteins with a distinctive array of zinc-finger motifs, including two to four C(3)H zinc-fingers, an unusual Cys/His arrangement that may represent a novel zinc-finger structure, and a highly conserved RING zinc-finger. To date, we have identified nine MKRN family loci distributed throughout the human genome. The human and mouse MKRN1 loci map to a conserved syntenic group near the T-cell receptor beta cluster (TCRB) in chromosome 7q34-q35 and chromosome 6A, respectively. MKRN1 is widely transcribed in mammals, with high levels in murine embryonic nervous system and adult testis. The ancient origin of MKRN1, high degree of conservation, and expression pattern suggest important developmental and functional roles for this gene and its expressed family members.     

The Evolution of Mammalian Olfactory Receptor Genes

We performed a comparative study of four subfamilies of olfactory receptor genes first identified in the dog to assess changes in the gene family during mammalian evolution, and to begin linking the dog genetic map to that of humans. The human subfamilies were localized to chromosomes 7, 11, and 19. The two subfamilies that were tightly linked in the dog genome were also tightly linked in the human genome. The four subfamilies were compared in human (primate), horse (perissodactyl), and a variety of artiodactyls and carnivores. Some changes in gene number were detected, but overall subfamily size appeared to have been established before the divergence of these mammals 60-100 million years ago.     

Computer identification of snoRNA genes using a Mammalian Orthologous Intron Database

Based on comparative genomics, we created a bioinformatic package for computer prediction of small nucleolar RNA (snoRNA) genes in mammalian introns. The core of our approach was the use of the Mammalian Orthologous Intron Database (MOID), which contains all known introns within the human, mouse and rat genomes. Introns from orthologous genes from these three species, that have the same position relative to the reading frame, are grouped in a special orthologous intron table. Our program SNO.pl searches for conserved snoRNA motifs within MOID and reports all cases when characteristic snoRNA-like structures are present in all three orthologous introns of human, mouse and rat sequences. Here we report an example of the SNO.pl usage for searching a particular pattern of conserved C/D-box snoRNA motifs (canonical C- and D-boxes and the 6 nt long terminal stem). In this computer analysis, we detected 57 triplets of snoRNA-like structures in three mammals. Among them were 15 triplets that represented known C/D-box snoRNA genes. Six triplets represented snoRNA genes that had only been partially characterized in the mouse genome. One case represented a novel snoRNA gene, and another three cases, putative snoRNAs. Our programs are publicly available and can be easily adapted and/or modified for searching any conserved motifs within mammalian introns.     

Identification and developmental expression of the ets gene family in the sea urchin (Strongylocentrotus purpuratus)

A systematic search in the available scaffolds of the Strongylocentrotus purpuratus genome has revealed that this sea urchin has 11 members of the ets gene family. A phylogenetic analysis of these genes showed that almost all vertebrate ets subfamilies, with the exception of one, so far found only in mammals, are each represented by one orthologous sea urchin gene. The temporal and spatial expression of the identified ETS factors was also analyzed during embryogenesis. Five ets genes (Sp-Ets1/2, Sp-Tel, Sp-Pea, Sp-Ets4, Sp-Erf) are also maternally expressed. Three genes (Sp-Elk, Sp-Elf, Sp-Erf) are ubiquitously expressed during embryogenesis, while two others (Sp-Gabp, Sp-Pu.1) are not transcribed until late larval stages. Remarkably, five of the nine sea urchin ets genes expressed during embryogenesis are exclusively (Sp-Ets1/2, Sp-Erg, Sp-Ese) or additionally (Sp-Tel, Sp-Pea) expressed in mesenchyme cells and/or their progenitors. Functional analysis of Sp-Ets1/2 has previously demonstrated an essential role of this gene in the specification of the skeletogenic mesenchyme lineage. The dynamic, and in some cases overlapping and/or unique, developmental expression pattern of the latter five genes suggests a complex, non-redundant function for ETS factors in sea urchin mesenchyme formation and differentiation.     

A shared promoter region suggests a common ancestor for the human VCX/Y, SPANX, and CSAG gene families and the murine CYPT family

Many testis-specific genes from the sex chromosomes are subject to rapid evolution, which can make it difficult to identify murine genes in the human genome. The murine CYPT gene family includes 15 members, but orthologs were undetectable in the human genome. However, using refined homology search, sequences corresponding to the shared promoter region of the CYPT family were identified at 39 loci. Most loci were located immediately upstream of genes belonging to the VCX/Y, SPANX, or CSAG gene families. Sequence comparison of the loci revealed a conserved CYPT promoter-like (CPL) element featuring TATA and CCAAT boxes. The expression of members of the three families harboring the CPL resembled the murine expression of the CYPT family, with weak expression in late pachytene spermatocytes and predominant expression in spermatids, but some genes were also weakly expressed in somatic cells and in other germ cell types. The genomic regions harboring the gene families were rich in direct and inverted segmental duplications (SD), which may facilitate gene conversion and rapid evolution. The conserved CPL and the common expression profiles suggest that the human VCX/Y, SPANX, and CSAG2 gene families together with the murine SPANX gene and the CYPT family may share a common ancestor. Finally, we present evidence that VCX/Y and SPANX may be paralogs with a similar protein structure consisting of C terminal acidic repeats of variable lengths.     

Generation and characterization of a transgenic mouse with a functional human TSPY

To generate an animal model that is suitable for the analysis of regulation and expression of human testis-specific protein, Y-encoded TSPY, a transgenic mouse line, TgTSPY9, harboring a complete structural human TSPY gene was generated. Fluorescence in situ hybridization and Southern analyses show that approximately 50 copies of the human TSPY transgene are integrated at a single chromosomal site that maps to the distal long arm of the Y chromosome. The transgene is correctly transcribed and spliced according to the human pattern and is mainly expressed in testicular tissue, with spermatogonia and early primary spermatocytes (leptotene and zygotene) as expressing germ cells. TSPY transgenic mice are phenotypically normal, and spermatogenesis is neither impaired nor enhanced by the human transgene. The present study shows that a human TSPY gene integrated into the mouse genome follows the human expression pattern although murine tspy had lost its function in rodent evolution millions of years ago.     

Molecular analysis of the Aedes aegypti carboxypeptidase gene family

To gain a better understanding of coordinate regulation of protease gene expression in the mosquito midgut, we undertook a comprehensive molecular study of digestive carboxypeptidases in Aedes aegypti. Through a combination of cDNA cloning using degenerate PCR primers, and database mining of the recently completed A. aegypti genome, we cloned and characterized 18 A. aegypti carboxypeptidase genes. Bioinformatic analysis revealed that 11 of these genes belong to the carboxypeptidase A family (AaCPA-I through AaCPA-XI), and seven to the carboxypeptidase B gene family (AaCPB-I through AaCPB-VII). Phylogenetic analysis of 32 mosquito carboxypeptidases from five different species indicated that most of the sequence divergence in the carboxypeptidase gene family occurred prior to the separation of Aedes and Anopheles mosquito lineages. Unlike the CPA genes that are scattered throughout the A. aegypti genome, six of seven CPB genes were found to be located within a single 120 kb genome contig, suggesting that they most likely arose from multiple gene duplication events. Quantitative expression analysis revealed that 11 of the A. aegypti carboxypeptidase genes were induced up to 40-fold in the midgut in response to blood meal feeding, with peak expression times ranging from 3 to 36 h post-feeding depending on the gene.