Sequence-Based Analysis of Structural Organization and Composition of the Cultivated Sunflower (Helianthus annuus L.) Genome

The Sunflower WRINKLED1 Transcription Factor Regulates Fatty Acid Biosynthesis Genes through an AW Box Binding Sequence with a Particular Base Bias

Sunflower is an important oilseed crop in which the biochemical pathways leading to seed oil synthesis and accumulation have been widely studied. However, how these pathways are regulated is less well understood. The WRINKLED1 (WRI1) transcription factor is considered a key regulator in the control of triacylglycerol biosynthesis, acting through the AW box binding element (CNTNG(N)₇CG). Here, we identified the sunflower WRI1 gene and characterized its activity in electrophoretic mobility shift assays. We studied its role as a co-regulator of sunflower genes involved in plastidial fatty acid synthesis. Sunflower WRI1-targets included genes encoding the pyruvate dehydrogenase complex, the α-CT and BCCP genes, genes encoding ACPs and the fatty acid synthase complex, together with the FATA1 gene. As such, sunflower WRI1 regulates genes involved in seed plastidial fatty acid biosynthesis in a coordinated manner, establishing a WRI1 push and pull strategy that drives oleic acid synthesis for its export into the cytosol. We also determined the base bias at the N positions in the active sunflower AW box motif. The sunflower AW box is sequence-sensitive at the non-conserved positions, enabling WRI1-binding. Moreover, sunflower WRI1 could bind to a non-canonical AW-box motif, opening the possibility of searching for new target genes.

High genetic and epigenetic variation of transposable elements: Potential drivers to rapid adaptive evolution for the noxious invasive weed Mikania micrantha

Why invasive species can rapidly adapt to novel environments is a puzzling question known as the genetic paradox of invasive species. This paradox is explainable in terms of transposable elements (TEs) activity, which are theorized to be powerful mutational forces to create genetic variation. Mikania micrantha, a noxious invasive weed, in this sense provides an excellent opportunity to test the explanation. The genetic and epigenetic variation of 21 invasive populations of M. micrantha in southern China have been examined by using transposon display (TD) and transposon methylation display (TMD) techniques to survey 12 TE superfamilies. Our results showed that M. micrantha populations maintained an almost equally high level of TE-based genetic and epigenetic variation and they have been differentiated into subpopulations genetically and epigenetically. A similar positive spatial genetic and epigenetic structure pattern was observed within 300 m. Six and seven TE superfamilies presented significant genetic and epigenetic isolation by distance (IBD) pattern. In total, 59 genetic and 86 epigenetic adaptive TE loci were identified. Of them, 51 genetic and 44 epigenetic loci were found to correlate with 25 environmental variables (including precipitation, temperature, vegetation coverage, and soil metals). Twenty-five transposon-inserted genes were sequenced and homology-based annotated, which are found to be involved in a variety of molecular and cellular functions. Our research consolidates the importance of TE-associated genetic and epigenetic variation in the rapid adaptation and invasion of M. micrantha.

Diversity and evolution of the repetitive genomic content in Cannabis sativa

Background

The repetitive content of the genome, once considered to be “junk DNA”, is in fact an essential component of genomic architecture and evolution. In this study, we used the genomes of three varieties of Cannabis sativa, three varieties of Humulus lupulus and one genotype of Morus notabilis to explore their repetitive content using a graph-based clustering method, designed to explore and compare repeat content in genomes that have not been fully assembled.

Results

The repetitive content in the C. sativa genome is mainly composed of the retrotransposons LTR/Copia and LTR/Gypsy (14% and 14.8%, respectively), ribosomal DNA (2%), and low-complexity sequences (29%). We observed a recent copy number expansion in some transposable element families. Simple repeats and low complexity regions of the genome show higher intra and inter species variation.

Conclusions

As with other sequenced genomes, the repetitive content of C. sativa’s genome exhibits a wide range of evolutionary patterns. Some repeat types have patterns of diversity consistent with expansions followed by losses in copy number, while others may have expanded more slowly and reached a steady state. Still, other repetitive sequences, particularly ribosomal DNA (rDNA), show signs of concerted evolution playing a major role in homogenizing sequence variation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4494-3) contains supplementary material, which is available to authorized users.

Genomic abundance and transcriptional activity of diverse gypsy and copia long terminal repeat retrotransposons in three wild sunflower species

BACKGROUND

