Transposons played a major role in the diversification between the closely related almond and peach genomes: results from the almond genome sequence.

We sequenced the genome of the highly heterozygous almond Prunus dulcis cv. Texas combining short- and long-read sequencing. We obtained a genome assembly totaling 227.6 Mb of the estimated almond genome size of 238 Mb, of which 91% is anchored to eight pseudomolecules corresponding to its haploid chromosome complement, and annotated 27 969 protein-coding genes and 6747 non-coding transcripts. By phylogenomic comparison with the genomes of 16 additional close and distant species we estimated that almond and peach (Prunus persica) diverged around 5.88 million years ago. These two genomes are highly syntenic and show a high degree of sequence conservation (20 nucleotide substitutions per kb). However, they also exhibit a high number of presence/absence variants, many attributable to the movement of transposable elements (TEs). Transposable elements have generated an important number of presence/absence variants between almond and peach, and we show that the recent history of TE movement seems markedly different between them. Transposable elements may also be at the origin of important phenotypic differences between both species, and in particular for the sweet kernel phenotype, a key agronomic and domestication character for almond. Here we show that in sweet almond cultivars, highly methylated TE insertions surround a gene involved in the biosynthesis of amygdalin, whose reduced expression has been correlated with the sweet almond phenotype. Altogether, our results suggest a key role of TEs in the recent history and diversification of almond and its close relative peach.

Variability among Cucurbitaceae species (melon, cucumber and watermelon) in a genomic region containing a cluster of NBS-LRR genes.

BACKGROUND: Cucurbitaceae species contain a significantly lower number of genes coding for proteins with similarity to plant resistance genes belonging to the NBS-LRR family than other plant species of similar genome size. A large proportion of these genes are organized in clusters that appear to be hotspots of variability. The genomes of the Cucurbitaceae species measured until now are intermediate in size (between 350 and 450 Mb) and they apparently have not undergone any genome duplications beside those at the origin of eudicots. The cluster containing the largest number of NBS-LRR genes has previously been analyzed in melon and related species and showed a high degree of interspecific and intraspecific variability. It was of interest to study whether similar behavior occurred in other cluster of the same family of genes. RESULTS: The cluster of NBS-LRR genes located in melon chromosome 9 was analyzed and compared with the syntenic regions in other cucurbit genomes. This is the second cluster in number within this species and it contains nine sequences with a NBS-LRR annotation including two genes, Fom1 and Prv, providing resistance against Fusarium and Ppapaya ring-spot virus (PRSV). The variability within the melon species appears to consist essentially of single nucleotide polymorphisms. Clusters of similar genes are present in the syntenic regions of the two species of Cucurbitaceae that were sequenced, cucumber and watermelon. Most of the genes in the syntenic clusters can be aligned between species and a hypothesis of generation of the cluster is proposed. The number of genes in the watermelon cluster is similar to that in melon while a higher number of genes (12) is present in cucumber, a species with a smaller genome than melon. After comparing genome resequencing data of 115 cucumber varieties, deletion of a group of genes is observed in a group of varieties of Indian origin. CONCLUSIONS: Clusters of genes coding for NBS-LRR proteins in cucurbits appear to have specific variability in different regions of the genome and between different species. This observation is in favour of considering that the adaptation of plant species to changing environments is based upon the variability that may occur at any location in the genome and that has been produced by specific mechanisms of sequence variation acting on plant genomes. This information could be useful both to understand the evolution of species and for plant breeding.

The genome of melon (Cucumis melo L.).

We report the genome sequence of melon, an important horticultural crop worldwide. We assembled 375 Mb of the double-haploid line DHL92, representing 83.3% of the estimated melon genome. We predicted 27,427 protein-coding genes, which we analyzed by reconstructing 22,218 phylogenetic trees, allowing mapping of the orthology and paralogy relationships of sequenced plant genomes. We observed the absence of recent whole-genome duplications in the melon lineage since the ancient eudicot triplication, and our data suggest that transposon amplification may in part explain the increased size of the melon genome compared with the close relative cucumber. A low number of nucleotide-binding site-leucine-rich repeat disease resistance genes were annotated, suggesting the existence of specific defense mechanisms in this species. The DHL92 genome was compared with that of its parental lines allowing the quantification of sequence variability in the species. The use of the genome sequence in future investigations will facilitate the understanding of evolution of cucurbits and the improvement of breeding strategies.

