Whole-genome Shotgun Contigs Research Articles

First report of Erwinia pyrifoliae causing flower blight and fruit rot on greenhouse-grown Fragaria × ananassa in the United States.

Erwinia pyrifoliae causes disease of pear (Pyrus spp.), apple (Malus spp.), and strawberry (Fragaria × ananassa) (Wenneker and Bergsma-Vlami 2015), which are economically important commodities in the US. Disease symptoms on pear and apple are indistinguishable from those caused by the non-quarantine fire blight pathogen, E. amylovora (Kim et al. 1999), which also causes disease on strawberries (Atanasova et al. 2005). Samples of greenhouse-grown strawberry 'Albion' from Ohio were submitted to the Purdue Plant and Pest Diagnostic Lab in December 2023. Fruits were stunted with brown lesions, while sepals and pedicels had brown-black water-soaked lesions. Cut fruit exuded bacterial ooze from the main vascular bundle. Bacterial streaming was observed microscopically from symptomatic tissue which tested positive with the E. amylovora ImmunoStrip® (Agdia Inc., Elkhart, IN); reported by the manufacturer to cross-react with E. pyrifoliae. Isolation from symptomatic tissue produced pure cultures on sucrose peptone agar and Kings medium B after incubation at 27°C for 48 hr, and colonies appeared circular and white/opaque. Crude DNA extractions were prepared by boiling colony suspensions in Tris-EDTA buffer. Two independent real-time PCR tests specific for E. pyrifoliae (Lehman et al. 2008; Yasuhara-Bell et al. 2024) produced positive results. Conventional PCR using an E. pyrifoliae-specific primer set targeting a divergent region between pstS and glmS genes (Wensing et al. 2011) also produced positive results. The amplicon was Sanger-sequenced and deposited into NCBI GenBank (Accession PP757383). BLASTn analysis using the Nucleotide collection and Whole-genome shotgun contigs revealed top matches (100% query coverage; 97.5% identity to type strain DSM 12163) with E. pyrifoliae only; next closest match was E. amylovora (53% query coverage). To confirm Koch's postulates, immature fruit of six healthy strawberry 'Albion' plants were wounded with a sterile pipette tip and then submersed in a bacterial suspension in sterile deionized water (DI H2O) (3.1×107 cells/ml). Fruit of six additional plants were mock inoculated using sterile DI H2O. Plants were placed in plastic bags for 48 hr at room temperature with a 12-hr photoperiod. Symptoms were first observed on inoculated plants 1.5 days post-inoculation (DPI). Brown discoloration was observed within fruit and as spreading lesions on fruit pedicels by 4 DPI; mock-inoculated plants remained asymptomatic. Bacterial streaming from symptomatic tissue allowed successful re-isolation of the bacterium. Molecular testing confirmed isolates to be E. pyrifoliae, thus completing Koch's postulates. Following initial confirmation, additional samples of infected strawberry ('Albion' and 'GB96') from the same greenhouse were confirmed positive for E. pyrifoliae by molecular testing and sequencing. To our knowledge this is the first time Erwinia pyrifoliae was detected in the US. There are many known pathways of introduction from Asia and Europe; however, pstS-glmS sequence comparison with strawberry isolates from the Netherlands (sequences provided by M.J.C. Pel) suggests this US strawberry strain is unique, but most closely related to Japanese strains (98.5% identity). Potential origin of this strain is unknown, but comparative genomics studies to investigate relatedness among strains are planned.

Presence of p25alpha-Domain in Seed Plants (Spermatophyta): Microbial/Animal Contaminations and/or Orthologs.

Genome and transcriptome assembly data often contain DNA and RNA contaminations from external organisms, introduced during nucleotide extraction or sequencing. In this study, contamination of seed plant (Spermatophyta) transcriptomes/genomes with p25alpha domain encoding RNA/DNA was systematically investigated. This domain only occurs in organisms possessing a eukaryotic flagellum (cilium), which seed plants usually do not have. Nucleotide sequences available at the National Center for Biotechnology Information website, including transcriptome shotgun assemblies (TSAs), whole-genome shotgun contigs (WGSs), and expressed sequence tags (ESTs), were searched for sequences containing a p25alpha domain in Spermatophyta. Despite the lack of proteins containing the p25alpha domain, such fragments or complete mRNAs in some EST and TSA databases were found. A phylogenetic analysis showed that these were contaminations whose possible sources were microorganisms (flagellated fungi, protists) and arthropods/worms; however, there were cases where it cannot be excluded that the sequences found were genuine hits and not of external origin.

