The effect of transposable elements on phenotypic variation: insights from plants to humans

Forces driving transposable element load variation during Arabidopsis range expansion

Genetic load refers to the accumulated and potentially life-threatening deleterious mutations in populations. Understanding the mechanisms underlying genetic load variation of transposable element (TE) insertion, a major large-effect mutation, during range expansion is an intriguing question in biology. Here, we used 1,115 global natural accessions of Arabidopsis (Arabidopsis thaliana) to study the driving forces of TE load variation during its range expansion. TE load increased with range expansion, especially in the recently established Yangtze River basin population. Effective population size, which explains 62.0% of the variance in TE load, high transposition rate, and selective sweeps contributed to TE accumulation in the expanded populations. We genetically mapped and identified multiple candidate causal genes and TEs, and revealed the genetic architecture of TE load variation. Overall, this study reveals the variation in TE genetic load during Arabidopsis expansion and highlights the causes of TE load variation from the perspectives of both population genetics and quantitative genetics.

Profiling the selected hotspots for ear traits in two maize-teosinte populations

KEY MESSAGE

Eight selected hotspots related to ear traits were identified from two maize-teosinte populations. Throughout the history of maize cultivation, ear-related traits have been selected. However, little is known about the specific genes involved in shaping these traits from their origins in the wild progenitor, teosinte, to the characteristics observed in modern maize. In this study, five ear traits (kernel row number [KRN], ear length [EL], kernel number per row [KNR], cob diameter [CD], and ear diameter [ED]) were investigated, and eight quantitative trait loci (QTL) hotspots were identified in two maize-teosinte populations. Notably, our findings revealed a significant enrichment of genes showing a selection signature and expressed in the ear in qbdCD1.1, qbdCD5.1, qbpCD2.1, qbdED1.1, qbpEL1.1, qbpEL5.1, qbdKNR1.1, and qbdKNR10.1, suggesting that these eight QTL are selected hotspots involved in shaping the maize ear. By combining the results of the QTL analysis with data from previous genome-wide association study (GWAS) involving two natural panels, we identified eight candidate selected genes related to KRN, KNR, CD, and ED. Among these, considering their expression pattern and sequence variation, Zm00001d025111, encoding a WD40/YVTN protein, was proposed as a positive regulator of KNR. This study presents a framework for understanding the genomic distribution of selected loci crucial in determining ear-related traits.

Toward Transgene-Free Transposon-Mediated Biological Mutagenesis for Plant Breeding

Genetic diversity is a key factor for plant breeding. The birth of novel genic and genomic variants is also crucial for plant adaptation in nature. Therefore, the genomes of almost all living organisms possess natural mutagenic mechanisms. Transposable elements (TEs) are a major mutagenic force driving genetic diversity in wild plants and modern crops. The relatively rare TE transposition activity during the thousand-year crop domestication process has led to the phenotypic diversity of many cultivated species. The utilization of TE mutagenesis by artificial and transient acceleration of their activity in a controlled mode is an attractive foundation for a novel type of mutagenesis called TE-mediated biological mutagenesis. Here, I focus on TEs as mutagenic sources for plant breeding and discuss existing and emerging transgene-free approaches for TE activation in plants. Furthermore, I also review the non-randomness of TE insertions in a plant genome and the molecular and epigenetic factors involved in shaping TE insertion preferences. Additionally, I discuss the molecular mechanisms that prevent TE transpositions in germline plant cells (e.g., meiocytes, pollen, egg and embryo cells, and shoot apical meristem), thereby reducing the chances of TE insertion inheritance. Knowledge of these mechanisms can expand the TE activation toolbox using novel gene targeting approaches. Finally, the challenges and future perspectives of plant populations with induced novel TE insertions (iTE plant collections) are discussed.

Comparative genome and metabolome analyses uncover the evolution and flavonoid biosynthesis between Apocynum venetum and Apocynum hendersonii

Apocynum species have great application prospects in textile and phytoremediation of saline soil, are rich in flavonoids, and possess medicinal significance. Here, we report the draft genomes of Apocynum venetum and Apocynum hendersonii, and elucidate their evolutionary relationship. The high synteny and collinearity between the two suggested that they have experienced the same WGD event. Comparative analysis revealed that flavone 3-hydroxylase (ApF3H) and differentially evolved flavonoid 3-O-glucosyl transferase (ApUFGT) genes are critical for determining natural variation in flavonoid biosynthesis between the species. Overexpression of ApF3H-1 enhanced the total flavonoid content and promoted the antioxidant capacity of transformed plants compared to the wild-type. ApUFGT5 and 6 explained the diversification of flavonoids or their derivatives. These data provide biochemical insight and knowledge on the genetic regulation of flavonoid biosynthesis, supporting the adoption of these genes in breeding programs aimed at the multipurpose utilization of the plants.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

The rubber tree kinome: Genome-wide characterization and insights into coexpression patterns associated with abiotic stress responses

The protein kinase (PK) superfamily constitutes one of the largest and most conserved protein families in eukaryotic genomes, comprising core components of signaling pathways in cell regulation. Despite its remarkable relevance, only a few kinase families have been studied in Hevea brasiliensis. A comprehensive characterization and global expression analysis of the PK superfamily, however, is currently lacking. In this study, with the aim of providing novel inferences about the mechanisms associated with the stress response developed by PKs and retained throughout evolution, we identified and characterized the entire set of PKs, also known as the kinome, present in the Hevea genome. Different RNA-sequencing datasets were employed to identify tissue-specific expression patterns and potential correspondences between different rubber tree genotypes. In addition, coexpression networks under several abiotic stress conditions, such as cold, drought and latex overexploitation, were employed to elucidate associations between families and tissues/stresses. A total of 1,809 PK genes were identified using the current reference genome assembly at the scaffold level, and 1,379 PK genes were identified using the latest chromosome-level assembly and combined into a single set of 2,842 PKs. These proteins were further classified into 20 different groups and 122 families, exhibiting high compositional similarities among family members and with two phylogenetically close species Manihot esculenta and Ricinus communis. Through the joint investigation of tandemly duplicated kinases, transposable elements, gene expression patterns, and coexpression events, we provided insights into the understanding of the cell regulation mechanisms in response to several conditions, which can often lead to a significant reduction in rubber yield.

