A Rickettsia genome overrun by mobile genetic elements provides insight into the acquisition of genes characteristic of an obligate intracellular lifestyle

We present the draft genome for the Rickettsia endosymbiont of Ixodes scapularis (REIS), a symbiont of the deer tick vector of Lyme disease in North America. Among Rickettsia species (Alphaproteobacteria: Rickettsiales), REIS has the largest genome sequenced to date (>2 Mb) and contains 2,309 genes across the chromosome and four plasmids (pREIS1 to pREIS4). The most remarkable finding within the REIS genome is the extraordinary proliferation of mobile genetic elements (MGEs), which contributes to a limited synteny with other Rickettsia genomes. In particular, an integrative conjugative element named RAGE (for Rickettsiales amplified genetic element), previously identified in scrub typhus rickettsiae (Orientia tsutsugamushi) genomes, is present on both the REIS chromosome and plasmids. Unlike the pseudogene-laden RAGEs of O. tsutsugamushi, REIS encodes nine conserved RAGEs that include F-like type IV secretion systems similar to that of the tra genes encoded in the Rickettsia bellii and R. massiliae genomes. An unparalleled abundance of encoded transposases (>650) relative to genome size, together with the RAGEs and other MGEs, comprise ~35% of the total genome, making REIS one of the most plastic and repetitive bacterial genomes sequenced to date. We present evidence that conserved rickettsial genes associated with an intracellular lifestyle were acquired via MGEs, especially the RAGE, through a continuum of genomic invasions. Robust phylogeny estimation suggests REIS is ancestral to the virulent spotted fever group of rickettsiae. As REIS is not known to invade vertebrate cells and has no known pathogenic effects on I. scapularis, its genome sequence provides insight on the origin of mechanisms of rickettsial pathogenicity.     

桃WRKY基因家族全基因组鉴定和表达分析

WRKY转录因子是植物中最大的转录调控家族之一,在生物和非生物胁迫以及植物生长和发育过程中起着重要调控作用.本文利用HMMER 3.0软件,使用WRKY保守域全蛋白序列(Pfam数据库编号:PF03106)鉴定桃(Prunus persica L.)基因组中的WRKY基因;利用DNAMAN 5.0,WebLogo 3,MEGA5.1,MapInspect和MEME等软件对其蛋白序列进行生物信息学分析.本文共鉴定得到61个桃WRKY基因.进化树分析结果显示,桃WRKY蛋白分为Ⅰ,Ⅱ和Ⅲ类型,类型Ⅰ分为Ⅰ-C亚组和Ⅰ-N亚组,类型Ⅱ分为Ⅱ-a,II-b,II-c,II-d和II-e亚组.WRKY结构域分析显示,WRKY结构域高度保守,绝大多数都含有WRKYGQK七肽和锌指结构.染色体定位分析显示,桃WRKY基因分布于8条染色体中,呈不均匀分布.内含子和外显子结构分析表明,WRKY基因结构进化高度保守.保守元件分析表明,桃WRKY基因家族包含5个保守元件,元件1,2和3为WRKY盒,元件4,5为未知盒.桃WRKY基因家族都包含有WRKY盒,类型Ⅰ中含有2个WRKY盒;II-d亚组中含有未知元件5.半定量和荧光定量PCR结果显示,16个WRKY基因均在桃的根,茎,叶,花和果中表达,但其相对表达水平不同.     

Rapid genome evolution revealed by comparative sequence analysis of orthologous regions from four triticeae genomes

Bread wheat (Triticum aestivum) is an allohexaploid species, consisting of three subgenomes (A, B, and D). To study the molecular evolution of these closely related genomes, we compared the sequence of a 307-kb physical contig covering the high molecular weight (HMW)-glutenin locus from the A genome of durum wheat (Triticum turgidum, AABB) with the orthologous regions from the B genome of the same wheat and the D genome of the diploid wheat Aegilops tauschii (Anderson et al., 2003; Kong et al., 2004). Although gene colinearity appears to be retained, four out of six genes including the two paralogous HMW-glutenin genes are disrupted in the orthologous region of the A genome. Mechanisms involved in gene disruption in the A genome include retroelement insertions, sequence deletions, and mutations causing in-frame stop codons in the coding sequences. Comparative sequence analysis also revealed that sequences in the colinear intergenic regions of these different genomes were generally not conserved. The rapid genome evolution in these regions is attributable mainly to the large number of retrotransposon insertions that occurred after the divergence of the three wheat genomes. Our comparative studies indicate that the B genome diverged prior to the separation of the A and D genomes. Furthermore, sequence comparison of two distinct types of allelic variations at the HMW-glutenin loci in the A genomes of different hexaploid wheat cultivars with the A genome locus of durum wheat indicates that hexaploid wheat may have more than one tetraploid ancestor.     

