Cyanobacteria toxins in the Salton Sea

Background

Sequenced archaeal genomes contain a variety of bacterial and eukaryotic DNA repair gene homologs, but relatively little is known about how these microorganisms actually perform DNA repair. At least some archaea, including the extreme halophile Halobacterium sp. NRC-1, are able to repair ultraviolet light (UV) induced DNA damage in the absence of light-dependent photoreactivation but this 'dark' repair capacity remains largely uncharacterized. Halobacterium sp. NRC-1 possesses homologs of the bacterial uvrA, uvrB, and uvrC nucleotide excision repair genes as well as several eukaryotic repair genes and it has been thought that multiple DNA repair pathways may account for the high UV resistance and dark repair capacity of this model halophilic archaeon. We have carried out a functional analysis, measuring repair capability in uvrA, uvrB and uvrC deletion mutants.

Results

Deletion mutants lacking functional uvrA, uvrB or uvrC genes, including a uvrA uvrC double mutant, are hypersensitive to UV and are unable to remove cyclobutane pyrimidine dimers or 6–4 photoproducts from their DNA after irradiation with 150 J/m² of 254 nm UV-C. The UV sensitivity of the uvr mutants is greatly attenuated following incubation under visible light, emphasizing that photoreactivation is highly efficient in this organism. Phylogenetic analysis of the Halobacterium uvr genes indicates a complex ancestry.

Conclusion

Our results demonstrate that homologs of the bacterial nucleotide excision repair genes uvrA, uvrB, and uvrC are required for the removal of UV damage in the absence of photoreactivating light in Halobacterium sp. NRC-1. Deletion of these genes renders cells hypersensitive to UV and abolishes their ability to remove cyclobutane pyrimidine dimers and 6–4 photoproducts in the absence of photoreactivating light. In spite of this inability to repair UV damaged DNA, uvrA, uvrB and uvrC deletion mutants are substantially less UV sensitive than excision repair mutants of E. coli or yeast. This may be due to efficient damage tolerance mechanisms such as recombinational lesion bypass, bypass DNA polymerase(s) and the existence of multiple genomes in Halobacterium. Phylogenetic analysis provides no clear evidence for lateral transfer of these genes from bacteria to archaea.     

Fine mapping and chromosome walking towards the Ror1 locus in barley (Hordeum vulgare L.)

Key message

The Ror1 gene was fine-mapped to the pericentric region of barley chromosome 1HL.

Abstract

Recessively inherited loss-of-function alleles of the barley (Hordeum vulgare) Mildew resistance locus o (Mlo) gene confer durable broad-spectrum disease resistance against the obligate biotrophic fungal powdery mildew pathogen Blumeria graminis f.sp. hordei. Previous genetic analyses revealed two barley genes, Ror1 and Ror2, that are Required for mlo-specified resistance and basal defence. While Ror2 was cloned and shown to encode a t-SNARE protein (syntaxin), the molecular nature or Ror1 remained elusive. Ror1 was previously mapped to the centromeric region of the long arm of barley chromosome 1H. Here, we narrowed the barley Ror1 interval to 0.18 cM and initiated a chromosome walk using barley yeast artificial chromosome (YAC) clones, next-generation DNA sequencing and fluorescence in situ hybridization. Two non-overlapping YAC contigs containing Ror1 flanking genes were identified. Despite a high degree of synteny observed between barley and the sequenced genomes of the grasses rice (Oryza sativa), Brachypodium distachyon and Sorghum bicolor across the wider chromosomal area, the genes in the YAC contigs showed extensive interspecific rearrangements in orientation and order. Consequently, the position of a Ror1 homolog in these species could not be precisely predicted, nor was a barley gene co-segregating with Ror1 identified. These factors have prevented the molecular identification of the Ror1 gene for the time being.     

Distribution,functional impact,and origin mechanisms of copy number variation in the barley genome

Background

There is growing evidence for the prevalence of copy number variation (CNV) and its role in phenotypic variation in many eukaryotic species. Here we use array comparative genomic hybridization to explore the extent of this type of structural variation in domesticated barley cultivars and wild barleys.

Results

