Genome sequence of Haloarcula marismortui: a halophilic archaeon from the Dead Sea

We report the complete sequence of the 4,274,642-bp genome of Haloarcula marismortui, a halophilic archaeal isolate from the Dead Sea. The genome is organized into nine circular replicons of varying G+C compositions ranging from 54% to 62%. Comparison of the genome architectures of Halobacterium sp. NRC-1 and H. marismortui suggests a common ancestor for the two organisms and a genome of significantly reduced size in the former. Both of these halophilic archaea use the same strategy of high surface negative charge of folded proteins as means to circumvent the salting-out phenomenon in a hypersaline cytoplasm. A multitiered annotation approach, including primary sequence similarities, protein family signatures, structure prediction, and a protein function association network, has assigned putative functions for at least 58% of the 4242 predicted proteins, a far larger number than is usually achieved in most newly sequenced microorganisms. Among these assigned functions were genes encoding six opsins, 19 MCP and/or HAMP domain signal transducers, and an unusually large number of environmental response regulators-nearly five times as many as those encoded in Halobacterium sp. NRC-1--suggesting H. marismortui is significantly more physiologically capable of exploiting diverse environments. In comparing the physiologies of the two halophilic archaea, in addition to the expected extensive similarity, we discovered several differences in their metabolic strategies and physiological responses such as distinct pathways for arginine breakdown in each halophile. Finally, as expected from the larger genome, H. marismortui encodes many more functions and seems to have fewer nutritional requirements for survival than does Halobacterium sp. NRC-1.     

Directed networks reveal genomic barriers and DNA repair bypasses to lateral gene transfer among prokaryotes

Lateral gene transfer (LGT) plays a major role in prokaryote evolution with only a few genes that are resistant to it; yet the nature and magnitude of barriers to lateral transfer are still debated. Here, we implement directed networks to investigate donor-recipient events of recent lateral gene transfer among 657 sequenced prokaryote genomes. For 2,129,548 genes investigated, we detected 446,854 recent lateral gene transfer events through nucleotide pattern analysis. Among these, donor-recipient relationships could be specified through phylogenetic reconstruction for 7% of the pairs, yielding 32,028 polarized recent gene acquisition events, which constitute the edges of our directed networks. We find that the frequency of recent LGT is linearly correlated both with genome sequence similarity and with proteome similarity of donor-recipient pairs. Genome sequence similarity accounts for 25% of the variation in gene-transfer frequency, with proteome similarity adding only 1% to the variability explained. The range of donor-recipient GC content similarity within the network is extremely narrow, with 86% of the LGTs occurring between donor-recipient pairs having ≤5% difference in GC content. Hence, genome sequence similarity and GC content similarity are strong barriers to LGT in prokaryotes. But they are not insurmountable, as we detected 1530 recent transfers between distantly related genomes. The directed network revealed that recipient genomes of distant transfers encode proteins of nonhomologous end-joining (NHEJ; a DNA repair mechanism) far more frequently than the recipient lacking that mechanism. This implicates NHEJ in genes spread across distantly related prokaryotes through bypassing the donor-recipient sequence similarity barrier.     

A comparative proteome analysis of human metaphase chromosomes isolated from two different cell lines reveals a set of conserved chromosome-associated proteins

A comparative proteome analysis of human metaphase chromosomes between a typical epithelial-like cell, HeLa S3, and a lymphoma-type cell, BALL-1, was performed. One-dimensional (1-D) SDS-PAGE and radical-free and highly reducing two-dimensional electrophoresis (RFHR 2-DE) detected more than 200 proteins from chromosomes isolated from HeLa S3 cells, among which 189 proteins were identified by mass spectrometry (MS). Consistent with our recent four-layer structural model of a metaphase chromosome, all the identified proteins were grouped into four distinct levels of abundance. Both HeLa S3 and BALL-1 chromosomes contained specific sets of abundant chromosome structural and peripheral proteins in addition to less abundant chromosome coating proteins (CCPs). Furthermore, titin array analysis and a proteome analysis of the ultra-high molecular mass region indicated an absence of titin with their molecular weight (MW) more than 1000 kDa. Consequently, the present proteome analyses together with previous information on chromosome proteins provide the comprehensive list of proteins essential for the metaphase chromosome architecture.     

Sequence analysis of infectious pancreatic necrosis virus genome segment B and its encoded VP1 protein: a putative RNA-dependent RNA polymerase lacking the Gly-Asp-Asp motif

