Massive haplotypes underlie ecotypic differentiation in sunflowers.

Species often include multiple ecotypes that are adapted to different environments1. However, it is unclear how ecotypes arise and how their distinctive combinations of adaptive alleles are maintained despite hybridization with non-adapted populations2-4. Here, by resequencing 1,506 wild sunflowers from 3 species (Helianthus annuus, Helianthus petiolaris and Helianthus argophyllus), we identify 37 large (1-100 Mbp in size), non-recombining haplotype blocks that are associated with numerous ecologically relevant traits, as well as soil and climate characteristics. Limited recombination in these haplotype blocks keeps adaptive alleles together, and these regions differentiate sunflower ecotypes. For example, haplotype blocks control a 77-day difference in flowering between ecotypes of the silverleaf sunflower H. argophyllus (probably through deletion of a homologue of FLOWERING LOCUS T (FT)), and are associated with seed size, flowering time and soil fertility in dune-adapted sunflowers. These haplotypes are highly divergent, frequently associated with structural variants and often appear to represent introgressions from other-possibly now-extinct-congeners. These results highlight a pervasive role of structural variation in ecotypic adaptation.

The genomics of linkage drag in inbred lines of sunflower.

Crop wild relatives represent valuable sources of alleles for crop improvement, including adaptation to climate change and emerging diseases. However, introgressions from wild relatives might have deleterious effects on desirable traits, including yield, due to linkage drag. Here, we analyzed the genomic and phenotypic impacts of wild introgressions in inbred lines of cultivated sunflower to estimate the impacts of linkage drag. First, we generated reference sequences for seven cultivated and one wild sunflower genotype, as well as improved assemblies for two additional cultivars. Next, relying on previously generated sequences from wild donor species, we identified introgressions in the cultivated reference sequences, as well as the sequence and structural variants they contain. We then used a ridge-regression best linear unbiased prediction (BLUP) model to test the effects of the introgressions on phenotypic traits in the cultivated sunflower association mapping population. We found that introgression has introduced substantial sequence and structural variation into the cultivated sunflower gene pool, including >3,000 new genes. While introgressions reduced genetic load at protein-coding sequences, they mostly had negative impacts on yield and quality traits. Introgressions found at high frequency in the cultivated gene pool had larger effects than low-frequency introgressions, suggesting that the former likely were targeted by artificial selection. Also, introgressions from more distantly related species were more likely to be maladaptive than those from the wild progenitor of cultivated sunflower. Thus, breeding efforts should focus, as far as possible, on closely related and fully compatible wild relatives.

Re-evaluating Homoploid Reticulate Evolution in Helianthus Sunflowers.

Sunflowers of the genus Helianthus are models for hybridization research and contain three of the best-studied examples of homoploid hybrid speciation. To understand a broader picture of hybridization within the annual sunflowers, we used whole-genome resequencing to conduct a phylogenomic analysis and test for gene flow between lineages. We find that all annual sunflower species tested have evidence of admixture, suggesting hybridization was common during the radiation of the genus. Support for the major species tree decreases with increasing recombination rate, consistent with hybridization and introgression contributing to discordant topologies. Admixture graphs found hybridization to be associated with the origins of the three putative hybrid species (Helianthus anomalus, Helianthus deserticola, and Helianthus paradoxus). However, the hybridization events are more ancient than suggested by previous work. Furthermore, H. anomalus and H. deserticola appear to have arisen from a single hybridization event involving an unexpected donor, rather than through multiple independent events as previously proposed. This means our results are consistent with, but not definitive proof of, two ancient independent homoploid hybrid speciation events in the genus. Using a broader data set that covers the whole Helianthus genus, including perennial species, we find that signals of introgression span the genus and beyond, suggesting highly divergent introgression and/or the sorting of ancient haplotypes. Thus, Helianthus can be viewed as a syngameon in which largely reproductively isolated species are linked together by occasional or frequent gene flow.