Insertion of a TRIM-like sequence in MdFLS2-1 promoter is associated with its allele-specific expression in response to Alternaria alternata in apple

Alternaria blotch disease, caused by Alternaria alternata apple pathotype (AAAP), is one of the major fungal diseases in apple. Early field observations revealed, the anther-derived homozygote Hanfu line (HFTH1) was highly susceptible to AAAP, whereas Hanfu (HF) exhibited resistance to AAAP. To understand the molecular mechanisms underlying the difference in sensitivity of HF and HFTH1 to AAAP, we performed allele-specific expression (ASE) analysis and comparative transcriptomic analysis before and after AAAP inoculation. We reported an important immune gene, namely, MdFLS2, which displayed strong ASE in HF with much lower expression levels of HFTH1-derived alleles. Transient overexpression of the dominant allele of MdFLS2-1 from HF in GL-3 apple leaves could enhance resistance to AAAP and induce expression of genes related to salicylic acid pathway. In addition, MdFLS2-1 was identified with an insertion of an 85-bp terminal-repeat retrotransposon in miniature (TRIM) element-like sequence in the upstream region of the nonreference allele. In contrast, only one terminal direct repeat (TDR) from TRIM-like sequence was present in the upstream region of the HFTH1-derived allele MdFLS2-2. Furthermore, the results of luciferase and β-glucuronidase reporter assays demonstrated that the intact TRIM-like sequence has enhancer activity. This suggested that insertion of the TRIM-like sequence regulates the expression level of the allele of MdFLS2, in turn, affecting the sensitivity of HF and HFTH1 to AAAP.

Inpactor2: a software based on deep learning to identify and classify LTR-retrotransposons in plant genomes

LTR-retrotransposons are the most abundant repeat sequences in plant genomes and play an important role in evolution and biodiversity. Their characterization is of great importance to understand their dynamics. However, the identification and classification of these elements remains a challenge today. Moreover, current software can be relatively slow (from hours to days), sometimes involve a lot of manual work and do not reach satisfactory levels in terms of precision and sensitivity. Here we present Inpactor2, an accurate and fast application that creates LTR-retrotransposon reference libraries in a very short time. Inpactor2 takes an assembled genome as input and follows a hybrid approach (deep learning and structure-based) to detect elements, filter partial sequences and finally classify intact sequences into superfamilies and, as very few tools do, into lineages. This tool takes advantage of multi-core and GPU architectures to decrease execution times. Using the rice genome, Inpactor2 showed a run time of 5 minutes (faster than other tools) and has the best accuracy and F1-Score of the tools tested here, also having the second best accuracy and specificity only surpassed by EDTA, but achieving 28% higher sensitivity. For large genomes, Inpactor2 is up to seven times faster than other available bioinformatics tools.

Legume Pangenome: Status and Scope for Crop Improvement

In the last decade, legume genomics research has seen a paradigm shift due to advances in genome sequencing technologies, assembly algorithms, and computational genomics that enabled the construction of high-quality reference genome assemblies of major legume crops. These advances have certainly facilitated the identification of novel genetic variants underlying the traits of agronomic importance in many legume crops. Furthermore, these robust sequencing technologies have allowed us to study structural variations across the whole genome in multiple individuals and at the species level using 'pangenome analysis.' This review updates the progress of constructing pangenome assemblies for various legume crops and discusses the prospects for these pangenomes and how to harness the information to improve various traits of economic importance through molecular breeding to increase genetic gain in legumes and tackle the increasing global food crisis.

An LTR retrotransposon in the promoter of a PsMYB10.2 gene associated with the regulation of fruit flesh color in Japanese plum

Japanese plums exhibit wide diversity of fruit coloration. The red to black hues are caused by the accumulation of anthocyanins, while their absence results in yellow, orange or green fruits. In Prunus, MYB10 genes are determinants for anthocyanin accumulation. In peach, QTLs for red plant organ traits map in an LG3 region with three MYB10 copies (PpMYB10.1, PpMYB10.2 and PpMYB10.3). In Japanese plum the gene copy number in this region differs with respect to peach: there are at least three copies of PsMYB10.1, with the expression of one of them (PsMYB10.1a) correlating with fruit skin color. The objective of this study was to determine a possible role of LG3-PsMYB10 genes in the natural variability of the flesh color trait and to develop a molecular marker for marker-assisted selection (MAS). We explored the variability within the LG3-PsMYB10 region using long-range sequences obtained in previous studies through CRISPR-Cas9 enrichment sequencing. We found that the PsMYB10.2 gene was only expressed in red flesh fruits. Its role in promoting anthocyanin biosynthesis was validated by transient overexpression in Japanese plum fruits. The analysis of long-range sequences identified an LTR retrotransposon in the promoter of the expressed PsMYB10.2 gene that explained the trait in 93.1% of the 145 individuals analyzed. We hypothesize that the LTR retrotransposon may promote the PsMYB10.2 expression and activate the anthocyanin biosynthesis pathway. We propose for the first time the PsMYB10.2 gene as candidate for the flesh color natural variation in Japanese plum and provide a molecular marker for MAS.

Impact of LTR-Retrotransposons on Genome Structure, Evolution, and Function in Curcurbitaceae Species

Long terminal repeat (LTR)-retrotransposons (LTR-RTs) comprise a major portion of many plant genomes and may exert a profound impact on genome structure, function, and evolution. Although many studies have focused on these elements in an individual species, their dynamics on a family level remains elusive. Here, we investigated the abundance, evolutionary dynamics, and impact on associated genes of LTR-RTs in 16 species in an economically important plant family, Cucurbitaceae. Results showed that full-length LTR-RT numbers and LTR-RT content varied greatly among different species, and they were highly correlated with genome size. Most of the full-length LTR-RTs were amplified after the speciation event, reflecting the ongoing rapid evolution of these genomes. LTR-RTs highly contributed to genome size variation via species-specific distinct proliferations. The Angela and Tekay lineages with a greater evolutionary age were amplified in Trichosanthes anguina, whereas a recent activity burst of Reina and another ancient round of Tekay activity burst were examined in Sechium edule. In addition, Tekay and Retand lineages belonging to the Gypsy superfamily underwent a recent burst in Gynostemma pentaphyllum. Detailed investigation of genes with intronic and promoter LTR-RT insertion showed diverse functions, but the term of metabolism was enriched in most species. Further gene expression analysis in G.pentaphyllum revealed that the LTR-RTs within introns suppress the corresponding gene expression, whereas the LTR-RTs within promoters exert a complex influence on the downstream gene expression, with the main function of promoting gene expression. This study provides novel insights into the organization, evolution, and function of LTR-RTs in Cucurbitaceae genomes.

