A Network of Genes Antagonistic to the LIN-35 Retinoblastoma Protein of Caenorhabditis elegans

The Caenorhabditis elegans pRb ortholog, LIN-35, functions in a wide range of cellular and developmental processes. This includes a role of LIN-35 in nutrient utilization by the intestine, which it carries out redundantly with SLR-2, a zinc-finger protein. This and other redundant functions of LIN-35 were identified in genetic screens for mutations that display synthetic phenotypes in conjunction with loss of lin-35. To explore the intestinal role of LIN-35, we conducted a genome-wide RNA-interference-feeding screen for suppressors of lin-35; slr-2 early larval arrest. Of the 26 suppressors identified, 17 fall into three functional classes: (1) ribosome biogenesis genes, (2) mitochondrial prohibitins, and (3) chromatin regulators. Further characterization indicates that different categories of suppressors act through distinct molecular mechanisms. We also tested lin-35; slr-2 suppressors, as well as suppressors of the synthetic multivulval phenotype, to determine the spectrum of lin-35-synthetic phenotypes that could be suppressed following inhibition of these genes. We identified 19 genes, most of which are evolutionarily conserved, that can suppress multiple unrelated lin-35-synthetic phenotypes. Our study reveals a network of genes broadly antagonistic to LIN-35 as well as genes specific to the role of LIN-35 in intestinal and vulval development. Suppressors of multiple lin-35 phenotypes may be candidate targets for anticancer therapies. Moreover, screening for suppressors of phenotypically distinct synthetic interactions, which share a common altered gene, may prove to be a novel and effective approach for identifying genes whose activities are most directly relevant to the core functions of the shared gene.     

Editor's choice: AnABlast: a new in silico strategy for the genome-wide search of novel genes and fossil regions

Genome annotation, assisted by computer programs, is one of the great advances in modern biology. Nevertheless, the in silico identification of small and complex coding sequences is still challenging. We observed that amino acid sequences inferred from coding—but rarely from non-coding—DNA sequences accumulated alignments in low-stringency BLAST searches, suggesting that this alignments accumulation could be used to highlight coding regions in sequenced DNA. To investigate this possibility, we developed a computer program (AnABlast) that generates profiles of accumulated alignments in query amino acid sequences using a low-stringency BLAST strategy. To validate this approach, all six-frame translations of DNA sequences between every two annotated exons of the fission yeast genome were analysed with AnABlast. AnABlast-generated profiles identified three new copies of known genes, and four new genes supported by experimental evidence. New pseudogenes, ancestral carboxyl- and amino-terminal subtractions, complex gene rearrangements, and ancient fragments of mitDNA and of bacterial origin, were also inferred. Thus, this novel in silico approach provides a powerful tool to uncover new genes, as well as fossil-coding sequences, thus providing insight into the evolutionary history of annotated genomes.     

Annotation of hypothetical proteins orthologous in Pongo abelii and Sus scrofa

A hypothetical protein is predicted to be expressed from an open reading frame without known experimental evidence of translation. They constitute a substantial fraction of proteomes. Domain extraction from these hypothetical sequences helps to search for protein coding genes for protein structural and functional annotation. We describe the analysis of prediction data in a sequence dataset of hypothetical protein orthologs of Pongo abelii (orangutan) and Sus scrofa (pig). It should be noted that these orangutan-pig orthologs are also non-homologous to human proteins. These predicted data find application in the genome wide annotation of proteins in poorly understood genomes.

Abbreviations

PDB - Protein Data Bank, DEG - Database of Essential Genes, CDD - Conserved Domain Database, IUCN - International Union for Conservation of Nature.     

Differences in DBA/1J and DBA/2J reveal lipid QTL genes

Recent advances in mouse genomics have revealed considerable variation in the form of single-nucleotide polymorphisms (SNPs) among common inbred strains. This has made it possible to characterize closely related strains and to identify genes that differ; such genes may be causal for quantitative phenotypes. The mouse strains DBA/1J and DBA/2J differ by just 5.6% at the SNP level. These strains exhibit differences in a number of metabolic and lipid phenotypes, such as plasma levels of triglycerides (TGs) and HDL. A cross between these strains revealed multiple quantitative trait loci (QTLs) in 294 progeny. We identified significant TG QTLs on chromosomes (Chrs) 1, 2, 3, 4, 8, 9, 10, 11, 12, 13, 14, 16, and 19, and significant HDL QTLs on Chrs 3, 9, and 16. Some QTLs mapped to chromosomes with limited variability between the two strains, thus facilitating the identification of candidate genes. We suggest that Tshr is the QTL gene for Chr 12 TG and HDL levels and that Ihh may account for the TG QTL on Chr 1. This cross highlights the advantage of crossing closely related strains for subsequent identification of QTL genes.     

