PipMaker: a World Wide Web server for genomic sequence alignments

Bacteriophage Taxonomy: A Continually Evolving Discipline

While taxonomy is an often underappreciated branch of science, it serves very important roles. Bacteriophage taxonomy has evolved from a discipline based mainly on morphology, characterized by the work of David Bradley and Hans-Wolfgang Ackermann, to the sequence-based approach that is taken today. The Bacterial Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV) takes a holistic approach to classifying prokaryote viruses by measuring overall DNA and protein similarity and phylogeny before making decisions about the taxonomic position of a new virus. The huge number of complete genomes being deposited with the National Center for Biotechnology Information (NCBI) and other public databases has resulted in a reassessment of the taxonomy of many viruses, and the future will see the introduction of new viral families and higher orders.

An evolutionarily ancient mechanism for regulation of hemoglobin expression in vertebrate red cells

The oxygen transport function of hemoglobin (HB) is thought to have arisen ∼500 million years ago, roughly coinciding with the divergence between jawless (Agnatha) and jawed (Gnathostomata) vertebrates. Intriguingly, extant HBs of jawless and jawed vertebrates were shown to have evolved twice, and independently, from different ancestral globin proteins. This raises the question of whether erythroid-specific expression of HB also evolved twice independently. In all jawed vertebrates studied to date, one of the HB gene clusters is linked to the widely expressed NPRL3 gene. Here we show that the nprl3-linked hb locus of a jawless vertebrate, the river lamprey (Lampetra fluviatilis), shares a range of structural and functional properties with the equivalent jawed vertebrate HB locus. Functional analysis demonstrates that an erythroid-specific enhancer is located in intron 7 of lamprey nprl3, which corresponds to the NPRL3 intron 7 MCS-R1 enhancer of jawed vertebrates. Collectively, our findings signify the presence of an nprl3-linked multiglobin gene locus, which contains a remote enhancer that drives globin expression in erythroid cells, before the divergence of jawless and jawed vertebrates. Different globin genes from this ancestral cluster evolved in the current NPRL3-linked HB genes in jawless and jawed vertebrates. This provides an explanation of the enigma of how, in different species, globin genes linked to the same adjacent gene could undergo convergent evolution.

Natural selection supports escape from concerted evolution of a recently duplicated CEACAM1 paralog in the ruminant CEA gene family

Concerted evolution is often observed in multigene families such as the CEA gene family. As a result, sequence similarity of paralogous genes is significantly higher than expected from their evolutionary distance. Gene conversion, a “copy paste” DNA repair mechanism that transfers sequences from one gene to another and homologous recombination are drivers of concerted evolution. Nevertheless, some gene family members escape concerted evolution and acquire sufficient sequence differences that orthologous genes can be assigned in descendant species. Reasons why some gene family members can escape while others are captured by concerted evolution are poorly understood. By analyzing the entire CEA gene family in cattle (Bos taurus) we identified a member (CEACAM32) that was created by gene duplication and cooption of a unique transmembrane domain exon in the most recent ancestor of ruminants. CEACAM32 shows a unique, testis-specific expression pattern. Phylogenetic analysis indicated that CEACAM32 is not involved in concerted evolution of CEACAM1 paralogs in ruminants. However, analysis of gene conversion events revealed that CEACAM32 is subject to gene conversion but remarkably, these events are found in the leader exon and intron sequences but not in exons coding for the Ig-like domains. These findings suggest that natural selection hinders gene conversion affecting protein sequences of the mature protein and thereby support escape of CEACAM32 from concerted evolution.

Nested retrotransposition in the East Asian mouse genome causes the classical nonagouti mutation

Black coat color (nonagouti) is a widespread classical mutation in laboratory mouse strains. The intronic insertion of endogenous retrovirus VL30 in the nonagouti (a) allele of agouti gene was previously reported as the cause of the nonagouti phenotype. Here, we report agouti mouse strains from East Asia that carry the VL30 insertion, indicating that VL30 alone does not cause the nonagouti phenotype. We find that a rare type of endogenous retrovirus, β4, was integrated into the VL30 region at the a allele through nested retrotransposition, causing abnormal splicing. Targeted complete deletion of the β4 element restores agouti gene expression and agouti coat color, whereas deletion of β4 except for a single long terminal repeat results in black-and-tan coat color. Phylogenetic analyses show that the a allele and the β4 retrovirus originated from an East Asian mouse lineage most likely related to Japanese fancy mice. These findings reveal the causal mechanism and historic origin of the classical nonagouti mutation.