The genome segment B sequence of infectious pancreatic necrosis virus was determined for both the Jasper and Sp serotypes. The sequences are 2784 and 2630 bp long, respectively, and contain a single large open reading frame encoding the VP1 protein, the putative RNA-dependent RNA polymerase (RdRp) of IPNV. The proteins exhibit an 88% homology with each other, but only 41% with infectious bursal disease virus (IBDV) VP1, another member of the Birnaviridae. Despite the low overall homology between the IPNV and IBDV VP1 proteins, homologous regions were detected within the central portion of the proteins. The carboxy-proximal regions of the VP1, which contain very low amino acid homology, displayed evidence of conservation in structural features such as a hydrophilic, highly basic domain. Consensus sequences associated with GTP-binding proteins and RdRps were also detected in VP1. However, unlike the RdRps associated with single-stranded plus RNA viruses, the birnavirus RdRp lacks the Gly-Asp-Asp motif characteristic of this enzyme family.     

Exceptionally high levels of recombination across the honey bee genome

The first draft of the honey bee genome sequence and improved genetic maps are utilized to analyze a genome displaying 10 times higher levels of recombination (19 cM/Mb) than previously analyzed genomes of higher eukaryotes. The exceptionally high recombination rate is distributed genome-wide, but varies by two orders of magnitude. Analysis of chromosome, sequence, and gene parameters with respect to recombination showed that local recombination rate is associated with distance to the telomere, GC content, and the number of simple repeats as described for low-recombining genomes. Recombination rate does not decrease with chromosome size. On average 5.7 recombination events per chromosome pair per meiosis are found in the honey bee genome. This contrasts with a wide range of taxa that have a uniform recombination frequency of about 1.6 per chromosome pair. The excess of recombination activity does not support a mechanistic role of recombination in stabilizing pairs of homologous chromosome during chromosome pairing. Recombination rate is associated with gene size, suggesting that introns are larger in regions of low recombination and may improve the efficacy of selection in these regions. Very few transposons and no retrotransposons are present in the high-recombining genome. We propose evolutionary explanations for the exceptionally high genome-wide recombination rate.     

The common fragile site FRA16D and its associated gene WWOX are highly conserved in the mouse at Fra8E1

Recently, several common fragile sites (CFSs) have been cloned and characterized, including the two most frequently observed in the human population, FRA3B and FRA16D. In addition to their high frequency of breakage, FRA3B and FRA16D colocalize with genes crossing large regions of breakage. At FRA3B, the fragile histidine triad (FHIT) gene spans more than 1 Mb, and at FRA16D, the WWOX gene spans more than 750 kb. It has also been shown that in Mus musculus, a CFS Fra14A2 and the mouse Fhit gene are conserved in the orthologous region of the genome. In this study, we positioned the ortholog to WWOX (Wox1) at chromosome band 8E1 in the mouse genome. To determine whether, like Fra14A2 and Fhit, Fra8E1 and Wox1 colocalized in the mouse, we prepared bacterial and yeast artificial chromosome probes, and we hybridized them to aphidicolin-treated mouse metaphase chromosomes. Our data demonstrate that Wox1 colocalizes with Fra8E1. Furthermore, the sequence from this region, including introns, is highly conserved over at least a 100-kb region. This evolutionary conservation suggests that the two most active CFSs share many features, and that CFSs and their associated genes may be necessary for cell survival.     

Systems level insights into the stress response to UV radiation in the halophilic archaeon Halobacterium NRC-1

We report a remarkably high UV-radiation resistance in the extremely halophilic archaeon Halobacterium NRC-1 withstanding up to 110 J/m2 with no loss of viability. Gene knockout analysis in two putative photolyase-like genes (phr1 and phr2) implicated only phr2 in photoreactivation. The UV-response was further characterized by analyzing simultaneously, along with gene function and protein interactions inferred through comparative genomics approaches, mRNA changes for all 2400 genes during light and dark repair. In addition to photoreactivation, three other putative repair mechanisms were identified including d(CTAG) methylation-directed mismatch repair, four oxidative damage repair enzymes, and two proteases for eliminating damaged proteins. Moreover, a UV-induced down-regulation of many important metabolic functions was observed during light repair and seems to be a phenomenon shared by all three domains of life. The systems analysis has facilitated the assignment of putative functions to 26 of 33 key proteins in the UV response through sequence-based methods and/or similarities of their predicted three-dimensional structures to known structures in the PDB. Finally, the systems analysis has raised, through the integration of experimentally determined and computationally inferred data, many experimentally testable hypotheses that describe the metabolic and regulatory networks of Halobacterium NRC-1.     