Isolation,Cloning and Characterization of Resistance Gene Analogues in Pearl Millet Based on Conserved Nucleotide‐binding Sites

Sequence analysis of plant disease resistance genes shows similarity among themselves, with the presence of conserved motifs common to the nucleotide‐binding site (NBS). Oligonucleotide degenerate primers designed from the conserved NBS motifs encoded by several plant disease resistance genes were used to amplify resistance gene analogues (RGAs) corresponding to the NBS sequences from the genomic DNA of various plant species. Using specific primers designed from the conserved NBS regions, 22 RGAs were cloned and sequenced from pearl millet (Pennisetum glaucum L. Br.). Phylogenetic analysis of the predicted amino acid sequences grouped the RGAs into nine distinct classes. GenBank database searches with the consensus protein sequences of each of the nine classes revealed their conserved NBS domains and similarity to other known R genes of various crop species. One RGA 213 was mapped onto LG1 and LG7 in the pearl millet linkage map. This is the first report of the isolation and characterization of RGAs from pearl millet, which will facilitate the improvement of marker‐assisted breeding strategies.     

Insights into the completely annotated genome of Lactobacillus buchneri CD034, a strain isolated from stable grass silage

Lactobacillus buchneri belongs to the group of heterofermentative lactic acid bacteria and is a common member of the silage microbiome. Here we report the completely annotated genomic sequence of L. buchneri CD034, a strain isolated from stable grass silage. The whole genome of L. buchneri CD034 was sequenced on the Roche Genome Sequencer FLX platform. It was found to consist of four replicons, a circular chromosome, and three plasmids. The circular chromosome was predicted to encode 2319 proteins and contains a genomic island and two prophages which significantly differ in G+C-content from the remaining chromosome. It possesses all genes for enzymes of a complete phosphoketolase pathway, whereas two enzymes necessary for glycolysis are lacking. This confirms the classification of L. buchneri CD034 as an obligate heterofermentative lactic acid bacterium. A set of genes considered to be involved in the lactate degradation pathway and genes putatively involved in the breakdown of plant cell wall polymers were identified. Moreover, several genes encoding putative S-layer proteins and two CRISPR systems, belonging to the subclasses I-E and II-A, are located on the chromosome. The largest plasmid pCD034-3 was predicted to encode 57 genes, including a putative polysaccharide synthesis gene cluster, whereas the functions of the two smaller plasmids, pCD034-1 and pCD034-2, remain cryptic. Phylogenetic analysis based on sequence comparison of the conserved marker gene rpoA reveals that L. buchneri CD034 is more closely related to Lactobacillus hilgardii strains than to Lactobacillus brevis and Lactobacillus plantarum strains. Comparison of the L. buchneri CD034 core genome to other fully sequenced and closely related members of the genus Lactobacillus disclosed a high degree of conservation between L. buchneri CD034 and the recently sequenced L. buchneri strain NRRL B-30929 and a more distant relationship to L. buchneri ATCC 11577 and L. brevis ssp. gravesensis ATCC 27305, which cluster together with L. hilgardii type strain ATCC 8290. L. buchneri CD034 genome information will certainly provide the basis for further postgenome studies with the objective to optimize application of the strain in silage production.     

Plastid translation and transcription genes in a non-photosynthetic plant: intact, missing and pseudo genes

The non-photosynthetic, parasitic flowering plant Epifagus virginiana has recently been shown to contain a grossly reduced plastid genome that has lost many photosynthetic and chloro-respiratory genes. We have cloned and sequenced a 3.9 kb domain of plastid DNA from Epifagus to investigate the patterns of evolutionary change in such a reduced genome and to determine which genes are still present and likely to be functional. This 3.9 kb domain is colinear with a 35.4 kb region of tobacco chloroplast DNA, differing from it by a minimum of 11 large deletions varying in length from 354 bp to 11.5 kb, as well as by a number of small deletions and insertions. The nine genes retained in Epifagus encode seven tRNAs and two ribosomal proteins and are coextensive and highly conserved in sequence with homologs in photosynthetic plants. This suggests that these genes are functional in Epifagus and, together with evidence that the Epifagus plastid genome is transcribed, implies that plastid gene products play a role in processes other than photosynthesis and gene expression. Genes that are completely absent include not only photosynthetic genes, but surprisingly, genes encoding three subunits of RNA polymerase, four tRNAs and one ribosomal protein. In addition, only pseudogenes are found for two other tRNAs. Despite these defunct tRNA genes, codon and amino acid usage in Epifagus protein genes is normal. We therefore hypothesize that the expression of plastid genes in Epifagus relies on the import of nuclear encoded tRNAs and RNA polymerase from the cytoplasm.     