New features of Asian Crassostrea oyster mitochondrial genomes: A novel alloacceptor tRNA gene recruitment and two novel ORFs

A feasible way to perform evolutionary analyses is to compare characters divergent enough to observe significant differences, but sufficiently similar to exclude saturation of the differences that occurred. Thus, comparisons of invertebrate mitochondrial (mt) genomes at low taxonomic levels can be extremely helpful in investigating patterns of variation and evolutionary dynamics of genomes, as intermediate stages of the process may be identified. Fortunately, in this study, we newly sequenced the mt genome of the eighth member of Asian Crassostrea oysters which can provide necessary intermediate characters for us to believe that the variation of Crassostrea mt genomes is considerably greater than previously acknowledged. Several new features of Asian Crassostrea oyster mitochondrial genomes were revealed, and our results are particularly significant as they 1) suggest a novel model of alloacceptor tRNA gene recruitment, namely "vertical" tRNA gene recruitment, which can be successfully used to explain the origination of the unusually additional trnK and trnQ genes (annotated as trnK(2) and trnQ(2) respectively) in the mt genomes of the five Asian oysters, and we speculate that this recruitment progress may be a common phenomenon in the evolution of the tRNA multigene family; 2) reveal the existence of two additional, lineage-specific, mtDNA-encoded genes that may originate from duplication of nad2 followed by rapid evolutionary change. Each of these two genes encodes a unique amino terminal signal peptide, thus each might possess an unknown function; and 3) identify for the first time the atp8 gene in oysters. The present study thus gives further credence to the comparison of congeneric bivalves as a meaningful strategy to investigate mt genomic evolutionary trends in genome organization, tRNA multigene family, and gene loss and/or duplication that are difficult to undertake at higher taxonomic levels. In particular, our study provides new evidence for the identification and characterization of ORFs in the "non-coding region" of animal mt genomes.     

Trypanosomatid comparative genomics: Contributions to the study of parasite biology and different parasitic diseases

In 2005, draft sequences of the genomes of Trypanosoma brucei, Trypanosoma cruzi and Leishmania major, also known as the Tri-Tryp genomes, were published. These protozoan parasites are the causative agents of three distinct insect-borne diseases, namely sleeping sickness, Chagas disease and leishmaniasis, all with a worldwide distribution. Despite the large estimated evolutionary distance among them, a conserved core of ~6,200 trypanosomatid genes was found among the Tri-Tryp genomes. Extensive analysis of these genomic sequences has greatly increased our understanding of the biology of these parasites and their host-parasite interactions. In this article, we review the recent advances in the comparative genomics of these three species. This analysis also includes data on additional sequences derived from other trypanosmatid species, as well as recent data on gene expression and functional genomics. In addition to facilitating the identification of key parasite molecules that may provide a better understanding of these complex diseases, genome studies offer a rich source of new information that can be used to define potential new drug targets and vaccine candidates for controlling these parasitic infections.     

Characterization of the complete mitochondrial genomes of Cnaphalocrocis medinalis and Chilo suppressalis (Lepidoptera: Pyralidae)

The complete mitochondrial genomes (mitogenomes) of Cnaphalocrocis medinalis and Chilo suppressalis (Lepidoptera: Pyralidae) were determined and analyzed. The circular genomes were 15,388 bp long for C. medinalis and 15,395 bp long for C. suppressalis. Both mitogenomes contained 37 genes, with gene order similar to that of other lepidopterans. Notably, 12 protein-coding genes (PCGs) utilized the standard ATN, but the cox1 gene used CGA as the initiation codon; the cox1, cox2, and nad4 genes in the two mitogenomes had the truncated termination codons T, T, and TA, respectively, but the nad5 gene was found to use T as the termination codon only in the C. medinalis mitogenome. Additionally, the codon distribution and Relative Synonymous Codon Usage of the 13 PCGs in the C. medinalis mitogenome were very different from those in other pyralid moth mitogenomes. Most of the tRNA genes had typical cloverleaf secondary structures. However, the dihydrouridine (DHU) arm of the trnS1(AGN) gene did not form a stable stem-loop structure. Forty-nine helices in six domains, and 33 helices in three domains were present in the secondary structures of the rrnL and rrnS genes of the two mitogenomes, respectively. There were four major intergenic spacers, except for the A+T-rich region, spanning at least 12 bp in the two mitogenomes. The A+T-rich region contained an 'ATAGT(A)'-like motif followed by a poly-T stretch in the two mitogenomes. In addition, there were a potential stem-loop structure, a duplicated 25-bp repeat element, and a microsatellite '(TA)(13)' observed in the A+T-rich region of the C. medinalis mitogenome. A poly-T motif, a duplicated 31-bp repeat element, and a 19-bp triplication were found in the C. suppressalis mitogenome. However, there are many differences in the A+T-rich regions between the C. suppressalis mitogenome sequence in the present study and previous reports. Finally, the phylogenetic relationships of these insects were reconstructed based on amino acid sequences of mitochondrial 13 PCGs using Bayesian inference and maximum likelihood methods. These molecular-based phylogenies support the traditional morphologically based view of relationships within the Pyralidae.     