Long-range sequence composition mirrors linkage disequilibrium pattern in a 1.13 Mb region of human chromosome 22.

Association studies, the most powerful tool for the identification of genes underlying complex traits, depend on the observation of linkage disequilibrium (LD) between marker alleles and the trait. The LD pattern of the human genome which determines the regional density of required markers is non-uniform, with regions of long-range LD over several hundred kilobases and regions where LD extends only over a few kilobases. Studying LD in the NF1 gene region we encountered a transition from long-range to short-range LD which coincides with a switch in the isochore pattern. This observation prompted us to investigate the regional variation in the extent of LD more systematically and we selected an isochore transition within the MN1/PITPNB gene region on chromosome 22q12.1. Long-range LD characterizes the GC-poor (40% GC) parts of the sequences. No LD can be observed between closely spaced markers throughout the whole range of the GC-rich (50% GC) parts. In both cases, the NF1 and the MN1/PITPNB gene region, a clear-cut transition of the long-range GC content precisely coincides with a change in the extent of observable LD. The results can be explained by a 72-fold lower recombination frequency in the GC-poor, compared to the GC-rich isochores. Although recombination is not the only factor governing LD, our findings can help to predict levels of LD and marker densities required for association studies on the basis of regional GC content.     

尿道致病性大肠杆菌新基因R049特征的研究

目的 明确从国内分离的尿道致病性大肠杆菌(UPEC) 132中发现的新基因R049的特征及其表达蛋白在菌体内的定位。方法 提取UPEC132染色体DNA,鸟枪法构建文库,高通量焦磷酸法测序,Newbler软件拼接,进行生物信息学分析,确定R049相关片段的特征。提取UPEC132的内膜和外膜蛋白,与其全菌裂解物一起进行SDS-PAGE,用R049重组蛋白抗血清进行Western blot,确定R049表达蛋白的菌体内位置。结果 成功构建UPEC 132染色体文库,获取R049-contig169022 bp,与UPEC536染色体第233074至451502位碱基序列同源性较高,其中含有R049基因的20773 bp片段替代相当于UPEC536染色体上PAIⅢ536的位置,G+C含量为46.97％,两端具有正向重复序列,并与插入元件整合酶基因和thrW tRNA基因相邻,包含25个ORF,命名为R049 -GI。R049基因位于R049-GI的第13个ORF。SDS-PAGE和Western blot分析显示R049基因表达相对分子质量为47.0× 103蛋白是UPEC132的外膜蛋白。结论 R049基因编码细菌的外膜蛋白,是UPEC132通过基因水平转移获得的染色体基因组岛的组成部分。     

Analysis of Arabidopsis genome-wide variations before and after meiosis and meiotic recombination by resequencing Landsberg erecta and all four products of a single meiosis

Meiotic recombination, including crossovers (COs) and gene conversions (GCs), impacts natural variation and is an important evolutionary force. COs increase genetic diversity by redistributing existing variation, whereas GCs can alter allelic frequency. Here, we sequenced Arabidopsis Landsberg erecta (Ler) and two sets of all four meiotic products from a Columbia (Col)/Ler hybrid to investigate genome-wide variation and meiotic recombination at nucleotide resolution. Comparing Ler and Col sequences uncovered 349,171 Single Nucleotide Polymorphisms (SNPs), 58,085 small and 2315 large insertions/deletions (indels), with highly correlated genome-wide distributions of SNPs, and small indels. A total of 443 genes have at least 10 nonsynonymous substitutions in protein-coding regions, with enrichment for disease-resistance genes. Another 316 genes are affected by large indels, including 130 genes with complete deletion of coding regions in Ler. Using the Arabidopsis qrt1 mutant, two sets of four meiotic products were generated and analyzed by sequencing for meiotic recombination, representing the first tetrad analysis with whole-genome sequencing in a nonfungal species. We detected 18 COs, six of which had an associated GC event, and four GCs without COs (NCOs), and revealed that Arabidopsis GCs are likely fewer and with shorter tracts than those in yeast. Meiotic recombination and chromosome assortment events dramatically redistributed genome variation in meiotic products, contributing to population diversity. In particular, meiosis provides a rapid mechanism to generate copy-number variation (CNV) of sequences that have different chromosomal positions in Col and Ler.     

Whole proteome analysis of post-translational modifications: applications of mass-spectrometry for proteogenomic annotation