A collection of 14 barley genotypes including eight cultivars and six wild barleys were used for comparative genomic hybridization. CNV affects 14.9% of all the sequences that were assessed. Higher levels of CNV diversity are present in the wild accessions relative to cultivated barley. CNVs are enriched near the ends of all chromosomes except 4H, which exhibits the lowest frequency of CNVs. CNV affects 9.5% of the coding sequences represented on the array and the genes affected by CNV are enriched for sequences annotated as disease-resistance proteins and protein kinases. Sequence-based comparisons of CNV between cultivars Barke and Morex provided evidence that DNA repair mechanisms of double-strand breaks via single-stranded annealing and synthesis-dependent strand annealing play an important role in the origin of CNV in barley.

Conclusions

We present the first catalog of CNVs in a diploid Triticeae species, which opens the door for future genome diversity research in a tribe that comprises the economically important cereal species wheat, barley, and rye. Our findings constitute a valuable resource for the identification of CNV affecting genes of agronomic importance. We also identify potential mechanisms that can generate variation in copy number in plant genomes.     

Novel AlkB dioxygenases--alternative models for in silico and in vivo studies

Background

ALKBH proteins, the homologs of Escherichia coli AlkB dioxygenase, constitute a direct, single-protein repair system, protecting cellular DNA and RNA against the cytotoxic and mutagenic activity of alkylating agents, chemicals significantly contributing to tumor formation and used in cancer therapy. In silico analysis and in vivo studies have shown the existence of AlkB homologs in almost all organisms. Nine AlkB homologs (ALKBH1–8 and FTO) have been identified in humans. High ALKBH levels have been found to encourage tumor development, questioning the use of alkylating agents in chemotherapy. The aim of this work was to assign biological significance to multiple AlkB homologs by characterizing their activity in the repair of nucleic acids in prokaryotes and their subcellular localization in eukaryotes.

Methodology and Findings

Bioinformatic analysis of protein sequence databases identified 1943 AlkB sequences with eight new AlkB subfamilies. Since Cyanobacteria and Arabidopsis thaliana contain multiple AlkB homologs, they were selected as model organisms for in vivo research. Using E. coli alkB ⁻ mutant and plasmids expressing cyanobacterial AlkBs, we studied the repair of methyl methanesulfonate (MMS) and chloroacetaldehyde (CAA) induced lesions in ssDNA, ssRNA, and genomic DNA. On the basis of GFP fusions, we investigated the subcellular localization of ALKBHs in A. thaliana and established its mostly nucleo-cytoplasmic distribution. Some of the ALKBH proteins were found to change their localization upon MMS treatment.

Conclusions

Our in vivo studies showed highly specific activity of cyanobacterial AlkB proteins towards lesions and nucleic acid type. Subcellular localization and translocation of ALKBHs in A. thaliana indicates a possible role for these proteins in the repair of alkyl lesions. We hypothesize that the multiplicity of ALKBHs is due to their involvement in the metabolism of nucleo-protein complexes; we find their repair by ALKBH proteins to be economical and effective alternative to degradation and de novo synthesis.     

Comparative genomics of mutualistic viruses of Glyptapanteles parasitic wasps

Background

Polydnaviruses, double-stranded DNA viruses with segmented genomes, have evolved as obligate endosymbionts of parasitoid wasps. Virus particles are replication deficient and produced by female wasps from proviral sequences integrated into the wasp genome. These particles are co-injected with eggs into caterpillar hosts, where viral gene expression facilitates parasitoid survival and, thereby, survival of proviral DNA. Here we characterize and compare the encapsidated viral genome sequences of bracoviruses in the family Polydnaviridae associated with Glyptapanteles gypsy moth parasitoids, along with near complete proviral sequences from which both viral genomes are derived.

Results

The encapsidated Glyptapanteles indiensis and Glyptapanteles flavicoxis bracoviral genomes, each composed of 29 different size segments, total approximately 517 and 594 kbp, respectively. They are generated from a minimum of seven distinct loci in the wasp genome. Annotation of these sequences revealed numerous novel features for polydnaviruses, including insect-like sugar transporter genes and transposable elements. Evolutionary analyses suggest that positive selection is widespread among bracoviral genes.

Conclusions