Interspecific and intraspecific gene variability in a 1-Mb region containing the highest density of NBS-LRR genes found in the melon genome.

BACKGROUND: Plant NBS-LRR -resistance genes tend to be found in clusters, which have been shown to be hot spots of genome variability. In melon, half of the 81 predicted NBS-LRR genes group in nine clusters, and a 1 Mb region on linkage group V contains the highest density of R-genes and presence/absence gene polymorphisms found in the melon genome. This region is known to contain the locus of Vat, an agronomically important gene that confers resistance to aphids. However, the presence of duplications makes the sequencing and annotation of R-gene clusters difficult, usually resulting in multi-gapped sequences with higher than average errors. RESULTS: A 1-Mb sequence that contains the largest NBS-LRR gene cluster found in melon was improved using a strategy that combines Illumina paired-end mapping and PCR-based gap closing. Unknown sequence was decreased by 70% while about 3,000 SNPs and small indels were corrected. As a result, the annotations of 18 of a total of 23 NBS-LRR genes found in this region were modified, including additional coding sequences, amino acid changes, correction of splicing boundaries, or fussion of ORFs in common transcription units. A phylogeny analysis of the R-genes and their comparison with syntenic sequences in other cucurbits point to a pattern of local gene amplifications since the diversification of cucurbits from other families, and through speciation within the family. A candidate Vat gene is proposed based on the sequence similarity between a reported Vat gene from a Korean melon cultivar and a sequence fragment previously absent in the unrefined sequence. CONCLUSIONS: A sequence refinement strategy allowed substantial improvement of a 1 Mb fragment of the melon genome and the re-annotation of the largest cluster of NBS-LRR gene homologues found in melon. Analysis of the cluster revealed that resistance genes have been produced by sequence duplication in adjacent genome locations since the divergence of cucurbits from other close families, and through the process of speciation within the family a candidate Vat gene was also identified using sequence previously unavailable, which demonstrates the advantages of genome assembly refinements when analyzing complex regions such as those containing clusters of highly similar genes.

High presence/absence gene variability in defense-related gene clusters of Cucumis melo.

BACKGROUND: Changes in the copy number of DNA sequences are one of the main mechanisms generating genome variability in eukaryotes. These changes are often related to phenotypic effects such as genetic disorders or novel pathogen resistance. The increasing availability of genome sequences through the application of next-generation massive sequencing technologies has allowed the study of genomic polymorphisms at both the interspecific and intraspecific levels, thus helping to understand how species adapt to changing environments through genome variability. RESULTS: Data on gene presence/absence variation (PAV) in melon was obtained by resequencing a cultivated accession and an old-relative melon variety, and using previously obtained resequencing data from three other melon cultivars, among them DHL92, on which the current draft melon genome sequence is based. A total of 1,697 PAV events were detected, involving 4.4% of the predicted melon gene complement. In all, an average 1.5% of genes were absent from each analyzed cultivar as compared to the DHL92 reference genome. The most populated functional category among the 304 PAV genes of known function was that of stress response proteins (30% of all classified PAVs). Our results suggest that genes from multi-copy families are five times more likely to be affected by PAV than singleton genes. Also, the chance of genes present in the genome in tandem arrays being affected by PAV is double that of isolated genes, with PAV genes tending to be in longer clusters. The highest concentration of PAV events detected in the melon genome was found in a 1.1 Mb region of linkage group V, which also shows the highest density of melon stress-response genes. In particular, this region contains the longest continuous gene-containing PAV sequence so far identified in melon. CONCLUSIONS: The first genome-wide report of PAV variation among several melon cultivars is presented here. Multi-copy and clustered genes, especially those with putative stress-response functions, were found to be particularly affected by PAV polymorphisms. As cucurbits are known to possess a significantly lower number of defense-related genes compared to other plant species, PAV variation may play an important role in generating new pathogen resistances at the subspecies level. In addition, these results show the limitations of single reference genome sequences as the only basis for characterization and cloning of resistance genes.

Possible EU futures for CRISPR-edited plants: Little margin for optimism?