Bacteriophage Taxonomy: An Evolving Discipline

While taxonomy is an often-unappreciated branch of science it serves very important roles. Bacteriophage taxonomy has evolved from a mainly morphology-based discipline, characterized by the work of David Bradley and Hans-Wolfgang Ackermann, to the holistic approach that is taken today. The Bacterial and Archaeal Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV) takes a comprehensive approach to classifying prokaryote viruses measuring overall DNA and protein identity and phylogeny before making decisions about the taxonomic position of a new virus. The huge number of complete genomes being deposited with NCBI and other public databases has resulted in a reassessment of the taxonomy of many viruses, and the future will see the introduction of new viral families and higher orders.

Identification of Candidate Genes for Reactivity in Guzerat (Bos indicus) Cattle: A Genome-Wide Association Study

Temperament is fundamental to animal production due to its direct influence on the animal-herdsman relationship. When compared to calm animals, the aggressive, anxious or fearful ones exhibit less weight gain, lower reproductive efficiency, decreased milk production and higher herd maintenance costs, all of which contribute to reduced profits. However, temperament is a trait that is complex and difficult to assess. Recently, a new quantitative system, REATEST^®, for assessing reactivity, a phenotype of temperament, was developed. Herein, we describe the results of a Genome-wide association study for reactivity, assessed using REATEST^® with a sample of 754 females from five dual-purpose (milk and meat production) Guzerat (Bos indicus) herds. Genotyping was performed using a 50k SNP chip and a two-step mixed model approach (Grammar-Gamma) with a one-by-one marker regression was used to identify QTLs. QTLs for reactivity were identified on chromosomes BTA1, BTA5, BTA14, and BTA25. Five intronic and two intergenic markers were significantly associated with reactivity. POU1F1, DRD3, VWA3A, ZBTB20, EPHA6, SNRPF and NTN4 were identified as candidate genes. Previous QTL reports for temperament traits, covering areas surrounding the SNPs/genes identified here, further corroborate these associations. The seven genes identified in the present study explain 20.5% of reactivity variance and give a better understanding of temperament biology.

An Analysis of Adenovirus Genomes Using Whole Genome Software Tools

The evolution of sequencing technology has lead to an enormous increase in the number of genomes that have been sequenced. This is especially true in the field of virus genomics. In order to extract meaningful biological information from these genomes, whole genome data mining software tools must be utilized. Hundreds of tools have been developed to analyze biological sequence data. However, only some of these tools are user-friendly to biologists. Several of these tools that have been successfully used to analyze adenovirus genomes are described here. These include Artemis, EMBOSS, pDRAW, zPicture, CoreGenes, GeneOrder, and PipMaker. These tools provide functionalities such as visualization, restriction enzyme analysis, alignment, and proteome comparisons that are extremely useful in the bioinformatics analysis of adenovirus genomes.

SOX10 regulates expression of the SH3-domain kinase binding protein 1 (Sh3kbp1) locus in Schwann cells via an alternative promoter

The transcription factor SOX10 has essential roles in neural crest-derived cell populations, including myelinating Schwann cells-specialized glial cells responsible for ensheathing axons in the peripheral nervous system. Importantly, SOX10 directly regulates the expression of genes essential for proper myelin function. To date, only a handful of SOX10 target loci have been characterized in Schwann cells. Addressing this lack of knowledge will provide a better understanding of Schwann cell biology and candidate loci for relevant diseases such as demyelinating peripheral neuropathies. We have identified a highly-conserved SOX10 binding site within an alternative promoter at the SH3-domain kinase binding protein 1 (Sh3kbp1) locus. The genomic segment identified at Sh3kbp1 binds to SOX10 and displays strong promoter activity in Schwann cells in vitro and in vivo. Mutation of the SOX10 binding site ablates promoter activity, and ectopic expression of SOX10 in SOX10-negative cells promotes the expression of endogenous Sh3kbp1. Combined, these data reveal Sh3kbp1 as a novel target of SOX10 and raise important questions regarding the function of SH3KBP1 isoforms in Schwann cells.

