Exploiting a wheat EST database to assess genetic diversity

Expressed sequence tag (EST) markers have been used to assess variety and genetic diversity in wheat (Triticum aestivum). In this study, 1549 ESTs from wheat infested with yellow rust were used to examine the genetic diversity of six susceptible and resistant wheat cultivars. The aim of using these cultivars was to improve the competitiveness of public wheat breeding programs through the intensive use of modern, particularly marker-assisted, selection technologies. The F₂ individuals derived from cultivar crosses were screened for resistance to yellow rust at the seedling stage in greenhouses and adult stage in the field to identify DNA markers genetically linked to resistance. Five hundred and sixty ESTs were assembled into 136 contigs and 989 singletons. BlastX search results showed that 39 (29%) contigs and 96 (10%) singletons were homologous to wheat genes. The database-matched contigs and singletons were assigned to eight functional groups related to protein synthesis, photosynthesis, metabolism and energy, stress proteins, transporter proteins, protein breakdown and recycling, cell growth and division and reactive oxygen scavengers. PCR analyses with primers based on the contigs and singletons showed that the most polymorphic functional categories were photosynthesis (contigs) and metabolism and energy (singletons). EST analysis revealed considerable genetic variability among the Turkish wheat cultivars resistant and susceptible to yellow rust disease and allowed calculation of the mean genetic distance between cultivars, with the greatest similarity (0.725) being between Harmankaya99 and Sönmez2001, and the lowest (0.622) between Aytin98 and Izgi01.     

Permanent draft genome sequence of Dethiosulfovibrio peptidovorans type strain (SEBR 4207)

Dethiosulfovibrio peptidovorans Magot et al. 1997 is the type species of the genus Dethiosulfovibrio of the family Synergistaceae in the recently created phylum Synergistetes. The strictly anaerobic, vibriod, thiosulfate-reducing bacterium utilizes peptides and amino acids, but neither sugars nor fatty acids. It was isolated from an offshore oil well where it was been reported to be involved in pitting corrosion of mild steel. Initially, this bacterium was described as a distant relative of the genus Thermoanaerobacter, but was not assigned to a genus, it was subsequently placed into the novel phylum Synergistetes. A large number of repeats in the genome sequence prevented an economically justifiable closure of the last gaps. This is only the third published genome from a member of the phylum Synergistetes. The 2,576,359 bp long genome consists of three contigs with 2,458 protein-coding and 59 RNA genes and is part of the Genomic Encyclopedia of Bacteria and Archaea project.     

Molecular variation and fingerprinting of Leucadendron cultivars (Proteaceae) by ISSR markers

BACKGROUND AND AIMS: There are more than 80 species of Leucadendron and most are used as cut flowers. Currently, more than 100 cultivars are used by industry and many of them are interspecific hybrids. The origin of most cultivars is unclear and their genetic diversity and relationships have not been studied. This investigation was carried out to evaluate the genetic variation and relationships among 30 Leucadendron cultivars. METHODS: ISSR markers were applied to determine the genetic variation and to discriminate Leucadendron cultivars. Sixty-four ISSR primers were screened and 25 primers were selected for their ability to produce clear and reproducible patterns of multiple bands. KEY RESULTS: A total of 584 bands of 305-2400 bp were amplified, of which 97 % were polymorphic. A dendrogram generated using the Unweighted Pair Group Method with Arithmetic Average based on a distance measure of total character difference showed that the Leucadendron cultivars clustered into two main groups. Twenty-four of the 30 cultivars can be unequivocally differentiated, but identical profiles were observed for three cultivar pairs, 'Katie's Blush' and 'Silvan Red', 'Highlights' and 'Maui Sunset', and 'Yellow Crest' and 'Yellow Devil'. CONCLUSIONS: ISSR profiling is a powerful method for the identification and molecular classification of Leucadendron cultivars. A fingerprinting key was generated based on the banding patterns produced using two ISSR primers (UBC856 and UBC857). In addition cultivar-specific ISSR bands were obtained for 17 of the 30 Leucadendron cultivars tested.     

Correlation between grass carp (Ctenopharyngodon idella) resistance to grass carp reovirus and the genetic insert-deletion polymorphisms in promoter and intron of RIG-I gene

