Biology of <Emphasis Type="Italic">B. sorokiniana</Emphasis> (syn. <Emphasis Type="Italic">Cochliobolus sativus</Emphasis>) in genomics era

Bipolaris sorokiniana (Sacc.) Shoemaker is a hemi-biotrophic fungal pathogen, which is an anamorph (teleomorph Cochlibolus sativus). It causes spot blotch, root rot and leaf spot diseases in a number of cereals including wheat, barley and other small grain cereals. In the genomics era, the fungus has been subjected to a variety of studies using molecular approaches. Correct chromosome number was determined and molecular karyotypes were prepared using contour-clamped homogeneous electric field. Molecular maps were prepared using markers like RFLPs, SSRs, RAPDs and SNPs. For this purpose, segregating progenies derived from crosses between diverse isolates of the pathogen were used. Whole genome sequencing (WGS) data was collected not only for B. sorokiniana isolates, but also for several species of Cochliobolus. Genes involved in secondary metabolism and virulence were identified from genome sequences. The WGS data has also been utilized for comparative genomics giving useful information about evolutionary trends. A brief account of this information is presented in this review.     

Comparative genomics of the social amoebae <Emphasis Type="Italic">Dictyostelium discoideum</Emphasis> and <Emphasis Type="Italic">Dictyostelium purpureum</Emphasis>

Background

The social amoebae (Dictyostelia) are a diverse group of Amoebozoa that achieve multicellularity by aggregation and undergo morphogenesis into fruiting bodies with terminally differentiated spores and stalk cells. There are four groups of dictyostelids, with the most derived being a group that contains the model species Dictyostelium discoideum.

Results

We have produced a draft genome sequence of another group dictyostelid, Dictyostelium purpureum, and compare it to the D. discoideum genome. The assembly (8.41 × coverage) comprises 799 scaffolds totaling 33.0 Mb, comparable to the D. discoideum genome size. Sequence comparisons suggest that these two dictyostelids shared a common ancestor approximately 400 million years ago. In spite of this divergence, most orthologs reside in small clusters of conserved synteny. Comparative analyses revealed a core set of orthologous genes that illuminate dictyostelid physiology, as well as differences in gene family content. Interesting patterns of gene conservation and divergence are also evident, suggesting function differences; some protein families, such as the histidine kinases, have undergone little functional change, whereas others, such as the polyketide synthases, have undergone extensive diversification. The abundant amino acid homopolymers encoded in both genomes are generally not found in homologous positions within proteins, so they are unlikely to derive from ancestral DNA triplet repeats. Genes involved in the social stage evolved more rapidly than others, consistent with either relaxed selection or accelerated evolution due to social conflict.

Conclusions

The findings from this new genome sequence and comparative analysis shed light on the biology and evolution of the Dictyostelia.

     

Genome-wide identification and evolution of <Emphasis Type="Italic">TC1</Emphasis>/<Emphasis Type="Italic">Mariner</Emphasis> in the silkworm (<Emphasis Type="Italic">Bombyx mori</Emphasis>) genome

TC1/Mariner transposons belong to class II transposable elements (TEs) that use DNA-mediated “cut and paste” mechanism to transpose, and they have been identified in almost all organisms. Although silkworm (Bombyx mori) has a large amount of TC1/Mariner elements, the genome wide information of this superfamily in the silkworm is unknown. In this study, we have identified 2670 TC1/Mariner (Bmmar) elements in the silkworm genome. All the TEs were classified into 22 families by means of fgclust, a tool of repetitive sequence classification, seven of which was first reported in this study. Phylogenetic and structure analyses based on the catalytic domain (DDxD/E) of transposase sequences indicated that all members of TC1/Mariner were grouped into five subgroups: Mariner, Tc1, maT, DD40D and DD41D/E. Of these five subgroups, maT rather than Mariner possessed most members of TC1/Mariner (51.23%) in the silkworm genome. In particular, phylogenetic analysis and structure analysis revealed that Bmmar15 (DD40D) formed a new basal subgroup of TC1/Mariner element in insects, which was referred to as bmori. Furthermore, we concluded that DD40D appeared to intermediate between mariner and Tc1. Finally, we estimated the insertion time for each copy of TC1/Mariner in the silkworm and found that most of members were dramatically amplified during a period from 0 to 1 mya. Moreover, the detailed functional data analysis showed that Bmmar1, Bmmar6 and Bmmar9 had EST evidence and intact transposases. These implied that TC1/Mariner might have potential transpositional activity. In conclusion, this study provides some new insights into the landscape, origin and evolution of TC1/Mariner in the insect genomes.     