Genome sequencing of a single cell of the widely distributed marine subsurface Dehalococcoidia,phylum Chloroflexi

Bacteria of the class Dehalococcoidia (DEH), phylum Chloroflexi, are widely distributed in the marine subsurface, yet metabolic properties of the many uncultivated lineages are completely unknown. This study therefore analysed genomic content from a single DEH cell designated ‘DEH-J10'' obtained from the sediments of Aarhus Bay, Denmark. Real-time PCR showed the DEH-J10 phylotype was abundant in upper sediments but was absent below 160 cm below sea floor. A 1.44 Mbp assembly was obtained and was estimated to represent up to 60.8% of the full genome. The predicted genome is much larger than genomes of cultivated DEH and appears to confer metabolic versatility. Numerous genes encoding enzymes of core and auxiliary beta-oxidation pathways were identified, suggesting that this organism is capable of oxidising various fatty acids and/or structurally related substrates. Additional substrate versatility was indicated by genes, which may enable the bacterium to oxidise aromatic compounds. Genes encoding enzymes of the reductive acetyl-CoA pathway were identified, which may also enable the fixation of CO₂ or oxidation of organics completely to CO₂. Genes encoding a putative dimethylsulphoxide reductase were the only evidence for a respiratory terminal reductase. No evidence for reductive dehalogenase genes was found. Genetic evidence also suggests that the organism could synthesise ATP by converting acetyl-CoA to acetate by substrate-level phosphorylation. Other encoded enzymes putatively conferring marine adaptations such as salt tolerance and organo-sulphate sulfohydrolysis were identified. Together, these analyses provide the first insights into the potential metabolic traits that may enable members of the DEH to occupy an ecological niche in marine sediments.     

Spliceosomal snRNAs in the unicellular eukaryote Trichomonas vaginalis are structurally conserved but lack a 5'-cap structure

Few genes in the divergent eukaryote Trichomonas vaginalis have introns, despite the unusually large gene repertoire of this human-infective parasite. These introns are characterized by extended conserved regulatory motifs at the 5' and 3' boundaries, a feature shared with another divergent eukaryote, Giardia lamblia, but not with metazoan introns. This unusual characteristic of T. vaginalis introns led us to examine spliceosomal small nuclear RNAs (snRNAs) predicted to mediate splicing reactions via interaction with intron motifs. Here we identify T. vaginalis U1, U2, U4, U5, and U6 snRNAs, present predictions of their secondary structures, and provide evidence for interaction between the U2/U6 snRNA complex and a T. vaginalis intron. Structural models predict that T. vaginalis snRNAs contain conserved sequences and motifs similar to those found in other examined eukaryotes. These data indicate that mechanisms of intron recognition as well as coordination of the two catalytic steps of splicing have been conserved throughout eukaryotic evolution. Unexpectedly, we found that T. vaginalis spliceosomal snRNAs lack the 5' trimethylguanosine cap typical of snRNAs and appear to possess unmodified 5' ends. Despite the lack of a cap structure, U1, U2, U4, and U5 genes are transcribed by RNA polymerase II, whereas the U6 gene is transcribed by RNA polymerase III.     

The glypican 3-hosted murine mir717 gene: sequence conservation, seed region polymorphisms and putative targets

Mir717 (mmu-mir-717) was first reported in mouse and resides in the intron 3 of glypican 3 (Gpc3) gene. Our present study revealed that this microRNA (miRNA) gene is conserved among 26 mammalian species and harbors polymorphic sites within the mature seed region in mice. Our finding represents a rare four layer genomic overlap consisting of growth associated quantitative trait locus (QTL), body mass associated Gpc3 gene, highly conserved miRNA gene and mature miRNA seed single nucleotide polymorphism (SNP) identified in the lean mouse strain 129/Sv. Additionally, genes potentially targeted by Mir717 include 91 genes associated with obesity and related phenotypes in mammals. Our analysis provides a basis for further experiments to causally connect the identified SNP and Mir717 gene itself to obesity regulation. Furthermore, our bioinformatics analysis now enables functional annotation of Mir717 orthologs in other species, thus determining the effect of its target genes on fat-related traits.     

The mitochondrial genomes of Amphiascoides atopus and Schizopera knabeni (Harpacticoida: Miraciidae) reveal similarities between the copepod orders Harpacticoida and Poecilostomatoida