RIG-I (retinoic acid-inducible gene I) is an essential cytosolic pathogen recognition receptor that binds to a variety of viral RNA or DNA to induce type I interferons. In the present study, insert–deletion polymorphisms in promoter and introns of CiRIG-I (Ctenopharyngodon idella RIG-I) were explored, their associations with resistance/susceptibility to grass carp reovirus (GCRV) were analyzed. To this end, genomic sequence of CiRIG-I gene was obtained, and twenty pairs of primers were prepared for the detection of insert–deletion polymorphisms. Five insert–deletion mutations were found, a 2-bp mutation and an 8-bp mutation existed in the promoter and other three sizes in 74 bp, 146 bp and 53 bp were sited in the intron 8. After a challenge experiment, only the genotype and allele of − 740 insert–deletion mutation in the promoter and allele of 6804 insert–deletion mutation were significantly associated with resistance/susceptibility to GCRV among the five mutations (P < 0.05). To further identify this correlation, another independent challenge test was carried out. The result revealed that the cumulative mortality in ins/ins genotype individuals (43.75%) at − 740 insert–deletion mutation was significantly lower than that in ins/del (72.09%) and del/del (74.19%) genotypes (P < 0.05). Linkage disequilibrium and haplotype analysis showed 6610 insert–deletion mutation and 6804 insert–deletion mutation were linkage disequilibrium. The haplotype ins–ins (6610ins–6804ins) was significantly susceptible to GCRV, and ins–del (6610ins–6804del) was significantly resistant to GCRV (P < 0.05). Those could be potential gene markers for the future molecular selection of strains that are resistant to GCRV.     

Draft Genome Sequencing and Comparative Analysis of Aspergillus sojae NBRC4239

We conducted genome sequencing of the filamentous fungus Aspergillus sojae NBRC4239 isolated from the koji used to prepare Japanese soy sauce. We used the 454 pyrosequencing technology and investigated the genome with respect to enzymes and secondary metabolites in comparison with other Aspergilli sequenced. Assembly of 454 reads generated a non-redundant sequence of 39.5-Mb possessing 13 033 putative genes and 65 scaffolds composed of 557 contigs. Of the 2847 open reading frames with Pfam domain scores of >150 found in A. sojae NBRC4239, 81.7% had a high degree of similarity with the genes of A. oryzae. Comparative analysis identified serine carboxypeptidase and aspartic protease genes unique to A. sojae NBRC4239. While A. oryzae possessed three copies of α-amyalse gene, A. sojae NBRC4239 possessed only a single copy. Comparison of 56 gene clusters for secondary metabolites between A. sojae NBRC4239 and A. oryzae revealed that 24 clusters were conserved, whereas 32 clusters differed between them that included a deletion of 18 508 bp containing mfs1, mao1, dmaT, and pks-nrps for the cyclopiazonic acid (CPA) biosynthesis, explaining the no productivity of CPA in A. sojae. The A. sojae NBRC4239 genome data will be useful to characterize functional features of the koji moulds used in Japanese industries.     

High quality SNPs/Indels mining and characterization in ginger from ESTs data base

Ginger (Zingiber officinale Rosc.) is an important herb of the family Zingiberaceae. It is accepted as a universal cure for a multitude of diseases in Indian systems of medicine and its rhizomes are equally popular as a spice ingredient throughout Asia. SNPs, the definitive genetic markers, representing the finest resolution of a DNA sequence, are abundantly found in populations having a lower rate of mutation and are used for genomic analysis. The public ESTs sequences mostly lack quality files, making high quality SNPs detection more difficult since it is exclusively based on sequence comparisons. In the present study, current dbESTs of NCBI was mined and 38115 ginger ESTs sequences were obtained and assembled into contigs using CAP3 program. In this analysis, recent software tool QualitySNP was used to detect 11523 potential SNPs sites, 8810 high quality SNPs and 1008 indels polymorphisms with a frequency of 1.61 SNPs / 10 kbp. Of ESTs libraries generated from three ginger tissues together, rhizomes had a frequency of 0.32 SNPs and 0.03 indels per 10 kbp whereas the leaves had a frequency of 2.51 SNPs and 0.23 indels per 10 kbp and root is showing relative frequency of 0.76/10 kbp SNPs and 0.02/10 kbp indels. The present analysis provides additional information about the tissue wise presence of haplotypes (222), distribution of high quality exonic (2355) and intronic (6455) SNPs and information about singletons (7538) in addition to contigs transitions and transversions ratio (0.57). Among all tissue detected SNPs, transversions number is higher in comparison to the number of transitions. Quality SNPs detected in this work can be used as markers for further ginger genetic experiments.     

Differentiation of Species and Populations of Aphelenchoides and of Ditylenchus angustus Using a Fragment of Ribosomal DNA

The polymerase chain reaction (PCR) was used to amplify a fragment of the ribosomal DNA (rDNA) from species and undescribed populations of Aphelenchoides and Ditylenchus angustus. The PCR primers used were based on conserved sequences in the 18S and 26S ribosomal RNA genes of Caenorhabditis elegans. In C. elegans, these primers amplify a 1,292 base pair (bp) fragment, which consists of the two internal transcribed spacers and the entire 5.8S gene. Amplification products from crude DNA preparations of 12 species and populations of Aphelenchoides and from D. angustus ranged in size from approximately 860-1,100bp. Southern blots probed with a cloned ribosomal repeat from C. elegans confirmed the identity of these amplified bands as ribosomal fragments. In addition to the differing sizes of the amplified rDNA fragments, the relative intensity of hybridization with the C. elegans probe indicated varying degrees of sequence divergence between species and populations. In some cases, amplified rDNA from the fungal host was evident. Storage of A. composticola at - 45 C for 2 years did not affect the ability to obtain appropriate amplified products from crude DNA preparations. Amplified rDNA fragments were cut with six restriction enzymes, and the restriction fragments produced revealed useful diagnostic differences between species and some undescribed populations. These results were consistent with previous studies based on morphology and isoenzymes. Three undescribed populations of Aphelenchoides were found to be different from all the species examined and from each other.     