Genetic variation and virulence of Autographa californica multiple nucleopolyhedrovirus and Trichoplusia ni single nucleopolyhedrovirus isolates

To determine the genetic diversity within the baculovirus species Autographa calfornica multiple nucleopolyhedrovirus (AcMNPV; Baculoviridae: Alphabaculovirus), a PCR-based method was used to identify and classify baculoviruses found in virus samples from the lepidopteran host species A. californica, Autographa gamma, Trichoplusia ni, Rachiplusia ou, Anagrapha falcifera, Galleria mellonella, and Heliothis virescens. Alignment and phylogenetic inference from partial nucleotide sequences of three highly conserved genes (lef-8, lef-9, and polh) indicated that 45 of 74 samples contained isolates of AcMNPV, while six samples contained isolates of Rachiplusia ou multiple nucleopolyhedrovirus strain R1 (RoMNPV-R1) and 25 samples contained isolates of the species Trichoplusia ni single nucleopolyhedrovirus (TnSNPV; Alphabaculovirus). One sample from A. californica contained a previously undescribed NPV related to alphabaculoviruses of the armyworm genus Spodoptera. Data from PCR and sequence analysis of the ie-2 gene and a region containing ORF ac86 in samples from the AcMNPV and RoMNPV clades indicated a distinct group of viruses, mostly from G. mellonella, that are characterized by an unusual ie-2 gene previously found in the strain Plutella xylostella multiple nucleopolyhedrovirus CL3 (PlxyMNPV-CL3) and a large deletion within ac86 previously described in the AcMNPV isolate 1.2 and PlxyMNPV-CL3. PCR and sequence analysis of baculovirus repeated ORF (bro) genes revealed that the bro gene ac2 was split into two separate bro genes in some samples from the AcMNPV clade. Comparison of sequences in this region suggests that ac2 was formed by a deletion that fused the two novel bro genes together. In bioassays of a selection of isolates against T. ni, significant differences were observed in the insecticidal properties of individual isolates, but no trends were observed among the AcMNPV, TnSNPV, or RoMNPV groups of isolates. This study expands on what we know about the variation of AcMNPV, AcMNPV-like and TnSNPV viruses, provides novel information on the distinct groups in which AcMNPV isolates occur, and contributes to data useful for the registration, evaluation, and improvement of AcMNPV, AcMNPV-like, and TnSNPV isolates as biological control agents.     

Insertions of mini-Tn10 transposon T-POP in Salmonella enterica sv. typhi

We have mutagenized a clinical strain of Salmonella enterica sv. typhi with mini-transposon Tn10dTet (T-POP) to obtain conditional lethal (tetracycline-dependent) mutants with T-POP insertions upstream of essential genes. Generalized transducing phage P22 was used to introduce T-POP from a S. typhimurium donor into a S. typhi recipient. Chromosomal DNA was purified from the mutagenized donor strains, fragmented, and then electroporated into S. typhi to backcross the original T-POP insertions. Four tetracycline-dependent mutants with two distinct terminal phenotypes were found among 1700 mutants with T-POP insertions. When grown in the absence of tetracycline, two of the four tetracycline-dependent mutants arrest at a late stage in the cell cycle, can be rescued by outgrowth in media with tetracycline, and define a reversible checkpoint late in the cell cycle. One of these insertions creates an operon fusion with a gene, yqgF, that is conserved among gram-negative bacteria and likely encodes an essential Holliday junction resolvase. T-POP insertions can be used not only to identify essential S. typhi genes but also to reveal novel phenotypes resulting from the depletion of their products.     

A chromosome‐level genome assembly of the woolly apple aphid,Eriosoma lanigerum Hausmann (Hemiptera: Aphididae)

Woolly apple aphid (WAA, Eriosoma lanigerum Hausmann) (Hemiptera: Aphididae) is a major pest of apple trees (Malus domestica, order Rosales) and is critical to the economics of the apple industry in most parts of the world. Here, we generated a chromosome‐level genome assembly of WAA—representing the first genome sequence from the aphid subfamily Eriosomatinae—using a combination of 10X Genomics linked‐reads and in vivo Hi‐C data. The final genome assembly is 327 Mb, with 91% of the assembled sequences anchored into six chromosomes. The contig and scaffold N50 values are 158 kb and 71 Mb, respectively, and we predicted a total of 28,186 protein‐coding genes. The assembly is highly complete, including 97% of conserved arthropod single‐copy orthologues based on Benchmarking Universal Single‐Copy Orthologs (busco ) analysis. Phylogenomic analysis of WAA and nine previously published aphid genomes, spanning four aphid tribes and three subfamilies, reveals that the tribe Eriosomatini (represented by WAA) is recovered as a sister group to Aphidini + Macrosiphini (subfamily Aphidinae). We identified syntenic blocks of genes between our WAA assembly and the genomes of other aphid species and find that two WAA chromosomes (El5 and El6) map to the conserved Macrosiphini and Aphidini X chromosome. Our high‐quality WAA genome assembly and annotation provides a valuable resource for research in a broad range of areas such as comparative and population genomics, insect–plant interactions and pest resistance management.     

