Genetics of alternative splicing evolution during sunflower domestication.

Alternative splicing enables organisms to produce the diversity of proteins necessary for multicellular life by using relatively few protein-coding genes. Although differences in splicing have been identified among divergent taxa, the shorter-term evolution of splicing is understudied. The origins of novel splice forms, and the contributions of alternative splicing to major evolutionary transitions, are largely unknown. This study used transcriptomes of wild and domesticated sunflowers to examine splice differentiation and regulation during domestication. We identified substantial splicing divergence between wild and domesticated sunflowers, mainly in the form of intron retention. Transcripts with divergent splicing were enriched for seed-development functions, suggesting that artificial selection impacted splicing patterns. Mapping of quantitative trait loci (QTLs) associated with 144 differential splicing cases revealed primarily trans-acting variation affecting splicing patterns. A large proportion of identified QTLs contain known spliceosome proteins and are associated with splicing variation in multiple genes. Examining a broader set of wild and domesticated sunflower genotypes revealed that most differential splicing patterns in domesticated sunflowers likely arose from standing variation in wild Helianthus annuus and gained frequency during the domestication process. However, several domesticate-associated splicing patterns appear to be introgressed from other Helianthus species. These results suggest that sunflower domestication involved selection on pleiotropic regulatory alleles. More generally, our findings indicate that substantial differences in isoform abundances arose rapidly during a recent evolutionary transition and appear to contribute to adaptation and population divergence.

Genomics of homoploid hybrid speciation: diversity and transcriptional activity of long terminal repeat retrotransposons in hybrid sunflowers.

Hybridization is thought to play an important role in plant evolution by introducing novel genetic combinations and promoting genome restructuring. However, surprisingly little is known about the impact of hybridization on transposable element (TE) proliferation and the genomic response to TE activity. In this paper, we first review the mechanisms by which homoploid hybrid species may arise in nature. We then present hybrid sunflowers as a case study to examine transcriptional activity of long terminal repeat retrotransposons in the annual sunflowers Helianthus annuus, Helianthus petiolaris and their homoploid hybrid derivatives (H. paradoxus, H. anomalus and H. deserticola) using high-throughput transcriptome sequencing technologies (RNAseq). Sampling homoploid hybrid sunflower taxa revealed abundant variation in TE transcript accumulation. In addition, genetic diversity for several candidate genes hypothesized to regulate TE activity was characterized. Specifically, we highlight one candidate chromatin remodelling factor gene with a direct role in repressing TE activity in a hybrid species. This paper shows that TE amplification in hybrid lineages is more idiosyncratic than previously believed and provides a first step towards identifying the mechanisms responsible for regulating and repressing TE expansions.

Genetics of cryptic speciation within an Arctic mustard, Draba nivalis.

Crossing experiments indicate that hybrid sterility barriers frequently have developed within diploid, circumpolar plant species of the genus Draba. To gain insight into the rapid evolution of postzygotic reproductive isolation in this system, we augmented the linkage map of one of these species, D. nivalis, and searched for quantitative trait loci (QTLs) associated with reproductive isolation. The map adds 63 new dominant markers to a previously published dataset of 31 co-dominant microsatellites. These markers include 52 amplified fragment length polymorphisms (AFLPs) and 11 sequence-specific amplified polymorphisms (SSAPs) based on retrotransposon sequence. 22 markers displaying transmission ratio distortion were further included in the map. We resolved eight linkage groups with a total map length of 894 cM. Significant genotype-trait associations, or quantitative trait loci (QTL), were detected for reproductive phenotypes including pollen fertility (4 QTLs), seed set (3 QTLs), flowering time (3 QTLs) and number of flowers (4 QTLs). Observed patterns of inheritance were consistent with the influence of both nuclear-nuclear interactions and chromosomal changes on these traits. All seed set QTLs and one pollen fertility QTL displayed underdominant effects suggestive of the involvement of chromosomal rearrangements in hybrid sterility. Interestingly, D. nivalis is predominantly self-fertilizing, which may facilitate the establishment of underdominant loci and contribute to reproductive isolation.