Complete chloroplast genome sequence of <Emphasis Type="Italic">Capsicum</Emphasis><Emphasis Type="Italic">baccatum</Emphasis> var. <Emphasis Type="Italic">baccatum</Emphasis>

Molecular markers derived from the complete chloroplast genome can provide effective tools for species identification and phylogenetic resolution. Complete chloroplast (cp) genome sequences of Capsicum species have been reported. We herein report the complete chloroplast genome sequence of Capsicum baccatum var. baccatum, a wild Capsicum species. The total length of the chloroplast genome is 157,145 bp with 37.7 % overall GC content. One pair of inverted repeats, 25,910 bp in length, was separated by a small single-copy region (17,974 bp) and large single-copy region (87,351 bp). This region contains 86 protein-coding genes, 30 tRNA genes, 4 rRNA genes, and 11 genes contain one or two introns. Pair-wise alignments of chloroplast genome were performed for genome-wide comparison. Analysis revealed a total of 134 simple sequence repeat (SSR) motifs and 282 insertions or deletions variants in the C. baccatum var. baccatum cp genome. The types and abundances of repeat units in Capsicum species were relatively conserved, and these loci could be used in future studies to investigate and conserve the genetic diversity of the Capsicum species.     

Allele-specific CAPS marker in a <Emphasis Type="Italic">Ve1</Emphasis> homolog of <Emphasis Type="Italic">Capsicum annuum</Emphasis> for improved selection of <Emphasis Type="Italic">Verticillium dahliae</Emphasis> resistance

Verticillium wilt (Verticillium dahliae) is an economically important disease for many high-value crops. The pathogen is difficult to manage due to the long viability of its resting structures, wide host range, and the inability of fungicides to affect the pathogen once in the plant vascular system. In chile pepper (Capsicum annuum), breeding for resistance to Verticillium wilt is especially challenging due to the limited resistance sources. The dominant Ve locus in tomato (Solanum lycopersicum) contains two closely linked and inversely oriented genes, Ve1 and Ve2. Homologs of Ve1 have been characterized in diverse plant species, and interfamily transfer of Ve1 confers race-specific resistance. Queries in the chile pepper WGS database in NCBI with Ve1 and Ve2 sequences identified one open reading frame (ORF) with homology to the tomato Ve genes. Comparison of the candidate CaVe (Capsicum annuum Ve) gene sequences from susceptible and resistant accessions revealed 16 single nucleotide polymorphisms (SNPs) and several haplotypes. A homozygous haplotype was identified for the susceptible accessions and for resistant accessions. We developed a cleaved amplified polymorphic sequence (CAPS) molecular marker within the coding region of CaVe and screened diverse germplasm that has been previously reported as being resistant to Verticillium wilt in other regions. Based on our phenotyping using the New Mexico V. dahliae isolate, the marker could select resistance accessions with 48% accuracy. This molecular marker is a promising tool towards marker-assisted selection for Verticillium wilt resistance and has the potential to improve the efficacy of chile pepper breeding programs, but does not eliminate the need for a bioassay. Furthermore, this work provides a basis for future research in this important pathosystem.     

Horizontal transfer of the C-termini of <Emphasis Type="Italic">tccC</Emphasis> genes in <Emphasis Type="Italic">Photorhabdus</Emphasis> and <Emphasis Type="Italic">Xenorhabdus</Emphasis>