While bacterial genome annotations have significantly improved in recent years, techniques for bacterial proteome annotation (including post-translational chemical modifications, signal peptides, proteolytic events, etc.) are still in their infancy. At the same time, the number of sequenced bacterial genomes is rising sharply, far outpacing our ability to validate the predicted genes, let alone annotate bacterial proteomes. In this study, we use tandem mass spectrometry (MS/MS) to annotate the proteome of Shewanella oneidensis MR-1, an important microbe for bioremediation. In particular, we provide the first comprehensive map of post-translational modifications in a bacterial genome, including a large number of chemical modifications, signal peptide cleavages, and cleavages of N-terminal methionine residues. We also detect multiple genes that were missed or assigned incorrect start positions by gene prediction programs, and suggest corrections to improve the gene annotation. This study demonstrates that complementing every genome sequencing project by an MS/MS project would significantly improve both genome and proteome annotations for a reasonable cost.     

Analysis of the complete DNA sequence of the temperate bacteriophage TP901-1: evolution, structure, and genome organization of lactococcal bacteriophages

A complete analysis of the entire genome of the temperate lactococcal bacteriophage TP901-1 has been performed and the function of 21 of 56 TP901-1-encoded ORFs has been assigned. This knowledge has been used to propose 10 functional modules each responsible for specific functions during bacteriophage TP901-1 proliferation. Short regions of microhomology in intergenic regions present in several lactococcal bacteriophages and chromosomal fragments of Lactococcus lactis are suggested to be points of exchange of genetic material through homologous recombination. Our results indicate that TP901-1 may have evolved by homologous recombination between the host chromosome and a mother phage and support the observation that phage remnants as well as prophages located in the Lactococcus chromosome contribute significantly to bacteriophage evolution. Some proteins encoded in the early transcribed region of the TP901-1 genome were more homologous to proteins encoded by phages infecting gram-positive hosts other than L. lactis. This protein homology argues for the occurrence of horizontal genetic exchange among these bacteriophages and indicates that they have access to a common gene pool.     

Genomic structure of chromosome 17 deletions in BRCA1-associated ovarian cancers

The aim of this study was to determine the genomic structure of the deletions on chromosome 17 in ovarian carcinomas from women with inherited BRCA1 mutations. Normal and tumor DNA from 14 ovarian tumors associated with inherited BRCA1 mutations were extracted and tested for loss of heterozygosity (LOH) at microsatellite markers along chromosome 17. Finer mapping using more microsatellite markers and single nucleotide polymorphisms helped further define the LOH margins. The genomic repeated elements within the LOH breakpoint regions were identified using the University of California Santa Cruz Genome Database, and the frequencies were compared to regions of equal GC percentages across the genome. Of the 14 ovarian tumors, 12 showed LOH of the entire chromosome 17. The other two tumors lost the distal end of the 17q arm. The breakpoints of these two tumors occurred in regions with significantly high frequencies of short interspersed nuclear elements (SINE), specifically Alu elements. Ovarian tumors of high grade and stage have large regions of LOH along chromosome 17, with most tumors showing loss of the entire chromosome. In those tumors with retention of part of chromosome 17, LOH margins suggest that a high Alu content may have a role in the deletions.     

Classical and Molecular Cytogenetics of the Pufferfish Tetraodon Nigroviridis

Because of its highly compact genome, the pufferfish has become an important animal model in genome research. Although the small chromosome size renders chromosome analysis difficult, we have established both classical and molecular cytogenetics in the freshwater pufferfish Tetraodon nigroviridis (TNI). The karyotype of T. nigroviridis consists of 2n = 42 biarmed chromosomes, in contrast to the known 2n = 44 chromosomes of the Japanese pufferfish Fugu rubripes (FRU). RBA banding can identify homologous chromosomes in both species. TNI 1 corresponds to two smaller FRU chromosomes, explaining the difference in chromosome number. TNI 2 is homologous to FRU 1. Fluorescence in-situ hybridization (FISH) allows one to map single-copy sequences, i.e. the Huntingtin gene, on chromosomes of the species of origin and also on chromosomes of the heterologous pufferfish species. Hybridization of total genomic DNA shows large blocks of (species-specific) repetitive sequences in the pericentromeric region of all TNI and FRU chromosomes. Hybridization with cloned human rDNA and classical silver staining reveal two large and actively transcribed rRNA gene clusters. Similar to the situation in mammals, the highly compact pufferfish genome is endowed with considerable amounts of localized repeat DNAs.