Genome-scale transcriptional analyses of first-generation interspecific sunflower hybrids reveals broad regulatory compatibility.

BACKGROUND: Interspecific hybridization creates individuals harboring diverged genomes. The interaction of these genomes can generate successful evolutionary novelty or disadvantageous genomic conflict. Annual sunflowers Helianthus annuus and H. petiolaris have a rich history of hybridization in natural populations. Although first-generation hybrids generally have low fertility, hybrid swarms that include later generation and fully fertile backcross plants have been identified, as well as at least three independently-originated stable hybrid taxa. We examine patterns of transcript accumulation in the earliest stages of hybridization of these species via analyses of transcriptome sequences from laboratory-derived F1 offspring of an inbred H. annuus cultivar and a wild H. petiolaris accession. RESULTS: While nearly 14% of the reference transcriptome showed significant accumulation differences between parental accessions, total F1 transcript levels showed little evidence of dominance, as midparent transcript levels were highly predictive of transcript accumulation in F1 plants. Allelic bias in F1 transcript accumulation was detected in 20% of transcripts containing sufficient polymorphism to distinguish parental alleles; however the magnitude of these biases were generally smaller than differences among parental accessions. CONCLUSIONS: While analyses of allelic bias suggest that cis regulatory differences between H. annuus and H. petiolaris are common, their effect on transcript levels may be more subtle than trans-acting regulatory differences. Overall, these analyses found little evidence of regulatory incompatibility or dominance interactions between parental genomes within F1 hybrid individuals, although it is unclear whether this is a legacy or an enabler of introgression between species.

Genomic resources notes accepted 1 February 2013-31 March 2013.

This article documents the public availability of three transcriptome sequences. These genomic resources were developed for the following species: Helianthus anomalus, Helianthus deserticola and Helianthus paradoxus.

Response of sunflower (Helianthus annuus L.) leaf surface defenses to exogenous methyl jasmonate.

Helianthus annuus, the common sunflower, produces a complex array of secondary compounds that are secreted into glandular trichomes, specialized structures found on leaf surfaces and anther appendages of flowers. The primary components of these trichome secretions are sesquiterpene lactones (STL), a diverse class of compounds produced abundantly by the plant family Compositae and believed to contribute to plant defense against herbivory. We treated wild and cultivated H. annuus accessions with exogenous methyl jasmonate, a plant hormone that mediates plant defense against insect herbivores and certain classes of fungal pathogens. The wild sunflower produced a higher density of glandular trichomes on its leaves than the cultivar. Comparison of the profiles of glandular trichome extracts obtained by liquid chromatography-mass spectroscopy (LC-MS) showed that wild and cultivated H. annuus were qualitatively similar in surface chemistry, although differing in the relative size and proportion of various compounds detected. Despite observing consistent transcriptional responses to methyl jasmonate treatment, we detected no significant effect on glandular trichome density or LC-MS profile in cultivated or wild sunflower, with wild sunflower exhibiting a declining trend in overall STL production and foliar glandular trichome density of jasmonate-treated plants. These results suggest that glandular trichomes and associated compounds may act as constitutive defenses or require greater levels of stimulus for induction than the observed transcriptional responses to exogenous jasmonate. Reduced defense investment in domesticated lines is consistent with predicted tradeoffs caused by selection for increased yield; future research will focus on the development of genetic resources to explicitly test the ecological roles of glandular trichomes and associated effects on plant growth and fitness.

Combining genome-wide association mapping and transcriptional networks to identify novel genes controlling glucosinolates in Arabidopsis thaliana.