Analysis of the formation of flower shapes in wild species and cultivars of tree peony using the MADS-box subfamily gene

Tree peony (Paeonia suffricotisa) cultivars have a unique character compared with wild species; the stamen petalody results in increased whorls of petals and generates different flower forms, which are one of the most important traits for cultivar classification. In order to investigate how petaloid stamens are formed, we obtained the coding sequence (666 bp) and genomic DNA sequence of the PsTM6 genes (belongs to B subfamily of MADS-box gene family) from 23 tree peony samples, Five introns and six exons consisted of the genomic DNA sequence. The analysis of cis-acting regulatory elements in the third and fourth intron indicated that they were highly conserved in all samples. Partial putative amino acids were analyzed and the results suggested that functional differentiation of PsTM6 paralogs apparently affected stamen petalody and flower shape formation due to due to amino acid substitution caused by differences in polarity and electronic charge. Sliding window analysis indicated that the different regions of PsTM6 were subjected to different selection forces, especially in the K domain. This is the first attempt to investigate genetic control of the stamen petalody based on the PsTM6 sequence. This will provide a basis for understanding the evolution of PsTM6 and its the function of in determining stamen morphology of tree peony.     

Performance comparison of second- and third-generation sequencers using a bacterial genome with two chromosomes

Background

The availability of diverse second- and third-generation sequencing technologies enables the rapid determination of the sequences of bacterial genomes. However, identifying the sequencing technology most suitable for producing a finished genome with multiple chromosomes remains a challenge. We evaluated the abilities of the following three second-generation sequencers: Roche 454 GS Junior (GS Jr), Life Technologies Ion PGM (Ion PGM), and Illumina MiSeq (MiSeq) and a third-generation sequencer, the Pacific Biosciences RS sequencer (PacBio), by sequencing and assembling the genome of Vibrio parahaemolyticus, which consists of a 5-Mb genome comprising two circular chromosomes.

Results

We sequenced the genome of V. parahaemolyticus with GS Jr, Ion PGM, MiSeq, and PacBio and performed de novo assembly with several genome assemblers. Although GS Jr generated the longest mean read length of 418 bp among the second-generation sequencers, the maximum contig length of the best assembly from GS Jr was 165 kbp, and the number of contigs was 309. Single runs of Ion PGM and MiSeq produced data of considerably greater sequencing coverage, 279× and 1,927×, respectively. The optimized result for Ion PGM contained 61 contigs assembled from reads of 77× coverage, and the longest contig was 895 kbp in size. Those for MiSeq were 34 contigs, 58× coverage, and 733 kbp, respectively. These results suggest that higher coverage depth is unnecessary for a better assembly result. We observed that multiple rRNA coding regions were fragmented in the assemblies from the second-generation sequencers, whereas PacBio generated two exceptionally long contigs of 3,288,561 and 1,875,537 bps, each of which was from a single chromosome, with 73× coverage and mean read length 3,119 bp, allowing us to determine the absolute positions of all rRNA operons.

Conclusions

PacBio outperformed the other sequencers in terms of the length of contigs and reconstructed the greatest portion of the genome, achieving a genome assembly of “finished grade” because of its long reads. It showed the potential to assemble more complex genomes with multiple chromosomes containing more repetitive sequences.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-699) contains supplementary material, which is available to authorized users.     

The complete mitochondrial genome sequence and gene organization of the mud crab (Scylla paramamosain) with phylogenetic consideration

The complete mitochondrial genome is of great importance for better understanding the genome-level characteristics and phylogenetic relationships among related species. In the present study, we determined the complete mitochondrial genome DNA sequence of the mud crab (Scylla paramamosain) by 454 deep sequencing and Sanger sequencing approaches. The complete genome DNA was 15,824 bp in length and contained a typical set of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a putative control region (CR). Of 37 genes, twenty-three were encoded by the heavy strand (H-strand), while the other ones were encoded by light strand (L-strand). The gene order in the mitochondrial genome was largely identical to those obtained in most arthropods, although the relative position of gene tRNA^His differed from other arthropods. Among 13 protein-coding genes, three (ATPase subunit 6 (ATP6), NADH dehydrogenase subunits 1 (ND1) and ND3) started with a rare start codon ATT, whereas, one gene cytochrome c oxidase subunit I (COI) ended with the incomplete stop codon TA. All 22 tRNAs could fold into a typical clover-leaf secondary structure, with the gene sizes ranging from 63 to 73 bp. The phylogenetic analysis based on 12 concatenated protein-coding genes showed that the molecular genetic relationship of 19 species of 11 genera was identical to the traditional taxonomy.