Organization and epigenomic control of RNA polymerase III-transcribed genes in plants

The genetic information linearly scripted in chromosomes is wrapped in a ribonucleoprotein complex called chromatin. The adaptation of its compaction level and spatiotemporal organization refines gene expression in response to developmental and environmental cues. RNA polymerase III (RNAPIII) is responsible for the biogenesis of elementary non-coding RNAs. Their genes are subjected to high duplication and mutational rates, and invade nuclear genomes. Their insertion into different epigenomic environments raises the question of how chromatin packing affects their individual transcription. In this review, we provide a unique perspective to this issue in plants. In addition, we discuss how the genomic organization of RNAPIII-transcribed loci, combined with epigenetic differences, might participate to plant trait variations.

Alternative splicing triggered by the insertion of a CACTA transposon attenuates LsGLK and leads to the development of pale-green leaves in lettuce

Lettuce (Lactuca sativa) is one of the most important vegetable crops in the world. As a leafy vegetable, the polymorphism of lettuce leaves from dark to pale green is an important trait. However, the genetic and molecular mechanisms underlying such variations remain poorly understood. In this study, one major locus controlling the polymorphism of dark- and pale-green leaves in lettuce was identified using genome-wide association studies (GWAS). This locus was then fine mapped to an interval of 5375 bp on chromosome 4 using a segregating population containing 2480 progeny. Only one gene, homologous to the GLK genes in Arabidopsis and other plants, is present in the candidate region. A complementation test confirmed that the candidate gene, LsGLK, contributes to the variation of dark- and pale-green leaves. Sequence analysis showed that a CACTA transposon of 7434 bp was inserted 10 bp downstream of the stop codon of LsGLK, followed by a duplication of a 1826-bp fragment covering exons 3-6 of the LsGLK gene. The transposon insertion did not change the expression level of the LsGLK gene. However, because of alternative splicing, only 6% of the transcripts produced from the transposon insertion were wild-type transcripts, which led to the production of pale-green leaves. An evolutionary analysis revealed that the insertion of the CACTA transposon occurred in cultivated lettuce and might have been selected in particular cultivars to satisfy the diverse demands of consumers. In this study, we demonstrated that a transposon insertion near a gene may affect its splicing and consequently generate phenotypic variations.

Plant transfer RNA-derived fragments: Biogenesis and functions

Processing of mature transfer RNAs (tRNAs) produces complex populations of tRNA-derived fragments (tRFs). Emerging evidence shows that tRFs have important functions in bacteria, animals, and plants. Here, we review recent advances in understanding plant tRFs, focusing on their biological and cellular functions, such as regulating stress responses, mediating plant-pathogen interactions, and modulating post-transcriptional gene silencing and translation. We also review sequencing strategies and bioinformatics resources for studying tRFs in plants. Finally, we discuss future directions for plant tRF research, which will expand our knowledge of plant non-coding RNAs.

Impact of Repetitive DNA Elements on Snake Genome Biology and Evolution

The distinctive biology and unique evolutionary features of snakes make them fascinating model systems to elucidate how genomes evolve and how variation at the genomic level is interlinked with phenotypic-level evolution. Similar to other eukaryotic genomes, large proportions of snake genomes contain repetitive DNA, including transposable elements (TEs) and satellite repeats. The importance of repetitive DNA and its structural and functional role in the snake genome, remain unclear. This review highlights the major types of repeats and their proportions in snake genomes, reflecting the high diversity and composition of snake repeats. We present snakes as an emerging and important model system for the study of repetitive DNA under the impact of sex and microchromosome evolution. We assemble evidence to show that certain repetitive elements in snakes are transcriptionally active and demonstrate highly dynamic lineage-specific patterns as repeat sequences. We hypothesize that particular TEs can trigger different genomic mechanisms that might contribute to driving adaptive evolution in snakes. Finally, we review emerging approaches that may be used to study the expression of repetitive elements in complex genomes, such as snakes. The specific aspects presented here will stimulate further discussion on the role of genomic repeats in shaping snake evolution.

Androgenic-Induced Transposable Elements Dependent Sequence Variation in Barley

A plant genome usually encompasses different families of transposable elements (TEs) that may constitute up to 85% of nuclear DNA. Under stressful conditions, some of them may activate, leading to sequence variation. In vitro plant regeneration may induce either phenotypic or genetic and epigenetic changes. While DNA methylation alternations might be related, i.e., to the Yang cycle problems, DNA pattern changes, especially DNA demethylation, may activate TEs that could result in point mutations in DNA sequence changes. Thus, TEs have the highest input into sequence variation (SV). A set of barley regenerants were derived via in vitro anther culture. High Performance Liquid Chromatography (RP-HPLC), used to study the global DNA methylation of donor plants and their regenerants, showed that the level of DNA methylation increased in regenerants by 1.45% compared to the donors. The Methyl-Sensitive Transposon Display (MSTD) based on methylation-sensitive Amplified Fragment Length Polymorphism (metAFLP) approach demonstrated that, depending on the selected elements belonging to the TEs family analyzed, varying levels of sequence variation were evaluated. DNA sequence contexts may have a different impact on SV generated by distinct mobile elements belonged to various TE families. Based on the presented study, some of the selected mobile elements contribute differently to TE-related SV. The surrounding context of the TEs DNA sequence is possibly important here, and the study explained some part of SV related to those contexts.

Identifying transposable element expression dynamics and heterogeneity during development at the single-cell level with a processing pipeline scTE