BACKGROUND: Genome-wide association (GWA) is gaining popularity as a means to study the architecture of complex quantitative traits, partially due to the improvement of high-throughput low-cost genotyping and phenotyping technologies. Glucosinolate (GSL) secondary metabolites within Arabidopsis spp. can serve as a model system to understand the genomic architecture of adaptive quantitative traits. GSL are key anti-herbivory defenses that impart adaptive advantages within field trials. While little is known about how variation in the external or internal environment of an organism may influence the efficiency of GWA, GSL variation is known to be highly dependent upon the external stresses and developmental processes of the plant lending it to be an excellent model for studying conditional GWA. METHODOLOGY/PRINCIPAL FINDINGS: To understand how development and environment can influence GWA, we conducted a study using 96 Arabidopsis thaliana accessions, >40 GSL phenotypes across three conditions (one developmental comparison and one environmental comparison) and ∼230,000 SNPs. Developmental stage had dramatic effects on the outcome of GWA, with each stage identifying different loci associated with GSL traits. Further, while the molecular bases of numerous quantitative trait loci (QTL) controlling GSL traits have been identified, there is currently no estimate of how many additional genes may control natural variation in these traits. We developed a novel co-expression network approach to prioritize the thousands of GWA candidates and successfully validated a large number of these genes as influencing GSL accumulation within A. thaliana using single gene isogenic lines. CONCLUSIONS/SIGNIFICANCE: Together, these results suggest that complex traits imparting environmentally contingent adaptive advantages are likely influenced by up to thousands of loci that are sensitive to fluctuations in the environment or developmental state of the organism. Additionally, while GWA is highly conditional upon genetics, the use of additional genomic information can rapidly identify causal loci en masse.

The complex genetic architecture of the metabolome.

Discovering links between the genotype of an organism and its metabolite levels can increase our understanding of metabolism, its controls, and the indirect effects of metabolism on other quantitative traits. Recent technological advances in both DNA sequencing and metabolite profiling allow the use of broad-spectrum, untargeted metabolite profiling to generate phenotypic data for genome-wide association studies that investigate quantitative genetic control of metabolism within species. We conducted a genome-wide association study of natural variation in plant metabolism using the results of untargeted metabolite analyses performed on a collection of wild Arabidopsis thaliana accessions. Testing 327 metabolites against >200,000 single nucleotide polymorphisms identified numerous genotype-metabolite associations distributed non-randomly within the genome. These clusters of genotype-metabolite associations (hotspots) included regions of the A. thaliana genome previously identified as subject to recent strong positive selection (selective sweeps) and regions showing trans-linkage to these putative sweeps, suggesting that these selective forces have impacted genome-wide control of A. thaliana metabolism. Comparing the metabolic variation detected within this collection of wild accessions to a laboratory-derived population of recombinant inbred lines (derived from two of the accessions used in this study) showed that the higher level of genetic variation present within the wild accessions did not correspond to higher variance in metabolic phenotypes, suggesting that evolutionary constraints limit metabolic variation. While a major goal of genome-wide association studies is to develop catalogues of intraspecific variation, the results of multiple independent experiments performed for this study showed that the genotype-metabolite associations identified are sensitive to environmental fluctuations. Thus, studies of intraspecific variation conducted via genome-wide association will require analyses of genotype by environment interaction. Interestingly, the network structure of metabolite linkages was also sensitive to environmental differences, suggesting that key aspects of network architecture are malleable.

Deficiencies in jasmonate-mediated plant defense reveal quantitative variation in Botrytis cinerea pathogenesis.

Despite the described central role of jasmonate signaling in plant defense against necrotrophic pathogens, the existence of intraspecific variation in pathogen capacity to activate or evade plant jasmonate-mediated defenses is rarely considered. Experimental infection of jasmonate-deficient and jasmonate-insensitive Arabidopsis thaliana with diverse isolates of the necrotrophic fungal pathogen Botrytis cinerea revealed pathogen variation for virulence inhibition by jasmonate-mediated plant defenses and induction of plant defense metabolites. Comparison of the transcriptional effects of infection by two distinct B. cinerea isolates showed only minor differences in transcriptional responses of wild-type plants, but notable isolate-specific transcript differences in jasmonate-insensitive plants. These transcriptional differences suggest B. cinerea activation of plant defenses that require plant jasmonate signaling for activity in response to only one of the two B. cinerea isolates tested. Thus, similar infection phenotypes observed in wild-type plants result from different signaling interactions with the plant that are likely integrated by jasmonate signaling.