The high molecular weight insecticidal toxin complexes (Tcs), including four toxin-complex loci (tca, tcb, tcc and tcd), were first identified in Photorhabdus luminescens W14. Each member of tca, tcb or tcc is required for oral toxicity of Tcs. However, the sequence sources of the C-termini of tccC3, tccC4, tccC6 and tccC7 are unknown. Here, we performed a whole genome survey to identify the orthologs of Tc genes, and found 165 such genes in 14 bacterial genomes, including 40 genes homologous to tccC1-7 in P. luminescens TT01. The sequence sources of the C-termini of tccC2-6 were determined by sequence analysis. Further phylogenetic investigations suggested that the C-termini of 6 tccC genes experienced horizontal gene transfer events.     

Structural organization characteristics of the <Emphasis Type="Italic">DIP1</Emphasis> gene in <Emphasis Type="Italic">Drosophila melanogaster</Emphasis> strains mutant for the <Emphasis Type="Italic">flamenco</Emphasis> gene

Molecular cloning of the DIP1 gene located in the 20A4-5 region has been performed from the following strains with the flamenco phenotype: flam ^SS (SS) and flam ^MS (MS) characterized by a high transposition rate of retrotransposon gypsy (mdg4), flam ^{py + (P)} carrying the insertion of a construction based on the P element into the region of the flamenco gene, and flamenco ⁺. The results of restriction analysis and sequencing cloned DNA fragments has shown that strains flam ^SS , flam ^MS considerably differ from flam ^{py + (P)}, and flamenco ⁺ in the structure of DIP1. Strains flam ^SS and flam ^MS have no DraI restriction site at position 1765 in the coding region of the gene, specifically, in the domain determining the signal of the nuclear localization of the DIP1 protein. This mutation has been found to consist in a nucleotide substitution in the recognition site of DraI restriction endonuclease, which is transformed from TTTAAA into TTTAAG and, hence, is not recognized by the enzyme. This substitution changes codon AAA into AAG and is translationally insignificant, because both triplets encode the same amino acid, lysine. The DIP1 gene of strains flam ^SS and flam ^MS has been found to contain a 182-bp insertion denoted IdSS (insertion in DIP1 strain SS); it is located in the second intron of the gene. The IdSS sequence is part of the open reading frame encoding the putative transposase of the mobile genetic element HB1 belonging to the Tc1/mariner family. This insertion is presumed to disturb the conformations of DNA and the chromosome, in particular, by forming loops, which alters the expression of DIP1 and, probably, neighboring genes. In strains flamenco ⁺ and flam ^{py + (P)}, the IdSS insertion within the HB1 sequence is deleted. The deletion encompasses five C-terminal amino acid residues of the conserved domain and the entire C-terminal region of the putative HB1 transposase. The obtained data suggest that DIP1 is involved in the control of gypsy transpositions either directly or through interaction with other elements of the genome.     

Evaluation of the Diversity of Probiotic <Emphasis Type="Italic">Bacillus</Emphasis>, <Emphasis Type="Italic">Clostridium</Emphasis>, and <Emphasis Type="Italic">Bifidobacterium</Emphasis> Using the Illumina-Based Sequencing Method

Bacterial species of Bacillus, Lactobacillus, and Bifidobacterium in the intestinal tract have been used as probiotics. Selections for probiotic candidates by the culture-based approaches are time-consuming and labor-consuming. The aim of this study was to develop a new method based on sequencing strategies to select the probiotic Bacillus, Lactobacillus, and Bifidobacterium. The Illumina-based sequencing strategies with different specific primers for Bacillus, Clostridium, and Bifidobacterium were applied to analyze diversity of the genera in goat feces. The average number of different Bacillus, Clostridium, and Bifidobacterium OTUs (operational taxonomic units) at the 97% similarity level ranged from 1922 to 63172. The coverage index values of Bacillus, Clostridium, and Bifidobacterium calculated from the bacterial OTUs were 0.89, 0.99, and 1.00, respectively. The most genera of Bacillus (37.9%), Clostridium (53%), and Bifidobacterium (99%) were detected in goat feces by the Illumina-based sequencing with the specific primers of the genera, respectively. Higher phylogenetic resolutions of the genera in goat feces were successfully established. The results suggest that the selection for probiotic Bacillus, Clostridium, and Bifidobacterium based on the Illumina sequencing with their specific primers is reliable and feasible, and the core Bacillus, Clostridium, and Bifidobacterium species of healthy goats possess the potentials as probiotic microbial consortia.     