This article addresses the scenarios that may be encountered by the first application for pre-market approval of a CRISPR-edited plant in the EU. Two alternative scenarios are considered in the short and medium term. One of these possible EU futures depends on the final drafting and approval of EU legislation on certain New Genomic Techniques, which was started in 2021 and is due to be quite advanced before the next European Parliament elections in 2024. Since the proposed legislation excludes plants with foreign DNA, two different approval processes for CRISPR-edited plants will coexist if the legislation enters into force: one for plants whose genome has been altered, resulting in mutagenesis, cisgenesis and intragenesis; and the second for plants whose alterations result in transgenesis in general. In the event that this legislative process does not succeed, CRISPR-edited plants in the EU could face a regulatory scenario whose foundations were laid in the 1990s: the regulatory framework that applies to GM crops, food and feed. In this review, an ad hoc analytical framework has been built that considers in depth the two possible futures for CRISPR-edited plants in the EU. This framework emphasises the way in which the European Union and the Member States (MS), with their respective national interests, have historically shaped the regulatory framework for plant breeding in the EU. On the basis of the analyses carried out on the two possible futures for CRISPR-edited plants and of their potential with respect to plant breeding, the main conclusions are the following. Firstly, that the regulatory review that started in 2021 is not in itself "good enough" for plant breeding and CRISPR-edited plants. Secondly, that compared to its alternative, the regulatory review currently underway contains at least some promising improvements in the short term. Hence, thirdly, in addition to adopting the current regulation, the MS need to continue to work towards a substantial improvement in the legal status of plant breeding in the EU in the medium term.

ZmMYB31 directly represses maize lignin genes and redirects the phenylpropanoid metabolic flux.

Few regulators of phenylpropanoids have been identified in monocots having potential as biofuel crops. Here we demonstrate the role of the maize (Zea mays) R2R3-MYB factor ZmMYB31 in the control of the phenylpropanoid pathway. We determined its in vitro consensus DNA-binding sequence as ACC(T)/(A) ACC, and chromatin immunoprecipitation (ChIP) established that it interacts with two lignin gene promoters in vivo. To explore the potential of ZmMYB31 as a regulator of phenylpropanoids in other plants, its role in the regulation of the phenylpropanoid pathway was further investigated in Arabidopsis thaliana. ZmMYB31 downregulates several genes involved in the synthesis of monolignols and transgenic plants are dwarf and show a significantly reduced lignin content with unaltered polymer composition. We demonstrate that these changes increase cell wall degradability of the transgenic plants. In addition, ZmMYB31 represses the synthesis of sinapoylmalate, resulting in plants that are more sensitive to UV irradiation, and induces several stress-related proteins. Our results suggest that, as an indirect effect of repression of lignin biosynthesis, transgenic plants redirect carbon flux towards the biosynthesis of anthocyanins. Thus, ZmMYB31 can be considered a good candidate for the manipulation of lignin biosynthesis in biotechnological applications.

Determination of the melon chloroplast and mitochondrial genome sequences reveals that the largest reported mitochondrial genome in plants contains a significant amount of DNA having a nuclear origin.

BACKGROUND: The melon belongs to the Cucurbitaceae family, whose economic importance among vegetable crops is second only to Solanaceae. The melon has a small genome size (454 Mb), which makes it suitable for molecular and genetic studies. Despite similar nuclear and chloroplast genome sizes, cucurbits show great variation when their mitochondrial genomes are compared. The melon possesses the largest plant mitochondrial genome, as much as eight times larger than that of other cucurbits. RESULTS: The nucleotide sequences of the melon chloroplast and mitochondrial genomes were determined. The chloroplast genome (156,017 bp) included 132 genes, with 98 single-copy genes dispersed between the small (SSC) and large (LSC) single-copy regions and 17 duplicated genes in the inverted repeat regions (IRa and IRb). A comparison of the cucumber and melon chloroplast genomes showed differences in only approximately 5% of nucleotides, mainly due to short indels and SNPs. Additionally, 2.74 Mb of mitochondrial sequence, accounting for 95% of the estimated mitochondrial genome size, were assembled into five scaffolds and four additional unscaffolded contigs. An 84% of the mitochondrial genome is contained in a single scaffold. The gene-coding region accounted for 1.7% (45,926 bp) of the total sequence, including 51 protein-coding genes, 4 conserved ORFs, 3 rRNA genes and 24 tRNA genes. Despite the differences observed in the mitochondrial genome sizes of cucurbit species, Citrullus lanatus (379 kb), Cucurbita pepo (983 kb) and Cucumis melo (2,740 kb) share 120 kb of sequence, including the predicted protein-coding regions. Nevertheless, melon contained a high number of repetitive sequences and a high content of DNA of nuclear origin, which represented 42% and 47% of the total sequence, respectively. CONCLUSIONS: Whereas the size and gene organisation of chloroplast genomes are similar among the cucurbit species, mitochondrial genomes show a wide variety of sizes, with a non-conserved structure both in gene number and organisation, as well as in the features of the noncoding DNA. The transfer of nuclear DNA to the melon mitochondrial genome and the high proportion of repetitive DNA appear to explain the size of the largest mitochondrial genome reported so far.