Transposable elements (TEs) make up a majority of a typical eukaryote’s genome, and contribute to cell heterogeneity in unclear ways. Single-cell sequencing technologies are powerful tools to explore cells, however analysis is typically gene-centric and TE expression has not been addressed. Here, we develop a single-cell TE processing pipeline, scTE, and report the expression of TEs in single cells in a range of biological contexts. Specific TE types are expressed in subpopulations of embryonic stem cells and are dynamically regulated during pluripotency reprogramming, differentiation, and embryogenesis. Unexpectedly, TEs are expressed in somatic cells, including human disease-specific TEs that are undetectable in bulk analyses. Finally, we apply scTE to single-cell ATAC-seq data, and demonstrate that scTE can discriminate cell type using chromatin accessibly of TEs alone. Overall, our results classify the dynamic patterns of TEs in single cells and their contributions to cell heterogeneity.

How transposable elements (TE) contribute to cell fate changes is unclear. Here, the authors generate a pipeline to quantify TE expression from single cell data. They show the dynamic expression of TEs from gastrulation to somatic cell reprogramming and human disease

Dynamics of DNA Methylation and Its Functions in Plant Growth and Development

Epigenetic modifications in DNA bases and histone proteins play important roles in the regulation of gene expression and genome stability. Chemical modification of DNA base (e.g., addition of a methyl group at the fifth carbon of cytosine residue) switches on/off the gene expression during developmental process and environmental stresses. The dynamics of DNA base methylation depends mainly on the activities of the writer/eraser guided by non-coding RNA (ncRNA) and regulated by the developmental/environmental cues. De novo DNA methylation and active demethylation activities control the methylation level and regulate the gene expression. Identification of ncRNA involved in de novo DNA methylation, increased DNA methylation proteins guiding DNA demethylase, and methylation monitoring sequence that helps maintaining a balance between DNA methylation and demethylation is the recent developments that may resolve some of the enigmas. Such discoveries provide a better understanding of the dynamics/functions of DNA base methylation and epigenetic regulation of growth, development, and stress tolerance in crop plants. Identification of epigenetic pathways in animals, their existence/orthologs in plants, and functional validation might improve future strategies for epigenome editing toward climate-resilient, sustainable agriculture in this era of global climate change. The present review discusses the dynamics of DNA methylation (cytosine/adenine) in plants, its functions in regulating gene expression under abiotic/biotic stresses, developmental processes, and genome stability.

Structural and Functional Annotation of Transposable Elements Revealed a Potential Regulation of Genes Involved in Rubber Biosynthesis by TE-Derived siRNA Interference in Hevea brasiliensis

The natural rubber biosynthetic pathway is well described in Hevea, although the final stages of rubber elongation are still poorly understood. Small Rubber Particle Proteins and Rubber Elongation Factors (SRPPs and REFs) are proteins with major function in rubber particle formation and stabilization. Their corresponding genes are clustered on a scaffold1222 of the reference genomic sequence of the Hevea brasiliensis genome. Apart from gene expression by transcriptomic analyses, to date, no deep analyses have been carried out for the genomic environment of SRPPs and REFs loci. By integrative analyses on transposable element annotation, small RNAs production and gene expression, we analysed their role in the control of the transcription of rubber biosynthetic genes. The first in-depth annotation of TEs (Transposable Elements) and their capacity to produce TE-derived siRNAs (small interfering RNAs) is presented, only possible in the Hevea brasiliensis clone PB 260 for which all data are available. We observed that 11% of genes are located near TEs and their presence may interfere in their transcription at both genetic and epigenetic level. We hypothesized that the genomic environment of rubber biosynthesis genes has been shaped by TE and TE-derived siRNAs with possible transcriptional interference on their gene expression. We discussed possible functionalization of TEs as enhancers and as donors of alternative transcription start sites in promoter sequences, possibly through the modelling of genetic and epigenetic landscapes.

Complex Evolutionary History of Mboumar, a Mariner Element Widely Represented in Ant Genomes

Mboumar-9 is an active mariner-transposable element previously isolated in the ant Messor bouvieri. In this work, a mariner-like element, Mboumar, isolated from 22 species of ants, is analyzed. These species belong to nine different subfamilies, including Leptanillinae, the most primitive ant subfamily, and Myrmicinae and Formicidae, the most derived ones. Consequently, Mboumar-like elements seem to be well-represented in ant genomes. The phylogenetic tree drawn for mariner elements is highly inconsistent with the phylogeny of host ants, with almost identical elements found in clearly distant species and, on the contrary, more variable elements in closely related species. The inconsistency between the two phylogenetic trees indicates that these transposable elements have evolved independently from the speciation events of the ants that host them. Besides, we found closer genetic relationships among elements than among their host ants. We also found potential coding copies with an uninterrupted open reading frame of 345 aa in 11 species. The putative transposase codified by them showed a high sequence identity with the active Mboumar-9 transposase. The results of selection tests suggest the intervention of purifying selection in the evolution of these elements. Overall, our study suggests a complex evolutionary history of the Mboumar-like mariner in ants, with important participation of horizontal transfer events. We also suggest that the evolutionary dynamics of Mboumar-like elements can be influenced by the genetic system of their host ants, which are eusocial insects with a haplodiploid genetic system.

Dynamic patterns of circular and linear RNAs in maize hybrid and parental lines

Hybrid vigor, also known as heterosis, has been widely utilized in agronomic production of maize (Zea mays L.) and other crops. However, the molecular mechanisms underlying heterosis are still not fully understood. To provide a more complete understanding of the transcriptomic dynamics associated with heterosis, we collected a comprehensive set of sequence data on linear mRNA transcripts and circular RNAs (circRNAs) from seedling leaves of two widely used maize inbred lines and their F₁ hybrid at the V4 growth stage. We detected over 25,000 expressed genes with more than 1200 circRNAs that showed dramatic and distinct variations in expression level across the three genotypes. Although most linear and circular transcripts exhibited additive expression in the hybrid, the expression of circRNAs was more likely to be nonadditive. Interestingly, the levels of linear transcripts and their corresponding circRNAs from the same loci showed a significant relationship and coordinated expression mode across all three genotypes. Notably, in the hybrid, allele-specific expression of linear transcripts was significantly associated with the expression of circRNAs from the same locus, suggesting potential regulatory cross talk between linear and circular transcripts. Our study provides a deeper understanding of dynamic variations for both the linear and circular transcriptome in a classical hybrid triplet of maize.

The Past, Present, and Future of Maize Improvement: Domestication, Genomics, and Functional Genomic Routes toward Crop Enhancement