Sex determination in the first-described sexual fungus

The original report of sex in fungi dates 2 centuries ago to the species Syzygites megalocarpus (Mucoromycotina). The organism was subsequently used in 1904 to represent self-fertile homothallic species when the concepts of heterothallism and homothallism were developed for the fungal kingdom. In this study, two putative sex/MAT loci were identified in individual strains of S. megalocarpus, accounting for its homothallic behavior. The strains encode both of the high-mobility-group domain-containing proteins, SexM and SexP, flanked by RNA helicase and glutathione oxidoreductase genes that are found adjacent to the mating-type loci in other Mucoromycotina species. The presence of pseudogenes and the arrangement of genes suggest that the origin of homothallism in this species is from a heterothallic relative, obtained via a chromosomal rearrangement to switch two alleles into two separated loci within a single genetic background. Similar events have given rise to homothallic species from heterothallic species in ascomycete fungi, demonstrating that conserved forces shape the evolution of sex determination and speciation in highly diverged fungi.

A rare myelin protein zero (MPZ) variant alters enhancer activity in vitro and in vivo

Background

Myelin protein zero (MPZ) is a critical structural component of myelin in the peripheral nervous system. The MPZ gene is regulated, in part, by the transcription factors SOX10 and EGR2. Mutations in MPZ, SOX10, and EGR2 have been implicated in demyelinating peripheral neuropathies, suggesting that components of this transcriptional network are candidates for harboring disease-causing mutations (or otherwise functional variants) that affect MPZ expression.

Methodology

We utilized a combination of multi-species sequence comparisons, transcription factor-binding site predictions, targeted human DNA re-sequencing, and in vitro and in vivo enhancer assays to study human non-coding MPZ variants.

Principal Findings

Our efforts revealed a variant within the first intron of MPZ that resides within a previously described SOX10 binding site is associated with decreased enhancer activity, and alters binding of nuclear proteins. Additionally, the genomic segment harboring this variant directs tissue-relevant reporter gene expression in zebrafish.

Conclusions

This is the first reported MPZ variant within a cis-acting transcriptional regulatory element. While we were unable to implicate this variant in disease onset, our data suggests that similar non-coding sequences should be screened for mutations in patients with neurological disease. Furthermore, our multi-faceted approach for examining the functional significance of non-coding variants can be readily generalized to study other loci important for myelin structure and function.

Comparative genomic analysis of the false killer whale (Pseudorca crassidens) LMBR1 locus

The sequencing and comparative genomic analysis of LMBR1 loci in mammals or other species, including human, would be very important in understanding evolutionary genetic changes underlying the evolution of limb development. In this regard, comparative genomic annotation of the false killer whale LMBR1 locus could shed new light on the evolution of limb development. We sequenced two false killer whale BAC clones, corresponding to 156 kb and 144 kb, respectively, harboring the tightly linked RNF32, LMBR1, and NOM1 genes. Our annotation of the false killer whale LMBR1 gene showed that it consists of 17 exons (1473 bp), in contrast to 18 exons (1596 bp) in human, and it displays 93.1% and 95.6% nucleotide and amino acid sequence similarity, respectively, compared with the human gene. In particular, we discovered that exon 10, deleted in the false killer whale LMBR1 gene, is present only in primates, and this fact strongly implies that exon 10 might be crucial in determining primate-specific limb development. ZRS and TFBS sequences have been well conserved across 11 species, suggesting that these regions could be involved in an important function of limb development and limb patterning. The neighboring gene RNF32 showed several lineage-conserved exons, such as exons 2 through 9 conserved in eutherian mammals, exons 3 through 9 conserved in mammals, and exons 5 through 9 conserved in vertebrates. The other neighboring gene, NOM1, had undergone a substitution (ATG→GTA) at the start codon, giving rise to a 36 bp shorter N-terminal sequence compared with the human sequence. Our comparative analysis of the false killer whale LMBR1 genomic locus provides important clues regarding the genetic regions that may play crucial roles in limb development and patterning.

Identification of a distant cis-regulatory element controlling pharyngeal arch-specific expression of zebrafish gdf6a/radar