The structure and organization of G. indiensis and G. flavicoxis bracovirus proviral segments as multiple loci containing one to many viral segments, flanked and separated by wasp gene-encoding DNA, is confirmed. Rapid evolution of bracovirus genes supports the hypothesis of bracovirus genes in an 'arms race' between bracovirus and caterpillar. Phylogenetic analyses of the bracoviral genes encoding sugar transporters provides the first robust evidence of a wasp origin for some polydnavirus genes. We hypothesize transposable elements, such as those described here, could facilitate transfer of genes between proviral segments and host DNA.     

MutS and MutL Are Dispensable for Maintenance of the Genomic Mutation Rate in the Halophilic Archaeon Halobacterium salinarum NRC-1

Background

The genome of the halophilic archaeon Halobacterium salinarum NRC-1 encodes for homologs of MutS and MutL, which are key proteins of a DNA mismatch repair pathway conserved in Bacteria and Eukarya. Mismatch repair is essential for retaining the fidelity of genetic information and defects in this pathway result in the deleterious accumulation of mutations and in hereditary diseases in humans.

Methodology/Principal Findings

We calculated the spontaneous genomic mutation rate of H. salinarum NRC-1 using fluctuation tests targeting genes of the uracil monophosphate biosynthesis pathway. We found that H. salinarum NRC-1 has a low incidence of mutation suggesting the presence of active mechanisms to control spontaneous mutations during replication. The spectrum of mutational changes found in H. salinarum NRC-1, and in other archaea, appears to be unique to this domain of life and might be a consequence of their adaption to extreme environmental conditions. In-frame targeted gene deletions of H. salinarum NRC-1 mismatch repair genes and phenotypic characterization of the mutants demonstrated that the mutS and mutL genes are not required for maintenance of the observed mutation rate.

Conclusions/Significance

We established that H. salinarum NRC-1 mutS and mutL genes are redundant to an alternative system that limits spontaneous mutation in this organism. This finding leads to the puzzling question of what mechanism is responsible for maintenance of the low genomic mutation rates observed in the Archaea, which for the most part do not have MutS and MutL homologs.     

Genome-wide identification and characterization of replication origins by deep sequencing

Background

DNA replication initiates at distinct origins in eukaryotic genomes, but the genomic features that define these sites are not well understood.

Results

We have taken a combined experimental and bioinformatic approach to identify and characterize origins of replication in three distantly related fission yeasts: Schizosaccharomyces pombe, Schizosaccharomyces octosporus and Schizosaccharomyces japonicus. Using single-molecule deep sequencing to construct amplification-free high-resolution replication profiles, we located origins and identified sequence motifs that predict origin function. We then mapped nucleosome occupancy by deep sequencing of mononucleosomal DNA from the corresponding species, finding that origins tend to occupy nucleosome-depleted regions.

Conclusions

The sequences that specify origins are evolutionarily plastic, with low complexity nucleosome-excluding sequences functioning in S. pombe and S. octosporus, and binding sites for trans-acting nucleosome-excluding proteins functioning in S. japonicus. Furthermore, chromosome-scale variation in replication timing is conserved independently of origin location and via a mechanism distinct from known heterochromatic effects on origin function. These results are consistent with a model in which origins are simply the nucleosome-depleted regions of the genome with the highest affinity for the origin recognition complex. This approach provides a general strategy for understanding the mechanisms that define DNA replication origins in eukaryotes.     

AtRAD5A is a DNA translocase harboring a HIRAN domain which confers binding to branched DNA structures and is required for DNA repair in vivo

DNA lesions such as crosslinks represent obstacles for the replication machinery. Nonetheless, replication can proceed via the DNA damage tolerance pathway also known as postreplicative repair pathway. SNF2 ATPase Rad5 homologs, such as RAD5A of the model plant Arabidopsis thaliana, are important for the error‐free mode of this pathway. We able to demonstrate before, that RAD5A is a key factor in the repair of DNA crosslinks in Arabidopsis. Here, we show by in vitro analysis that AtRAD5A protein is a DNA translocase able to catalyse fork regression. Interestingly, replication forks with a gap in the leading strand are processed best, in line with its suggested function. Furthermore AtRAD5A catalyses branch migration of a Holliday junction and is furthermore not impaired by the DNA binding of a model protein, which is indicative of its ability to displace other proteins. Rad5 homologs possess HIRAN (Hip116, Rad5; N‐terminal) domains. By biochemical analysis we were able to demonstrate that the HIRAN domain variant from Arabidopsis RAD5A mediates structure selective DNA binding without the necessity for a free 3′OH group as has been shown to be required for binding of HIRAN domains in a mammalian RAD5 homolog. The biological importance of the HIRAN domain in AtRAD5A is demonstrated by our result that it is required for its function in DNA crosslink repair in vivo.     