After being domesticated from teosinte, cultivated maize (Zea mays ssp. mays) spread worldwide and now is one of the most important staple crops. Due to its tremendous phenotypic and genotypic diversity, maize also becomes to be one of the most widely used model plant species for fundamental research, with many important discoveries reported by maize researchers. Here, we provide an overview of the history of maize domestication and key genes controlling major domestication-related traits, review the currently available resources for functional genomics studies in maize, and discuss the functions of most of the maize genes that have been positionally cloned and can be used for crop improvement. Finally, we provide some perspectives on future directions regarding functional genomics research and the breeding of maize and other crops.

Influence of CNV on transcript levels of HvCBF genes at Fr-H2 locus revealed by resequencing in resistant barley cv. 'Nure' and expression analysis

Resequencing in resistant cultivar 'Nure' and structural comparison with the same region of susceptible 'Morex' was performed in order to gain a better insight into barley Frost-resistance-H2 locus. Accurate annotation showed copy number variation (CNV) in the proximal part of the locus. In 'Nure', two exact copies of the HvCBF4-HvCBF2A region and one of the HvCBF4-HvCBF2B segment were observed, while in 'Morex' the corresponding region harboured a single HvCBF4-HvCBF2A (22 kb) segment. Abundance and diversity of repetitive element classes, gene function gain/losses, regulatory motifs and SNPs in gene sequences were identified. An expression study of key HvCBFs with/without CNV on selected genotypes contrasting for frost resistance and estimated HvCBF4-HvCBF2B copy number (2-10 copies) was also performed. Under light stimulus at warm temperature (23 °C), CNV of HvCBF2A and HvCBF4 correlated with their expression levels and reported frost resistance of genotypes; moreover, expression levels of HvCBF2A and HvCBF14 were strongly correlated (r = 0.908, p < 0.01). On the other hand, frost resistance correlated to HvCBF14 expression (r = 0.871, p < 0.01) only after cold induction (6°C) in the dark. A complex interplay of HvCBFs expression levels under different light/temperature stimuli is discussed in light of CNV and presence/number of regulatory elements that integrate different signal transduction pathways.

Positional differences of intronic transposons in pAMT affect the pungency level in chili pepper through altered splicing efficiency

Capsaicinoids are unique compounds that give chili pepper fruits their pungent taste. Capsaicinoid levels vary widely among pungent cultivars, which range from low pungency to extremely pungent. However, the molecular mechanisms underlying this quantitative variation have not been elucidated. Our previous study identified various loss-of-function alleles of the pAMT gene which led to low pungency. The mutations in these alleles are commonly defined by Tcc transposon insertion and its footprint. In this study, we identified two leaky pamt alleles (pamt^L1 and pamt^L2 ) with different levels of putative aminotransferase (pAMT) activity. Notably, both alleles had a Tcc transposon insertion in intron 3, but the locations of the insertions within the intron were different. Genetic analysis revealed that pamt^L1 , pamt^L2 and a loss-of-function pamt allele reduced capsaicinoid levels to about 50%, 10% and less than 1%, respectively. pamt^L1 and pamt^L2 encoded functional pAMT proteins, but they exhibited lower transcript levels than the functional type. RNA sequencing analysis showed that intronic transposons disrupted splicing in intron 3, which resulted in simultaneous expression of functional pAMT mRNA and non-functional splice variants containing partial sequences of Tcc. The non-functional splice variants were more dominant in pamt^L2 than in pamt^L1 . This suggested that the difference in position of the intronic transposons could alter splicing efficiency, leading to different pAMT activities and reducing capsaicinoid content to different levels. Our results provide a striking example of allelic variations caused by intronic transposons; these variations contribute to quantitative differences in secondary metabolite contents.

Circular RNAs exhibit extensive intraspecific variation in maize

Comprehensive transcriptome profiling uncovers extensive intraspecific variation of circular RNAs in maize, shedding light on genomic and phenotypic variation among maize inbred lines. Circular RNAs (circRNAs) are single-strand, covalently closed transcripts. A substantial number of circRNAs have been identified and shown to be associated with phenotypic variation in various species. However, little is known about the intraspecific variation of circRNAs in maize (Zea mays L.). Here, we collected a large transcriptomic dataset (by circRNA-seq and mRNA-seq) from seedling leaves of the reference maize inbred lines B73 and Mo17. We identified over 1500 circRNAs in these lines using two circRNA detection methods, CIRCexplorer2 and CIRI. Notably, a substantial proportion of circRNAs varied in terms of sequence or expression level between lines, pointing to extensive intraspecific variation of circRNAs in maize. GO and KEGG analyses showed that genes producing circRNAs with intraspecific variation were more likely to be enriched in multiple functional groups, compared with those that did not produce circRNAs. These findings suggest that circRNAs could be utilized as an indicator of genomic and phenotypic variation among maize inbred lines. Ribosomal profiling revealed that several circRNAs might have translational capacity in maize. These results uncover the extensive intraspecific variation of circRNAs and pave the way for further understanding the molecular mechanisms underlying phenotypic variation at the circRNA level in maize.

Structural variants in 3000 rice genomes

Investigation of large structural variants (SVs) is a challenging yet important task in understanding trait differences in highly repetitive genomes. Combining different bioinformatic approaches for SV detection, we analyzed whole-genome sequencing data from 3000 rice genomes and identified 63 million individual SV calls that grouped into 1.5 million allelic variants. We found enrichment of long SVs in promoters and an excess of shorter variants in 5′ UTRs. Across the rice genomes, we identified regions of high SV frequency enriched in stress response genes. We demonstrated how SVs may help in finding causative variants in genome-wide association analysis. These new insights into rice genome biology are valuable for understanding the effects SVs have on gene function, with the prospect of identifying novel agronomically important alleles that can be utilized to improve cultivated rice.

A high-quality apple genome assembly reveals the association of a retrotransposon and red fruit colour