Proteomic analysis of MON810 and comparable non-GM maize varieties grown in agricultural fields.

Worldwide maize is the second major agricultural commodity and around one-fourth is currently biotech, with significant application of the insect resistant event MON810 particularly in the European Union. Grains are the major commercialized part of the plant, and can be harvested after maturity (for food and feed purposes) or at late milky-starchy stage (for forage uses, with the whole plant). We assessed possible proteomic unintended effects of the MON810 transgene using two-dimensional gel electrophoresis coupled to mass spectrometry. To keep in a realistic scenario we used plants grown in agricultural fields in a region where ~50% of maize was MON810, and analyzed grains at milky-starchy stage. In maize, differential transcripts and metabolites between GM and comparable non-GM varieties tend to be variety specific. Thus, we analyzed two variety pairs, DKC6575/Tietar and PR33P67/PR33P66 which are considered representative of Food and Agriculture Organization 700 and 600 varieties commercially grown in the region. MON810 and non-GM milky-starchy grains had virtually identical proteomic patterns, with a very small number of spots showing fold-variations in the 1-1.8 range. They were all variety specific and had divergent identities and functions. Although 2DE allows the analysis of a limited dataset our results support substantial equivalence between MON810 and comparable non-GM varieties.

Generation of a BAC-based physical map of the melon genome.

BACKGROUND: Cucumis melo (melon) belongs to the Cucurbitaceae family, whose economic importance among horticulture crops is second only to Solanaceae. Melon has high intra-specific genetic variation, morphologic diversity and a small genome size (450 Mb), which make this species suitable for a great variety of molecular and genetic studies that can lead to the development of tools for breeding varieties of the species. A number of genetic and genomic resources have already been developed, such as several genetic maps and BAC genomic libraries. These tools are essential for the construction of a physical map, a valuable resource for map-based cloning, comparative genomics and assembly of whole genome sequencing data. However, no physical map of any Cucurbitaceae has yet been developed. A project has recently been started to sequence the complete melon genome following a whole-genome shotgun strategy, which makes use of massive sequencing data. A BAC-based melon physical map will be a useful tool to help assemble and refine the draft genome data that is being produced. RESULTS: A melon physical map was constructed using a 5.7 x BAC library and a genetic map previously developed in our laboratories. High-information-content fingerprinting (HICF) was carried out on 23,040 BAC clones, digesting with five restriction enzymes and SNaPshot labeling, followed by contig assembly with FPC software. The physical map has 1,355 contigs and 441 singletons, with an estimated physical length of 407 Mb (0.9 x coverage of the genome) and the longest contig being 3.2 Mb. The anchoring of 845 BAC clones to 178 genetic markers (100 RFLPs, 76 SNPs and 2 SSRs) also allowed the genetic positioning of 183 physical map contigs/singletons, representing 55 Mb (12%) of the melon genome, to individual chromosomal loci. The melon FPC database is available for download at http://melonomics.upv.es/static/files/public/physical_map/. CONCLUSIONS: Here we report the construction of the first physical map of a Cucurbitaceae species described so far. The physical map was integrated with the genetic map so that a number of physical contigs, representing 12% of the melon genome, could be anchored to known genetic positions. The data presented is already helping to improve the quality of the melon genomic sequence available as a result of a project currently being carried out in Spain, adopting a whole genome shotgun approach based on 454 sequencing data.

Genome-wide BAC-end sequencing of Cucumis melo using two BAC libraries.

