Peter A. Schweitzer

Transcriptome Analysis of Pseudomonas syringae Identifies New Genes, Noncoding RNAs, and Antisense Activity

To fully understand how bacteria respond to their environment, it is essential to assess genome-wide transcriptional activity. New high-throughput sequencing technologies make it possible to query the transcriptome of an organism in an efficient unbiased manner. We applied a strand-specific method to sequence bacterial transcripts using Illuminas high-throughput sequencing technology. The resulting sequences were used to construct genome-wide transcriptional profiles. Novel bioinformatics analyses were developed and used in combination with proteomics data for the qualitative classification of transcriptional activity in defined regions. As expected, most transcriptional activity was consistent with predictions from the genome annotation. Importantly, we identified and confirmed transcriptional activity in areas of the genome inconsistent with the annotation and in unannotated regions. Further analyses revealed potential RpoN-dependent promoter sequences upstream of several noncoding RNAs (ncRNAs), suggesting a role for these ncRNAs in RpoN-dependent phenotypes. We were also able to validate a number of transcriptional start sites, many of which were consistent with predicted promoter motifs. Overall, our approach provides an efficient way to survey global transcriptional activity in bacteria and enables rapid discovery of specific areas in the genome that merit further investigation.

Multi-platform assessment of transcriptome profiling using RNA-seq in the ABRF next-generation sequencing study.

High-throughput RNA sequencing (RNA-seq) greatly expands the potential for genomics discoveries, but the wide variety of platforms, protocols and performance capabilitites has created the need for comprehensive reference data. Here we describe the Association of Biomolecular Resource Facilities next-generation sequencing (ABRF-NGS) study on RNA-seq. We carried out replicate experiments across 15 laboratory sites using reference RNA standards to test four protocols (poly-A–selected, ribo-depleted, size-selected and degraded) on five sequencing platforms (Illumina HiSeq, Life Technologies PGM and Proton, Pacific Biosciences RS and Roche 454). The results show high intraplatform (Spearman rank R > 0.86) and inter-platform (R > 0.83) concordance for expression measures across the deep-count platforms, but highly variable efficiency and cost for splice junction and variant detection between all platforms. For intact RNA, gene expression profiles from rRNA-depletion and poly-A enrichment are similar. In addition, rRNA depletion enables effective analysis of degraded RNA samples. This study provides a broad foundation for cross-platform standardization, evaluation and improvement of RNA-seq.

Genome-Wide Identification of Transcriptional Start Sites in the Plant Pathogen Pseudomonas syringae pv. tomato str. DC3000

RNA-Seq has provided valuable insights into global gene expression in a wide variety of organisms. Using a modified RNA-Seq approach and Illuminas high-throughput sequencing technology, we globally identified 5′-ends of transcripts for the plant pathogen Pseudomonas syringae pv. tomato str. DC3000. A substantial fraction of 5′-ends obtained by this method were consistent with results obtained using global RNA-Seq and 5′RACE. As expected, many 5′-ends were positioned a short distance upstream of annotated genes. We also captured 5′-ends within intergenic regions, providing evidence for the expression of un-annotated genes and non-coding RNAs, and detected numerous examples of antisense transcription, suggesting additional levels of complexity in gene regulation in DC3000. Importantly, targeted searches for sequence patterns in the vicinity of 5′-ends revealed over 1200 putative promoters and other regulatory motifs, establishing a broad foundation for future investigations of regulation at the genomic and single gene levels.

Evaluation of a fibrillin 2 gene haplotype associated with hip dysplasia and incipient osteoarthritis in dogs

OBJECTIVE To determine whether a mutation in the fibrillin 2 gene (FBN2) is associated with canine hip dysplasia (CHD) and osteoarthritis in dogs. ANIMALS 1,551 dogs. Procedures-Hip conformation was measured radiographically. The FBN2 was sequenced from genomic DNA of 21 Labrador Retrievers and 2 Greyhounds, and a haplotype in intron 30 of FBN2 was sequenced in 90 additional Labrador Retrievers and 143 dogs of 6 other breeds. Steady-state values of FBN2 mRNA and control genes were measured in hip joint tissues of fourteen 8-month-old Labrador Retriever-Greyhound crossbreeds. RESULTS The Labrador Retrievers homozygous for a 10-bp deletion haplotype in intron 30 of FBN2 had significantly worse CHD as measured via higher distraction index and extended-hip joint radiograph score and a lower Norberg angle and dorsolateral subluxation score. Among 143 dogs of 6 other breeds, those homozygous for the same deletion haplotype also had significantly worse radiographic CHD. Among the 14 crossbred dogs, as the dorsolateral subluxation score decreased, the capsular FBN2 mRNA increased significantly. Those dogs with incipient hip joint osteoarthritis had significantly increased capsular FBN2 mRNA, compared with those dogs without osteoarthritis. Dogs homozygous for the FBN2 deletion haplotype had significantly less FBN2 mRNA in their femoral head articular cartilage. CONCLUSIONS AND CLINICAL RELEVANCE The FBN2 deletion haplotype was associated with CHD. Capsular gene expression of FBN2 was confounded by incipient secondary osteoarthritis in dysplastic hip joints. Genes influencing complex traits in dogs can be identified by genome-wide screening, fine mapping, and candidate gene screening.