Characterization of a BAC Library from Channel Catfish <Emphasis Type="Italic">Ictalurus punctatus</Emphasis>: Indications of High Levels of Chromosomal Reshuffling Among Teleost Genomes

The CHORI-212 bacterial artificial chromosome (BAC) library was constructed by cloning EcoRI/EcoRI partially digested DNA into the pTARBAC2.1 vector. The library has an average insert size of 161 kb, and provides 10.6-fold coverage of the channel catfish haploid genome. Screening of 32 genes using overgo or cDNA probes indicated that this library had a good representation of the genome as all tested genes existed in the library. We previously reported sequencing of approximately 25,000 BAC ends that generated 20,366 high-quality BAC end sequences (BES) and identified a large number of sequences similar to known genes using BLASTX searches. In this work, particular attention was given to identification of BAC mate pairs with known genes from both ends. When identified, comparative genome analysis was conducted to determine syntenic regions of the catfish genome with the genomes of zebrafish and Tetraodon. Of the 141 mate pairs with known genes from channel catfish, conserved syntenies were identified in 34 (24.1%), with 30 conserved in the zebrafish genome and 14 conserved in the Tetraodon genome. Additional analysis of three of the 34 conserved syntenic groups by direct sequencing indicated conserved gene contents in all three species. This indicates that comparative genome analysis may provide shortcuts to genome analysis in catfish, especially for short genomic regions once the conserved syntenies are identified. Shaolin Wang and Peng Xu contributed equally to the article.     

Comparative genome mapping of Pseudomonas aeruginosa PAO with P. aeruginosa C, which belongs to a major clone in cystic fibrosis patients and aquatic habitats.

A physical and genetic map was constructed for Pseudomonas aeruginosa C. Mainly, two-dimensional methods were used to place 47 SpeI, 8 PacI, 5 SwaI, and 4 I-CeuI sites onto the 6.5-Mb circular chromosome. A total of 21 genes, including the rrn operons and the origin of replication, were located on the physical map. Comparison of the physical and genetic map of strain C with that of the almost 600-kb-smaller genome of P. aeruginosa reference strain PAO revealed conservation of gene order between the two strains. A large-scale mosaic structure which was due to insertions of blocks of new genetic elements which had sizes of 23 to 155 kb and contained new SpeI sites was detected in the strain C chromosome. Most of these insertions were concentrated in three locations: two are congruent with the ends of the region rich in biosynthetic genes, and the third is located in the proposed region of the replication terminus. In addition, three insertions were scattered in the region rich in biosynthetic genes. The arrangement of the rrn operons around the origin of replication was conserved in C, PAO, and nine other examined independent strains.     

Comparative genomics of the T4-Like Escherichia coli phage JS98: implications for the evolution of T4 phages

About 130 kb of sequence information was obtained from the coliphage JS98 isolated from the stool of a pediatric diarrhea patient in Bangladesh. The DNA shared up to 81% base pair identity with phage T4. The most conserved regions between JS98 and T4 were the structural genes, but their degree of conservation was not uniform. The head genes showed the highest sequence conservation, followed by the tail, baseplate, and tail fiber genes. Many tail fiber genes shared only protein sequence identity. Except for the insertion of endonuclease genes in T4 and gene 24 duplication in JS98, the structural gene maps of the two phages were colinear. The receptor-recognizing tail fiber proteins gp37 and gp38 were only distantly related to T4, but shared up to 83% amino acid identity to other T6-like phages, suggesting lateral gene transfer. A greater degree of variability was seen between JS98 and T4 over DNA replication and DNA transaction genes. While most of these genes came in the same order and shared up to 76% protein sequence identity, a few rearrangements, insertions, and replacements of genes were observed. Many putative gene insertions in the DNA replication module of T4 were flanked by intron-related endonuclease genes, suggesting mobile DNA elements. A hotspot of genome diversification was located downstream of the DNA polymerase gene 43 and the DNA binding gene 32. Comparative genomics of 100-kb genome sequence revealed that T4-like phages diversify more by the accumulation of point mutations and occasional gene duplication events than by modular exchanges.