BACKGROUND: Although melon (Cucumis melo L.) is an economically important fruit crop, no genome-wide sequence information is openly available at the current time. We therefore sequenced BAC-ends representing a total of 33,024 clones, half of them from a previously described melon BAC library generated with restriction endonucleases and the remainder from a new random-shear BAC library. RESULTS: We generated a total of 47,140 high-quality BAC-end sequences (BES), 91.7% of which were paired-BES. Both libraries were assembled independently and then cross-assembled to obtain a final set of 33,372 non-redundant, high-quality sequences. These were grouped into 6,411 contigs (4.5 Mb) and 26,961 non-assembled BES (14.4 Mb), representing ~4.2% of the melon genome. The sequences were used to screen genomic databases, identifying 7,198 simple sequence repeats (corresponding to one microsatellite every 2.6 kb) and 2,484 additional repeats of which 95.9% represented transposable elements. The sequences were also used to screen expressed sequence tag (EST) databases, revealing 11,372 BES that were homologous to ESTs. This suggests that ~30% of the melon genome consists of coding DNA. We observed regions of microsynteny between melon paired-BES and six other dicotyledonous plant genomes. CONCLUSION: The analysis of nearly 50,000 BES from two complementary genomic libraries covered ~4.2% of the melon genome, providing insight into properties such as microsatellite and transposable element distribution, and the percentage of coding DNA. The observed synteny between melon paired-BES and six other plant genomes showed that useful comparative genomic data can be derived through large scale BAC-end sequencing by anchoring a small proportion of the melon genome to other sequenced genomes.

Stability of the MON 810 transgene in maize.

We analysed the DNA variability of the transgene insert and its flanking regions in maize MON 810 commercial varieties. Southern analysis demonstrates that breeding, since the initial transformation event more than 10 years ago, has not resulted in any rearrangements. A detailed analysis on the DNA variability at the nucleotide level, using DNA mismatch endonuclease assays, showed the lack of polymorphisms in the transgene insert. We conclude that the mutation rate of the transgene is not significantly different from that observed in the maize endogenous genes. Six SNPs were observed in the 5'flanking region, corresponding to a Zeon1 retrotransposon long terminal repeat. All six SNPs are more than 500 bp upstream of the point of insertion of the transgene and do not affect the reliability of the established PCR-based transgene detection and quantification methods. The mutation rate of the flanking region is similar to that expected for a maize repetitive sequence. We detected low levels of cytosine methylation in leaves of different transgenic varieties, with no significant differences on comparing different transgenic varieties, and minor differences in cytosine methylation when comparing leaves at different developmental stages. There was also a reduction in cryIAb mRNA accumulation during leaf development.

Sequencing of 6.7 Mb of the melon genome using a BAC pooling strategy.

BACKGROUND: Cucumis melo (melon) belongs to the Cucurbitaceae family, whose economic importance among horticulture crops is second only to Solanaceae. Melon has a high intra-specific genetic variation, morphologic diversity and a small genome size (454 Mb), which make it suitable for a great variety of molecular and genetic studies. A number of genetic and genomic resources have already been developed, such as several genetic maps, BAC genomic libraries, a BAC-based physical map and EST collections. Sequence information would be invaluable to complete the picture of the melon genomic landscape, furthering our understanding of this species' evolution from its relatives and providing an important genetic tool. However, to this day there is little sequence data available, only a few melon genes and genomic regions are deposited in public databases. The development of massively parallel sequencing methods allows envisaging new strategies to obtain long fragments of genomic sequence at higher speed and lower cost than previous Sanger-based methods. RESULTS: In order to gain insight into the structure of a significant portion of the melon genome we set out to perform massive sequencing of pools of BAC clones. For this, a set of 57 BAC clones from a double haploid line was sequenced in two pools with the 454 system using both shotgun and paired-end approaches. The final assembly consists of an estimated 95% of the actual size of the melon BAC clones, with most likely complete sequences for 50 of the BACs, and a total sequence coverage of 39x. The accuracy of the assembly was assessed by comparing the previously available Sanger sequence of one of the BACs against its 454 sequence, and the polymorphisms found involved only 1.7 differences every 10,000 bp that were localized in 15 homopolymeric regions and two dinucleotide tandem repeats. Overall, the study provides approximately 6.7 Mb or 1.5% of the melon genome. The analysis of this new data has allowed us to gain further insight into characteristics of the melon genome such as gene density, average protein length, or microsatellite and transposon content. The annotation of the BAC sequences revealed a high degree of collinearity and protein sequence identity between melon and its close relative Cucumis sativus (cucumber). Transposon content analysis of the syntenic regions suggests that transposition activity after the split of both cucurbit species has been low in cucumber but very high in melon. CONCLUSIONS: The results presented here show that the strategy followed, which combines shotgun and BAC-end sequencing together with anchored marker information, is an excellent method for sequencing specific genomic regions, especially from relatively compact genomes such as that of melon. However, in agreement with other results, this map-based, BAC approach is confirmed to be an expensive way of sequencing a whole plant genome. Our results also provide a partial description of the melon genome's structure. Namely, our analysis shows that the melon genome is highly collinear with the smaller one of cucumber, the size difference being mainly due to the expansion of intergenic regions and proliferation of transposable elements.