<Emphasis Type="Italic">HB</Emphasis> mobile element in the <Emphasis Type="Italic">Drosophila melanogaster</Emphasis> genome: Structural and functional analyses

The structure was analyzed for 60 annotated copies of the mobile genetic element (MGE) HB from the Drosophila melanogaster genome. The genomic distribution of HB copies was studied, and preferential insertion sites (hot spots) were identified, which presumably amount to several kilobases. Structural analysis of the open reading frame (ORF) and terminal repeats of HB was performed. All 26 HB copies retaining the ORF sequence have a stop codon in the same position. Consequently, the HB ORF proved indeed to code for an enzyme of 148 amino acid residues, relatively small for Tc1-family transposases. The ORF consensus sequence was established. HB{}1185 was identified as the only HB copy potentially coding for a functional protein. All 37 repeat-containing HB copies were analyzed. Of these, only four had functional terminal sequences, lacking, however, a functional transposase gene. A new 7762-bp copy of MGE roo was found in the D. melanogaster genome; the copy was earlier unavailable from databases and represents an insert in the HB{}1605 sequence.     

New resources for functional analysis of omics data for the genus <Emphasis Type="Italic">Aspergillus</Emphasis>

Background

Detailed and comprehensive genome annotation can be considered a prerequisite for effective analysis and interpretation of omics data. As such, Gene Ontology (GO) annotation has become a well accepted framework for functional annotation. The genus Aspergillus comprises fungal species that are important model organisms, plant and human pathogens as well as industrial workhorses. However, GO annotation based on both computational predictions and extended manual curation has so far only been available for one of its species, namely A. nidulans.

Results

Based on protein homology, we mapped 97% of the 3,498 GO annotated A. nidulans genes to at least one of seven other Aspergillus species: A. niger, A. fumigatus, A. flavus, A. clavatus, A. terreus, A. oryzae and Neosartorya fischeri. GO annotation files compatible with diverse publicly available tools have been generated and deposited online. To further improve their accessibility, we developed a web application for GO enrichment analysis named FetGOat and integrated GO annotations for all Aspergillus species with public genome sequences. Both the annotation files and the web application FetGOat are accessible via the Broad Institute's website (http://www.broadinstitute.org/fetgoat/index.html). To demonstrate the value of those new resources for functional analysis of omics data for the genus Aspergillus, we performed two case studies analyzing microarray data recently published for A. nidulans, A. niger and A. oryzae.

Conclusions

We mapped A. nidulans GO annotation to seven other Aspergilli. By depositing the newly mapped GO annotation online as well as integrating it into the web tool FetGOat, we provide new, valuable and easily accessible resources for omics data analysis and interpretation for the genus Aspergillus. Furthermore, we have given a general example of how a well annotated genome can help improving GO annotation of related species to subsequently facilitate the interpretation of omics data.

     

High-throughput sequencing of <Emphasis Type="Italic">Prunus ferganensis</Emphasis> indicates that it is a geographical population of <Emphasis Type="Italic">P. persica</Emphasis>

Peach belongs to the genus Prunus, which includes Prunus persica and its relative species, P. mira, P. davidiana, P. kansuensis, and P. ferganensis. Of these, P. ferganensis have been classified as a species, subspecies, or geographical population of P. persica. To explore the genetic difference between P. ferganensis and P. persica, high-throughput sequencing was used in different peach accessions belonging to different species. First, low-depth sequencing data of peach accessions belonging to four categories revealed that similarity between P. ferganensis and P. persica was similar to that between P. persica accessions from different geographical populations. Then, to further detect the genomic variation in P. ferganensis, the P. ferganensis accession “Xinjiang Pan Tao 1” and the P. persica accession “Xia Miao 1” were sequenced with high depth, and sequence reads were assembled. The results showed that the collinearity of “Xinjiang Pan Tao 1” with the reference genome “Lovell” was higher than that of “Xia Miao 1” and “Lovell” peach. Additionally, the number of genetic variants, including single nucleotide polymorphisms (SNPs), structural variations (SVs), and the specific genes annotated from unmapped sequence in “Xia Miao 1” was higher than that in “Xinjiang Pan Tao 1” peach. The data showed that there was a close distance between “Xinjiang Pan Tao 1” (P. ferganensis) and reference genome which belong to P. persica, comparing “Xia Miao 1” (P. persica) and reference ones. The results accompany with phylogenetic tree and structure analysis confirmed that P. ferganensis should be considered as a geographic population of P. persica rather than a subspecies or a distinct species. Furthermore, gene ontology analysis was performed using the gene comprising large-effect variation to understand the phenotypic difference between two accessions. The result revealed that the pathways of gene function affected by SVs but SNPs and insertion-deletions markedly differed between the two peach accessions.     