Identification and characterization of epicuticular proteins of nematodes sharing motifs with cuticular proteins of arthropods.

Specific collagens and insoluble proteins called cuticlins are major constituents of the nematode cuticles. The epicuticle, which forms the outermost electron-dense layer of the cuticle, is composed of another category of insoluble proteins called epicuticlins. It is distinct from the insoluble cuticlins localized in the cortical layer and the fibrous ribbon underneath lateral alae. Our objective was to identify and characterize genes and their encoded proteins forming the epicuticle. The combination between previously obtained laboratory results and recently made available data through the whole-genome shotgun contigs (WGS) and the transcriptome Shotgun Assembly (TSA) sequencing projects of Ascaris suum allowed us to identify the first epicuticlin gene, Asu-epic-1, on the chromosome VI. This gene is formed of exon1 (55 bp) and exon2 (1067 bp), separated by an intron of 1593 bp. Exon 2 is formed of tandem repeats (TR) whose number varies in different cDNA and genomic clones of Asu-epic-1. These variations could be due to slippage of the polymerases during DNA replication and RNA transcription leading to insertions and deletions (Indels). The deduced protein, Asu-EPIC-1, consists of a signal peptide of 20 amino acids followed by 353 amino acids composed of seven TR of 49 or 51 amino acids each. Three highly conserved tyrosine motifs characterize each repeat. The GYR motif is the Pfam motif PF02756 present in several cuticular proteins of arthropods. Asu-EPIC-1 is an intrinsically disordered protein (IDP) containing seven predicted molecular recognition features (MoRFs). This type of protein undergoes a disorder-to-order transition upon binding protein partners. Three epicuticular sequences have been identified in A. suum, Ascaris lumbricoides, and Toxocara canis. Homologous epicuticular proteins were identified in over 50 other nematode species. The potential of this new category of proteins in forming the nematode cuticle through covalent interactions with other cuticular components, particularly with collagens, is discussed. Their localization in the outermost layer of the nematode body and their unique structure render them crucial candidates for biochemical and molecular interaction studies and targets for new biotechnological and biomedical applications.

Structural analysis of the α subunit of Na(+)/K(+) ATPase genes in invertebrates.

The Na(+)/K(+) ATPase is a ubiquitous pump coordinating the transport of Na(+) and K(+) across the membrane of cells and its role is fundamental to cellular functions. It is heteromer in eukaryotes including two or three subunits (α, β and γ which is specific to the vertebrates). The catalytic functions of the enzyme have been attributed to the α subunit. Several complete α protein sequences are available, but only few gene structures were characterized. We identified the genomic sequences coding the α-subunit of the Na(+)/K(+) ATPase, from the whole-genome shotgun contigs (WGS), NCBI Genomes (chromosome), Genomic Survey Sequences (GSS) and High Throughput Genomic Sequences (HTGS) databases across distinct phyla. One copy of the α subunit gene was found in Annelida, Arthropoda, Cnidaria, Echinodermata, Hemichordata, Mollusca, Placozoa, Porifera, Platyhelminthes, Urochordata, but the nematodes seem to possess 2 to 4 copies. The number of introns varied from 0 (Platyhelminthes) to 26 (Porifera); and their localization and length are also highly variable. Molecular phylogenies (Maximum Likelihood and Maximum Parsimony methods) showed some clusters constituted by (Chordata/(Echinodermata/Hemichordata)) or (Plathelminthes/(Annelida/Mollusca)) and a basal position for Porifera. These structural analyses increase our knowledge about the evolutionary events of the α subunit genes in the invertebrates.

Phylogenetic analysis and molecular evolution patterns in the MIR482-MIR1448 polycistron of Populus L.