Single nucleotide polymorphisms refine QTL intervals for hip joint laxity in dogs

Hip laxity is one characteristic of canine hip dysplasia (CHD), an inheritable disease that leads to hip osteoarthritis. Using a genome-wide screen with 250 microsatellites in a crossbreed pedigree of 159 dysplastic Labrador retrievers and unaffected greyhounds, we previously identified putative (P < 0.01) QTL on canine chromosomes 11 and 29 (CFA11 and CFA29). To refine these QTL locations, we have genotyped 257 dogs including 105 Labrador retrievers, seven greyhounds, four generations of their crossbreed offspring and three German shepherds for 111 and 171 SNPs on CFA11 and CFA29 respectively. The distraction index (DI, a measure of maximum hip laxity) was used as an intermediate phenotype that predicts whether a hip joint will or will not develop osteoarthritis. Using a multipoint linkage analysis, significant evidence (95% posterior probability) was found for QTL contributing to hip laxity in the 16.2-21 cM region on CFA11 that explained 15-18% of the total variance in DI. Evidence for an independent QTL on CFA29 was weaker than that on CFA11. Identification of the causative mutation(s) will lead to better understanding of biochemical pathways in both dogs and humans with hip laxity and dysplasia.

A next-generation marker genotyping platform (AmpSeq) in heterozygous crops: a case study for marker-assisted selection in grapevine

Marker-assisted selection (MAS) is often employed in crop breeding programs to accelerate and enhance cultivar development, via selection during the juvenile phase and parental selection prior to crossing. Next-generation sequencing and its derivative technologies have been used for genome-wide molecular marker discovery. To bridge the gap between marker development and MAS implementation, this study developed a novel practical strategy with a semi-automated pipeline that incorporates trait-associated single nucleotide polymorphism marker discovery, low-cost genotyping through amplicon sequencing (AmpSeq) and decision making. The results document the development of a MAS package derived from genotyping-by-sequencing using three traits (flower sex, disease resistance and acylated anthocyanins) in grapevine breeding. The vast majority of sequence reads (⩾99%) were from the targeted regions. Across 380 individuals and up to 31 amplicons sequenced in each lane of MiSeq data, most amplicons (83 to 87%) had <10% missing data, and read depth had a median of 220–244×. Several strengths of the AmpSeq platform that make this approach of broad interest in diverse crop species include accuracy, flexibility, speed, high-throughput, low-cost and easily automated analysis.

Evaluation of quantitative trait loci for hip dysplasia in Labrador Retrievers

OBJECTIVE To identify the quantitative trait loci (QTL) that contribute to hip dysplasia in dogs. ANIMALS 192 Labrador Retrievers. PROCEDURES Hip dysplasia was measured by use of the Norberg angle (NA), dorsolateral subluxation (DLS) score, and distraction index (DI). Genome-wide screening was conducted by use of 276 unique microsatellites. Linkage analysis was performed with a variance-based linear model. Logarithm of the odds (LOD) scores were reported when values were > 2.0. RESULTS Canis familiaris autosomes (CFAs) 01, 02, 10, 20, 22, and 32 harbored significant QTL at LOD scores > 2.0. Among the 6 QTL, the QTL on CFA02 had not been reported to harbor QTL for hip dysplasia. The highest LOD score of 3.32 on CFA20 contributed to the second principal component of the DLS score and NA of the right hip joint. The QTL that was mapped on CFA01 (LOD score of 3.13 at 55 centimorgans) was located on the same chromosome reported to harbor a QTL for hip dysplasia in Portuguese Water Dogs and German Shepherd Dogs. In this study, CFAs 10, 20, 22, and 32 harbored QTL for hip dysplasia that have been identified in a Labrador Retriever-Greyhound pedigree and in German Shepherd Dogs. CONCLUSIONS AND CLINICAL RELEVANCE Multiple QTL were clearly involved with hip dysplasia. Identification of these QTL will enable fine-resolution mapping and subsequent assessment of candidate genes within the refined intervals to enable researchers to develop genetic screening tests and preventative and novel therapeutic regimens.