A complex interplay of three R2R3 MYB transcription factors determines the profile of aliphatic glucosinolates in Arabidopsis.

While R2R3 MYB transcription factors are a large gene family of transcription factors within plants, comprehensive functional data in planta are still scarce. A model for studying R2R3 MYB control of metabolic networks is the glucosinolates (GLSs), secondary metabolites that control plant resistance against insects and pathogens and carry cancer-preventive properties. Three related members of the R2R3 MYB transcription factor family within Arabidopsis (Arabidopsis thaliana), MYB28, MYB29, and MYB76, are the commonly defined regulators of aliphatic GLS biosynthesis. We utilized new genotypes and systems analysis techniques to test the existing regulatory model in which MYB28 is the dominant regulator, MYB29 plays a minor rheostat role, and MYB76 is largely uninvolved. We unequivocally show that MYB76 is not dependent on MYB28 and MYB29 for induction of aliphatic GLSs and that MYB76 plays a role in determining the spatial distribution of aliphatic GLSs within the leaf, pointing at a potential role of MYB76 in transport regulation. Transcriptional profiling of knockout mutants revealed that GLS metabolite levels are uncoupled from the level of transcript accumulation for aliphatic GLS biosynthetic genes. This uncoupling of chemotypes from biosynthetic transcripts suggests revising our view of the regulation of GLS metabolism from a simple linear transcription factor-promoter model to a more modular system in which transcription factors cause similar chemotypes via nonoverlapping regulatory patterns. Similar regulatory networks might exist in other secondary pathways.

Understanding the evolution of defense metabolites in Arabidopsis thaliana using genome-wide association mapping.

With the improvement and decline in cost of high-throughput genotyping and phenotyping technologies, genome-wide association (GWA) studies are fast becoming a preferred approach for dissecting complex quantitative traits. Glucosinolate (GSL) secondary metabolites within Arabidopsis spp. can serve as a model system to understand the genomic architecture of quantitative traits. GSLs are key defenses against insects in the wild and the relatively large number of cloned quantitative trait locus (QTL) controlling GSL traits allows comparison of GWA to previous QTL analyses. To better understand the specieswide genomic architecture controlling plant-insect interactions and the relative strengths of GWA and QTL studies, we conducted a GWA mapping study using 96 A. thaliana accessions, 43 GSL phenotypes, and approximately 230,000 SNPs. Our GWA analysis identified the two major polymorphic loci controlling GSL variation (AOP and MAM) in natural populations within large blocks of positive associations encompassing dozens of genes. These blocks of positive associations showed extended linkage disequilibrium (LD) that we hypothesize to have arisen from balancing or fluctuating selective sweeps at both the AOP and MAM loci. These potential sweep blocks are likely linked with the formation of new defensive chemistries that alter plant fitness in natural environments. Interestingly, this GWA analysis did not identify the majority of previously identified QTL even though these polymorphisms were present in the GWA population. This may be partly explained by a nonrandom distribution of phenotypic variation across population subgroups that links population structure and GSL variation, suggesting that natural selection can hinder the detection of phenotype-genotype associations in natural populations.

The chromatin remodeler SPLAYED regulates specific stress signaling pathways.