Members of subclass Copepoda are abundant, diverse, and—as a result of their variety of ecological roles in marine and freshwater environments—important, but their phylogenetic interrelationships are unclear. Recent studies of arthropods have used gene arrangements in the mitochondrial (mt) genome to infer phylogenies, but for copepods, only seven complete mt genomes have been published. These data revealed several within-order and few among-order similarities. To increase the data available for comparisons, we sequenced the complete mt genome (13,831 base pairs) of Amphiascoides atopus and 10,649 base pairs of the mt genome of Schizopera knabeni (both in the family Miraciidae of the order Harpacticoida). Comparison of our data to those for Tigriopus japonicus (family Harpacticidae, order Harpacticoida) revealed similarities in gene arrangement among these three species that were consistent with those found within and among families of other copepod orders. Comparison of the mt genomes of our species with those known from other copepod orders revealed the arrangement of mt genes of our Harpacticoida species to be more similar to that of Sinergasilus polycolpus (order Poecilostomatoida) than to that of T. japonicus. The similarities between S. polycolpus and our species are the first to be noted across the boundaries of copepod orders and support the possibility that mt-gene arrangement might be used to infer copepod phylogenies. We also found that our two species had extremely truncated transfer RNAs and that gene overlaps occurred much more frequently than has been reported for other copepod mt genomes.     

A unique tRNA gene family and a novel,highly expressed ORF in the mitochondrial genome of the silver-lip pearl oyster, Pinctada maxima (Bivalvia: Pteriidae)

Characteristics of mitochondrial (mt) DNA such as gene content and arrangement, as well as mt tRNA secondary structure, are frequently used in comparative genomic analyses because they provide valuable phylogenetic information. However, most analyses do not characterize the relationship of tRNA genes from the same mt genome and, in some cases, analyses overlook possible novel open reading frames (ORFs) when the 13 expected protein-coding genes are already annotated. In this study, we describe the sequence and characterization of the complete mt genome of the silver-lip pearl oyster, Pinctada maxima. The 16,994-bp mt genome contains the same 13 protein-coding genes (PCGs) and two ribosomal RNA genes typical of metazoans. The gene arrangement, however, is completely distinct from that of all other available bivalve mt genomes, and a unique tRNA gene family is observed in this genome. The unique tRNA gene family includes two trnS^− AGY and trnQ genes, a trnM isomerism, but it lacks trnS^− CUN. We also report the first clear evidence of alloacceptor tRNA gene recruitment (trnP → trnS^− AGY) in mollusks. In addition, a novel ORF (orfUR1) expressed at high levels is present in the mt genome of this pearl oyster. This gene contains a conserved domain, “Oxidored_q1_N”, which is a member of Complex I and thus may play an important role in key biological functions. Because orfUR1 has a very similar nucleotide composition and codon bias to that of other genes in this genome, we hypothesize that this gene may have been moved to the mt genome via gene transfer from the nuclear genome at an early stage of speciation of P. maxima, or it may have evolved as a result of gene duplication, followed by rapid sequence divergence. Lastly, a 319-bp region was identified as the possible control region (CR) even though it does not correspond to the longest non-coding region in the genome. Unlike other studies of mt genomes, this study compares the evolutionary patterns of all available bivalve mt tRNA and atp8 genes.     

Effector-Triggered Immune Response in Arabidopsis thaliana Is a Quantitative Trait

We identified loci responsible for natural variation in Arabidopsis thaliana (Arabidopsis) responses to a bacterial pathogen virulence factor, HopAM1. HopAM1 is a type III effector protein secreted by the virulent Pseudomonas syringae strain Pto DC3000. Delivery of HopAM1 from disarmed Pseudomonas strains leads to local cell death, meristem chlorosis, or both, with varying intensities in different Arabidopsis accessions. These phenotypes are not associated with differences in bacterial growth restriction. We treated the two phenotypes as quantitative traits to identify host loci controlling responses to HopAM1. Genome-wide association (GWA) of 64 Arabidopsis accessions identified independent variants highly correlated with response to each phenotype. Quantitative trait locus (QTL) mapping in a recombinant inbred population between Bur-0 and Col-0 accessions revealed genetic linkage to regions distinct from the top GWA hits. Two major QTL associated with HopAM1-induced cell death were also associated with HopAM1-induced chlorosis. HopAM1-induced changes in Arabidopsis gene expression showed that rapid HopAM1-dependent cell death in Bur-0 is correlated with effector-triggered immune responses. Studies of the effect of mutations in known plant immune system genes showed, surprisingly, that both cell death and chlorosis phenotypes are enhanced by loss of EDS1, a regulatory hub in the plant immune-signaling network. Our results reveal complex genetic architecture for response to this particular type III virulence effector, in contrast to the typical monogenic control of cell death and disease resistance triggered by most type III effectors.