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.

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.

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.

An improved assembly and annotation of the melon (Cucumis melo L.) reference genome.

We report an improved assembly (v3.6.1) of the melon (Cucumis melo L.) genome and a new genome annotation (v4.0). The optical mapping approach allowed correcting the order and the orientation of 21 previous scaffolds and permitted to correctly define the gap-size extension along the 12 pseudomolecules. A new comprehensive annotation was also built in order to update the previous annotation v3.5.1, released more than six years ago. Using an integrative annotation pipeline, based on exhaustive RNA-Seq collections and ad-hoc transposable element annotation, we identified 29,980 protein-coding loci. Compared to the previous version, the v4.0 annotation improved gene models in terms of completeness of gene structure, UTR regions definition, intron-exon junctions and reduction of fragmented genes. More than 8,000 new genes were identified, one third of them being well supported by RNA-Seq data. To make all the new resources easily exploitable and completely available for the scientific community, a redesigned Melonomics genomic platform was released at http://melonomics.net . The resources produced in this work considerably increase the reliability of the melon genome assembly and resolution of the gene models paving the way for further studies in melon and related species.

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.

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.

The use of massive sequencing to detect differences between immature embryos of MON810 and a comparable non-GM maize variety.

The insect resistant maize YieldGard MON810 was studied to assess the extent to which introduction of a transgene may putatively alter the expression of endogenous genes by comparison of various GM lines vs. their non-transgenic counterparts. To assess the extent to which introduction of a transgene may putatively alter the expression of endogenous genes, GM lines of the insect resistant maize YieldGard MON810 were compared with non-transgenic counterparts. For a more in-depth study, high-throughput deep sequencing together with microarrays were used to compare the transcriptomes of immature embryos of the MON810 variety DKC6575, with a cryIA(b) transgene, and its near-isogenic variety Tietar, grown under controlled environmental conditions. This technique also allows characterisation of the transgenic mRNAs produced. 3'UTR-anchored mRNA-seq produced 1,802,571 sequences from DKC6575 and 1,170,973 from Tietar, which mapped to 14,712 and 14,854 unigenes, respectively. Up to 32 reads from the transgenic embryos matched to the synthetic cry1A(b) sequence, similar to medium-abundant mRNAs. Gene expression analysis, using the R-bioconductor packages EdgeR and DEseq, revealed 140 differentially expressed genes mainly involved in carbohydrate metabolism, protein metabolism and chromatin organisation. Comparison of the expression of 30 selected genes in two additional MON810 and near-isogenic variety pairs showed that most of them were differentially expressed in the three pairs of varieties analysed. Analysis of functional annotation and the precise moment of expression of the differentially expressed genes and physiological data obtained suggest a slight but significant delay in seed and plant maturation of MON810 plants. However, these transcriptomic changes were not associated to undesirable changes in the global phenotype or plant behaviour. Moreover, while most expression changes in MON810 immature embryos were maintained in other transgenic varieties, some gene expression was found to be modulated by the genetic background in which the transgene was introduced through conventional breeding programs.

Dealing with scientific integrity issues: the Spanish experience.

Integrity has been an important matter of concern for the scientific community as it affects the basis of its activities. Most countries having a significant scientific activity have dealt with this problem by different means, including drafting specific legal or soft law regulations and the appointment of stable or ad hoc committees that take care of these questions. This has also been the case in Spain. After the period of transition between dictatorship to a democratic regime, and, particularly, after the entrance in the European Union, scientific activity has increased in the country. As it could be expected, problems of misconduct have appeared and different institutions have been dealing with these matters. One of the best examples is that of Consejo Superior de Investigaciones Cientificas (CSIC), the largest institution devoted to scientific research belonging to the Spanish Government. The experience of the CSIC's Ethics Committee in dealing with conflicts related to scientific practices is discussed here.

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.

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.

Altered lignin biosynthesis improves cellulosic bioethanol production in transgenic maize plants down-regulated for cinnamyl alcohol dehydrogenase.

Cinnamyl alcohol dehydrogenase (CAD) is a key enzyme involved in the last step of monolignol biosynthesis. The effect of CAD down-regulation on lignin production was investigated through a transgenic approach in maize. Transgenic CAD-RNAi plants show a different degree of enzymatic reduction depending on the analyzed tissue and show alterations in cell wall composition. Cell walls of CAD-RNAi stems contain a lignin polymer with a slight reduction in the S-to-G ratio without affecting the total lignin content. In addition, these cell walls accumulate higher levels of cellulose and arabinoxylans. In contrast, cell walls of CAD-RNAi midribs present a reduction in the total lignin content and of cell wall polysaccharides. In vitro degradability assays showed that, although to a different extent, the changes induced by the repression of CAD activity produced midribs and stems more degradable than wild-type plants. CAD-RNAi plants grown in the field presented a wild-type phenotype and produced higher amounts of dry biomass. Cellulosic bioethanol assays revealed that CAD-RNAi biomass produced higher levels of ethanol compared to wild-type, making CAD a good target to improve both the nutritional and energetic values of maize lignocellulosic biomass.

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.

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.

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.

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.