Organisms are continuously exposed to a myriad of environmental stresses. Central to an organism's survival is the ability to mount a robust transcriptional response to the imposed stress. An emerging mechanism of transcriptional control involves dynamic changes in chromatin structure. Alterations in chromatin structure are brought about by a number of different mechanisms, including chromatin modifications, which covalently modify histone proteins; incorporation of histone variants; and chromatin remodeling, which utilizes ATP hydrolysis to alter histone-DNA contacts. While considerable insight into the mechanisms of chromatin remodeling has been gained, the biological role of chromatin remodeling complexes beyond their function as regulators of cellular differentiation and development has remained poorly understood. Here, we provide genetic, biochemical, and biological evidence for the critical role of chromatin remodeling in mediating plant defense against specific biotic stresses. We found that the Arabidopsis SWI/SNF class chromatin remodeling ATPase SPLAYED (SYD) is required for the expression of selected genes downstream of the jasmonate (JA) and ethylene (ET) signaling pathways. SYD is also directly recruited to the promoters of several of these genes. Furthermore, we show that SYD is required for resistance against the necrotrophic pathogen Botrytis cinerea but not the biotrophic pathogen Pseudomonas syringae. These findings demonstrate not only that chromatin remodeling is required for selective pathogen resistance, but also that chromatin remodelers such as SYD can regulate specific pathways within biotic stress signaling networks.

Complex genetics control natural variation in Arabidopsis thaliana resistance to Botrytis cinerea.

The genetic architecture of plant defense against microbial pathogens may be influenced by pathogen lifestyle. While plant interactions with biotrophic pathogens are frequently controlled by the action of large-effect resistance genes that follow classic Mendelian inheritance, our study suggests that plant defense against the necrotrophic pathogen Botrytis cinerea is primarily quantitative and genetically complex. Few studies of quantitative resistance to necrotrophic pathogens have used large plant mapping populations to dissect the genetic structure of resistance. Using a large structured mapping population of Arabidopsis thaliana, we identified quantitative trait loci influencing plant response to B. cinerea, measured as expansion of necrotic lesions on leaves and accumulation of the antimicrobial compound camalexin. Testing multiple B. cinerea isolates, we identified 23 separate QTL in this population, ranging in isolate-specificity from being identified with a single isolate to controlling resistance against all isolates tested. We identified a set of QTL controlling accumulation of camalexin in response to pathogen infection that largely colocalized with lesion QTL. The identified resistance QTL appear to function in epistatic networks involving three or more loci. Detection of multilocus connections suggests that natural variation in specific signaling or response networks may control A. thaliana-B. cinerea interaction in this population.

Biochemical networks and epistasis shape the Arabidopsis thaliana metabolome.

Genomic approaches have accelerated the study of the quantitative genetics that underlie phenotypic variation. These approaches associate genome-scale analyses such as transcript profiling with targeted phenotypes such as measurements of specific metabolites. Additionally, these approaches can help identify uncharacterized networks or pathways. However, little is known about the genomic architecture underlying data sets such as metabolomics or the potential of such data sets to reveal networks. To describe the genetic regulation of variation in the Arabidopsis thaliana metabolome and test our ability to integrate unknown metabolites into biochemical networks, we conducted a replicated metabolomic analysis on 210 lines of an Arabidopsis population that was previously used for targeted metabolite quantitative trait locus (QTL) and global expression QTL analysis. Metabolic traits were less heritable than the average transcript trait, suggesting that there are differences in the power to detect QTLs between transcript and metabolite traits. We used statistical analysis to identify a large number of metabolite QTLs with moderate phenotypic effects and found frequent epistatic interactions controlling a majority of the variation. The distribution of metabolite QTLs across the genome included 11 QTL clusters; 8 of these clusters were associated in an epistatic network that regulated plant central metabolism. We also generated two de novo biochemical network models from the available data, one of unknown function and the other associated with central plant metabolism.

Linking metabolic QTLs with network and cis-eQTLs controlling biosynthetic pathways.

