Dina A. St. Clair

Global eQTL Mapping Reveals the Complex Genetic Architecture of Transcript-Level Variation in Arabidopsis

The genetic architecture of transcript-level variation is largely unknown. The genetic determinants of transcript-level variation were characterized in a recombinant inbred line (RIL) population (n = 211) of Arabidopsis thaliana using whole-genome microarray analysis and expression quantitative trait loci (eQTL) mapping of transcript levels as expression traits (e-traits). Genetic control of transcription was highly complex: one-third of the quantitatively controlled transcripts/e-traits were regulated by cis-eQTL, and many trans-eQTL mapped to hotspots that regulated hundreds to thousands of e-traits. Several thousand eQTL of large phenotypic effect were detected, but almost all (93%) of the 36,871 eQTL were associated with small phenotypic effects (R2 < 0.3). Many transcripts/e-traits were controlled by multiple eQTL with opposite allelic effects and exhibited higher heritability in the RILs than their parents, suggesting nonadditive genetic variation. To our knowledge, this is the first large-scale global eQTL study in a relatively large plant mapping population. It reveals that the genetic control of transcript level is highly variable and multifaceted and that this complexity may be a general characteristic of eukaryotes.

Quantitative Disease Resistance and Quantitative Resistance Loci in Breeding

Quantitative disease resistance (QDR) has been observed within many crop plants but is not as well understood as qualitative (monogenic) disease resistance and has not been used as extensively in breeding. Mapping quantitative trait loci (QTLs) is a powerful tool for genetic dissection of QDR. DNA markers tightly linked to quantitative resistance loci (QRLs) controlling QDR can be used for marker-assisted selection (MAS) to incorporate these valuable traits. QDR confers a reduction, rather than lack, of disease and has diverse biological and molecular bases as revealed by cloning of QRLs and identification of the candidate gene(s) underlying QRLs. Increasing our biological knowledge of QDR and QRLs will enhance understanding of how QDR differs from qualitative resistance and provide the necessary information to better deploy these resources in breeding. Application of MAS for QRLs in breeding for QDR to diverse pathogens is illustrated by examples from wheat, barley, common bean, tomato, and pepper. Strategies for optimum deployment of QRLs require research to understand effects of QDR on pathogen populations over time.

Global expression analysis of nucleotide binding site-leucine rich repeat-encoding and related genes in Arabidopsis

BackgroundNucleotide binding site-leucine rich repeat (NBS-LRR)-encoding genes comprise the largest class of plant disease resistance genes. The 149 NBS-LRR-encoding genes and the 58 related genes that do not encode LRRs represent approximately 0.8% of all ORFs so far annotated in Arabidopsis ecotype Col-0. Despite their prevalence in the genome and functional importance, there was little information regarding expression of these genes.ResultsWe analyzed the expression patterns of ~170 NBS-LRR-encoding and related genes in Arabidopsis Col-0 using multiple analytical approaches: expressed sequenced tag (EST) representation, massively parallel signature sequencing (MPSS), microarray analysis, rapid amplification of cDNA ends (RACE) PCR, and gene trap lines. Most of these genes were expressed at low levels with a variety of tissue specificities. Expression was detected by at least one approach for all but 10 of these genes. The expression of some but not the majority of NBS-LRR-encoding and related genes was affected by salicylic acid (SA) treatment; the response to SA varied among different accessions. An analysis of previously published microarray data indicated that ten NBS-LRR-encoding and related genes exhibited increased expression in wild-type Landsberg erecta (Ler) after flagellin treatment. Several of these ten genes also showed altered expression after SA treatment, consistent with the regulation of R gene expression during defense responses and overlap between the basal defense response and salicylic acid signaling pathways. Enhancer trap analysis indicated that neither jasmonic acid nor benzothiadiazole (BTH), a salicylic acid analog, induced detectable expression of the five NBS-LRR-encoding genes and one TIR-NBS-encoding gene tested; however, BTH did induce detectable expression of the other TIR-NBS-encoding gene analyzed. Evidence for alternative mRNA polyadenylation sites was observed for many of the tested genes. Evidence for alternative splicing was found for at least 12 genes, 11 of which encode TIR-NBS-LRR proteins. There was no obvious correlation between expression pattern, phylogenetic relationship or genomic location of the NBS-LRR-encoding and related genes studied.ConclusionTranscripts of many NBS-LRR-encoding and related genes were defined. Most were present at low levels and exhibited tissue-specific expression patterns. Expression data are consistent with most Arabidopsis NBS-LRR-encoding and related genes functioning in plant defense responses but do not preclude other biological roles.