The influence of gene-environment interactions on GHR and IGF-1 expression and their association with growth in brook charr, Salvelinus fontinalis (Mitchill)

Background

Information transfer systems in Archaea, including many components of the DNA replication machinery, are similar to those found in eukaryotes. Functional assignments of archaeal DNA replication genes have been primarily based upon sequence homology and biochemical studies of replisome components, but few genetic studies have been conducted thus far. We have developed a tractable genetic system for knockout analysis of genes in the model halophilic archaeon, Halobacterium sp. NRC-1, and used it to determine which DNA replication genes are essential.

Results

Using a directed in-frame gene knockout method in Halobacterium sp. NRC-1, we examined nineteen genes predicted to be involved in DNA replication. Preliminary bioinformatic analysis of the large haloarchaeal Orc/Cdc6 family, related to eukaryotic Orc1 and Cdc6, showed five distinct clades of Orc/Cdc6 proteins conserved in all sequenced haloarchaea. Of ten orc/cdc6 genes in Halobacterium sp. NRC-1, only two were found to be essential, orc10, on the large chromosome, and orc2, on the minichromosome, pNRC200. Of the three replicative-type DNA polymerase genes, two were essential: the chromosomally encoded B family, polB1, and the chromosomally encoded euryarchaeal-specific D family, polD1/D2 (formerly called polA1/polA2 in the Halobacterium sp. NRC-1 genome sequence). The pNRC200-encoded B family polymerase, polB2, was non-essential. Accessory genes for DNA replication initiation and elongation factors, including the putative replicative helicase, mcm, the eukaryotic-type DNA primase, pri1/pri2, the DNA polymerase sliding clamp, pcn, and the flap endonuclease, rad2, were all essential. Targeted genes were classified as non-essential if knockouts were obtained and essential based on statistical analysis and/or by demonstrating the inability to isolate chromosomal knockouts except in the presence of a complementing plasmid copy of the gene.

Conclusion

The results showed that ten out of nineteen eukaryotic-type DNA replication genes are essential for Halobacterium sp. NRC-1, consistent with their requirement for DNA replication. The essential genes code for two of ten Orc/Cdc6 proteins, two out of three DNA polymerases, the MCM helicase, two DNA primase subunits, the DNA polymerase sliding clamp, and the flap endonuclease.     

Feeling manipulated: cytomegalovirus immune manipulation

Background

ORF59 DNA polymerase processivity factor of the human rhadinovirus, Kaposi's sarcoma-associated herpesvirus (KSHV), is required for efficient copying of the genome during virus replication. KSHV ORF59 is antigenic in the infected host and is used as a marker for virus activation and replication.

Results

We cloned, sequenced and expressed the genes encoding related ORF59 proteins from the RV1 rhadinovirus homologs of KSHV from chimpanzee (PtrRV1) and three species of macaques (RFHVMm, RFHVMn and RFHVMf), and have compared them with ORF59 proteins obtained from members of the more distantly-related RV2 rhadinovirus lineage infecting the same non-human primate species (PtrRV2, RRV, MneRV2, and MfaRV2, respectively). We found that ORF59 homologs of the RV1 and RV2 Old World primate rhadinoviruses are highly conserved with distinct phylogenetic clustering of the two rhadinovirus lineages. RV1 and RV2 ORF59 C-terminal domains exhibit a strong lineage-specific conservation. Rabbit antiserum was developed against a C-terminal polypeptide that is highly conserved between the macaque RV2 ORF59 sequences. This anti-serum showed strong reactivity towards ORF59 encoded by the macaque RV2 rhadinoviruses, RRV (rhesus) and MneRV2 (pig-tail), with no cross reaction to human or macaque RV1 ORF59 proteins. Using this antiserum and RT-qPCR, we determined that RRV ORF59 is expressed early after permissive infection of both rhesus primary fetal fibroblasts and African green monkey kidney epithelial cells (Vero) in vitro. RRV- and MneRV2-infected foci showed strong nuclear expression of ORF59 that correlated with production of infectious progeny virus. Immunohistochemical studies of an MneRV2-infected macaque revealed strong nuclear expression of ORF59 in infected cells within the differentiating layer of epidermis corroborating previous observations that differentiated epithelial cells are permissive for replication of KSHV-like rhadinoviruses.