A complete and accurate genome sequence provides a fundamental tool for functional genomics and DNA-informed breeding. Here, we assemble a high-quality genome (contig N50 of 6.99 Mb) of the apple anther-derived homozygous line HFTH1, including 22 telomere sequences, using a combination of PacBio single-molecule real-time (SMRT) sequencing, chromosome conformation capture (Hi-C) sequencing, and optical mapping. In comparison to the Golden Delicious reference genome, we identify 18,047 deletions, 12,101 insertions and 14 large inversions. We reveal that these extensive genomic variations are largely attributable to activity of transposable elements. Interestingly, we find that a long terminal repeat (LTR) retrotransposon insertion upstream of MdMYB1, a core transcriptional activator of anthocyanin biosynthesis, is associated with red-skinned phenotype. This finding provides insights into the molecular mechanisms underlying red fruit coloration, and highlights the utility of this high-quality genome assembly in deciphering agriculturally important trait in apple.

Transposable elements drive rapid phenotypic variation in Capsella rubella

The mechanisms underlying rapid adaptation to changing environments in species with reduced genetic variation, referred to as the “genetic paradox of invasion,” are unknown. We report that transposable elements (TEs) are highly enriched in the gene promoter regions of Capsella rubella compared with its outcrossing sister species Capsella grandiflora. We also show that a number of polymorphic TEs in C. rubella are associated with changes in gene expression. Frequent TE insertions at FLOWERING LOCUS C of C. rubella affect flowering-time variation, an important life history trait correlated with fitness. These results indicate that TE insertions drive rapid phenotypic variation, which could potentially help adapting to novel environments in species with limited genetic variation.

Rapid phenotypic changes in traits of adaptive significance are crucial for organisms to thrive in changing environments. How such phenotypic variation is achieved rapidly, despite limited genetic variation in species that experience a genetic bottleneck is unknown. Capsella rubella, an annual and inbreeding forb (Brassicaceae), is a great system for studying this basic question. Its distribution is wider than those of its congeneric species, despite an extreme genetic bottleneck event that severely diminished its genetic variation. Here, we demonstrate that transposable elements (TEs) are an important source of genetic variation that could account for its high phenotypic diversity. TEs are (i) highly enriched in C. rubella compared with its outcrossing sister species Capsella grandiflora, and (ii) 4.2% of polymorphic TEs in C. rubella are associated with variation in the expression levels of their adjacent genes. Furthermore, we show that frequent TE insertions at FLOWERING LOCUS C (FLC) in natural populations of C. rubella could explain 12.5% of the natural variation in flowering time, a key life history trait correlated with fitness and adaptation. In particular, we show that a recent TE insertion at the 3′ UTR of FLC affects mRNA stability, which results in reducing its steady-state expression levels, to promote the onset of flowering. Our results highlight that TE insertions can drive rapid phenotypic variation, which could potentially help with adaptation to changing environments in a species with limited standing genetic variation.

Comparative effect of allopolyploidy on transposable element composition and gene expression between Gossypium hirsutum and its two diploid progenitors

An allopolyploidization event formed allotetraploid Gossypium species from an A-genome diploid species and a D-genome diploid species. To explore the responses of transposable elements (TEs) to allopolyploidy, we assembled parallel TE datasets from G. hirsutum, G. arboreum and G. raimondii and analyzed the TE types and the effects of TEs on orthologous gene expression in the three Gossypium genomes. Gypsy was the most abundant TE type and most TEs were located ∼500 bp from genes in all three genomes. In G. hirsutum, 35.6% of genes harbored TE insertions, whereas insertions were more frequent in G. arboreum and G. raimondii. G. hirsutum had the highest proportion of uniquely matching 24-nt small interfering RNAs (siRNAs) that targeted TEs. TEs, particularly those targeted by 24-nt siRNAs, were associated with reduced gene expression, but the effect of TEs on orthologous gene expression varied substantially among species. Orthologous gene expression levels in G. hirsutum were intermediate between those of G. arboreum and G. raimondii, which did not experience TE expansion or reduction resulting from allopolyploidization. This study underscores the diversity of TEs co-opted by host genes and provides insights into the roles of TEs in regulating gene expression in Gossypium.

Partial maintenance of organ-specific epigenetic marks during plant asexual reproduction leads to heritable phenotypic variation

While clonally propagated individuals should share identical genomes, there is often substantial phenotypic variation among them. Both genetic and epigenetic modifications induced during regeneration have been associated with this phenomenon. Here we investigated the fate of the epigenome after asexual propagation by generating clonal individuals from differentiated somatic cells through the manipulation of a zygotic transcription factor. We found that phenotypic novelty in clonal progeny was linked to epigenetic imprints that reflect the organ used for regeneration. Some of these organ-specific imprints can be maintained during the cloning process and subsequent rounds of meiosis. Our findings are fundamental for understanding the significance of epigenetic variability arising from asexual reproduction and have significant implications for future biotechnological applications.

Plants differ from animals in their capability to easily regenerate fertile adult individuals from terminally differentiated cells. This unique developmental plasticity is commonly observed in nature, where many species can reproduce asexually through the ectopic initiation of organogenic or embryogenic developmental programs. While organ-specific epigenetic marks are not passed on during sexual reproduction, the fate of epigenetic marks during asexual reproduction and the implications for clonal progeny remain unclear. Here we report that organ-specific epigenetic imprints in Arabidopsis thaliana can be partially maintained during asexual propagation from somatic cells in which a zygotic program is artificially induced. The altered marks are inherited even over multiple rounds of sexual reproduction, becoming fixed in hybrids and resulting in heritable molecular and physiological phenotypes that depend on the identity of the founder tissue. Consequently, clonal plants display distinct interactions with beneficial and pathogenic microorganisms. Our results demonstrate how novel phenotypic variation in plants can be unlocked through altered inheritance of epigenetic marks upon asexual propagation.

Transposable Elements Contribute to the Adaptation of Arabidopsis thaliana

Transposable elements (TEs) are mobile genetic elements with very high mutation rates that play important roles in shaping genome architecture and regulating phenotypic variation. However, the extent to which TEs influence the adaptation of organisms in their natural habitats is largely unknown. Here, we scanned 201 representative resequenced genomes from the model plant Arabidopsis thaliana and identified 2,311 polymorphic TEs from noncentromeric regions. We found expansion and contraction of different types of TEs in different A. thaliana populations. More importantly, we identified two TE insertions that are likely candidates to play a role in adaptive evolution. Our results highlight the importance of variations in TEs for the adaptation of plants in general in the context of rapid global climate change.

De novo assembly of a Chinese soybean genome