Genomic Survey of Gene Expression Diversity in Arabidopsis thaliana

Differential gene expression controls variation in numerous plant traits, such as flowering time and plant/pest interactions, but little is known about the genomic distribution of the determinants of transcript levels and their associated variation. Affymetrix ATH1 GeneChip microarrays representing 22,810 genes were used to survey the transcriptome of seven Arabidopsis thaliana accessions in the presence and absence of exogenously applied salicylic acid (SA). These accessions encompassed ∼80% of the moderate- to high-frequency nucleotide polymorphisms in Arabidopsis. A factorial design, consisting of three biological replicates per accession for the two treatments at three time points (4, 28, and 52 hr post-treatment), and a total of 126 microarrays were used. Between any pair of Arabidopsis accessions, we detected on average 2234 genes (ranging from 1428 to 3334) that were significantly differentially expressed under the conditions of this experiment, using a split-plot analysis of variance. Upward of 6433 genes were differentially expressed between at least one pair of accessions. These results suggest that analysis of additional genetic, developmental, and environmental conditions may show that a significant fraction of the Arabidopsis genome is differentially expressed. Examination of sequence diversity demonstrated a significant positive association with diversity in gene expression.

Identification of QTLs controlling gene expression networks defined a priori

BackgroundGene expression microarrays allow the quantification of transcript accumulation for many or all genes in a genome. This technology has been utilized for a range of investigations, from assessments of gene regulation in response to genetic or environmental fluctuation to global expression QTL (eQTL) analyses of natural variation. Current analysis techniques facilitate the statistical querying of individual genes to evaluate the significance of a change in response, also known as differential expression. Since genes are also known to respond as groups due to their membership in networks, effective approaches are needed to investigate transcriptome variation as related to gene network responses.ResultsWe describe a statistical approach that is capable of assessing higher-order a priori defined gene network response, as measured by microarrays. This analysis detected significant network variation between two Arabidopsis thaliana accessions, Bay-0 and Shahdara. By extending this approach, we were able to identify eQTLs controlling network responses for 18 out of 20 a priori-defined gene networks in a recombinant inbred line population derived from accessions Bay-0 and Shahdara.ConclusionThis approach has the potential to be expanded to facilitate direct tests of the relationship between phenotypic trait and transcript genetic architecture. The use of a priori definitions for network eQTL identification has enormous potential for providing direction toward future eQTL analyses.

Natural Variation among Arabidopsis thaliana Accessions for Transcriptome Response to Exogenous Salicylic Acid

Little is known about how gene expression variation within a given species controls phenotypic variation under different treatments or environments. Here, we surveyed the transcriptome response of seven diverse Arabidopsis thaliana accessions in response to two treatments: the presence and absence of exogenously applied salicylic acid (SA), an important signaling molecule in plant defense. A factorial experiment was conducted with three biological replicates per accession with and without applications of SA and sampled at three time points posttreatment. Transcript level data from Affymetrix ATH1 microarrays were analyzed on both per-gene and gene-network levels to detect expression level polymorphisms associated with SA response. Significant variation in transcript levels for response to SA was detected among the accessions, with relatively few genes responding similarly across all accessions and time points. Twenty-five of 54 defined gene networks identified from other microarray studies (pathogen-challenged Columbia [Col-0]) showed a significant response to SA in one or more accessions. A comparison of gene-network relationships in our data to the pathogen-challenged Col-0 data demonstrated a higher-order conservation of linkages between defense response gene networks. Cvi-1 and Mt-0 appeared to have globally different SA responsiveness in comparison to the other five accessions. Expression level polymorphisms for SA response were abundant at both individual gene and gene-network levels in the seven accessions, suggesting that natural variation for SA response is prevalent in Arabidopsis.

Genetic diversity of Alternaria alternata isolated from tomato in California assessed using RAPDs

Black mould lesions were caused by Alternaria alternata in 76% of 228 tomato fruit with characteristic sunken black lesions collected from fields of processing tomatoes in California. Analysis of 29 RAPD primers revealed a high level of genetic diversity among the 69 isolates tested. Two major phenetic groups (Group 1 with 55 isolates and Group 2 with 14) were identified independently by PCA and by UPGMA of Jaccard similarity coefficients. Only 34 of 137 RAPD markers were monomorphic for all isolates and the genetic similarity between the two groups was 50%. Co-infection of black mould lesions by genetically distinct strains of A. alternata occurred in two of 10 isolates tested. There was no evidence for geographic clustering of isolates with high levels of genetic similarity, suggesting that isolates are widely dispersed across California. Only one isolate was identified which also caused stem canker disease on a susceptible tomato cv., suggesting that these strains play a minor role in causing black mould on processing tomatoes in California. This isolate and two other known stem canker isolates were clustered together with 11 other isolates in Group 2. Group 2-specific bands were also identified in a survey of seven additional isolates known to produce host-specific toxins.