Conclusion

The ORF59 DNA polymerase processivity factor homologs of the Old World primate RV1 and RV2 rhadinovirus lineages are phylogenetically distinct yet demonstrate similar expression and localization characteristics that correlate with their use as lineage-specific markers for permissive infection and virus replication. These studies will aid in the characterization of virus activation from latency to the replicative state, an important step for understanding the biology and transmission of rhadinoviruses, such as KSHV.     

Whole proteome identification of plant candidate G-protein coupled receptors in Arabidopsis, rice, and poplar: computational prediction and in-vivo protein coupling

Background

The classic paradigm of heterotrimeric G-protein signaling describes a heptahelical, membrane-spanning G-protein coupled receptor that physically interacts with an intracellular Gα subunit of the G-protein heterotrimer to transduce signals. G-protein coupled receptors comprise the largest protein superfamily in metazoa and are physiologically important as they sense highly diverse stimuli and play key roles in human disease. The heterotrimeric G-protein signaling mechanism is conserved across metazoa, and also readily identifiable in plants, but the low sequence conservation of G-protein coupled receptors hampers the identification of novel ones. Using diverse computational methods, we performed whole-proteome analyses of the three dominant model plant species, the herbaceous dicot Arabidopsis thaliana (mouse-eared cress), the monocot Oryza sativa (rice), and the woody dicot Populus trichocarpa (poplar), to identify plant protein sequences most likely to be GPCRs.

Results

Our stringent bioinformatic pipeline allowed the high confidence identification of candidate G-protein coupled receptors within the Arabidopsis, Oryza, and Populus proteomes. We extended these computational results through actual wet-bench experiments where we tested over half of our highest ranking Arabidopsis candidate G-protein coupled receptors for the ability to physically couple with GPA1, the sole Gα in Arabidopsis. We found that seven out of eight tested candidate G-protein coupled receptors do in fact interact with GPA1. We show through G-protein coupled receptor classification and molecular evolutionary analyses that both individual G-protein coupled receptor candidates and candidate G-protein coupled receptor families are conserved across plant species and that, in some cases, this conservation extends to metazoans.

Conclusion

Our computational and wet-bench results provide the first step toward understanding the diversity, conservation, and functional roles of plant candidate G-protein coupled receptors.     

A first survey of the rye (Secale cereale) genome composition through BAC end sequencing of the short arm of chromosome 1R

Background

Thellungiella halophila (also known as Thellungiella salsuginea) is a model halophyte with a small plant size, short life cycle, and small genome. It easily undergoes genetic transformation by the floral dipping method used with its close relative, Arabidopsis thaliana. Thellungiella genes exhibit high sequence identity (approximately 90% at the cDNA level) with Arabidopsis genes. Furthermore, Thellungiella not only shows tolerance to extreme salinity stress, but also to chilling, freezing, and ozone stress, supporting the use of Thellungiella as a good genomic resource in studies of abiotic stress tolerance.

Results

We constructed a full-length enriched Thellungiella (Shan Dong ecotype) cDNA library from various tissues and whole plants subjected to environmental stresses, including high salinity, chilling, freezing, and abscisic acid treatment. We randomly selected about 20 000 clones and sequenced them from both ends to obtain a total of 35 171 sequences. CAP3 software was used to assemble the sequences and cluster them into 9569 nonredundant cDNA groups. We named these cDNAs "RTFL" (RIKEN Thellungiella Full-Length) cDNAs. Information on functional domains and Gene Ontology (GO) terms for the RTFL cDNAs were obtained using InterPro. The 8289 genes assigned to InterPro IDs were classified according to the GO terms using Plant GO Slim. Categorical comparison between the whole Arabidopsis genome and Thellungiella genes showing low identity to Arabidopsis genes revealed that the population of Thellungiella transport genes is approximately 1.5 times the size of the corresponding Arabidopsis genes. This suggests that these genes regulate a unique ion transportation system in Thellungiella.