Soybean was domesticated in China and has become one of the most important oilseed crops. Due to bottlenecks in their introduction and dissemination, soybeans from different geographic areas exhibit extensive genetic diversity. Asia is the largest soybean market; therefore, a high-quality soybean reference genome from this area is critical for soybean research and breeding. Here, we report the de novo assembly and sequence analysis of a Chinese soybean genome for "Zhonghuang 13" by a combination of SMRT, Hi-C and optical mapping data. The assembled genome size is 1.025 Gb with a contig N50 of 3.46 Mb and a scaffold N50 of 51.87 Mb. Comparisons between this genome and the previously reported reference genome (cv. Williams 82) uncovered more than 250,000 structure variations. A total of 52,051 protein coding genes and 36,429 transposable elements were annotated for this genome, and a gene co-expression network including 39,967 genes was also established. This high quality Chinese soybean genome and its sequence analysis will provide valuable information for soybean improvement in the future.

Role of the Mycobacterium marinum ESX-1 Secretion System in Sliding Motility and Biofilm Formation

Mycobacterium marinum is a close relative of Mycobacterium tuberculosis that can cause systemic tuberculosis-like infections in ectotherms and skin infections in humans. Sliding motility correlates with biofilm formation and virulence in most bacteria. In this study, we used a sliding motility assay to screen 2,304 transposon mutants of M. marinum NTUH-M6885 and identified five transposon mutants with decreased sliding motility. Transposons that interrupted the type VII secretion system (T7SS) ESX-1-related genes, espE (mmar_5439), espF (mmar_5440), and eccA1 (mmar_5443), were present in 3 mutants. We performed reverse-transcription polymerase chain reaction to verify genes from mmar_5438 to mmar_5450, which were found to belong to a single transcriptional unit. Deletion mutants of espE, espF, espG (mmar_5441), and espH (mmar_5442) displayed significant attenuation regarding sliding motility and biofilm formation. M. marinum NTUH-M6885 possesses a functional ESX-1 secretion system. However, deletion of espG or espH resulted in slightly decreased secretion of EsxB (which is also known as CFP-10). Thus, the M. marinum ESX-1 secretion system mediates sliding motility and is crucial for biofilm formation. These data provide new insight into M. marinum biofilm formation.

Circular RNAs mediated by transposons are associated with transcriptomic and phenotypic variation in maize

Circular RNAs (circRNAs) are covalently closed RNA molecules. Recent studies have shown that circRNAs can arise from the transcripts of transposons. Given the prevalence of transposons in the maize genome and dramatic genomic variation driven by transposons, we hypothesize that transposons in maize may be involved in the formation of circRNAs and further modulate phenotypic variation. We performed circRNA-Seq on B73 seedling leaves and uncovered 2804 high-confidence maize circRNAs, which show distinct genomic features. Comprehensive analyses demonstrated that sequences related to LINE1-like elements (LLEs) and their Reverse Complementary Pairs (LLERCPs) are significantly enriched in the flanking regions of circRNAs. Interestingly, as the number of LLERCPs increase, the accumulation of circRNAs varies, whereas that of linear transcripts decreases. Furthermore, genes with LLERCP-mediated circRNAs are enriched among loci that are associated with phenotypic variation. These results suggest that circRNAs are likely to be involved in the modulation of phenotypic variation by LLERCPs. Further, we showed that the presence/absence variation of LLERCPs was associated with expression variation of circRNA-circ1690 and was related to ear height, potentially through the interplay between circRNAs and functional linear transcripts. Our first study of maize circRNAs uncovers a potential new way for transposons to modulate transcriptomic and phenotypic variations.

Detection of active transposable elements in Arabidopsis thaliana using Oxford Nanopore Sequencing technology

BACKGROUND

Transposables elements (TEs) contribute to both structural and functional dynamics of most eukaryotic genomes. Because of their propensity to densely populate plant and animal genomes, the precise estimation of the impact of transposition on genomic diversity has been considered as one of the main challenges of today's genomics. The recent development of NGS (next generation sequencing) technologies has open new perspectives in population genomics by providing new methods for high throughput detection of Transposable Elements-associated Structural Variants (TEASV). However, these have relied on Illumina platform that generates short reads (up to 350 nucleotides). This limitation in size of sequence reads can cause high false discovery rate (FDR) and therefore limit the power of detection of TEASVs, especially in the case of large, complex genomes. The newest sequencing technologies, such as Oxford Nanopore Technologies (ONT) can generate kilobases-long reads thus representing a promising tool for TEASV detection in plant and animals.

RESULTS

We present the results of a pilot experiment for TEASV detection on the model plant species Arabidopsis thaliana using ONT sequencing and show that it can be used efficiently to detect TE movements. We generated a ~0.8X genome coverage of a met1-derived epigenetic recombinant inbred line (epiRIL) using a MinIon device with R7 chemistry. We were able to detect nine new copies of the LTR-retrotransposon Evadé (EVD). We also evidenced the activity of the DNA transposon CACTA, CAC1.

CONCLUSIONS

Even at a low sequence coverage (0.8X), ONT sequencing allowed us to reliably detect several TE insertions in Arabidopsis thaliana genome. The long read length allowed a precise and un-ambiguous mapping of the structural variations caused by the activity of TEs. This suggests that the trade-off between read length and genome coverage for TEASV detection may be in favor of the former. Should the technology be further improved both in terms of lower error rate and operation costs, it could be efficiently used in diversity studies at population level.

Bioinformatics in the plant genomic and phenomic domain: The German contribution to resources, services and perspectives

Plant genetic resources are a substantial opportunity for plant breeding, preservation and maintenance of biological diversity. As part of the German Network for Bioinformatics Infrastructure (de.NBI) the German Crop BioGreenformatics Network (GCBN) focuses mainly on crop plants and provides both data and software infrastructure which are tailored to the needs of the plant research community. Our mission and key objectives include: (1) provision of transparent access to germplasm seeds, (2) the delivery of improved workflows for plant gene annotation, and (3) implementation of bioinformatics services that link genotypes and phenotypes. This review introduces the GCBN's spectrum of web-services and integrated data resources that address common research problems in the plant genomics community.

A New Reference Genome Assembly for the Microcrustacean Daphnia pulex