Phenotypic variation between individuals of a species is often under quantitative genetic control. Genomic analysis of gene expression polymorphisms between individuals is rapidly gaining popularity as a way to query the underlying mechanistic causes of variation between individuals. However, there is little direct evidence of a linkage between global gene expression polymorphisms and phenotypic consequences. In this report, we have mapped quantitative trait loci (QTLs)-controlling glucosinolate content in a population of 403 Arabidopsis Bay x Sha recombinant inbred lines, 211 of which were previously used to identify expression QTLs controlling the transcript levels of biosynthetic genes. In a comparative study, we have directly tested two plant biosynthetic pathways for association between polymorphisms controlling biosynthetic gene transcripts and the resulting metabolites within the Arabidopsis Bay x Sha recombinant inbred line population. In this analysis, all loci controlling expression variation also affected the accumulation of the resulting metabolites. In addition, epistasis was detected more frequently for metabolic traits compared to transcript traits, even when both traits showed similar distributions. An analysis of candidate genes for QTL-controlling networks of transcripts and metabolites suggested that the controlling factors are a mix of enzymes and regulatory factors. This analysis showed that regulatory connections can feedback from metabolism to transcripts. Surprisingly, the most likely major regulator of both transcript level for nearly the entire pathway and aliphatic glucosinolate accumulation is variation in the last enzyme in the biosynthetic pathway, AOP2. This suggests that natural variation in transcripts may significantly impact phenotypic variation, but that natural variation in metabolites or their enzymatic loci can feed back to affect the transcripts.

Elevated genetic variation within virulence-associated Botrytis cinerea polygalacturonase loci.

Botrytis cinerea, or gray mold, is a necrotrophic fungal pathogen of hundreds of plant species. The genetic diversity of B. cinerea may contribute to its broad host range; however, the level and structure of genetic variation at pathogenesis-associated loci has not been described. B. cinerea possesses six distinct cell-wall-degrading polygalacturonases (PGs), enzymes of demonstrated importance to pathogenesis and interaction with host plant defenses. Sequencing a collection of 34 B. cinerea isolates at three PG-encoding loci, BcPG1, BcPG2, and BcPG3, revealed limited evidence of host-mediated genetic subdivision within loci, yet suggested differences in the action of evolutionary forces among loci. BcPG1 and BcPG2 are highly polymorphic, particularly when compared with previously published data from nonpathogenicity loci, whereas BcPG3 is relatively conserved. Sequence variation at BcPG1 and BcPG2 did not appear to be associated with virulence on Arabidopsis leaves; however, BcPG2 variation showed a statistically significant association with growth rate on pectin. Rather than providing evidence for host-mediated genetic subdivision at individual PG loci, our data support specialization among PGs and the potential diversification of PGs interacting directly with host defenses.

Secondary metabolites influence Arabidopsis/Botrytis interactions: variation in host production and pathogen sensitivity.

Numerous studies have suggested that plant/pathogen interactions are partially mediated via plant secondary metabolite production and corresponding pathogen tolerance. However, there are inconsistent reports on the ability of particular compounds to provide resistance to a pathogen. Most of these studies have focused on individual isolates of a given pathogen, suggesting that pathogens vary in their sensitivity to plant-produced toxins. We tested variability in virulence among pathogen isolates, and the impact on this by plant production of, and pathogen tolerance to, secondary metabolites. Botrytis cinerea isolates showed differing sensitivity to purified camalexin, and camalexin-sensitive isolates produced larger lesions on camalexin-deficient Arabidopsis genotypes than on the wild type. In contrast, the camalexin-insensitive isolate produced lesions of similar size on wild-type and camalexin-deficient Arabidopsis. Additional analysis with Arabidopsis secondary metabolite biosynthetic mutants suggests that Botrytis also has variable sensitivity to phenylpropanoids and glucosinolates. Furthermore, Botrytis infection generates a gradient of secondary metabolite responses emanating from the developing lesion, with the Botrytis isolate used determining the accumulation pattern. Collectively, our results indicate that Arabidopsis/Botrytis interactions are influenced at the metabolic level by variations in toxin production in the host and sensitivity in the pathogen.