Long terminal repeat (LTR) retrotransposons are highly abundant in plant genomes and require transcriptional activity for their proliferative mode of replication. These sequences exist in plant genomes as diverse sublineages within the main element superfamilies (i.e., gypsy and copia). While transcriptional activity of these elements is increasingly recognized as a regular attribute of plant transcriptomes, it is currently unknown the extent to which different sublineages of these elements are transcriptionally active both within and across species. In the current report, we utilize next generation sequencing methods to examine genomic copy number abundance of diverse LTR retrotransposon sublineages and their corresponding levels of transcriptional activity in three diploid wild sunflower species, Helianthus agrestis, H. carnosus and H. porteri.

RESULTS

The diploid sunflower species under investigation differ in genome size 2.75-fold, with 2C values of 22.93 for H. agrestis, 12.31 for H. carnosus and 8.33 for H. porteri. The same diverse gypsy and copia sublineages of LTR retrotransposons were identified across species, but with gypsy sequences consistently more abundant than copia and with global gypsy sequence abundance positively correlated with nuclear genome size. Transcriptional activity was detected for multiple copia and gypsy sequences, with significantly higher activity levels detected for copia versus gypsy. Interestingly, of 11 elements identified as transcriptionally active, 5 exhibited detectable expression in all three species and 3 exhibited detectable expression in two species.

CONCLUSIONS

Combined analyses of LTR retrotransposon genomic abundance and transcriptional activity across three sunflower species provides novel insights into genome size evolution and transposable element dynamics in this group. Despite considerable variation in nuclear genome size among species, relatively conserved patterns of LTR retrotransposon transcriptional activity were observed, with a highly overlapping set of copia and gypsy sequences observed to be transcriptionally active across species. A higher proportion of copia versus gypsy elements were found to be transcriptionally active and these sequences also were expressed at higher levels.

Long Terminal Repeat Retrotransposon Content in Eight Diploid Sunflower Species Inferred from Next-Generation Sequence Data

The most abundant transposable elements (TEs) in plant genomes are Class I long terminal repeat (LTR) retrotransposons represented by superfamilies gypsy and copia. Amplification of these superfamilies directly impacts genome structure and contributes to differential patterns of genome size evolution among plant lineages. Utilizing short-read Illumina data and sequence information from a panel of Helianthus annuus (sunflower) full-length gypsy and copia elements, we explore the contribution of these sequences to genome size variation among eight diploid Helianthus species and an outgroup taxon, Phoebanthus tenuifolius. We also explore transcriptional dynamics of these elements in both leaf and bud tissue via RT-PCR. We demonstrate that most LTR retrotransposon sublineages (i.e., families) display patterns of similar genomic abundance across species. A small number of LTR retrotransposon sublineages exhibit lineage-specific amplification, particularly in the genomes of species with larger estimated nuclear DNA content. RT-PCR assays reveal that some LTR retrotransposon sublineages are transcriptionally active across all species and tissue types, whereas others display species-specific and tissue-specific expression. The species with the largest estimated genome size, H. agrestis, has experienced amplification of LTR retrotransposon sublineages, some of which have proliferated independently in other lineages in the Helianthus phylogeny.

Diversity analysis and genomic prediction of Sclerotinia resistance in sunflower using a new 25 K SNP genotyping array

KEY MESSAGE

We have developed a SNP array for sunflower containing more than 25 K markers, representing single loci mostly in or near transcribed regions of the genome. The array was successfully applied to genotype a diversity panel of lines, hybrids, and mapping populations and represented well the genetic diversity of cultivated sunflower. Results of PCoA and population substructure analysis underlined the complexity of the genetic composition of current elite breeding material. The performance of this genotyping platform for genome-based prediction of phenotypes and detection of QTL with improved resolution could be demonstrated based on the re-evaluation of a population segregating for resistance to Sclerotinia midstalk rot. Given our results, the newly developed 25 K SNP array is expected to be of great utility for the most important applications in genome-based sunflower breeding and research.

ABSTRACT