Skeletal formation is an essential and intricately regulated part of vertebrate development. Humans and mice deficient in growth and differentiation factor 6 (Gdf6) have numerous skeletal abnormalities, including joint fusions and cartilage reductions. The expression of Gdf6 is dynamic and in part regulated by distant evolutionarily conserved cis-regulatory elements. radar/gdf6a is a zebrafish ortholog of Gdf6 and has an essential role in embryonic patterning. Here, we show that radar is transcribed in the cells surrounding and between the developing cartilages of the ventral pharyngeal arches, similar to mouse Gdf6. A 312 bp evolutionarily conserved region (ECR5), 122 kilobases downstream, drives expression in a pharyngeal arch-specific manner similar to endogenous radar/gdf6a. Deletion analysis identified a 78 bp region within ECR5 that is essential for transgene activity. This work illustrates that radar is regulated in the pharyngeal arches by a distant conserved element and suggests radar has similar functions in skeletal development in fish and mammals.

The specification of glycinergic neurons and the role of glycinergic transmission in development

Glycine's role as an inhibitory neurotransmitter in the adult vertebrate nervous system has been well characterized in a number of different model organisms. However, a full understanding of glycinergic transmission requires a knowledge of how glycinergic synapses emerge and the role of glycinergic signaling during development. Recent literature has provided a detailed picture of the developmental expression of many of the molecular components that comprise the glycinergic phenotype, namely the glycine transporters and the glycine receptor subunits; the transcriptional networks leading to the expression of this important neurotransmitter phenotype are also being elucidated. An equally important focus of research has revealed the critical role of glycinergic signaling in sculpting many different aspects of neural development. This review examines the current literature detailing the expression patterns of the components of the glycinergic phenotype in various vertebrate model organisms over the course of development and the molecular mechanisms governing the expression of the glycinergic phenotype. The review then surveys the recent work on the role of glycinergic signaling in the developing nervous system and concludes with an overview of areas for further research.

Mouse Dfa is a repressor of TATA-box promoters and interacts with the Abt1 activator of basal transcription

Our study of the mouse Ate1 arginyltransferase, a component of the N-end rule pathway, has shown that Ate1 pre-mRNA is produced from a bidirectional promoter that also expresses, in the opposite direction, a previously uncharacterized gene (Hu, R. G., Brower, C. S., Wang, H., Davydov, I. V., Sheng, J., Zhou, J., Kwon, Y. T., and Varshavsky, A. (2006) J. Biol. Chem. 281, 32559-32573). In this work, we began analyzing this gene, termed Dfa (divergent from Ate1). Mouse Dfa was found to be transcribed from both the bidirectional P(Ate1/Dfa) promoter and other nearby promoters. The resulting transcripts are alternatively spliced, yielding a complex set of Dfa mRNAs that are present largely, although not exclusively, in the testis. A specific Dfa mRNA encodes, via its 3'-terminal exon, a 217-residue protein termed Dfa(A). Other Dfa mRNAs also contain this exon. Dfa(A) is sequelogous (similar in sequence) to a region of the human/mouse HTEX4 protein, whose physiological function is unknown. We produced an affinity-purified antibody to recombinant mouse Dfa(A) that detected a 35-kDa protein in the mouse testis and in several cell lines. Experiments in which RNA interference was used to down-regulate Dfa indicated that the 35-kDa protein was indeed Dfa(A). Furthermore, Dfa(A) was present in the interchromatin granule clusters and was also found to bind to the Ggnbp1 gametogenetin-binding protein-1 and to the Abt1 activator of basal transcription that interacts with the TATA-binding protein. Given these results, RNA interference was used to probe the influence of Dfa levels in luciferase reporter assays. We found that Dfa(A) acts as a repressor of TATA-box transcriptional promoters.

Cross-species sequence comparisons: a review of methods and available resources

With the availability of whole-genome sequences for an increasing number of species, we are now faced with the challenge of decoding the information contained within these DNA sequences. Comparative analysis of DNA sequences from multiple species at varying evolutionary distances is a powerful approach for identifying coding and functional noncoding sequences, as well as sequences that are unique for a given organism. In this review, we outline the strategy for choosing DNA sequences from different species for comparative analyses and describe the methods used and the resources publicly available for these studies.