Development and characterization of a new sunflower source of resistance to race G of Orobanche cumana Wallr. derived from Helianthus anomalus.

KEY MESSAGE: A new OrAnom1 gene introgressed in cultivated sunflower from wild Helianthus anomalus confers late post-attachment resistance to Orobanche cumana race G and maps to a target interval in Chromosome 4 where two receptor-like kinases (RLKs) have been identified in the H. anomalus genome as putative candidates. Sunflower broomrape is a parasitic weed that infects sunflower (Helianthus annuus L.) roots causing severe yield losses. Breeding for resistance is the most effective and sustainable control method. In this study, we report the identification, introgression, and genetic and physiological characterization of a new sunflower source of resistance to race G of broomrape developed from the wild annual sunflower H. anomalus (accession PI 468642). Crosses between PI 468642 and the susceptible line P21 were carried out, and the genetic study was conducted in BC1F1, BC1F2, and its derived BC1F3 populations. A BC1F5 germplasm named ANOM1 was developed through selection for race G resistance and resemblance to cultivated sunflower. The resistant trait showed monogenic and dominant inheritance. The gene, named OrAnom1, was mapped to Chromosome 4 within a 1.2 cM interval and co-segregated with 7 SNP markers. This interval corresponds to a 1.32 Mb region in the sunflower reference genome, housing a cluster of receptor-like kinase and receptor-like protein (RLK-RLP) genes. Notably, the analysis of the H. anomalus genome revealed the absence of RLPs in the OrAnom1 target region but featured two RLKs as possible OrAnom1 candidates. Rhizotron and histological studies showed that OrAnom1 determines a late post-attachment resistance mechanism. Broomrape can establish a vascular connection with the host, but parasite growth is stopped before tubercle development, showing phenolic compounds accumulation and tubercle necrosis. ANOM1 will contribute to broadening the genetic basis of broomrape resistance in the cultivated sunflower pool and to a better understanding of the molecular basis of the sunflower-broomrape interaction.

Mutation Load in Sunflower Inversions Is Negatively Correlated with Inversion Heterozygosity.

Recombination is critical both for accelerating adaptation and purging deleterious mutations. Chromosomal inversions can act as recombination modifiers that suppress local recombination in heterozygotes and thus, under some conditions, are predicted to accumulate such mutations. In this study, we investigated patterns of recombination, transposable element abundance, and coding sequence evolution across the genomes of 1,445 individuals from three sunflower species, as well as within nine inversions segregating within species. We also analyzed the effects of inversion genotypes on 87 phenotypic traits to test for overdominance. We found significant negative correlations of long terminal repeat retrotransposon abundance and deleterious mutations with recombination rates across the genome in all three species. However, we failed to detect an increase in these features in the inversions, except for a modest increase in the proportion of stop codon mutations in several very large or rare inversions. Consistent with this finding, there was little evidence of overdominance of inversions in phenotypes that may relate to fitness. On the other hand, significantly greater load was observed for inversions in populations polymorphic for a given inversion compared to populations monomorphic for one of the arrangements, suggesting that the local state of inversion polymorphism affects deleterious load. These seemingly contradictory results can be explained by the low frequency of inversion heterozygotes in wild sunflower populations, apparently due to divergent selection and associated geographic structure. Inversions contributing to local adaptation represent ideal recombination modifiers, acting to facilitate adaptive divergence with gene flow, while largely escaping the accumulation of deleterious mutations.

Wild sunflower goes viral: Citizen science and comparative genomics allow tracking the origin and establishment of invasive sunflower in the Levant.