The single-stranded DNA-binding protein of <Emphasis Type="Italic">Deinococcus radiodurans</Emphasis>

Background

Deinococcus radiodurans R1 is one of the most radiation-resistant organisms known and is able to repair an unusually large amount of DNA damage without induced mutation. Single-stranded DNA-binding (SSB) protein is an essential protein in all organisms and is involved in DNA replication, recombination and repair. The published genomic sequence from Deinococcus radiodurans includes a putative single-stranded DNA-binding protein gene (ssb; DR0100) requiring a translational frameshift for synthesis of a complete SSB protein. The apparently tripartite gene has inspired considerable speculation in the literature about potentially novel frameshifting or RNA editing mechanisms. Immediately upstream of the ssb gene is another gene (DR0099) given an ssb-like annotation, but left unexplored.

Results

A segment of the Deinococcus radiodurans strain R1 genome encompassing the ssb gene has been re-sequenced, and two errors involving omitted guanine nucleotides have been documented. The corrected sequence incorporates both of the open reading frames designated DR0099 and DR0100 into one contiguous ssb open reading frame (ORF). The corrected gene requires no translational frameshifts and contains two predicted oligonucleotide/oligosaccharide-binding (OB) folds. The protein has been purified and its sequence is closely related to the Thermus thermophilus and Thermus aquaticus SSB proteins. Like the Thermus SSB proteins, the SSB_Dr functions as a homodimer. The Deinococcus radiodurans SSB homodimer stimulates Deinococcus radiodurans RecA protein and Escherichia coli RecA protein-promoted DNA three-strand exchange reactions with at least the same efficiency as the Escherichia coli SSB homotetramer.

Conclusions

The correct Deinococcus radiodurans ssb gene is a contiguous open reading frame that codes for the largest bacterial SSB monomer identified to date. The Deinococcus radiodurans SSB protein includes two OB folds per monomer and functions as a homodimer. The Deinococcus radiodurans SSB protein efficiently stimulates Deinococcus radiodurans RecA and also Escherichia coli RecA protein-promoted DNA strand exchange reactions. The identification and purification of Deinococcus radiodurans SSB protein not only allows for greater understanding of the SSB protein family but provides an essential yet previously missing player in the current efforts to understand the extraordinary DNA repair capacity of Deinococcus radiodurans.

     

<Emphasis Type="Italic">BcMF11</Emphasis> and its homologous sequences may form a lncRNA family in <Emphasis Type="Italic">Brassica</Emphasis> diploids

BcMF11 is a long non-coding RNA that has been identified in Brassica rapa and shown to be involved in pollen development. Here, when re-cloned the gene sequence, multiple paralogous copies of BcMF11 were identified in B. rapa (A genome). Multiple paralogous copies of BcMF11 were also found in B. nigra (B genome) and Brassica oleracea (C genome), the other two primary diploids of Brassica U triangle. While in the early diverging Brassicaceae lineage including Arabidopsis thaliana, no BcMF11 homolog was found. Phylogenetic analysis showed that the BcMF11 homologous sequences cloned from A genome or C genome could be clustered into a separate branch, respectively. However, there was no distinct cluster defined for BcMF11 homologous sequences cloned from B genome. The expression of BcMF11 in B. rapa was investigated and revealed a different result in the previous study. In addition, 12 expressed sequence tags from B. napus and B. rapa showing high similarities with BcMF11 were identified in the NCBI database, which further verified that rather than the useless repeat fragments in the genome, the BcMF11 homologous genes could transcribe. It is possible that BcMF11 and its homologous sequences may form a large gene family which might be originated in the recent ancestral lineage of Brassica.     