Conclusion

As the number of Thellungiella halophila (Thellungiella salsuginea) expressed sequence tags (ESTs) was 9388 in July 2008, the number of ESTs has increased to approximately four times the original value as a result of this effort. Our sequences will thus contribute to correct future annotation of the Thellungiella genome sequence. The full-length enriched cDNA clones will enable the construction of overexpressing mutant plants by introduction of the cDNAs driven by a constitutive promoter, the complementation of Thellungiella mutants, and the determination of promoter regions in the Thellungiella genome.     

The Mitochondrial Genome of an Aquatic Plant,Spirodela polyrhiza

Background

Spirodela polyrhiza is a species of the order Alismatales, which represent the basal lineage of monocots with more ancestral features than the Poales. Its complete sequence of the mitochondrial (mt) genome could provide clues for the understanding of the evolution of mt genomes in plant.

Methods

Spirodela polyrhiza mt genome was sequenced from total genomic DNA without physical separation of chloroplast and nuclear DNA using the SOLiD platform. Using a genome copy number sensitive assembly algorithm, the mt genome was successfully assembled. Gap closure and accuracy was determined with PCR products sequenced with the dideoxy method.

Conclusions

This is the most compact monocot mitochondrial genome with 228,493 bp. A total of 57 genes encode 35 known proteins, 3 ribosomal RNAs, and 19 tRNAs that recognize 15 amino acids. There are about 600 RNA editing sites predicted and three lineage specific protein-coding-gene losses. The mitochondrial genes, pseudogenes, and other hypothetical genes (ORFs) cover 71,783 bp (31.0%) of the genome. Imported plastid DNA accounts for an additional 9,295 bp (4.1%) of the mitochondrial DNA. Absence of transposable element sequences suggests that very few nuclear sequences have migrated into Spirodela mtDNA. Phylogenetic analysis of conserved protein-coding genes suggests that Spirodela shares the common ancestor with other monocots, but there is no obvious synteny between Spirodela and rice mtDNAs. After eliminating genes, introns, ORFs, and plastid-derived DNA, nearly four-fifths of the Spirodela mitochondrial genome is of unknown origin and function. Although it contains a similar chloroplast DNA content and range of RNA editing as other monocots, it is void of nuclear insertions, active gene loss, and comprises large regions of sequences of unknown origin in non-coding regions. Moreover, the lack of synteny with known mitochondrial genomic sequences shed new light on the early evolution of monocot mitochondrial genomes.     

A diploid wheat TILLING resource for wheat functional genomics

Background

Cultivated peanut (Arachis hypogaea L.) is an important crop worldwide, valued for its edible oil and digestible protein. It has a very narrow genetic base that may well derive from a relatively recent single polyploidization event. Accordingly molecular markers have low levels of polymorphism and the number of polymorphic molecular markers available for cultivated peanut is still limiting.

Results

Here, we report a large set of BAC-end sequences (BES), use them for developing SSR (BES-SSR) markers, and apply them in genetic linkage mapping. The majority of BESs had no detectable homology to known genes (49.5%) followed by sequences with similarity to known genes (44.3%), and miscellaneous sequences (6.2%) such as transposable element, retroelement, and organelle sequences. A total of 1,424 SSRs were identified from 36,435 BESs. Among these identified SSRs, dinucleotide (47.4%) and trinucleotide (37.1%) SSRs were predominant. The new set of 1,152 SSRs as well as about 4,000 published or unpublished SSRs were screened against two parents of a mapping population, generating 385 polymorphic loci. A genetic linkage map was constructed, consisting of 318 loci onto 21 linkage groups and covering a total of 1,674.4 cM, with an average distance of 5.3 cM between adjacent loci. Two markers related to resistance gene homologs (RGH) were mapped to two different groups, thus anchoring 1 RGH-BAC contig and 1 singleton.

Conclusions

The SSRs mined from BESs will be of use in further molecular analysis of the peanut genome, providing a novel set of markers, genetically anchoring BAC clones, and incorporating gene sequences into a linkage map. This will aid in the identification of markers linked to genes of interest and map-based cloning.