Genomic Plasticity Enables Phenotypic Variation of Pseudomonas syringae pv. tomato DC3000

Whole genome sequencing revealed the presence of a genomic anomaly in the region of 4.7 to 4.9 Mb of the Pseudomonas syringae pv. tomato (Pst) DC3000 genome. The average read depth coverage of Pst DC3000 whole genome sequencing results suggested that a 165 kb segment of the chromosome had doubled in copy number. Further analysis confirmed the 165 kb duplication and that the two copies were arranged as a direct tandem repeat. Examination of the corresponding locus in Pst NCPPB1106, the parent strain of Pst DC3000, suggested that the 165 kb duplication most likely formed after the two strains diverged via transposition of an ISPsy5 insertion sequence (IS) followed by unequal crossing over between ISPsy5 elements at each end of the duplicated region. Deletion of one copy of the 165 kb region demonstrated that the duplication facilitated enhanced growth in some culture conditions, but did not affect pathogenic growth in host tomato plants. These types of chromosomal structures are predicted to be unstable and we have observed resolution of the 165 kb duplication to single copy and its subsequent re-duplication. These data demonstrate the role of IS elements in recombination events that facilitate genomic reorganization in P. syringae.

Antisense sequencing improves the accuracy and precision of A-to-I editing measurements using the peak height ratio method

BackgroundA-to-I RNA editing is found in all phyla of animals and contributes to transcript diversity that may have profound impacts on behavior and physiology. Many transcripts of genes involved in axonal conductance, synaptic transmission and modulation are the targets of A-to-I RNA editing. There are a number of methods to measure the extent of A-to-I RNA editing, but they are generally costly and time consuming. One way to determine the frequency of A-to-I RNA editing is the peak height ratio method, which compares the size of peaks on electropherograms that represent unedited and edited sites.FindingsSequencing of 4 editing sites of the Dα6 nicotinic acetylcholine receptor subunit with an antisense primer (which uses T/C peaks to measure unedited and edited sites, respectively) showed very accurate and precise measurements of A-to-I RNA editing. The accuracy and precision were excellent for all editing sites, including those edited with high or low frequencies. The frequency of A-to-I RNA editing was comparable to the editing frequency as measured by clone counting from the same sample. Sequencing these same sites with the sense primer (which uses A/G peaks) yielded inaccurate and imprecise measurements.ConclusionsWe have validated and improved the accuracy and precision of the peak height ratio method to measure the frequency of A-to-I RNA editing, and shown that results are primer specific. Thus, the correct sequencing primer must be utilized for the most dependable data. When compared to other methods used to measure the frequency of A-to-I RNA editing, the major benefits of the peak height ratio are that this method is inexpensive, fast, non-labor intensive and easily adaptable to many laboratory and field settings.

Illumina DNA test panel-based genotyping of whole genome amplified-DNA extracted from hair samples: performance and agreement with genotyping results from genomic DNA from buccal cells.

Abstract Background: Hair is a DNA source that can be collected easily and inexpensively from participants in epidemiological studies. However, there is concern about DNA quality and quantity. Therefore, we assessed genotyping performance of whole genome amplified (WGA)-DNA extracted from hair using the GenomePlex® method and evaluated its agreement with genotyping results of buccal cell DNA from the same individuals, using the Illumina GoldenGate platform. Methods: The Illumina DNA test panel includes 360 highly validated single nucleotide polymorphisms (SNPs) selected from the Linkage IV Panel that are distributed across the entire genome. DNA was extracted from both archived hair and buccal cell samples obtained from 44 randomly selected subjects participating in a large cohort study in Canada. Results: The genotyping success rate was 97.7% for 44 paired samples. However, WGA-DNA from hair failed more during genotyping in comparison to buccal cell DNA. Hair samples with a pre-WGA-DNA ≥1 ng/μL quantified using the PicoGreen® assay (n=33) showed an average genotyping completion rate of 98.8% and SNP concordance of 91.2% with genotyping performance of buccal cell DNA. In contrast, samples with a pre-WGA-DNA <1 ng/μL had lower genotyping completion rate (94%) and poor SNP concordance (49%). Conclusions: Results suggest that WGA-DNA obtained from hair can produce excellent genotyping call rates and show relatively good SNP concordance with results from buccal cell DNA using high-throughput technology. DNA quantity obtained from hair samples is a crucial determinant of genotyping performance. Larger studies are needed to examine the utility of hair DNA with different genotyping platforms. Clin Chem Lab Med 2009;47:516–22.