A <Emphasis Type="Italic">bean common mosaic virus</Emphasis> (BCMV)-resistance gene is fine-mapped to the same region as <Emphasis Type="Italic">Rsv1</Emphasis>-<Emphasis Type="Italic">h</Emphasis> in the soybean cultivar Suweon 97

Key message

In the soybean cultivar Suweon 97, BCMV-resistance gene was fine-mapped to a 58.1-kb region co-localizing with the Soybean mosaic virus (SMV)-resistance gene, Rsv1-h raising a possibility that the same gene is utilized against both viral pathogens.

Abstract

Certain soybean cultivars exhibit resistance against soybean mosaic virus (SMV) or bean common mosaic virus (BCMV). Although several SMV-resistance loci have been reported, the understanding of the mechanism underlying BCMV resistance in soybean is limited. Here, by crossing a resistant cultivar Suweon 97 with a susceptible cultivar Williams 82 and inoculating 220 F₂ individuals with a BCMV strain (HZZB011), we observed a 3:1 (resistant/susceptible) segregation ratio, suggesting that Suweon 97 possesses a single dominant resistance gene against BCMV. By performing bulked segregant analysis with 186 polymorphic simple sequence repeat (SSR) markers across the genome, the resistance gene was determined to be linked with marker BARSOYSSR_13_1109. Examining the genotypes of nearby SSR markers on all 220 F₂ individuals then narrowed down the gene between markers BARSOYSSR_13_1109 and BARSOYSSR_13_1122. Furthermore, 14 previously established F_2:3 lines showing crossovers between the two markers were assayed for their phenotypes upon BCMV inoculation. By developing six more SNP (single nucleotide polymorphism) markers, the resistance gene was finally delimited to a 58.1-kb interval flanked by BARSOYSSR_13_1114 and SNP-49. Five genes were annotated in this interval of the Williams 82 genome, including a characteristic coiled-coil nucleotide-binding site-leucine-rich repeat (CC-NBS-LRR, CNL)-type of resistance gene, Glyma13g184800. Coincidentally, the SMV-resistance allele Rsv1-h was previously mapped to almost the same region, thereby suggesting that soybean Suweon 97 likely relies on the same CNL-type R gene to resist both viral pathogens.

     

Map-based cloning of the dominant genic male sterile <Emphasis Type="Italic">Ms</Emphasis>-<Emphasis Type="Italic">cd1</Emphasis> gene in cabbage (<Emphasis Type="Italic">Brassica oleracea</Emphasis>)

Key message

Using map-based cloning, we delimited the Ms - cd1 gene responsible for the male sterile phenotype in B. oleracea to an approximately 39-kb fragment. Expression analysis suggests that a new predicted gene, a homolog of the Arabidopsis SIED1 gene, is a potential candidate gene.

Abstract

A dominant genic male sterile (DGMS) mutant 79-399-3 in Brassica oleracea (B. oleracea) is controlled by a single gene named Ms-cd1, which was genetically mapped on chromosome C09. The derived DGMS lines of 79-399-3 have been successfully applied in hybrid cabbage breeding and commercial hybrid seed production of several B. oleracea cultivars in China. However, the Ms-cd1 gene responsible for the DGMS has not been identified, and the molecular basis of the DGMS is unclear, which then limits its widespread application in hybrid cabbage seed production. In the present study, a large BC₉ population with 12,269 individuals was developed for map-based cloning of the Ms-cd1 gene, and Ms-cd1 was mapped to a 39.4-kb DNA fragment between two InDel markers, InDel14 and InDel24. Four genes were identified in this region, including two annotated genes based on the available B. oleracea annotation database and two new predicted open reading frames (ORFs). Finally, a newly predicted ORF designated Bol357N3 was identified as the candidate of the Ms-cd1 gene. These results will be useful to reveal the molecular mechanism of the DGMS and develop more practical DGMS lines with stable male sterility for hybrid seed production in cabbage.