Outcrossing in the homothallic oomycete, Pythium ultimum, detected with molecular markers

The oomycete Pythium ultimum is homothallic, thus a single isolate completes the sexual stage in pure culture. It has been generally assumed that homothallic oomycetes are predominantly inbreeding. In P. ultimum, antheridia occasionally develop from hyphae not directly connected to the oogonium and appear to participate in fertilization, suggesting a possible mechanism for outcrossing. We have used molecular markers to confirm that outcrossing can occur between isolates of P. ultimum. Genetic markers based on randomly amplified polymorphic DNA (RAPD) and restriction fragment length polymorphisms (RFLP) were used to distinguish isolates in a collection of P. ultimum. Two isolates displaying a high level of polymorphism were mixed and placed on media which allows the development of the sexual stage. RAPD markers were used to screen single oospore progeny to identify potential hybrids between the two parental isolates. Subsequent self-fertilization of one putative F1 yielded a F2 population which demonstrated segregation and independent assortment of RAPD and RFLP markers. A similar strategy was used to show that an isolate which is incapable of producing oospores in pure culture can outcross when mixed with a homothallic isolate. These results suggest that other homothallic oomycetes may be capable of outcrossing, and sexual reproduction may, therefore, play an important role in the generation of variation in homothallic oomycetes.

Resistance to bacterial canker in tomato (Lycopersicon hirsutum LA407) and its progeny derived from crosses to L. esculentum

Bacterial canker caused by Clavibacter michiganensis subsp. michiganensis causes significant yield losses on tomatoes grown in a humid environment. This study was conducted to identify a source of resistance that could be easily crossed to cultivated tomato and to study the inheritance of resistance. Diverse bacterial strains representative of the major DNA fingerprint classes endemic to North America were used to screen germ plasm and populations derived from wide crosses. Partial resistance to genetically characterized and distinct strains of C. michiganensis subsp. michiganensis was identified in a wild relative of cultivated tomato, Lycopersicon hirsutum Lycopersicon accession (LA)407. The level of resistance in LA407 was not significantly different from that of the resistant L. peruvianum control, LA2157. Resistance from LA407 was recovered in lines from a BC2S4 inbred backcross (IBC) population in both greenhouse and field trials. Linear correlations between field and greenhouse resistance scores were significant, though correlation coefficients tended to be low. Variance components for genetic and environmental variation in resistance were used to estimate broad-sense heritability in the IBC population. These estimates were moderate to high, ranging from 0.34 to 0.85. The number of genes contributing to resistance was estimated from four trials, with most estimates falling in the range of one to three loci. Two lines from the IBC population, IBL 2353 and IBL 2361, were identified as sources that retain resistance in a genetic background that has a theoretical L. esculentum genome content of 87.5%.

A comparison of microarray and MPSS technology platforms for expression analysis of Arabidopsis.

BackgroundSeveral high-throughput technologies can measure in parallel the abundance of many mRNA transcripts within a sample. These include the widely-used microarray as well as the more recently developed methods based on sequence tag abundances such as the Massively Parallel Signature Sequencing (MPSS) technology. A comparison of microarray and MPSS technologies can help to establish the metrics for data comparisons across these technology platforms and determine some of the factors affecting the measurement of mRNA abundances using different platforms.ResultsWe compared transcript abundance (gene expression) measurement data obtained using Affymetrix and Agilent microarrays with MPSS data. All three technologies were used to analyze the same set of mRNA samples; these samples were extracted from various wild type Arabidopsis thaliana tissues and floral mutants. We calculated correlations and used clustering methodology to compare the normalized expression data and expression ratios across samples and technologies. Abundance expression measurements were more similar between different samples measured by the same technology than between the same sample measured by different technologies. However, when expression ratios were employed, samples measured by different technologies were found to cluster together more frequently than with abundance expression levels.Furthermore, the two microarray technologies were more consistent with each other than with MPSS. We also investigated probe-position effects on Affymetrix data and tag-position effects in MPSS. We found a similar impact on Affymetrix and MPSS measurements, which suggests that these effects were more likely a characteristic of the RNA sample rather than technology-specific biases.ConclusionComparisons of transcript expression ratios showed greater consistency across platforms than measurements of transcript abundance. In addition, for measurements based on abundances, technology differences can mask the impact of biological differences between samples and tissues.