The microRNAs (miRNAs) miR482 and miR1448 are disease resistance-related miRNAs; the former is ubiquitously distributed in seed plants whereas the latter has only been reported in Populus trichocarpa. The precursor and mature sequences of poplar miR1448 are highly homologous to those of poplar miR482, and these two miRNAs are located in one transcript as a polycistron. Therefore, we hypothesized that the MIR1448 gene may have evolved from the MIR482 gene in poplar. However, the molecular evolution patterns of this process remain unclear. In this study, utilizing cloning and Blast analysis in NCBI ESTs and whole-genome shotgun contigs (WGS) dataset, we determined that the MIR482-MIR1448 polycistron is a family-specific clustered miRNA in Salicaceae. Moreover, phylogenetic analysis illustrated that MIR1448 is the product of a tandem duplication event from MIR482. Nucleotide substitution analysis revealed that both MIR482 and MIR1448 have more rapid evolution ratios than ribosomal DNA (rDNA) genes, and that compensatory mutations that occurred in the stem region of the secondary structure were the main mechanisms that drove the evolution of these MIRNA genes. Furthermore, by comparing the substitution patterns in the miRNA-target complexes of miR482 and miR1448, we inferred that co-evolution between miRNAs and their targets was the major force that drove the “duplicated MIR482” evolve to MIR1448. We propose a novel miRNA-target pairing pattern called the “frameshift targeted mechanism” to explain the gain of target genes by miR1448. The results also imply that the major role of miR482 was in resistance to disease or other stresses via NBS-LRR proteins, whereas the biological functions of miR1448 are more diverse.

Discovery, validation and characterization of 1039 cattle single nucleotide polymorphisms

We identified approximately 13 000 putative single nucleotide polymorphisms (SNPs) by comparison of repeat-masked BAC-end sequences from the cattle RPCI-42 BAC library with whole-genome shotgun contigs of cattle genome assembly Btau 1.0. Genotyping of a subset of these SNPs was performed on a panel containing 186 DNA samples from 18 cattle breeds including 43 trios. Of 1039 SNPs confirmed as polymorphic in the panel, 998 had minor allele frequency > or =0.25 among unrelated individuals of at least one breed. When Btau 4.0 became available, 974 of these validated SNPs were assigned in silico to known cattle chromosomes, while 41 SNPs were mapped to unassigned sequence scaffolds, yielding one SNP every approximately 3 Mbp on average. Twenty-four SNPs identified in Btau 1.0 were not mapped to Btau 4.0. Of the 1015 SNPs mapped to Btau 4.0, 959 SNPs had nucleotide bases identical in Btau 4.0 and Btau 1.0 contigs, whereas 56 bases were changed, resulting in the loss of the in silico SNP in Btau 4.0. Because these 1039 SNPs were all directly confirmed by genotyping on the multi-breed panel, it is likely that the original polymorphisms were correctly identified. The 1039 validated SNPs identified in this study represent a new and useful resource for genome-wide association studies and applications in animal breeding.

Annotation and BAC/PAC localization of nonredundant ESTs from drought-stressed seedlings of an indica rice.

To decipher the genes associated with drought stress response and to identify novel genes in rice, we utilized 1540 high-quality expressed sequence tags (ESTs) for functional annotation and mapping to rice genomic sequences. These ESTs were generated earlier by 3'-end single-pass sequencing of 2000 cDNA clones from normalized cDNA libraries constructed form drought-stressed seedlings of an indica rice. A rice UniGene set of 1025 transcripts was constructed from this collection through the BLASTN algorithm. Putative functions of 559 nonredundant ESTs were identified by BLAST similarity search against public databases. Putative functions were assigned at a stringency E value of 10(-6) in BLASTN and BLASTX algorithms. To understand the gene structure and function further, we have utilized the publicly available finished and unfinished rice BAC/PAC (BAC, bacterial artificial chromosome; PAC, P1 artificial chromosome) sequences for similarity search using the BLASTN algorithm. Further, 603 nonredundant ESTs have been mapped to BAC/PAC clones. BAC clones were assigned by a homology of above 95% identity along 90% of EST sequence length in the aligned region. In all, 700 ESTs showed rice EST hits in GenBank. Of the 325 novel ESTs, 128 were localized to BAC clones. In addition, 127 ESTs with identified putative functions but with no homology in IRGSP (International Rice Genome Sequencing Program) BAC/PAC sequences were mapped to the Chinese WGS (whole genome shotgun contigs) draft sequence of the rice genome. Functional annotation uncovered about a hundred candidate ESTs associated with abiotic stress in rice and Arabidopsis that were previously reported based on microarray analysis and other studies. This study is a major effort in identifying genes associated with drought stress response and will serve as a resource to rice geneticists and molecular biologists.