Genotyping with a large number of molecular markers is a prerequisite to conduct genome-based genetic analyses with high precision. Here, we report the design and performance of a 25 K SNP genotyping array for sunflower (Helianthus annuus L.). SNPs were discovered based on variant calling in de novo assembled, UniGene-based contigs of sunflower derived from whole genome sequencing and amplicon sequences originating from four and 48 inbred lines, respectively. After inclusion of publically available transcriptome-derived SNPs, in silico design of the Illumina(®) Infinium iSelect HD BeadChip yielded successful assays for 22,299 predominantly haplotype-specific SNPs. The array was validated in a sunflower diversity panel including inbred lines, open-pollinated varieties, introgression lines, landraces, recombinant inbred lines, and F2 populations. Validation provided 20,502 high-quality bi-allelic SNPs with stable cluster performance whereby each SNP marker represents a single locus mostly in or near transcribed regions of the sunflower genome. Analyses of population structure and quantitative resistance to Sclerotinia midstalk rot demonstrate that this array represents a significant improvement over currently available genomic tools for genetic diversity analyses, genome-wide marker-trait association studies, and genetic mapping in sunflower.

Phototropic solar tracking in sunflower plants: an integrative perspective

BACKGROUND

One of the best-known plant movements, phototropic solar tracking in sunflower (Helianthus annuus), has not yet been fully characterized. Two questions are still a matter of debate. (1) Is the adaptive significance solely an optimization of photosynthesis via the exposure of the leaves to the sun? (2) Is shade avoidance involved in this process? In this study, these concepts are discussed from a historical perspective and novel insights are provided.

SCOPE AND METHODS

Results from the primary literature on heliotropic growth movements led to the conclusion that these responses cease before anthesis, so that the flowering heads point to the East. Based on observations on 10-week-old plants, the diurnal East-West oscillations of the upper fifth of the growing stem and leaves in relation to the position of the sun (inclusive of nocturnal re-orientation) were documented, and photon fluence rates on the leaf surfaces on clear, cloudy and rainy days were determined. In addition, the light-response curve of net CO2 assimilation was determined on the upper leaves of the same batch of plants, and evidence for the occurrence of shade-avoidance responses in growing sunflower plants is summarized.

CONCLUSIONS

Only elongating, vegetative sunflower shoots and the upper leaves perform phototropic solar tracking. Photon fluence response and CO2 assimilation measurements cast doubt on the 'photosynthesis-optimization hypothesis' as the sole explanation for the evolution of these plant movements. We suggest that the shade-avoidance response, which maximizes light-driven CO2 assimilation, plays a major role in solar tracking populations of competing sunflower plants, and an integrative scheme of these growth movements is provided.

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.

Excisions of a defective transposable CACTA element (Tetu1) generate new alleles of a CYCLOIDEA-like gene of Helianthus annuus

Tubular ray flower (turf) is a sunflower mutant that caught attention because it bears actinomorphic ray flowers, due to the presence of an active, although non-autonomous CACTA transposon (Tetu1) in the TCP domain of a CYCLOIDEA-like gene, HaCYC2c, a major regulator of sunflower floral symmetry. Here, we analyzed its excision rates in F3 population deriving from independent crosses of turf with common sunflower accessions. Our results suggest that the excision rate, ranging from 1.21 to 6.29%, depends on genetic background; moreover, the absence of somatic sectors in inflorescences of revertant individuals analyzed (182) and genetic analyses suggests a tight developmental control of Tetu1 excision, likely restricted to germinal cells. We individuate events of Tetu1 excision through molecular analysis that restore the wild type (WT) HaCYC2c allele, but even transposon excisions during which footprints are left. All mutations we detected occurred at the TCP basic motif and cause a change in ray flower phenotype. In particular, we selected five mutants with a one-to-four amino acid change that influence the capacity of reproductive organ development and ray flower corolla shaping (MUT-1, -2, -3, -4, -5). Revertant alleles not affecting turf phenotype (i.e. reading frame mutations) have also been identified (MUT-6). In all mutants, Real-time quantitative PCR (qPCR) experiments revealed variations of the steady state level of HaCYC2c mRNA. MUT-1 and MUT-4 showed a significant HaCYC2c down-regulation with respect to WT. A large variation within the biological replicates of MUT-2, MUT-3 and MUT-5 was detected and not significant differences in transcription levels between mutants and WT were observed. We detected low steady state level of HaCYC2c mRNA both in turf as in MUT-6. A three dimensional (3D) structure prediction tool let us predict an incorrect folding of the TCP protein already after a single amino acid deletion. This in turn is detectable as the restore of traits that are not peculiar of WT ray flowers, such as male fertility. Our analysis of an active TE sheds light on the TCP motif of the HaCYC2c gene and suggests that Tetu1 may be useful to obtain new natural mutants and for transposon tagging in different inbred lines of sunflower.