Comparing genomes of closely related genotypes from populations with distinct demographic histories can help reveal the impact of effective population size on genome evolution. For this purpose, we present a high quality genome assembly of Daphnia pulex (PA42), and compare this with the first sequenced genome of this species (TCO), which was derived from an isolate from a population with >90% reduction in nucleotide diversity. PA42 has numerous similarities to TCO at the gene level, with an average amino acid sequence identity of 98.8 and >60% of orthologous proteins identical. Nonetheless, there is a highly elevated number of genes in the TCO genome annotation, with ∼7000 excess genes appearing to be false positives. This view is supported by the high GC content, lack of introns, and short length of these suspicious gene annotations. Consistent with the view that reduced effective population size can facilitate the accumulation of slightly deleterious genomic features, we observe more proliferation of transposable elements (TEs) and a higher frequency of gained introns in the TCO genome.

Transposon-mediated epigenetic regulation contributes to phenotypic diversity and environmental adaptation in rice

Transposable elements (TEs) have long been regarded as 'selfish DNA', and are generally silenced by epigenetic mechanisms. However, work in the past decade has identified positive roles for TEs in generating genomic novelty and diversity in plants. In particular, recent studies suggested that TE-induced epigenetic alterations and modification of gene expression contribute to phenotypic variation and adaptation to geography or stress. These findings have led many to regard TEs, not as junk DNA, but as sources of control elements and genomic diversity. As a staple food crop and model system for genomic research on monocot plants, rice (Oryza sativa) has a modest-sized genome that harbors massive numbers of DNA transposons (class II transposable elements) scattered across the genome, which may make TE regulation of genes more prevalent. In this review, we summarize recent progress in research on the functions of rice TEs in modulating gene expression and creating new genes. We also examine the contributions of TEs to phenotypic diversity and adaptation to environmental conditions.

An epiallele of rice AK1 affects photosynthetic capacity

Epigenetic gene variants, termed epialleles, can broaden genetic and phenotypic diversity in eukaryotes. Here, we identify a natural epiallele of OsAK1, which encodes a rice adenylate kinase. The Epi-ak1 plants show albino in young leaf and panicle with abnormal chloroplast structures. We found that no nucleotide sequence variation but hypermethylation at promoter region caused silencing of OsAK1 (Os08g01770) in Epi-ak1 plants. OsAK1 localizes to chloroplast and many genes associated with photosynthesis processes were downregulated in Epi-ak1. Thus, the work identified a novel rice epiallele caused by DNA methylation changes, shedding light on significant roles of DNA methylation on rice development.

Epigenetic regulation and epigenomic landscape in rice

Epigenetic regulation has been implicated in the control of complex agronomic traits in rice (Oryza sativa), a staple food crop and model monocot plant. Recent advances in high-throughput sequencing and the moderately complex genome of rice have made it possible to study epigenetic regulation in rice on a genome-wide scale. This review discusses recent advances in our understanding of epigenetic regulation in rice, with an emphasis on the roles of key epigenetic regulators, the epigenomic landscape, epigenetic variation, transposon repression, and plant development.

Recent advances in the research for the homolog of breast cancer associated gene AtROW1 in higher plants

BARD1 (BRCA1 associated RING domain protein 1), as an important animal tumor suppressor gene associated with many kinds of cancers, has been intensively studied for decades. Surprisingly, homolog of BARD1 was found in plants and it was renamed AtROW1 (repressor of Wuschel-1) according to its extremely important function with regard to plant stem cell homeostasis. Although great advances have been made in human BARD1, the function of this animal tumor-suppressor like gene in plant is not well studied and need to be further elucidated. Here, we review and summarize past and present work regarding this protein. Apart from its previously proposed role in DNA repair, recently it is found essential for shoot and root stem cell development and differentiation in plants. The study of AtROW1 in plant may provide an ideal model for further elucidating the functional mechanism of BARD1 in mammals.

Evolutionary direction of processed pseudogenes

While some pseudogenes have been reported to play important roles in gene regulation, little is known about the possible relationship between pseudogene functions and evolutionary process of pseudogenes, or about the forces responsible for the pseudogene evolution. In this study, we characterized human processed pseudogenes in terms of evolutionary dynamics. Our results show that pseudogenes tend to evolve toward: lower GC content, strong dinucleotide bias, reduced abundance of transcription factor binding motifs and short palindromes, and decreased ability to form nucleosomes. We explored possible evolutionary forces that shaped the evolution pattern of pseudogenes, and concluded that mutations in pseudogenes are likely determined, at least partially, by neighbor-dependent mutational bias and recombination-associated selection.

Structural organization of fatty acid desaturase loci in linseed lines with contrasting linolenic acid contents

Flax (Linum usitatissimum L.), the richest crop source of omega-3 fatty acids (FAs), is a diploid plant with an estimated genome size of ~370 Mb and is well suited for studying genomic organization of agronomically important traits. In this study, 12 bacterial artificial chromosome clones harbouring the six FA desaturase loci sad1, sad2, fad2a, fad2b, fad3a and fad3b from the conventional variety CDC Bethune and the high linolenic acid line M5791 were sequenced, analysed and compared to determine the structural organization of these loci and to gain insights into the genetic mechanisms underlying FA composition in flax. With one gene every 3.2-4.6 kb, the desaturase loci have a higher gene density than the genome's average of one gene per 7.8-8.2 kb. The gene order and orientation across the two genotypes were generally conserved with the exception of the sad1 locus that was predicted to have additional genes in CDC Bethune. High sequence conservation in both genic and intergenic regions of the sad and fad2b loci contrasted with the significant level of variation of the fad2a and fad3 loci, with SNPs being the most frequently observed mutation type. The fad2a locus had 297 SNPs and 36 indels over ~95 kb contrasting with the fad2b locus that had a mere seven SNPs and four indels in ~110 kb. Annotation of the gene-rich loci revealed other genes of known role in lipid or carbohydrate metabolic/catabolic pathways. The organization of the fad2b locus was particularly complex with seven copies of the fad2b gene in both genotypes. The presence of Gypsy, Copia, MITE, Mutator, hAT and other novel repeat elements at the desaturase loci was similar to that of the whole genome. This structural genomic analysis provided some insights into the genomic organization and composition of the main desaturase loci of linseed and of their complex evolution through both tandem and whole genome duplications.