Globalization and intensified volume of trade and transport around the world are accelerating the rate of biological invasions. It is therefore increasingly important to understand the processes through which invasive species colonize new habitats, often to the detriment of native flora. The initial steps of an invasion are particularly critical, as the introduced species relies on limited genetic diversity to adapt to a new environment. However, our understanding of this critical stage of the invasion is currently limited. We used a citizen science approach and social media to survey the distribution of invasive sunflower in Israel. We then sampled and sequenced a representative collection and compared it with available genomic data sets of North American wild sunflower, landraces and cultivars. We show that invasive wild sunflower is rapidly establishing throughout Israel, probably from a single, recent introduction from Texas, while maintaining high genetic diversity through ongoing gene flow. Since its introduction, invasive sunflower has spread quickly to most regions, and differentiation was detected despite extensive gene flow between clusters. Our findings suggest that rapid spread followed by continuous gene flow between diverging populations can serve as an efficient mechanism for maintaining sufficient genetic diversity at the early stages of invasion, promoting rapid adaptation and establishment in the new territory.

Standing variation rather than recent adaptive introgression probably underlies differentiation of the texanus subspecies of Helianthus annuus.

The origins of geographic races in wide-ranging species are poorly understood. In Texas, the texanus subspecies of Helianthus annuus has long been thought to have acquired its defining phenotypic traits via introgression from a local congener, H. debilis, but previous tests of this hypothesis were inconclusive. Here, we explore the origins of H. a. texanus using whole genome sequencing data from across the entire range of H. annuus and possible donor species, as well as phenotypic data from a common garden study. We found that although it is morphologically convergent with H. debilis, H. a. texanus has conflicting signals of introgression. Genome wide tests (Patterson's D and TreeMix) only found evidence of introgression from H. argophyllus (sister species to H. annuus and also sympatric), but not H. debilis, with the exception of one individual of 109 analysed. We further scanned the genome for localized signals of introgression using PCAdmix and found minimal but nonzero introgression from H. debilis and significant introgression from H. argophyllus in some populations. Given the paucity of introgression from H. debilis, we argue that the morphological convergence observed in Texas is probably from standing genetic variation. We also found that genomic differentiation in H. a. texanus is mostly driven by large segregating inversions, several of which have signatures of natural selection based on haplotype frequencies.

Parallel Evolution of Common Allelic Variants Confers Flowering Diversity in Capsella rubella.

Flowering time is an adaptive life history trait. Capsella rubella, a close relative of Arabidopsis thaliana and a young species, displays extensive variation for flowering time but low standing genetic variation due to an extreme bottleneck event, providing an excellent opportunity to understand how phenotypic diversity can occur with a limited initial gene pool. Here, we demonstrate that common allelic variation and parallel evolution at the FLC locus confer variation in flowering time in C. rubella. We show that two overlapping deletions in the 5' untranslated region (UTR) of C. rubella FLC, which are associated with local changes in chromatin conformation and histone modifications, reduce its expression levels and promote flowering. We further show that these two pervasive variants originated independently in natural C. rubella populations after speciation and spread to an intermediate frequency, suggesting a role of this parallel cis-regulatory change in adaptive evolution. Our results provide an example of how parallel mutations in the same 5' UTR region can shape phenotypic evolution in plants.

Modulation of ACD6 dependent hyperimmunity by natural alleles of an Arabidopsis thaliana NLR resistance gene.

Plants defend themselves against pathogens by activating an array of immune responses. Unfortunately, immunity programs may also cause unintended collateral damage to the plant itself. The quantitative disease resistance gene ACCELERATED CELL DEATH 6 (ACD6) serves to balance growth and pathogen resistance in natural populations of Arabidopsis thaliana. An autoimmune allele, ACD6-Est, which strongly reduces growth under specific laboratory conditions, is found in over 10% of wild strains. There is, however, extensive variation in the strength of the autoimmune phenotype expressed by strains with an ACD6-Est allele, indicative of genetic modifiers. Quantitative genetic analysis suggests that ACD6 activity can be modulated in diverse ways, with different strains often carrying different large-effect modifiers. One modifier is SUPPRESSOR OF NPR1-1, CONSTITUTIVE 1 (SNC1), located in a highly polymorphic cluster of nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes, which are prototypes for qualitative disease resistance genes. Allelic variation at SNC1 correlates with ACD6-Est activity in multiple accessions, and a common structural variant affecting the NL linker sequence can explain differences in SNC1 activity. Taken together, we find that an NLR gene can mask the activity of an ACD6 autoimmune allele in natural A. thaliana populations, thereby linking different arms of the plant immune system.

Novel allelic variants in ACD6 cause hybrid necrosis in local collection of Arabidopsis thaliana.

Hybrid necrosis is a common type of hybrid incompatibility in plants. This phenomenon is caused by deleterious epistatic interactions, resulting in spontaneous activation of plant defenses associated with leaf necrosis, stunted growth and reduced fertility in hybrids. Specific combinations of alleles of ACCELERATED CELL DEATH 6 (ACD6) have been shown to be a common cause of hybrid necrosis in Arabidopsis thaliana. Increased ACD6 activity confers broad-spectrum resistance against biotrophic pathogens but reduces biomass production. We generated 996 crosses among individuals derived from a single collection area around Tübingen (Germany) and screened them for hybrid necrosis. Necrotic hybrids were further investigated by genetic linkage, amiRNA silencing, genomic complementation and metabolic profiling. Restriction site associated DNA (RAD)-sequencing was used to understand genetic diversity in the collection sites containing necrosis-inducing alleles. Novel combinations of ACD6 alleles found in neighbouring stands were found to activate the A. thaliana immune system. In contrast to what we observed in controlled conditions, necrotic hybrids did not show reduced fitness in the field. Metabolic profiling revealed changes associated with the activation of the immune system in ACD6-dependent hybrid necrosis. This study expands our current understanding of the active role of ACD6 in mediating trade-offs between defense responses and growth in A.  thaliana.

The genetic dimension of pest pressure in the tropical rainforest.

Wet tropical forests are among the most diverse ecosystems on Earth and can host several hundreds of tree species per hectare. To maintain such diversity, the community must contain large numbers of relatively rare species rather than be dominated by a few very common trees, as is often the case in temperate forests. Explaining the mechanisms preventing dominance by common species has been a major task of tropical forest ecology. One of the most promising mechanisms is negative density dependence (NDD) of tree abundance driven by pests, including fungal diseases ('pest pressure'). NDD entails that the chance of survival of a sapling increases with the distance from a mature tree of the same species, thus preventing species from becoming locally dominant. Curiously, the strength of NDD is negatively correlated with abundance, meaning that tree species that are more common generally show weaker NDD (Comita et al. ). Interactions between plants and soil pathogens have been shown to play an important role in NDD (Klironomos ), and rare species are apparently more strongly affected (Mangan et al. ). However, the genetic mechanisms underlying this phenomenon have remained obscure. In this issue of Molecular Ecology, Marden et al. () suggest that reduced diversity of the genes involved in pathogen recognition (Resistance genes or R genes) could explain why NDD is stronger in locally rare species.

Activation of the Arabidopsis thaliana immune system by combinations of common ACD6 alleles.

A fundamental question in biology is how multicellular organisms distinguish self and non-self. The ability to make this distinction allows animals and plants to detect and respond to pathogens without triggering immune reactions directed against their own cells. In plants, inappropriate self-recognition results in the autonomous activation of the immune system, causing affected individuals to grow less well. These plants also suffer from spontaneous cell death, but are at the same time more resistant to pathogens. Known causes for such autonomous activation of the immune system are hyperactive alleles of immune regulators, or epistatic interactions between immune regulators and unlinked genes. We have discovered a third class, in which the Arabidopsis thaliana immune system is activated by interactions between natural alleles at a single locus, ACCELERATED CELL DEATH 6 (ACD6). There are two main types of these interacting alleles, one of which has evolved recently by partial resurrection of a pseudogene, and each type includes multiple functional variants. Most previously studies hybrid necrosis cases involve rare alleles found in geographically unrelated populations. These two types of ACD6 alleles instead occur at low frequency throughout the range of the species, and have risen to high frequency in the Northeast of Spain, suggesting a role in local adaptation. In addition, such hybrids occur in these populations in the wild. The extensive functional variation among ACD6 alleles points to a central role of this locus in fine-tuning pathogen defenses in natural populations.

Natural allelic variation underlying a major fitness trade-off in Arabidopsis thaliana.

Plants can defend themselves against a wide array of enemies, from microbes to large animals, yet there is great variability in the effectiveness of such defences, both within and between species. Some of this variation can be explained by conflicting pressures from pathogens with different modes of attack. A second explanation comes from an evolutionary 'tug of war', in which pathogens adapt to evade detection, until the plant has evolved new recognition capabilities for pathogen invasion. If selection is, however, sufficiently strong, susceptible hosts should remain rare. That this is not the case is best explained by costs incurred from constitutive defences in a pest-free environment. Using a combination of forward genetics and genome-wide association analyses, we demonstrate that allelic diversity at a single locus, ACCELERATED CELL DEATH 6 (ACD6), underpins marked pleiotropic differences in both vegetative growth and resistance to microbial infection and herbivory among natural Arabidopsis thaliana strains. A hyperactive ACD6 allele, compared to the reference allele, strongly enhances resistance to a broad range of pathogens from different phyla, but at the same time slows the production of new leaves and greatly reduces the biomass of mature leaves. This allele segregates at intermediate frequency both throughout the worldwide range of A. thaliana and within local populations, consistent with this allele providing substantial fitness benefits despite its marked impact on growth.

Genome-wide association study of 107 phenotypes in Arabidopsis thaliana inbred lines.

Although pioneered by human geneticists as a potential solution to the challenging problem of finding the genetic basis of common human diseases, genome-wide association (GWA) studies have, owing to advances in genotyping and sequencing technology, become an obvious general approach for studying the genetics of natural variation and traits of agricultural importance. They are particularly useful when inbred lines are available, because once these lines have been genotyped they can be phenotyped multiple times, making it possible (as well as extremely cost effective) to study many different traits in many different environments, while replicating the phenotypic measurements to reduce environmental noise. Here we demonstrate the power of this approach by carrying out a GWA study of 107 phenotypes in Arabidopsis thaliana, a widely distributed, predominantly self-fertilizing model plant known to harbour considerable genetic variation for many adaptively important traits. Our results are dramatically different from those of human GWA studies, in that we identify many common alleles of major effect, but they are also, in many cases, harder to interpret because confounding by complex genetics and population structure make it difficult to distinguish true associations from false. However, a-priori candidates are significantly over-represented among these associations as well, making many of them excellent candidates for follow-up experiments. Our study demonstrates the feasibility of GWA studies in A. thaliana and suggests that the approach will be appropriate for many other organisms.

Target mimicry provides a new mechanism for regulation of microRNA activity.

MicroRNAs (miRNA) regulate key aspects of development and physiology in animals and plants. These regulatory RNAs act as guides of effector complexes to recognize specific mRNA sequences based on sequence complementarity, resulting in translational repression or site-specific cleavage. In plants, most miRNA targets are cleaved and show almost perfect complementarity with the miRNAs around the cleavage site. Here, we examined the non-protein coding gene IPS1 (INDUCED BY PHOSPHATE STARVATION 1) from Arabidopsis thaliana. IPS1 contains a motif with sequence complementarity to the phosphate (Pi) starvation-induced miRNA miR-399, but the pairing is interrupted by a mismatched loop at the expected miRNA cleavage site. We show that IPS1 RNA is not cleaved but instead sequesters miR-399. Thus, IPS1 overexpression results in increased accumulation of the miR-399 target PHO2 mRNA and, concomitantly, in reduced shoot Pi content. Engineering of IPS1 to be cleavable abolishes its inhibitory activity on miR-399. We coin the term 'target mimicry' to define this mechanism of inhibition of miRNA activity. Target mimicry can be generalized beyond the control of Pi homeostasis, as demonstrated using artificial target mimics.

On reconciling the interactions between APETALA2, miR172 and AGAMOUS with the ABC model of flower development.

The ABC model of flower development explains how three classes of homeotic genes confer identity to the four types of floral organs. In Arabidopsis thaliana, APETALA2 (AP2) and AGAMOUS (AG) represent A- and C-class genes that act in an antagonistic fashion to specify perianth and reproductive organs, respectively. An apparent paradox was the finding that AP2 mRNA is supposedly uniformly distributed throughout young floral primordia. Although miR172 has a role in preventing AP2 protein accumulation, miR172 was reported to disappear from the periphery only several days after AG activation in the center of the flower. Here, we resolve the enigmatic behavior of AP2 and its negative regulator miR172 through careful expression analyses. We find that AP2 mRNA accumulates predominantly in the outer floral whorls, as expected for an A-class homeotic gene. Its pattern overlaps only transiently with that of miR172, which we find to be restricted to the center of young floral primordia from early stages on. MiR172 also accumulates in the shoot meristem upon floral induction, compatible with its known role in regulating AP2-related genes with a role in flowering. Furthermore, we show that AP2 can cause striking organ proliferation defects that are not limited to the center of the floral meristem, where its antagonist AG is required for terminating stem cell proliferation. Moreover, AP2 never expands uniformly into the center of ag mutant flowers, while miR172 is largely unaffected by loss of AG activity. We present a model in which the decision whether stamens or petals develop is based on the balance between AP2 and AG activities, rather than the two being mutually exclusive.

A collection of target mimics for comprehensive analysis of microRNA function in Arabidopsis thaliana.

Many targets of plant microRNAs (miRNAs) are thought to play important roles in plant physiology and development. However, because plant miRNAs are typically encoded by medium-size gene families, it has often been difficult to assess their precise function. We report the generation of a large-scale collection of knockdowns for Arabidopsis thaliana miRNA families; this has been achieved using artificial miRNA target mimics, a recently developed technique fashioned on an endogenous mechanism of miRNA regulation. Morphological defects in the aerial part were observed for approximately 20% of analyzed families, all of which are deeply conserved in land plants. In addition, we find that non-cleavable mimic sites can confer translational regulation in cis. Phenotypes of plants expressing target mimics directed against miRNAs involved in development were in several cases consistent with previous reports on plants expressing miRNA-resistant forms of individual target genes, indicating that a limited number of targets mediates most effects of these miRNAs. That less conserved miRNAs rarely had obvious effects on plant morphology suggests that most of them do not affect fundamental aspects of development. In addition to insight into modes of miRNA action, this study provides an important resource for the study of miRNA function in plants.

Repeatability of adaptation in sunflowers reveals that genomic regions harbouring inversions also drive adaptation in species lacking an inversion.

Local adaptation commonly involves alleles of large effect, which experience fitness advantages when in positive linkage disequilibrium (LD). Because segregating inversions suppress recombination and facilitate the maintenance of LD between locally adapted loci, they are also commonly found to be associated with adaptive divergence. However, it is unclear what fraction of an adaptive response can be attributed to inversions and alleles of large effect, and whether the loci within an inversion could still drive adaptation in the absence of its recombination-suppressing effect. Here, we use genome-wide association studies to explore patterns of local adaptation in three species of sunflower: Helianthus annuus, Helianthus argophyllus, and Helianthus petiolaris, which each harbour a large number of species-specific inversions. We find evidence of significant genome-wide repeatability in signatures of association to phenotypes and environments, which are particularly enriched within regions of the genome harbouring an inversion in one species. This shows that while inversions may facilitate local adaptation, at least some of the loci can still harbour mutations that make substantial contributions without the benefit of recombination suppression in species lacking a segregating inversion. While a large number of genomic regions show evidence of repeated adaptation, most of the strongest signatures of association still tend to be species-specific, indicating substantial genotypic redundancy for local adaptation in these species.

In plants, like in humans, DNA is arranged into sections known as genes that are in turn organised into structures called chromosomes. Mutations that modify the activity of these genes can help plant species to adapt to a new environment or to extreme conditions such as drought. However, successful adaptation often requires changes in many different genes. If these sets of genes are located close to each other on the same chromosome, any mutations will likely be passed onto the next generation together. If the genes are located further away, or even on different chromosomes, they may instead be inherited separately so that the next generation does not benefit as much from the adaptation. A chromosome inversion ­ when a segment of chromosome breaks off and reattaches the other way around ­ can increase the likelihood that sets of mutations on the same chromosome will be inherited together. Many previous studies have found that chromosome inversions tend to drive the ability of species to adapt to different environments by keeping together mutations that affect the same characteristics. However, it is not clear how inversions affect the repeatability of the adaptation, that is, if another group of closely related plants faced the same challenge in their environment would they evolve in the same way, or would they evolve a new response? To address this question, Soudi, Jahani et al. used a genetics approach known as a genome wide association study to explore how three closely related species of sunflower have adapted to their respective environments. Two of the species grow in various environments across the centre and west of the USA that are often hot and dry, whereas the third species is restricted to the more humid coastal plain of Texas, USA. The experiments found that a few key genes had changed in all three sunflower species. However, each species also had mutations in a larger set of unique genes that were not changed in the other species. Regions of chromosomes harbouring inversions in one of the species tended to have more of the key genes within them, compared to other genomic regions. This was also true for species that did not have inversions in those regions. This demonstrates that genes in regions affected by chromosome inversions can still help plants adapt to changes in the environment even in the absence of inversions. Sunflowers are widely grown for their edible oily seeds. In the future, some of the key genes identified in this work may be useful candidates for plant breeding to improve the resilience of sunflowers to drought, high temperatures and other environmental challenges.

Interplay of SLIM1 and miR395 in the regulation of sulfate assimilation in Arabidopsis.

MicroRNAs play a key role in the control of plant development and response to adverse environmental conditions. For example, microRNA395 (miR395), which targets three out of four isoforms of ATP sulfurylase, the first enzyme of sulfate assimilation, as well as a low-affinity sulfate transporter, SULTR2;1, is strongly induced by sulfate deficiency. However, other components of sulfate assimilation are induced by sulfate starvation, so that the role of miR395 is counterintuitive. Here, we describe the regulation of miR395 and its targets by sulfate starvation. We show that miR395 is important for the increased translocation of sulfate to the shoots during sulfate starvation. MiR395 together with the SULFUR LIMITATION 1 transcription factor maintain optimal levels of ATP sulfurylase transcripts to enable increased flux through the sulfate assimilation pathway in sulfate-deficient plants. Reduced expression of ATP sulfurylase (ATPS) alone affects both sulfate translocation and flux, but SULTR2;1 is important for the full rate of sulfate translocation to the shoots. Thus, miR395 is an integral part of the regulatory circuit controlling plant sulfate assimilation with a complex mechanism of action.

HeliantHOME, a public and centralized database of phenotypic sunflower data.

Genomic studies often attempt to link natural genetic variation with important phenotypic variation. To succeed, robust and reliable phenotypic data, as well as curated genomic assemblies, are required. Wild sunflowers, originally from North America, are adapted to diverse and often extreme environments and have historically been a widely used model plant system for the study of population genomics, adaptation, and speciation. Moreover, cultivated sunflower, domesticated from a wild relative (Helianthus annuus) is a global oil crop, ranking fourth in production of vegetable oils worldwide. Public availability of data resources both for the plant research community and for the associated agricultural sector, are extremely valuable. We have created HeliantHOME ( http://www.helianthome.org ), a curated, public, and interactive database of phenotypes including developmental, structural and environmental ones, obtained from a large collection of both wild and cultivated sunflower individuals. Additionally, the database is enriched with external genomic data and results of genome-wide association studies. Finally, being a community open-source platform, HeliantHOME is expected to expand as new knowledge and resources become available.