Characterization of polyadenylated cryIA(b) transcripts in maize MON810 commercial varieties.

The Zea mays L. event MON810 is one of the major commercialized genetically modified crops. The inserted expression cassette has a 3' truncation partially affecting the cryIA(b) coding sequence, resulting in the lack of the NOS terminator, with transcription of the transgene reported to read-through 3'-past the truncation site. Here, we demonstrate that the cryIA(b) transgene gives rise to a variety of polyadenylated transcripts of different sizes that extend to around 1 kbp downstream the truncation site. A Stop codon at position +7 downstream the truncation site indicates the production of a transgenic protein with two additional amino acids; which is compatible with the reported size of the CryIA(b) protein in MON810. There is no evidence of the existence of other translated products. Several main 3' transcription termination regions were detected close to the truncation site and in the transgene 3' flanking sequence. Next to these main termination sites, we identified some sequence motifs that could potentially act as 3'-end-processing elements and drive termination of the transgene transcripts. The MON810 transgene has been introduced into different commercial varieties through breeding programs. Here, we demonstrate that there are no significant differences among the levels of transgene mRNA accumulation, major transcript sizes and 3' termini profiles comparing a number of MON810 commercial varieties grown under similar environmental conditions. Commercial varieties of this event appear to be stable in terms of transgene expression.

An oligo-based microarray offers novel transcriptomic approaches for the analysis of pathogen resistance and fruit quality traits in melon (Cucumis melo L.).

BACKGROUND: Melon (Cucumis melo) is a horticultural specie of significant nutritional value, which belongs to the Cucurbitaceae family, whose economic importance is second only to the Solanaceae. Its small genome of approx. 450 Mb coupled to the high genetic diversity has prompted the development of genetic tools in the last decade. However, the unprecedented existence of a transcriptomic approaches in melon, highlight the importance of designing new tools for high-throughput analysis of gene expression. RESULTS: We report the construction of an oligo-based microarray using a total of 17,510 unigenes derived from 33,418 high-quality melon ESTs. This chip is particularly enriched with genes that are expressed in fruit and during interaction with pathogens. Hybridizations for three independent experiments allowed the characterization of global gene expression profiles during fruit ripening, as well as in response to viral and fungal infections in plant cotyledons and roots, respectively. Microarray construction, statistical analyses and validation together with functional-enrichment analysis are presented in this study. CONCLUSION: The platform validation and enrichment analyses shown in our study indicate that this oligo-based microarray is amenable for future genetic and functional genomic studies of a wide range of experimental conditions in melon.

The maize ZmMYB42 represses the phenylpropanoid pathway and affects the cell wall structure, composition and degradability in Arabidopsis thaliana.

The involvement of the maize ZmMYB42 R2R3-MYB factor in the phenylpropanoid pathway and cell wall structure and composition was investigated by overexpression in Arabidopsis thaliana. ZmMYB42 down-regulates several genes of the lignin pathway and this effect reduces the lignin content in all lignified tissues. In addition, ZmMYB42 plants generate a lignin polymer with a decreased S to G ratio through the enrichment in H and G subunits and depletion in S subunits. This transcription factor also regulates other genes involved in the synthesis of sinapate esters and flavonoids. Furthermore, ZmMYB42 affects the cell wall structure and degradability, and its polysaccharide composition. Together, these results suggest that ZmMYB42 may be part of the regulatory network controlling the phenylpropanoid biosynthetic pathway.

Use of omics analytical methods in the study of genetically modified maize varieties tested in 90 days feeding trials.

Genetically modified (GM) maize and their non-modified counterparts were compared using MON810 varieties, the only GMO event cultivated in Europe. The differences in grain samples were analysed by omics profiles, including transcriptomics, proteomics and metabolomics. Other cultivated maize varieties were analysed as a reference for the variability that will exist between cultivated varieties. The observed differences between modified and non-modified maize varieties do not exceed typical differences between non-modified varieties. The use of these advanced analytical approaches to analyse novel plant materials as compared to the results from animal feeding trials with whole foods is assessed. No indications were observed for changes in the GM varieties that warrant further investigations. Furthermore, it was shown that such indications will be obtained if maize samples of inferior quality are analysed similarly. Omics data provide detailed analytical information of the plant material, which facilitates a risk assessment procedure of new (GM) plant varieties.