# Chimeras in Merlot grapevine revealed by phased assembly

V. Sichel, G. Sarah, Nabil Girollet, Valérie Laucou, C. Roux, M. Roques, P. Mournet, L. Le Cunff, P.F. Bert, Patrice This, et al.

► To cite this version:

V. Sichel, G. Sarah, Nabil Girollet, Valérie Laucou, C. Roux, et al. Chimeras in Merlot grapevine revealed by phased assembly. BMC Genomics, 2023, 24 (1), pp.396. 10.1186/s12864-023-09453-8. hal-04187937

HAL Id: hal-04187937

https://hal.inrae.fr/hal-04187937

Submitted on 25 Aug 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

Open Access

<figure><img src="image_2.png" /></figure>

Check for
updates

# Chimeras in Merlot grapevine revealed by phased assembly

V. Sichel¹, G. Sarah¹,⁴, N. Girollet², V. Laucou¹,⁴, C. Roux¹,⁴, M. Roques³,⁴, P. Mournet¹,⁵, L. Le Cunff³,⁴, P.F. Bert², P. This¹,⁴*
and T. Lacombe¹,⁴

**Abstract**

Chimerism is the phenomenon when several genotypes coexist in a single individual. Used to understand plant ontogenesis they also have been valorised through new cultivar breeding. Viticulture has been taking economic advantage out of chimeras when the variant induced an important modification of wine type such as berry skin colour. Crucial agronomic characters may also be impacted by chimeras that aren't identified yet. Periclinal chimera where the variant has entirely colonised a cell layer is the most stable and can be propagated through cuttings. In grapevine, leaves are derived from both meristem layers, L1 and L2. However, lateral roots are formed from the L2 cell layer only. Thus, comparing DNA sequences of roots and leaves allows chimera detection. In this study we used new generation Hifi long reads sequencing, recent bioinformatics tools and trio-binning with parental sequences to detect periclinal chimeras on 'Merlot' grapevine cultivar. Sequencing of cv. 'Magdeleine Noire des Charentes' and 'Cabernet Franc', the parents of cv. 'Merlot', allowed haplotype resolved assembly. Pseudomolecules were built with a total of 33 to 47 contigs and in few occasions a unique contig for one chromosome. This high resolution allowed haplotype comparison. Annotation was transferred from PN40024 VCost.v3 to all pseudomolecules. After strong selection of variants, 51 and 53 'Merlot' specific periclinal chimeras were found on the Merlot-haplotype-CF and Merlot-haplotype-MG respectively, 9 and 7 been located in a coding region. A subset of positions was analysed using Molecular Inversion Probes (MIPseq) and 69% were unambiguously validated, 25% are doubtful because of technological noise or weak depth and 6% invalidated. These results open new perspectives on chimera detection as an important resource to improve cultivars through clonal selection or breeding.

**Keywords** Chimera, Hifi sequencing, Phased assembly, Whole genome, *Vitis vinifera*

**Background**

**Chimera phenomenon**

Individuals including cells with different genotypes are called chimeras or genetic mosaics. They are formed when a somatic genetic variation appears in a single cell in the meristem and is propagated through cell divisions. Occasionally the variation modifies a character of the plant and makes the chimera visible. This phenomenon called sporting was observed many centuries ago and has been fascinating scientists since then [1–3]. Chimeras can also induce variegation which is a very noticeable sporting type because it appears as a mosaic of colours in

<footnote>*Correspondence:
P. This
patrice.this@inrae.fr
¹UMR AGAP Institut, Univ Montpellier, CIRAD, INRAE, Institut Agro,
Montpellier F-34398, France
²EGFV, Université de Bordeaux, Bordeaux-Sciences Agro, INRAe, ISVV, 210
Chemin de Leysotte, F-33882 Villenave d'Ornon, France
³Institut Français de la Vigne et du Vin, Montpellier F-34398, France
⁴UMT Geno-Vigne®, IFV-INRAE-Institut Agro, Montpellier F-34398, France
⁵UMR AGAP Institut, CIRAD, Montpellier F-34398, France</footnote>
either leaves, flowers or fruits [4]. In 1907, Winkler was the first to use the term chimera while observing grafted plants [5]. Then, the investigation of variegation led Erwin Baur's on the path of non-Mendelian inheritance [6]. Later, colchicine treatment on Datura seeds revealed periclinal chimeras [7] which allowed the understanding of cell lineages and ontogenesis of plant organs [8–10]. Indeed, two main types exist: sectorial and periclinal chimeras. The differences between those are the complete (periclinal) or incomplete (sectorial) colonisation of a cell layer by the somatic variation [11–13]. Periclinal chimeras are the most stable and can propagate by vegetative multiplication from cuttings [14]. Chimeras have already highly contributed to plant ontogenesis comprehension [10, 15]. In most cases these mutations are silent although some may modify the plant's phenotype on important agronomic traits [16]. They can be used to study biosynthetic pathways [17] but should also be considered as an important source to improve current cultivars or breed new ones [18, 19].

Viticulture has already been taking economic advan-
tage of grapevine chimeras by propagating the new phe-
notype as a new cultivar. For instance, genetic mosaics
explain the origin and the evolution between cv. ‘Pinot
Blanc’ and ‘Pinot Gris’ with a modification of berry skin
colour [20, 21] and between several teinturier cultivars
[22]. A somatic mutation in cv. ‘Meunier’, derived from
‘Pinot Noir’, was used to produce a microvine which

strongly accelerates physiology, biology and genetics
studies [23].

For these reasons, grapevine is a good example to study
chimeras in woody plant species. While plant organs
generally originate from three meristematic cell layers,
grapevine organs come from only two functional cell
layers (L1 and L2) in the apical meristem [12]. Leaves
derived from L1 and L2 cell layers while other organs like
gametophytic tissues and lateral roots originate from L2
cell layer only [12, 24]. Indeed, lateral roots are formed
from the differentiation of the meristematic L2 cell layer
(Fig. 1). Comparing the whole DNA sequence obtained
through leaf samples (L1+L2) with lateral roots samples
(L2 only) could allow us to identify chimeras that may
not be visible but play an important role in intra-varietal
genetic diversity. In order to do this, a high quality whole
genome sequence to accurately validate a chimera against
sequencing mistakes is necessary. The assembly genome
also needs to be resolved per haplotype to distinguish
chimeras from grapevine heterozygosity.

Whole genome new generation sequencing

The first grapevine whole genome sequencing was pub-
lished in 2007 by the French-Italian consortium [25].
Since then, the ‘PN40024’ sequence obtained with a
nearly homozygous inbred ‘Pinot Noir’ plant has been
the reference genome for Vitis sp. The first version had
8X coverage and has been gradually updated through
the 12x.v0, the 12x.v2 [26] and the most recent one is

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Cellular layers in grapevine roots and leaves. Schematic representation of a grapevine plant. Leaf and lateral root cross sections are enlarged in order to present the different cell layers present in both organs. Leaves are derived from both L1 and L2 meristem cell layers, while lateral roots are only formed out of the L2 layer</figcaption></figure>
PN40024.v4 which has 40X coverage (European Nucleotide Archive - project PRJEB45423). Having a good quality whole genome reference has highly increased the understanding of *Vitis vinifera* genome but the nearly homozygous plant (PN40024 is homozygous on 93% of the genome) can hardly be considered as representative of cultivars used in grape production.

Short reads technology produces accurate genome sequencing, but because of the highly repetitive sequences of the grapevine genome they are difficult to assemble, therefore producing a whole reference genome with short reads can become very challenging [27]. Long reads technology has allowed massive improvement in genome assembly. The Falcon-Unzip phasing algorithm was very successful on Arabidopsis but had more difficulty with the grape cv. 'Cabernet-Sauvignon' sample because of the high rate of heterozygous position and the amount of repetitive sequences [28]. Purging treatment on haplotigs allowed to increase the assembly quality on the 'Chardonnay' sequence [29]. The 'Carmenère' phased assembly was also improved by optimising coverage, error correction, repeat masking methods and assembly parameters of FALCON-UNZIP [30]. Since then resolving haplotype phased assembly has become more

accessible and numerous cultivars of the *Vitis* genus have been sequenced (Table 1).

Third generation long reads sequencing with high accuracy brings us into a new perspective of whole genome sequencing. Bioinformatics engineering is adapting to these new sequencing technologies, long reads haplotype assembly is now possible for diploids through tools like Hi-Canu and Hifiasm [34]. To increase accuracy and have a better chance to resolve haplotype phasing, triobinning with parental sequences can be used to sort the child's reads in two groups [35]. Up to now different techniques have led to chimera detection: (i) random amplified polymorphic DNA [36]; (ii) comparing phenotypes of regenerated plants from different cell layers [14]; (iii) comparing microsatellites markers in wood or roots tissues (L2) against leaves (L1+L2) coming from the same plant [37]; (iv) flow cytometry measurements on pericarp and flesh fruit tissues in order to compare ploidy level between L1 and L2 [38]; (v) Real Time PCR on regenerated transgenic plants in order to evaluate the amount of chimeras and the uniformity of the transformation [39]; (vi) microsatellite (SSR) amplification by PCR when three alleles are found on one loci and confirmed by comparing different regenerated plants [40]; (vii) comparing

<table><caption>**Table 1** Whole grapevine genome sequences published until today:</caption><thead><tr><th>Authors</th><th>Year</th><th>Technology</th><th>Grapevine genotype</th><th>Size (mbp)</th><th>Haplotig size (mbp)</th><th>Coverage</th><th>N50 (kbp)</th><th>BUSCO</th></tr></thead><tbody><tr><td>[25] Jaillon et al.</td><td>2007</td><td>Sanger</td><td>PN 40,024 (12x.v1)</td><td>487</td><td></td><td>8X</td><td>65.9</td><td></td></tr><tr><td>[26] Canaguier et al.</td><td>2017</td><td>Sanger</td><td>PN 40,024 (12x.v2)</td><td></td><td></td><td>12X</td><td></td><td></td></tr><tr><td>[28] Chin et al.</td><td>2017</td><td>PacBio RS II</td><td>Cabernet-Sauvignon</td><td>591</td><td></td><td>140X</td><td>72</td><td>80%</td></tr><tr><td>[29] Roach et al.</td><td>2018</td><td>PacBio RS II</td><td>Chardonnay</td><td>490</td><td></td><td>115X</td><td>935.8</td><td>95%</td></tr><tr><td>[30] Minio et al.</td><td>2019</td><td>PacBio RS II</td><td>Carmenère</td><td>622</td><td></td><td>115X</td><td>1040</td><td>95%</td></tr><tr><td>[31] Girollet et al.</td><td>2019</td><td>PacBio RS II</td><td>Vitis riparia</td><td>500</td><td></td><td>225X</td><td>1000</td><td>95%</td></tr><tr><td rowspan="11">[32] Massonnet et al.</td><td rowspan="11">2020</td><td>PacBio RS II</td><td>Muscadinia rotundifolia</td><td>460</td><td>364</td><td>115X</td><td>4761</td><td>97%</td></tr><tr><td rowspan="10">PacBio Sequel II</td><td>Vitis arizonica b40-14</td><td>604</td><td>337</td><td>160X</td><td>1536</td><td>96%</td></tr><tr><td>Vitis vinifera subsp. sylvestris O34-16</td><td>678</td><td>252</td><td>85X</td><td>998</td><td>97%</td></tr><tr><td>Vitis vinifera subsp. sylvestris DVIT3603.16</td><td>667</td><td>314</td><td>82X</td><td>2661</td><td>97%</td></tr><tr><td>Vitis vinifera subsp. sylvestris DVIT3603.07</td><td>663</td><td>234</td><td>61X</td><td>1169</td><td>97%</td></tr><tr><td>Vitis vinifera subsp. sylvestris DVIT3351.27</td><td>670</td><td>311</td><td>77X</td><td>1778</td><td>97%</td></tr><tr><td>Zinfandel</td><td>591</td><td>306</td><td>94X</td><td>1062</td><td>97%</td></tr><tr><td>Merlot</td><td>606</td><td>244</td><td>64X</td><td>810</td><td>94%</td></tr><tr><td>Black Corinth 2.1</td><td>672</td><td>288</td><td>63X</td><td>1113</td><td>97%</td></tr><tr><td>Black Corinth seeded</td><td>650</td><td>357</td><td>75X</td><td>2309</td><td>97%</td></tr><tr><td>Cabernet-Sauvignon</td><td>449</td><td>444</td><td>147X</td><td>2161</td><td>97%</td></tr><tr><td rowspan="2">[33] Zou et al.</td><td rowspan="2">2021</td><td>PacBio Sequel I</td><td>Carmenère</td><td>623</td><td>420</td><td>112x</td><td>1039</td><td>97%</td></tr><tr><td>PacBio Sequel I</td><td>Riesling</td><td>742</td><td>323</td><td>118x</td><td>2970</td><td>98%</td></tr></tbody></table>

Each references, are sorted by authors and publication year. The sequencing technology used is specified as well as the genotype's name. The following data allows to compare sequencing quality by: the total size of the genome and the haplotig size, the average coverage, the N50 value which means that half of the genome is formed with contigs bigger than this size and the percentage of gene detected from BUSCO analysis
DNA sequences obtained from different tissues dissection of leaf or berry skin (L1+L2) against flesh or roots (L2 only) [20, 22, 41]. This last method has also been used in other species such as bananas by comparing DNA from leaf, stem, rhizome and roots [42]. Although all these experiments demonstrate the existence of chimeras in plants and sometimes their crucial impact on agronomical traits, genome wide chimera detection has yet not been possible. Validation of chimeras is also a challenge because we expect low alternative allele frequency since the variant will appear on only one haplotype and one cell layer representing a small proportion of leaf tissue. Sanger sequencing has been found to be limited when alternative allele frequency is under 15–20% [43] but other technologies such as Molecular Inversion Probes (MIP) [44] has proven to be efficient in this particular condition [45]. Widely used in medical programs to detect rare diseases [46–48] it has also been used in plants to detect pathogens [49] or assist in genomic selection [50]. Because of repetitive sequences, it is also difficult to design specific target sequences that only capture the identified SNVs, MIPs along with MIPGEN designing software [51] is efficient on this specific criteria, it also has the advantage of being performant with a small amount of DNA (200ng) which also makes it useful in forensic applications [52, 53]. According to this information MIPs should be an interesting technology for chimera validation.

## Importance of ‘Merlot’ grapevine cultivar

'Merlot', which is a cross between 'Cabernet Franc' and 'Magdeleine Noire des Charentes' [54], is the grape cultivar used in this study. It was first mentioned in southwestern France in the late 18th century and expended in Bordeaux area since the middle of 19th century; the impressive spreading of this cultivar in other French regions and worldwide only dates from the 1970s [54]. Today it is the fourth most planted cultivar in the world

<table><caption>**Table 2** Sequencing quality information for the 4 samples</caption><thead><tr><th></th><th>Number of sequences</th><th>Average length (bp)</th><th>Sequence length (bp)</th><th>Coverage</th></tr></thead><tbody><tr><td>Merlot leaf</td><td>1 724 709</td><td>13 500</td><td>46 to 48<br>499</td><td>47x</td></tr><tr><td>Merlot root</td><td>1 901 058</td><td>13 500</td><td>47 to 49<br>695</td><td>51x</td></tr><tr><td>Cabernet Franc</td><td>2 310 232</td><td>12 500</td><td>48 to 45<br>849</td><td>58x</td></tr><tr><td>Magdeleine Noire des Charentes</td><td>1 876 162</td><td>13 500</td><td>44 to 49<br>698</td><td>51x</td></tr></tbody></table>

Sequencing quality is diagnosed for each sample through: the total number of sequences, their average length, their minimal and maximal length and the overall estimated coverage. Merlot leaf and Merlot root refer respectively to the sequencing of DNA extracted from leaves or from lateral roots

for table and wine grapes and the second cultivar for wine, cultivated in at least thirty-seven countries on 266 000 ha [55]. It is also the most planted variety in France with 114 578 ha in 2018 [56]. The international success of cv. 'Merlot' is mainly explained by the high quality red wines produced in renowned Bordeaux vineyards [57]. This cultivar is also one of the earliest black varieties to be harvested and thus one of the most impacted by climate change. In some areas, cultivating 'Merlot' could become inappropriate to produce high quality red wines because of cooked aromas and too high alcohol content [58]. Therefore, exploring 'Merlot' genome could open new perspectives to better understand its genetic and physiologic functioning as well as its intravarietal diversity required for clonal preservation and selection. New knowledge on 'Merlot' genome and chimeras could also help future grape breeding in order to create improved varieties with a similar fruit phenotype.

Throughout this study we take advantage of the latest sequencing and bioinformatics technologies not only to obtain a whole phased assembly 'Merlot' genome but also to contribute to a better understanding of a complex biological phenomena. We used parental sequences of cv. 'Cabernet Franc' and 'Magdeleine Noire des Charentes' to bin 'Merlot' reads in two groups, assemble the reads per haplotype and build pseudo-molecules. We compared root and leaf sequence to detect periclinal chimeras on each haplotype. We transferred gene annotation from 'PN40024' Vcost.v3 [26] to our pseudo-molecules in order to have a functional interpretation of the chimera's location. Finally, a subset of the chimeras was analysed by MIP in order to validate them with an independent technology.

## Results

### Building pseudo-molecules

DNA samples from 'Merlot' lateral roots and leaves and from leaves only for 'Magdeleine Noire des Charentes' (maternal) and 'Cabernet Franc' (paternal) have been sequenced using Pacific Bioscience Sequel II technology. For each sample between 1.7 and 2.3 million HIFI reads are obtained with an average length of 13 kb and 99.9% phred score accuracy. Taking 500 Mb as *Vitis* genome size we estimate a mean coverage between 47x and 58x according to the sample (Table 2).

Hifiasm trio-binning resolved the raw assembly of each haplotype with high confidence in regards to the statistics presented below (Table 3). Using 500 Mb as the expected genome size of *Vitis vinifera*, the mean N90 of 7,39 Mb with the L90 of 27 allows us to consider these results as being very good quality.

After two successful alignments, first on the PN40024. v4 reference [43] then on the second haplotype (see Material and Methods section), each contig was assigned
<table><caption>**Table 3** Assembly quality after trio-binning</caption><thead><tr><th>Sample name</th><th>N50 (Mb)</th><th>L50</th><th>N90 (Mb)</th><th>L90</th><th>Number of contigs</th><th>Total length (Mb)</th></tr></thead><tbody><tr><td>Merlot-root-Hap-CF</td><td>18.67</td><td>10</td><td>7.50</td><td>26</td><td>509</td><td>530</td></tr><tr><td>Merlot-root-Hap-MG</td><td>14.95</td><td>12</td><td>7.67</td><td>30</td><td>1 849</td><td>586</td></tr><tr><td>Merlot-leaf-Hap-CF</td><td>21.66</td><td>10</td><td>8.24</td><td>23</td><td>206</td><td>519</td></tr><tr><td>Merlot-leaf-Hap-MG</td><td>15.83</td><td>12</td><td>6.18</td><td>31</td><td>1 254</td><td>550</td></tr><tr><td>Average for all samples</td><td>17.77</td><td>11</td><td>7.39</td><td>27</td><td>955</td><td>546</td></tr></tbody></table>

Assembly quality results are defined by N50, L50, N90, L90 statistics. They are completed with the number of contigs assembled and the total length of the assembly. Merlot-root-Hap-CF refers to the root haplotype of Merlot genome transmitted by Cabernet Franc, Merlot-root-Hap-MG, the one transmitted by Magdeleine noire des Charentes, and so on

to its chromosome, their order and orientation were found. In fine, a unique contig to a maximum of 5 were needed to shape chromosomes and between 33 and 47 were used for the whole genome (Table 4). Thus, “Merlot haplotype Cabernet-Franc” (Merlot-hap-CF) is set to 486–490 Mb and “Merlot haplotype Magdeleine Noire des Charentes” (Merlot-hap-MG) to 491 Mb. Because of the sequencing technology and the high performance of long read assembly, we obtain longer chromosomes compared to PN40024.v4. Chromosome lengths are very similar between leaves and roots but trio-binning which correctly phase the assembly, there is a slight difference between Merlot-hap-CF and Merlot-hap-MG (Fig. 2).

From 1,6% to 2,8% of the assembly can't be accurately placed in the pseudo-molecule. It is mainly small highly repetitive sequences that are found in different places throughout the genome. This proportion was not attributed to a specific location and couldn't be associated to a contig.

The results of Benchmarking Universal Single-Copy Orthologs (BUSCO) using embryophyta lineage-specific databases [59] are resumed in Table 5. Results show that up to 98.7% of genes searched are found in the pseudo-molecules and nothing is lost compared to the raw assembly. Duplicated genes are reduced to 1.2–2.2% while missing genes remain around 0.4–0.6%. Therefore, those unplaced contigs were confidently ignored in the further analysis.

The annotation ‘PN40024’ Vcost.v3 was transferred to each pseudo-molecule using Liftoff tool [60]. In average 95% of the 42 413 genes were positioned throughout the 19 chromosomes with difference of gene numbers between chromosomes as expected (Fig. 3). BUSCO analysis was also performed on protein sequences from ‘PN40024’ Vcost.v3 files and obtained a score of 97,3% complete genes against 95,2% for Merlot-leaf-hap-CF and 95,1% for Merlot-leaf-hap-MG.

### Haplotype comparison

‘Magdeleine Noire des Charentes’ haplotype (Hap-MG) is slightly longer than ‘Cabernet Franc’ haplotype (Hap-CF). Merlot-leaf-hap-MG has 111 more genes than

Merlot-leaf-hap-CF (Table 4). All six groups of reads (Cabernet Franc, Magdeleine Noire des Charentes, Merlot-root-hap-CF, Merlot-root-hap-MG, Merlot-leaf-hap-CF, Merlot-leaf-hap-MG) are aligned on Merlot-root-hap-CF pseudo-molecule and DeepVariant is used to perform variant calling (see Material and Methods) [61]. Mapping the reads from both Cabernet Franc haplotypes back on the Merlot-root-hap-CF consensus pseudo-molecule allows us to estimate the amount of potential errors in sequencing or assembling. Merlot-root-hap-CF and Merlot-leaf-hap-CF have respectively 3.2k and 3.7k variant sites, mostly in repeated sequences against the pseudo-molecule (Table 6) which are thus potential sequencing errors but could also mean that mosaic mutations appear more frequently in repetitive regions. Mapping the reads from Merlot-MG haplotype and from both ‘Cabernet Franc’ and ‘Magdeleine Noire des Charentes’ lets us accurately compare haplotypes. We identified about 3.5 millions of variants between Merlot-leaf-hap-MG or Merlot-root-hap-MG and Merlot-root-hap-CF pseudo-molecules. These variants are 89% Single Nucleotide Variants (SNV), mainly located in repeated sequences, around 30% are included in a gene region and 5% in a coding region (Table 6).

### Chimera detection

Periclinal chimeras were detected by variant calling from the alignment of Merlot-leaf-hap-CF reads on Merlot-root-hap-CF pseudo-molecule. As presented in Fig. 1, chimeras can be located on either L1 or L2 cell layers. The comparison of both haplotypes and parental sequences allows to distinguish one case from the other. When Merlot-root-hap-CF reads (L2) and some Merlot-leaf-hap-CF reads (L1+L2) carry the same allele but all other sequences carry an alternative one, also present in Merlot-leaf-hap-CF reads, we considered that the L2 cell layer has mutated (Fig. 4). On the contrary, when all sequences have the same allele but a variant is confidently detected on Merlot-leaf-hap-CF, we consider the variant allocated to the L1 cell layer. When these mutations are confirmed by all root reads and are present only in a subset but not in all leaf reads, it means that the
<table><caption>**Table 4** Pseudomolecules characteristics for the haplotypes from both roots and leaves genomes</caption><thead><tr><th rowspan="2"></th><th colspan="3">Merlot leaf CF</th><th colspan="3">Merlot leaf MG</th><th colspan="3">Merlot root CF</th><th colspan="3">Merlot root MG</th></tr><tr><th>Length (bp)</th><th>Contig number</th><th>Gene number</th><th>Length (bp)</th><th>Contig number</th><th>Gene number</th><th>Length (bp)</th><th>Contig number</th><th>Gene number</th><th>Length (bp)</th><th>Contig number</th><th>Gene number</th></tr></thead><tbody><tr><td>chr01</td><td>23,598,264</td><td>2</td><td>2091</td><td>25,285,196</td><td>3</td><td>2308</td><td>23,589,260</td><td>2</td><td>2090</td><td>25,812,159</td><td>3</td><td>2229</td></tr><tr><td>chr02</td><td>21,618,960</td><td>2</td><td>1711</td><td>20,834,059</td><td>2</td><td>1771</td><td>21,612,554</td><td>2</td><td>1710</td><td>20,815,648</td><td>2</td><td>1773</td></tr><tr><td>chr03</td><td>22,910,969</td><td>2</td><td>1896</td><td>23,738,413</td><td>3</td><td>1853</td><td>22,695,420</td><td>4</td><td>1900</td><td>22,863,606</td><td>2</td><td>1856</td></tr><tr><td>chr04</td><td>28,028,847</td><td>1</td><td>2171</td><td>28,067,391</td><td>1</td><td>2199</td><td>27,848,422</td><td>1</td><td>2172</td><td>28,047,555</td><td>1</td><td>2195</td></tr><tr><td>chr05</td><td>28,252,839</td><td>3</td><td>2354</td><td>27,194,848</td><td>4</td><td>2190</td><td>28,133,326</td><td>2</td><td>2359</td><td>27,245,740</td><td>5</td><td>2197</td></tr><tr><td>chr06</td><td>21,902,219</td><td>1</td><td>1997</td><td>25,066,429</td><td>2</td><td>1901</td><td>21,887,941</td><td>2</td><td>1999</td><td>25,064,422</td><td>2</td><td>1900</td></tr><tr><td>chr07</td><td>31,335,031</td><td>1</td><td>2892</td><td>30,285,517</td><td>2</td><td>2877</td><td>31,316,563</td><td>1</td><td>2889</td><td>30,272,463</td><td>2</td><td>2875</td></tr><tr><td>chr08</td><td>25,084,147</td><td>2</td><td>2260</td><td>25,133,383</td><td>2</td><td>2159</td><td>25,068,184</td><td>1</td><td>2264</td><td>25,106,757</td><td>2</td><td>2163</td></tr><tr><td>chr09</td><td>23,385,614</td><td>3</td><td>1774</td><td>25,283,510</td><td>3</td><td>1789</td><td>23,388,363</td><td>3</td><td>1772</td><td>25,079,599</td><td>3</td><td>1785</td></tr><tr><td>chr10</td><td>26,473,290</td><td>1</td><td>2156</td><td>25,356,578</td><td>3</td><td>2159</td><td>26,466,742</td><td>1</td><td>2159</td><td>25,373,432</td><td>3</td><td>2160</td></tr><tr><td>chr11</td><td>20,347,376</td><td>2</td><td>1571</td><td>20,416,243</td><td>1</td><td>1563</td><td>19,944,127</td><td>2</td><td>1572</td><td>20,402,067</td><td>1</td><td>1569</td></tr><tr><td>chr12</td><td>27,064,986</td><td>2</td><td>2447</td><td>24,119,349</td><td>4</td><td>2347</td><td>24,025,996</td><td>1</td><td>2400</td><td>24,096,300</td><td>3</td><td>2343</td></tr><tr><td>chr13</td><td>29,400,470</td><td>2</td><td>2277</td><td>29,181,387</td><td>2</td><td>2244</td><td>29,412,682</td><td>3</td><td>2277</td><td>29,182,803</td><td>3</td><td>2245</td></tr><tr><td>chr14</td><td>30,137,549</td><td>1</td><td>2577</td><td>31,130,771</td><td>2</td><td>2564</td><td>30,138,415</td><td>1</td><td>2570</td><td>31,044,170</td><td>1</td><td>2545</td></tr><tr><td>chr15</td><td>23,518,695</td><td>2</td><td>1552</td><td>23,434,760</td><td>2</td><td>1832</td><td>23,494,571</td><td>2</td><td>1545</td><td>23,422,104</td><td>2</td><td>1828</td></tr><tr><td>chr16</td><td>22,496,135</td><td>1</td><td>1813</td><td>22,861,382</td><td>3</td><td>1840</td><td>22,480,932</td><td>3</td><td>1815</td><td>22,716,031</td><td>2</td><td>1844</td></tr><tr><td>chr17</td><td>20,156,470</td><td>2</td><td>1552</td><td>20,678,794</td><td>5</td><td>1628</td><td>20,368,439</td><td>3</td><td>1554</td><td>20,784,459</td><td>4</td><td>1628</td></tr><tr><td>chr18</td><td>36,729,617</td><td>2</td><td>3197</td><td>37,876,454</td><td>1</td><td>3161</td><td>36,718,487</td><td>2</td><td>3192</td><td>37,875,940</td><td>1</td><td>3164</td></tr><tr><td>chr19</td><td>27,613,142</td><td>1</td><td>1978</td><td>25,782,670</td><td>2</td><td>1992</td><td>27,618,512</td><td>2</td><td>1979</td><td>25,770,078</td><td>2</td><td>1994</td></tr><tr><td>Total length</td><td>490,054,620</td><td>33</td><td>40,266</td><td>491,727,134</td><td>47</td><td>40,377</td><td>486,208,936</td><td>38</td><td>40,218</td><td>490,975,333</td><td>44</td><td>40,293</td></tr><tr><td>Total+UKN</td><td>499,510,259</td><td></td><td>499,510,259</td><td></td><td></td><td>500,119,474</td><td></td><td></td><td></td><td>503,747,540</td><td></td><td></td></tr><tr><td>chrUn</td><td>9,455,639</td><td></td><td>7,783,125</td><td></td><td></td><td>13,910,538</td><td></td><td></td><td></td><td>12,772,207</td><td></td><td></td></tr><tr><td>%chrUn in total length</td><td>1.90%</td><td></td><td>1.60%</td><td></td><td></td><td>2.80%</td><td></td><td></td><td></td><td>2.50%</td><td></td><td></td></tr></tbody></table>

For each pseudo-molecule and per chromosome the data shows: the total length, the number of contig necessary to cover the entire chromosome and the number of genes detected. The total length of the pseudo-molecule is found below without or with (total+UKN) the unknown chromosome; it is compared with the size of the unknown chromosome (chrUn) which contains all contigs that were not confidently placed on a chromosome
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Chromosome length per haplotype compared to PN40024_12X.v4 genome. Chromosome length per haplotype in Mbp for each pseudomolecule built (Merlot-Root-Hap-CF; Merlot-Leaf-Hap-CF; Merlot-Root-Hap-MG; Merlot-Leaf-Hap-MG) against PN40024.v4.</figcaption></figure>

<table><caption>**Table 5** Search for genome completion using BUSCO embryophyta odb 10</caption><thead><tr><th rowspan="3"></th><th colspan="4">Merlot Hap CF</th><th colspan="4">Merlot Hap MG</th></tr><tr><th colspan="2">Leaf</th><th colspan="2">Root</th><th colspan="2">Leaf</th><th colspan="2">Root</th></tr><tr><th>Raw</th><th>Pseudo</th><th>Raw</th><th>Pseudo</th><th>Raw</th><th>Pseudo</th><th>Raw</th><th>Pseudo</th></tr></thead><tbody><tr><td>Complete</td><td>98.7%</td><td>98.7%</td><td>98.6%</td><td>98.5%</td><td>98.2%</td><td>98.2%</td><td>98.2%</td><td>98.2%</td></tr><tr><td>Complete and single copy</td><td>93.2%</td><td>96.5%</td><td>94.2%</td><td>97.1%</td><td>96.0%</td><td>97.0%</td><td>95.1%</td><td>97.0%</td></tr><tr><td>Complete and duplicated</td><td>5.5%</td><td>2.2%</td><td>4.4%</td><td>1.4%</td><td>2.2%</td><td>1.2%</td><td>3.1%</td><td>1.2%</td></tr><tr><td>Fragmented</td><td>0.8%</td><td>0.9%</td><td>1.0%</td><td>0.9%</td><td>1.2%</td><td>1.2%</td><td>1.2%</td><td>1.2%</td></tr><tr><td>Missing</td><td>0.5%</td><td>0.4%</td><td>0.4%</td><td>0.6%</td><td>0.6%</td><td>0.6%</td><td>0.6%</td><td>0.6%</td></tr><tr><td>Total groups searched</td><td colspan="8">1614</td></tr></tbody></table>

BUSCO was first performed on raw reads of each haplotype, just after trio-binning (Raw), second BUSCO was done on the pseudomolecules (Pseudo). Complete means that the gene is entirely found, whereas fragmented is not. Single copy means the gene is only found once whereas duplicated means it is found several times. Missing is the percentage of genes expected but not found on the genome. Total groups searched is the number of genes expected

entire L2 meristem cell layer carries the mutation and can be called periclinal chimeras. To increase confidence in detection, we focused on haplotype specific chimeras when no variant is found in reads of the opposite haplotype. We also chose to only select 'Merlot' specific chimeras and therefore excluded variants if parental reads were heterozygous. Grapevine DNA also has a lot of repeated sequences known to evolve more rapidly (e.g. microsatellites, transposable elements). In this study we focused on periclinal chimera that are located in non-repeated sequences because they are more stable and less prone to mapping errors. Only SNVs were kept. The work was executed on each haplotype separately. In total, 51 positions match the requirements on Merlot-hap-CF, and 53 on Merlot-hap-MG (Table 7).

Respectively 37 and 36 are found on L1 cell layer for Merlot-Hap-CF and Merlot-Hap-MG and 14 and 17

respectively on L2 cell layer. A total of 19 and 16 on each haplotype are located in a gene region, 9 and 7 chimeras are in a coding region. The exact position of the chimeras in the genome, the nucleotide and the number of reads for each allele, the type of chimera and location in coding region are presented in Tables 8 and 9.

## Validation of chimeras by MIPs

Out of the 104 positions identified as chimeras, MIP-GEN software was able to design a ligation and extension probe for 95 target regions. MIPseq was performed on Merlot_leaf and Merlot_root samples. Amplification was obtained for 86 positions but only 32 had enough depth to compare both samples (Table 10 and Additional files 2 and 3). 22 positions have the expected alleles on each sample and validate with enough depth PacBio results. 8 positions have the expected alleles but also have a single
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Number of genes per chromosome, for each Merlot sample and each haplotype. Number of gene per chromosome for each pseudomolecule detected by transferred annotation from PN40024 Vcost.v3 using Liftoff</figcaption></figure>

<table><caption>**Table 6** Variant calling statistics when reads were aligned on Merlot-root-hap-CF pseudomolecule</caption><thead><tr><th>Reads aligned on<br>Merlot-root-Hap-CF</th><th>Filtered variant<br>sites</th><th>% of SNV in<br>filtered variants</th><th>% of filtered variants<br>included in repeated<br>Sequences</th><th>% of variant sites<br>included in genes</th><th>% of variant<br>sites included in coding<br>region</th></tr></thead><tbody><tr><td>Merlot root hap CF (himself)</td><td>3 207</td><td>71%</td><td>63%</td><td>28%</td><td>6%</td></tr><tr><td>Merlot leaf hap CF</td><td>3 731</td><td>67%</td><td>67%</td><td>24%</td><td>5%</td></tr><tr><td>Merlot leaf hap MG</td><td>3 561 942</td><td>89%</td><td>60%</td><td>29%</td><td>4%</td></tr><tr><td>Merlot root hap MG</td><td>3 490 079</td><td>89%</td><td>60%</td><td>29%</td><td>4%</td></tr><tr><td>Cabernet Franc (paternal)</td><td>2 730 140</td><td>86%</td><td>52%</td><td>36%</td><td>5%</td></tr><tr><td>Magdeleine N. C. (maternal)</td><td>4 987 693</td><td>89%</td><td>57%</td><td>31%</td><td>5%</td></tr></tbody></table>

This table contains the number of variants found by mapping each packet of reads to the Merlot-Root-Hap-CF pseudomolecule after processing quality filtering, the percentage of these variants that are Single Nucleotide Variant, the percentage of variants included in repeated sequences, the percentage included in genes and the percentage included in a coding region

read that is unexpected which makes them ambiguous. One position is invalid because an allele is missing on Merlot leaf sample (chr3-3419274 on MG haplotype). Finally, one position is classified as ambiguous (chr10-22333673 on MG haplotype) because it has both alleles on leaves and roots. This doesn't necessarily invalid the existence of the chimera but makes it unclear on which cell layer it is located on.

Results of molecular inversion probes sequencing is shown per haplotype; the location of the previously identified chimera is detailed. For both samples, first column contains the alleles found by PacBio sequencing, the second one the alleles revealed by MIP and the third one the number of reads that support each allele separated by "/" symbol. Valid positions have the expected alleles according to PacBio results and ambiguous or invalid have unexpected alleles.

## Discussion

The combination of long reads high quality sequencing, trio binning using parental sequences, and a long read assembler gives the opportunity to resolve accurately phased assembly. This new 'Merlot' genome encompasses a total length of about 500 Mbp, and was constructed with only 33 to 47 contigs. Some chromosomes were resolved with a single contig, while others needed up to 5. The number of contigs for each chromosome for the 4 haplotypes were however different, therefore not linked to a specific chromosome, but most probably to sequence depth. Chromosome lengths and gene number were slightly different between both haplotypes. The gene numbers may however be underestimated since no de novo annotation was done. In the objective towards the definition of grape pangenome, precise *de novo* annotation of 'Merlot' genome should be performed. According to Shumate and Salzberg [60], Liftoff can accurately transfer 99.9% of the genes when working intra-species. Here, PN40024 reference genome for *Vitis vinifera* has
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Schematic genetic interpretation of L2 periclinal chimeras in grape cv. 'Merlot': Different allele configurations expected for L2 periclinal SNV 'Merlot' specific chimeras; a L2 periclinal chimera should be identified in Merlot leaf and root of the same haplotype but not found in either the opposite haplotype nor in parental reads. Each cell layer is represented here as a "stick", two for leaves and only on for roots. The SNV is represented as either A or G. Cabernet Franc and Magdeleine noire des Charentes leaves are also represented</figcaption></figure>

<table><caption>**Table 7** Number of chimeras per haplotype and per cell layer</caption><thead><tr><th rowspan="2"></th><th colspan="3">Merlot Haplotype CF</th><th colspan="3">Merlot Haplotype MG</th></tr><tr><th>L1</th><th>L2</th><th>L1+L2</th><th>L1</th><th>L2</th><th>L1+L2</th></tr></thead><tbody><tr><td>SNV, Non-repeated sequences Merlot specific chimeras</td><td>37</td><td>14</td><td><strong>51</strong></td><td>36</td><td>17</td><td><strong>53</strong></td></tr><tr><td>Number of chimeras in non coding gene body</td><td>15</td><td>4</td><td><strong>19</strong></td><td>11</td><td>5</td><td><strong>16</strong></td></tr><tr><td>Number of chimeras in coding region</td><td>6</td><td>3</td><td><strong>9</strong></td><td>6</td><td>1</td><td><strong>7</strong></td></tr></tbody></table>

Chimeras are found by comparing root and leaves DNA. Results are shown per haplotype and per location on the cell layer: either L1 or L2, L1+L2 is the sum of both. First line gives the total number of position detected, second line is the number of position included in a gene body region but not in a coding region and the last line is the number included in a coding region inside the gene

been used to transfer genes on a *Vitis vinifera* cultivar. Because of this intra specific design it is expected that liftoff accurately transferred most of the genes, this is confirmed by the close results of BUSCO analysis performed on protein sequences, 97% of complete genes for PN40024 against 95% for Merlot pseudomolecules. Moreover, the program works by aligning the reference genes on the target sequence, so although it can't identify new genes the mapping correspondence can be considered accurate.

The pseudomolecules obtained in this study were compared to existent whole genome sequencing data

available in the literature, information is compiled in "Additional file 1". The closest genome found in terms of assembly quality is the 'Cabernet-Sauvignon' assembly [32]. However, taking advantage of the newest technologies the pseudo molecules obtained in this work have longer assemblies (~490 Mb), longer average sequence lengths (~26 Mb), longer maximum lengths (~37 Mb), longer N50 (~25 Mb) and higher Busco scores (~98%) with less missing genes (~0.5%). Moreover, these technologies are time saving as the complete assembly and pseudomolecule building can be done in a couple of days. Liftoff offers a fast tool to transfer annotation, 96%
<table><caption>**Table 8** SNVs Chimeras found on Merlot-hap-CF</caption><thead><tr><th>CELL LAYER</th><th>Chr</th><th>Pos</th><th>Merlot Leaf Hap CF</th><th>Merlot Root Hap CF</th><th>Merlot Leaf Hap MG</th><th>Merlot Root Hap MG</th><th>Cabernet</th><th>Magdeleine</th><th>Annotation of Merlot root hap CF pseudomolecule</th></tr></thead><tbody><tr><td>L1</td><td>chr01</td><td>622,811</td><td>G/T</td><td>16/7</td><td>G</td><td>22</td><td>G</td><td>3</td><td>G</td><td>G</td><td>59</td><td>No</td></tr><tr><td>L1</td><td>chr01</td><td>10,384,263</td><td>G/A</td><td>28/9</td><td>G</td><td>31</td><td>G</td><td>28</td><td>G</td><td>23</td><td>G</td><td>52</td><td>Gene</td></tr><tr><td>L1</td><td>chr02</td><td>7,380,049</td><td>T/C</td><td>22/6</td><td>T</td><td>23</td><td>T</td><td>14</td><td>T</td><td>19</td><td>T</td><td>43</td><td>CDS</td></tr><tr><td>L1</td><td>chr04</td><td>23,822,006</td><td>C/G</td><td>20/8</td><td>C</td><td>19</td><td>C</td><td>21</td><td>C</td><td>32</td><td>C</td><td>48</td><td>Gene</td></tr><tr><td>L1</td><td>chr04</td><td>25,819,229</td><td>A/G</td><td>12/5</td><td>A</td><td>26</td><td>A</td><td>29</td><td>A</td><td>43</td><td>A</td><td>65</td><td>No</td></tr><tr><td>L1</td><td>chr05</td><td>22,934,305</td><td>C/T</td><td>22/7</td><td>C</td><td>21</td><td>C</td><td>16</td><td>C</td><td>30</td><td>C</td><td>25</td><td>Gene</td></tr><tr><td>L1</td><td>chr05</td><td>23,010,268</td><td>A/G</td><td>7/31</td><td>A</td><td>33</td><td>A</td><td>29</td><td>A</td><td>19</td><td>A</td><td>45</td><td>No</td></tr><tr><td>L1</td><td>chr05</td><td>25,416,770</td><td>C/T</td><td>15/5</td><td>C</td><td>19</td><td>C</td><td>22</td><td>C</td><td>25</td><td>C</td><td>57</td><td>No</td></tr><tr><td>L1</td><td>chr06</td><td>3,621,220</td><td>T/C</td><td>10/8</td><td>T</td><td>32</td><td>T</td><td>2</td><td>T</td><td>4</td><td>T</td><td>65</td><td>No</td></tr><tr><td>L1</td><td>chr06</td><td>12,834,474</td><td>A/G</td><td>20/10</td><td>A</td><td>16</td><td>A</td><td>25</td><td>A</td><td>19</td><td>A</td><td>56</td><td>No</td></tr><tr><td>L1</td><td>chr07</td><td>8,599,845</td><td>C/T</td><td>22/8</td><td>C</td><td>30</td><td>C</td><td>26</td><td>C</td><td>33</td><td>C</td><td>44</td><td>Gene</td></tr><tr><td>L1</td><td>chr07</td><td>12,188,378</td><td>A/G</td><td>23/8</td><td>A</td><td>29</td><td>A</td><td>29</td><td>A</td><td>33</td><td>A</td><td>44</td><td>No</td></tr><tr><td>L1</td><td>chr08</td><td>20,988,006</td><td>C/T</td><td>17/11</td><td>C</td><td>29</td><td>C</td><td>24</td><td>C</td><td>35</td><td>C</td><td>43</td><td>CDS</td></tr><tr><td>L1</td><td>chr08</td><td>22,606,874</td><td>T/C</td><td>18/5</td><td>T</td><td>37</td><td>T</td><td>14</td><td>T</td><td>9</td><td>T</td><td>31</td><td>No</td></tr><tr><td>L1</td><td>chr10</td><td>16,298,623</td><td>A/G</td><td>22/13</td><td>A</td><td>30</td><td>A</td><td>21</td><td>A</td><td>28</td><td>A</td><td>51</td><td>Gene</td></tr><tr><td>L1</td><td>chr11</td><td>7,012,825</td><td>G/A</td><td>23/10</td><td>G</td><td>30</td><td>G</td><td>30</td><td>G</td><td>26</td><td>G</td><td>53</td><td>Gene</td></tr><tr><td>L1</td><td>chr12</td><td>8,907,103</td><td>C/A</td><td>17/8</td><td>C</td><td>23</td><td>C</td><td>23</td><td>C</td><td>17</td><td>C</td><td>28</td><td>CDS</td></tr><tr><td>L1</td><td>chr12</td><td>16,700,578</td><td>A/G</td><td>27/14</td><td>A</td><td>37</td><td>A</td><td>17</td><td>A</td><td>16</td><td>A</td><td>35</td><td>Gene</td></tr><tr><td>L1</td><td>chr13</td><td>1,078,880</td><td>A/G</td><td>22/10</td><td>A</td><td>29</td><td>A</td><td>22</td><td>A</td><td>25</td><td>A</td><td>53</td><td>No</td></tr><tr><td>L1</td><td>chr13</td><td>19,791,949</td><td>A/T</td><td>25/8</td><td>A</td><td>28</td><td>A</td><td>22</td><td>A</td><td>29</td><td>A</td><td>48</td><td>No</td></tr><tr><td>L1</td><td>chr13</td><td>23,444,751</td><td>A/T</td><td>23/9</td><td>A</td><td>30</td><td>A</td><td>28</td><td>A</td><td>29</td><td>A</td><td>63</td><td>Gene</td></tr><tr><td>L1</td><td>chr14</td><td>2,105,195</td><td>A/C</td><td>24/7</td><td>A</td><td>25</td><td>A</td><td>17</td><td>A</td><td>23</td><td>A</td><td>52</td><td>No</td></tr><tr><td>L1</td><td>chr14</td><td>8,620,581</td><td>T/C</td><td>24/7</td><td>T</td><td>24</td><td>T</td><td>21</td><td>T</td><td>21</td><td>T</td><td>42</td><td>Gene</td></tr><tr><td>L1</td><td>chr14</td><td>15,844,322</td><td>C/T</td><td>19/7</td><td>C</td><td>12</td><td>C</td><td>27</td><td>C</td><td>28</td><td>C</td><td>45</td><td>No</td></tr><tr><td>L1</td><td>chr15</td><td>19,297,087</td><td>C/T</td><td>27/8</td><td>C</td><td>33</td><td>C</td><td>26</td><td>C</td><td>33</td><td>C</td><td>53</td><td>CDS</td></tr><tr><td>L1</td><td>chr15</td><td>20,162,161</td><td>G/A</td><td>23/11</td><td>G</td><td>31</td><td>G</td><td>20</td><td>G</td><td>35</td><td>G</td><td>70</td><td>CDS</td></tr><tr><td>L1</td><td>chr15</td><td>20,298,163</td><td>C/A</td><td>16/8</td><td>C</td><td>22</td><td>C</td><td>16</td><td>C</td><td>24</td><td>C</td><td>38</td><td>Gene</td></tr><tr><td>L1</td><td>chr16</td><td>18,511,980</td><td>T/A</td><td>19/9</td><td>T</td><td>35</td><td>T</td><td>23</td><td>T</td><td>32</td><td>T</td><td>59</td><td>No</td></tr><tr><td>L1</td><td>chr17</td><td>8,087,306</td><td>C/A</td><td>14/7</td><td>C</td><td>23</td><td>C</td><td>13</td><td>C</td><td>15</td><td>C</td><td>47</td><td>No</td></tr><tr><td>L1</td><td>chr17</td><td>9,516,464</td><td>C/G</td><td>18/8</td><td>C</td><td>32</td><td>C</td><td>27</td><td>C</td><td>16</td><td>C</td><td>66</td><td>Gene</td></tr><tr><td>L1</td><td>chr17</td><td>12,223,344</td><td>A/C</td><td>18/10</td><td>A</td><td>32</td><td>A</td><td>9</td><td>A</td><td>2</td><td>A</td><td>22</td><td>Gene</td></tr><tr><td>L1</td><td>chr18</td><td>4,105,052</td><td>T/C</td><td>16/6</td><td>T</td><td>24</td><td>T</td><td>29</td><td>T</td><td>29</td><td>T</td><td>43</td><td>Gene</td></tr><tr><td>L1</td><td>chr18</td><td>5,367,698</td><td>C/A</td><td>23/6</td><td>C</td><td>34</td><td>C</td><td>33</td><td>C</td><td>41</td><td>C</td><td>63</td><td>Gene</td></tr><tr><td>L1</td><td>chr18</td><td>11,379,480</td><td>G/T</td><td>22/9</td><td>G</td><td>37</td><td>G</td><td>27</td><td>G</td><td>26</td><td>G</td><td>30</td><td>Gene</td></tr><tr><td>L1</td><td>chr18</td><td>12,312,986</td><td>C/A</td><td>12/4</td><td>C</td><td>27</td><td>C</td><td>26</td><td>C</td><td>32</td><td>C</td><td>26</td><td>No</td></tr><tr><td>L1</td><td>chr18</td><td>21,684,668</td><td>G/A</td><td>32/10</td><td>G</td><td>27</td><td>G</td><td>4</td><td>G</td><td>7</td><td>G</td><td>27</td><td>No</td></tr><tr><td>L1</td><td>chr19</td><td>6,420,354</td><td>C/T</td><td>9/7</td><td>C</td><td>22</td><td>C</td><td>34</td><td>C</td><td>30</td><td>C</td><td>66</td><td>CDS</td></tr></tbody></table>
<table><caption>Table 8 (continued)</caption><thead><tr><th>CELL LAYER</th><th>Chr</th><th>Pos</th><th>Merlot Leaf Hap CF</th><th>Merlot Root Hap CF</th><th>Merlot Leaf Hap MG</th><th>Merlot Root Hap MG</th><th>Cabernet</th><th>Magdeleine</th><th>Annotation of Merlot root hap CF pseudomolecule</th></tr></thead><tbody><tr><td>L2</td><td>chr04</td><td>26,224,776</td><td>T/C</td><td>21/7</td><td>T</td><td>34</td><td>C</td><td>29</td><td>C</td><td>56</td></tr><tr><td>L2</td><td>chr05</td><td>24,987,737</td><td>C/T</td><td>22/11</td><td>C</td><td>31</td><td>T</td><td>26</td><td>T</td><td>61</td></tr><tr><td>L2</td><td>chr06</td><td>17,243,628</td><td>T/A</td><td>20/10</td><td>T</td><td>44</td><td>A</td><td>20</td><td>A</td><td>58</td></tr><tr><td>L2</td><td>chr07</td><td>10,337,296</td><td>G/A</td><td>11/5</td><td>G</td><td>22</td><td>A</td><td>26</td><td>A</td><td>53</td></tr><tr><td>L2</td><td>chr10</td><td>8,445,101</td><td>T/A</td><td>17/12</td><td>T</td><td>24</td><td>A</td><td>26</td><td>A</td><td>47</td></tr><tr><td>L2</td><td>chr11</td><td>9,424,215</td><td>A/G</td><td>19/8</td><td>A</td><td>27</td><td>G</td><td>28</td><td>G</td><td>57</td></tr><tr><td>L2</td><td>chr12</td><td>22,268,147</td><td>A/T</td><td>22/9</td><td>A</td><td>21</td><td>T</td><td>11</td><td>T</td><td>27</td></tr><tr><td>L2</td><td>chr14</td><td>3,806,097</td><td>T/C</td><td>9/6</td><td>T</td><td>30</td><td>C</td><td>21</td><td>C</td><td>29</td></tr><tr><td>L2</td><td>chr15</td><td>13,526,354</td><td>A/G</td><td>23/8</td><td>A</td><td>29</td><td>G</td><td>26</td><td>G</td><td>60</td></tr><tr><td>L2</td><td>chr17</td><td>3,559,600</td><td>C/T</td><td>20/10</td><td>C</td><td>32</td><td>T</td><td>16</td><td>T</td><td>49</td></tr><tr><td>L2</td><td>chr17</td><td>7,610,259</td><td>C/T</td><td>15/5</td><td>C</td><td>24</td><td>T</td><td>30</td><td>T</td><td>46</td></tr><tr><td>L2</td><td>chr18</td><td>12,292,200</td><td>G/A</td><td>19/8</td><td>G</td><td>35</td><td>A</td><td>28</td><td>A</td><td>77</td></tr><tr><td>L2</td><td>chr19</td><td>468,364</td><td>G/T</td><td>25/10</td><td>G</td><td>32</td><td>T</td><td>32</td><td>T</td><td>66</td></tr><tr><td>L2</td><td>chr19</td><td>17,464,909</td><td>A/G</td><td>19/6</td><td>A</td><td>21</td><td>G</td><td>10</td><td>G</td><td>15</td></tr></tbody></table>

For each chimeria, the cell layer is identified, the chromosome and the exact location is given according to the reference Merlot-root-Hap-CF pseudo molecule. Then for each sample the allele(s) are completed by the number of reads confirming the nucleotide. If there is only one nucleotide, all reads of this haplotype converge, when two alleles are found the number of reads supporting each one is respectively given by the corresponding numbers. For the first position (chr 01, position 622 811), 16 reads have "G" nucleotide and 7 have "T" nucleotide. In the last column the position is either included in a coding region inside a gene (CDS), included in a gene but not in a coding region (gene) or not included in either (No).
of the genes from 'PN40024' Vcost.v3 were successfully replaced on all four pseudo-molecules allowing a functional interpretation of the results.

## Advantage of having parental reads

The originality of this study was to not only sequence the cultivar of interest but also its parents. Parental sequences have allowed to discriminate 'Merlot' reads from each haplotype. Each haplotype was then assembled independently, as if we had two homozygous individuals. This step increases confidence in haplotype comparison statistics. Haplotype differences in a single individual (~3.5 million of variant sites) are similar to what has already been reported [41]. The pseudomolecule used as reference for the variant calling is the Merlot-root-hap-CF which shares half of the 'Cabernet Franc's DNA. This explains why 2.7 million variants are detected for 'Cabernet Franc' (difference for only one haplotype) whereas 'Magdeleine Noire des Charentes' displayed 4.9 million variants, both haplotypes being different. In addition, being able to retrace parental origin could make it possible to know what agronomic character comes from each parent which increases possibilities in breeding and cultivar improvement. Considering the difference between both 'Merlot' haplotypes. ~60% variants are located on repeated sequences, ~30% are located in gene regions and ~6% in coding regions. These numbers align with the apportionment of each of one in the genome [25]. This suggests that the variants are not preferentially located in coding or repeated sequences. However, this doesn't fit with previous publications on 'Pinot Noir' [62] or 'Nebbiolo' [63] that found more variations in coding regions between both haplotypes. This could be a specificity of 'Merlot' haplotypes since differences in variation rates have already been noticed between 'Nebbiolo' and 'Zinfandel' or it could also be explained by trio binning technology that allows to rebuild each haplotype more accurately and therefore have a better appreciation of the comparison.

## Chimera identification and their impact on the phenotype

Until now, chimera detection was only possible with PCR sequencing when three alleles were found on the same locus or by dissecting tissues derived from different meristem cell layers as cited above. However, a genome wide screening of chimeras was not yet possible with these methods. Throughout this study we show that quantity and high quality sequencing, long reads, trio binning and organ comparison and strong selection open new doors in chimera detection. Around 3 000 variants were found when mapping Merlot-leaf-hap-CF and Merlot-root-hap-CF reads to the Merlot-root-hap-CF pseudomolecule, these are mainly SNVs (67–71%). The variants identified by mapping Merlot-root-hap-CF reads on the consensus

pseudomolecule most certainly correspond to sequencing, trio binning or assembly errors but it is also possible that some of these variants are sectorial chimeras and only located in a few cells. To detect chimeras, we remove these 3 000 positions. In addition, we focused on variants outside repeated sequences which are easier to map and more likely to be stable during evolution and less prone to errors. We also focused on variants that meet periclinal chimera definition because they are the most stable. They indeed meet very specific conditions but they are also the fewest. Nevertheless, the very selective criteria applied allow us to confidently identify these variants as being chimeras. It is not excluded that other types of chimeras exist but were not selected in this work. Indeed, a mutation can be present in a few cells of one or both cell layers and appear as a variant site but it would need extra experiments to truly validate them.

Similar amounts of SNV periclinal chimeras were found on each haplotype (51 and 53). These results seem to mean that they appear randomly and at the same frequency on both haplotypes. Among those, 70% correspond to mutations on L1 cell layer and 30% on L2 cell layer. Some sequencing errors detected on leaf samples and not on roots could explain this difference between L1 and L2 although such difference in frequency could also make sense because L1 cell layer is located on the surface of leaves and is more exposed to UV radiation. Moreover, L2 cell layer produces gametes and are probably more protected [64]. Validation on independent reads is needed the support this last theory.

The consequences of a chimera depend on its position on the genome, in our study 33% are located in a gene body region and 15% are located in a coding region and could modify the protein which can be perceived on the phenotype (Tables 7, 8, 9, 11 and 12). Although our data confirms this possibility, the phenomena appears to be a rare event.

MIP sequencing allowed to validate with confidence a subset of positions which makes chimera detection through hifi PacBio long reads and trio-binning reliable. However, MIP sequencing results overall did not have the depth expected compared to what is described in the literature leading to the loss of more than half of the positions tested. This means that either MIP target region design or laboratory protocol should be optimised. Having unexpected alleles only supported by a single read makes conclusion ambiguous but could be due to sequencing errors, mutation induced by PCR or it could be due to the higher sensitivity of the MIP technology to detect rare mutations. It appears that each technology has its own pros and cons and only a cross result between two sequencing technologies can bring a high confidence in the detection of the chimeras. However, PacBio technology seems trustworthy to detect SNVs on one hand
<table><caption>**Table 9** SNVs Chimeras found on Merlot-hap-MG</caption><thead><tr><th>Cell Layer</th><th>Chr</th><th>Position</th><th>Merlot Leaf Hap MG</th><th>Merlot Root Hap MG</th><th>Merlot Leaf Hap CF</th><th>Merlot Root Hap CF</th><th>Cabernet</th><th>Magdeleine</th><th>Annotation of Merlot root hap MG</th></tr></thead><tbody><tr><td>L1</td><td>chr02</td><td>3,007,441 A/G</td><td>26/7</td><td>A</td><td>33 A</td><td>24 A</td><td>27 A</td><td>24 A</td><td>32 No</td></tr><tr><td>L1</td><td>chr03</td><td>3,419,274 T/A</td><td>12/7</td><td>T</td><td>33 T</td><td>23 T</td><td>16 T</td><td>66 T</td><td>47 CDS</td></tr><tr><td>L1</td><td>chr04</td><td>27,055,780 A/T</td><td>18/7</td><td>A</td><td>27 A</td><td>22 A</td><td>35 A</td><td>64 A</td><td>38 Gene</td></tr><tr><td>L1</td><td>chr05</td><td>8,794,703 C/T</td><td>17/8</td><td>C</td><td>37 C</td><td>9 C</td><td>14 C</td><td>20 C</td><td>41 Gene</td></tr><tr><td>L1</td><td>chr05</td><td>8,797,586 C/A</td><td>23/9</td><td>C</td><td>31 C</td><td>8 C</td><td>3 C</td><td>19 C</td><td>27 Gene</td></tr><tr><td>L1</td><td>chr05</td><td>12,819,049 T/A</td><td>16/7</td><td>T</td><td>24 T</td><td>25 T</td><td>32 T</td><td>50 T</td><td>31 Gene</td></tr><tr><td>L1</td><td>chr05</td><td>18,874,107 C/T</td><td>22/10</td><td>C</td><td>23 C</td><td>24 C</td><td>26 C</td><td>47 C</td><td>52 Gene</td></tr><tr><td>L1</td><td>chr07</td><td>16,367,184 C/T</td><td>23/7</td><td>C</td><td>37 C</td><td>16 C</td><td>20 C</td><td>54 C</td><td>63 No</td></tr><tr><td>L1</td><td>chr08</td><td>13,911,900 C/A</td><td>16/9</td><td>C</td><td>33 C</td><td>15 C</td><td>22 C</td><td>44 C</td><td>63 CDS</td></tr><tr><td>L1</td><td>chr08</td><td>18,990,379 A/T</td><td>16/7</td><td>A</td><td>10 A</td><td>20 A</td><td>17 A</td><td>53 A</td><td>31 No</td></tr><tr><td>L1</td><td>chr09</td><td>2,285,675 G/T</td><td>23/7</td><td>G</td><td>34 G</td><td>22 G</td><td>24 G</td><td>61 G</td><td>51 No</td></tr><tr><td>L1</td><td>chr09</td><td>7,569,730 A/G</td><td>16/8</td><td>A</td><td>21 A</td><td>3 A</td><td>3 A</td><td>35 A</td><td>43 No</td></tr><tr><td>L1</td><td>chr09</td><td>11,451,290 G/A</td><td>20/8</td><td>G</td><td>30 G</td><td>28 G</td><td>23 G</td><td>30 G</td><td>32 No</td></tr><tr><td>L1</td><td>chr10</td><td>2,278,377 C/A</td><td>20/13</td><td>C</td><td>29 C</td><td>24 C</td><td>29 C</td><td>81 C</td><td>71 Gene</td></tr><tr><td>L1</td><td>chr10</td><td>5,337,832 A/G</td><td>16/10</td><td>A</td><td>24 A</td><td>11 A</td><td>12 A</td><td>44 A</td><td>48 CDS</td></tr><tr><td>L1</td><td>chr10</td><td>10,968,103 A/C</td><td>18/5</td><td>A</td><td>36 A</td><td>27 A</td><td>33 A</td><td>65 A</td><td>38 Gene</td></tr><tr><td>L1</td><td>chr10</td><td>22,333,673 T/A</td><td>18/6</td><td>T</td><td>28 T</td><td>28 T</td><td>14 T</td><td>85 T</td><td>63 Gene</td></tr><tr><td>L1</td><td>chr11</td><td>10,058,232 A/C</td><td>8/12</td><td>A</td><td>33 A</td><td>19 A</td><td>21 A</td><td>43 A</td><td>41 No</td></tr><tr><td>L1</td><td>chr11</td><td>10,357,099 T/C</td><td>15/5</td><td>T</td><td>27 T</td><td>16 T</td><td>13 T</td><td>42 T</td><td>49 No</td></tr><tr><td>L1</td><td>chr11</td><td>13,563,841 T/C</td><td>10/6</td><td>T</td><td>23 T</td><td>4 T</td><td>6 T</td><td>16 T</td><td>16 No</td></tr><tr><td>L1</td><td>chr11</td><td>13,732,865 T/C</td><td>21/4</td><td>T</td><td>33 T</td><td>16 T</td><td>27 T</td><td>61 T</td><td>68 CDS</td></tr><tr><td>L1</td><td>chr11</td><td>13,734,413 A/C</td><td>22/5</td><td>A</td><td>31 A</td><td>20 A</td><td>27 A</td><td>64 A</td><td>69 Gene</td></tr><tr><td>L1</td><td>chr13</td><td>15,741,877 A/T</td><td>14/7</td><td>A</td><td>18 A</td><td>27 A</td><td>24 A</td><td>86 A</td><td>74 No</td></tr><tr><td>L1</td><td>chr13</td><td>18,415,015 A/G</td><td>17/9</td><td>A</td><td>24 A</td><td>3 A</td><td>4 A</td><td>25 A</td><td>27 Gene</td></tr><tr><td>L1</td><td>chr13</td><td>22,993,757 T/A</td><td>16/14</td><td>T</td><td>33 T</td><td>31 T</td><td>37 T</td><td>53 T</td><td>56 No</td></tr><tr><td>L1</td><td>chr13</td><td>24,830,909 A/T</td><td>18/10</td><td>A</td><td>20 A</td><td>19 A</td><td>32 A</td><td>65 A</td><td>43 Gene</td></tr><tr><td>L1</td><td>chr14</td><td>5,457,647 T/A</td><td>22/6</td><td>T</td><td>19 T</td><td>3 T</td><td>6 T</td><td>17 T</td><td>28 No</td></tr><tr><td>L1</td><td>chr14</td><td>30,136,904 A/G</td><td>13/9</td><td>A</td><td>30 A</td><td>9 A</td><td>14 A</td><td>60 A</td><td>54 No</td></tr><tr><td>L1</td><td>chr15</td><td>2,385,711 C/T</td><td>10/6</td><td>C</td><td>16 C</td><td>2 C</td><td>2 C</td><td>26 C</td><td>21 No</td></tr><tr><td>L1</td><td>chr16</td><td>6,369,473 C/T</td><td>20/6</td><td>C</td><td>24 C</td><td>6 C</td><td>7 C</td><td>24 C</td><td>21 No</td></tr><tr><td>L1</td><td>chr16</td><td>15,506,287 C/T</td><td>22/7</td><td>C</td><td>29 C</td><td>14 C</td><td>22 C</td><td>62 C</td><td>27 No</td></tr><tr><td>L1</td><td>chr16</td><td>17,466,134 C/T</td><td>24/8</td><td>C</td><td>28 C</td><td>16 C</td><td>17 C</td><td>45 C</td><td>67 No</td></tr><tr><td>L1</td><td>chr17</td><td>7,800,769 G/A</td><td>13/5</td><td>G</td><td>22 G</td><td>20 G</td><td>14 G</td><td>49 G</td><td>21 CDS</td></tr><tr><td>L1</td><td>chr17</td><td>16,802,779 C/T</td><td>26/9</td><td>C</td><td>37 C</td><td>35 C</td><td>33 C</td><td>76 C</td><td>50 No</td></tr><tr><td>L1</td><td>chr17</td><td>20,480,355 C/T</td><td>10/4</td><td>C</td><td>32 C</td><td>25 C</td><td>28 C</td><td>49 C</td><td>39 CDS</td></tr><tr><td>L2</td><td>chr03</td><td>17,695,770 G/A</td><td>23/7</td><td>G</td><td>24 A</td><td>22 A</td><td>18 A</td><td>44 A</td><td>66 Gene</td></tr></tbody></table>
<table><caption>**Table 9 (continued)**</caption><thead><tr><th>Cell Layer</th><th>Chr</th><th>Position</th><th>Merlot Leaf Hap MG</th><th>Merlot Root Hap MG</th><th>Merlot Leaf Hap CF</th><th>Merlot Root Hap CF</th><th>Cabernet</th><th>Magdeleine</th><th>Annota-<br>tion of<br>Merlot root<br>hap MG</th></tr></thead><tbody><tr><td>L2</td><td>chr05</td><td>1,753,844 T/C</td><td>10/8</td><td>T</td><td>24 C</td><td>17 C</td><td>11 C</td><td>55 C</td><td>43 No</td></tr><tr><td>L2</td><td>chr05</td><td>5,349,008 A/G</td><td>19/9</td><td>A</td><td>23 G</td><td>28 G</td><td>35 G</td><td>63 G</td><td>75 Gene</td></tr><tr><td>L2</td><td>chr07</td><td>20,018,317 G/A</td><td>23/7</td><td>G</td><td>25 A</td><td>31 A</td><td>27 A</td><td>64 A</td><td>66 No</td></tr><tr><td>L2</td><td>chr08</td><td>5,971,148 A/G</td><td>20/10</td><td>A</td><td>30 G</td><td>26 G</td><td>17 G</td><td>71 G</td><td>56 Gene</td></tr><tr><td>L2</td><td>chr08</td><td>8,665,113 A/C</td><td>18/6</td><td>A</td><td>26 C</td><td>32 C</td><td>28 C</td><td>36 C</td><td>51 No</td></tr><tr><td>L2</td><td>chr08</td><td>10,405,377 C/T</td><td>17/4</td><td>C</td><td>27 T</td><td>28 T</td><td>30 T</td><td>71 T</td><td>77 No</td></tr><tr><td>L2</td><td>chr08</td><td>16,511,687 A/G</td><td>23/5</td><td>A</td><td>37 G</td><td>23 G</td><td>20 G</td><td>50 G</td><td>60 No</td></tr><tr><td>L2</td><td>chr08</td><td>18,095,503 C/T</td><td>24/8</td><td>C</td><td>33 T</td><td>26 T</td><td>17 T</td><td>60 T</td><td>52 No</td></tr><tr><td>L2</td><td>chr10</td><td>3,422,488 G/T</td><td>17/5</td><td>G</td><td>30 T</td><td>31 T</td><td>23 T</td><td>28 T</td><td>98 No</td></tr><tr><td>L2</td><td>chr11</td><td>1,451,292 C/T</td><td>22/8</td><td>C</td><td>32 T</td><td>32 T</td><td>43 T</td><td>80 T</td><td>63 No</td></tr><tr><td>L2</td><td>chr12</td><td>7,859,048 A/C</td><td>22/11</td><td>A</td><td>19 C</td><td>24 C</td><td>30 C</td><td>66 C</td><td>58 Gene</td></tr><tr><td>L2</td><td>chr12</td><td>14,034,710 T/G</td><td>15/8</td><td>T</td><td>16 G</td><td>13 G</td><td>25 G</td><td>50 G</td><td>24 No</td></tr><tr><td>L2</td><td>chr13</td><td>25,521,690 A/G</td><td>25/7</td><td>A</td><td>26 G</td><td>24 G</td><td>29 G</td><td>58 G</td><td>51 Gene</td></tr><tr><td>L2</td><td>chr14</td><td>16,511,432 T/A</td><td>10/8</td><td>T</td><td>21 A</td><td>22 A</td><td>22 A</td><td>42 A</td><td>51 No</td></tr><tr><td>L2</td><td>chr15</td><td>20,391,775 T/G</td><td>20/5</td><td>T</td><td>30 G</td><td>37 G</td><td>31 G</td><td>76 G</td><td>50 CDS</td></tr><tr><td>L2</td><td>chr17</td><td>17,822,477 C/T</td><td>20/7</td><td>C</td><td>25 T</td><td>12 T</td><td>19 T</td><td>47 T</td><td>62 No</td></tr></tbody></table>

For each chimera, the cell layer is identified, the chromosome and the exact location is given according to the reference Merlot-root-Hap-MG pseudo molecule. Then for each sample the allele(s) are completed by the number of reads confirming the nucleotide. If there is only one nucleotide, all reads of this haplotype converge, when two alleles are found the number of reads supporting each one is respectively given by the corresponding numbers. For the first position (chr 02, position 1 351 259), 15 reads have G nucleotide and 6 have T. In the last column the position is either included in a coding region inside a gene (CDS), included in a gene but not in a coding region (gene) or not included in either (No).
<table><caption>**Table 10** Chimera validation with molecular inversion probes sequencing</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Cell layer</th><th rowspan="2">Chro-mosome</th><th rowspan="2">Position</th><th colspan="2">Merlot 343 Leaf</th><th colspan="2">Merlot 343 Root</th><th rowspan="2">Validation conclusion</th></tr><tr><th>Expected alleles according to PacBio</th><th>Mip validation</th><th>Expected alleles according to PacBio</th><th>Mip validation</th></tr></thead><tbody><tr><td rowspan="15">Chimeras on Merlot haplotype 'Cabernet Franc'</td><td>L1</td><td>chr04</td><td>23 822 006</td><td>C/G</td><td>C/G</td><td>25/10</td><td>C</td><td>C</td><td>Validated</td></tr><tr><td>L1</td><td>chr05</td><td>25 416 770</td><td>C/T</td><td>C/T</td><td>125/24</td><td>C</td><td>C</td><td>Validated</td></tr><tr><td>L1</td><td>chr08</td><td>22 606 874</td><td>T/C</td><td>T/C</td><td>16/6</td><td>T</td><td>T</td><td>Validated</td></tr><tr><td>L1</td><td>chr10</td><td>16 298 623</td><td>A/G</td><td>A/G</td><td>101/10</td><td>A</td><td>A</td><td>Validated</td></tr><tr><td>L1</td><td>chr12</td><td>8 907 103</td><td>C/A</td><td>C/A</td><td>34/2</td><td>C</td><td>C</td><td>Validated</td></tr><tr><td>L1</td><td>chr12</td><td>16 700 578</td><td>A/G</td><td>A/G</td><td>47/15</td><td>A</td><td>A/G</td><td>1 unexpected read</td></tr><tr><td>L1</td><td>chr14</td><td>8 620 581</td><td>T/C</td><td>T/C</td><td>190/16</td><td>T</td><td>T/C</td><td>1 unexpected read</td></tr><tr><td>L1</td><td>chr14</td><td>15 844 322</td><td>C/T</td><td>C/T</td><td>71/4</td><td>C</td><td>C</td><td>Validated</td></tr><tr><td>L1</td><td>chr15</td><td>20 162 161</td><td>G/A</td><td>G/A</td><td>11/1</td><td>G</td><td>G/A</td><td>1 unexpected read</td></tr><tr><td>L1</td><td>chr18</td><td>4 105 052</td><td>T/C</td><td>T/C</td><td>102/6</td><td>T</td><td>T</td><td>Validated</td></tr><tr><td>L2</td><td>chr17</td><td>7 610 259</td><td>C/T</td><td>C/T</td><td>6/29</td><td>C/T</td><td>C/T</td><td>Validated</td></tr><tr><td>L2</td><td>chr18</td><td>12 292 200</td><td>G/A</td><td>G/A</td><td>56/67</td><td>G/A</td><td>G/A/C</td><td>1 unexpected read</td></tr><tr><td rowspan="10">Chimeras on Merlot Haplotype 'Magdeleine Noire des Charentes'</td><td>L1</td><td>chr03</td><td>3 419 274</td><td>T/A</td><td>T</td><td>75</td><td>T</td><td>T</td><td>Invalid</td></tr><tr><td>L1</td><td>chr05</td><td>8 797 586</td><td>C/A</td><td>C/A</td><td>35/5</td><td>C</td><td>C</td><td>Validated</td></tr><tr><td>L1</td><td>chr05</td><td>12 819 049</td><td>T/A</td><td>T/A</td><td>61/5</td><td>T</td><td>T</td><td>Validated</td></tr><tr><td>L1</td><td>chr07</td><td>16 367 184</td><td>C/T</td><td>C/T</td><td>172/16</td><td>C</td><td>C</td><td>Validated</td></tr><tr><td>L1</td><td>chr08</td><td>13 911 900</td><td>C/A</td><td>C/A/G</td><td>67/13/1</td><td>C</td><td>C</td><td>1 unexpected read</td></tr><tr><td>L1</td><td>chr09</td><td>7 569 730</td><td>A/G</td><td>A/G</td><td>20/1</td><td>A</td><td>A</td><td>Validated by only 1 read</td></tr><tr><td>L1</td><td>chr10</td><td>2 278 377</td><td>C/A</td><td>C/A</td><td>110/18</td><td>C</td><td>C</td><td>Validated</td></tr><tr><td>L1</td><td>chr10</td><td>22 333 673</td><td>T/A</td><td>T/A</td><td>53/26</td><td>T</td><td>T/A</td><td>Confusion on cell layer</td></tr><tr><td>L1</td><td>chr11</td><td>10 357 099</td><td>T/C</td><td>T/C</td><td>26/10</td><td>T</td><td>T</td><td>Validated</td></tr><tr><td>L1</td><td>chr13</td><td>22 993 757</td><td>T/A</td><td>T/A</td><td>17/2</td><td>T</td><td>T/C</td><td>1 unexpected read</td></tr></tbody></table>

and also makes it possible to determine on which haplotype and which cell layer chimeras are located on.

Throughout this study we have made a specific focus on single nucleotide variant because they are more stable. Yet some essential functions or characteristics of grapevine such as berry colour can be altered by structural variants [20], studying these types of variants would also be of interest.

Chimeras are rare but they can have a strong impact on phenotype. If they are identified and selected, they can lead to a new cultivar as it has already been reported with 'Pinot Gris'. In a less obvious evolution, perennial plants propagated over centuries only through cuttings, chimeras are most likely to accumulate over time and could slowly induce genetic diversity among the cultivar. By

continuously selecting the best plant to fit specific characteristics, breeders increase their chance to select and propagate useful chimeras. When chimeras are stable and conserved through several generations of cuttings, they could also be used to trace and identify clonal lineage. Since we have developed a tool for revealing chimeras, it would be interesting to analyse the presence of a subset of the chimeric mutations in different 'Merlot' clone 343 plants in order to check how stable these chimeras might be. For grapevine, clonal identification is an important issue because no low cost and rapid test can guarantee clonal origin, although it is the economic unit used today. Clonal lineage is only done by human traceability which can contain errors especially after a long period of time.
<table><caption>**Table 11** Description of the periclinal chimeras on Merlot-Root-Hap-CF located in a gene</caption><thead><tr><th>Chr</th><th>Position on Merlot Root Hap CF</th><th>Cell layer</th><th>Coding region</th><th>Gene orientation</th><th>REF Codon</th><th>REF Codon Amino Acid</th><th>ALT Codon Amino Acid</th><th>ALT Amino Acid</th><th>Gene name</th><th>LOC</th><th>NAME</th></tr></thead><tbody><tr><td>chr 01</td><td>10,384,263 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi01g00875</td><td>LOC100263854</td><td>disease resistance protein TAO1</td></tr><tr><td>chr 02</td><td>7,380,049 L1</td><td>CDS</td><td></td><td>reverse</td><td>AAC</td><td>N</td><td>AGC</td><td>S</td><td>Vitvi02g00607.CDS16</td><td>LOC100260794</td><td>uncharacterized</td></tr><tr><td>chr 04</td><td>23,822,006 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi04g02171.t01</td><td>LOC104879118</td><td>protein FREE1</td></tr><tr><td>chr 05</td><td>22,934,305 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi05g01344.t01</td><td>LOC100248036</td><td>eukaryotic translation initiation factor 3 subunit H</td></tr><tr><td>chr 07</td><td>8,599,845 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi07g00711</td><td>LOC100251541</td><td>WD repeat-containing protein 26</td></tr><tr><td>chr 08</td><td>20,988,006 L1</td><td>CDS</td><td></td><td>reverse</td><td>AAG</td><td>K</td><td>AAA</td><td>K</td><td>Vitvi08g01573.CDS1</td><td>LOC100259003</td><td>caffeic acid 3-O-methyltransferase</td></tr><tr><td>chr 10</td><td>8,445,101 L2</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi10g01813.t01</td><td>LOC104880497</td><td>uncharacterized</td></tr><tr><td>chr 10</td><td>16,298,623 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi10g01120.t01</td><td>LOC100252938</td><td>lysine-specific demethylase JMJ18</td></tr><tr><td>chr 11</td><td>7,012,825 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi19g00471.t01</td><td></td><td></td></tr><tr><td>chr 11</td><td>9,424,215 L2</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi11g00777.t01</td><td>LOC100267388</td><td>autophagy-related protein 11</td></tr><tr><td>chr 12</td><td>8,907,103 L1</td><td>CDS</td><td></td><td>forward</td><td>CTC</td><td>L</td><td>CTA</td><td>L</td><td>Vitvi12g02494.CDS1</td><td>LOC100251126</td><td>myelin transcription factor 1-like protein</td></tr><tr><td>chr 12</td><td>16,700,578 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi12g01666.t01</td><td>LOC100245109</td><td>protein FAM135B</td></tr><tr><td>chr 12</td><td>22,268,147 L2</td><td>CDS</td><td></td><td>reverse</td><td>TAG</td><td>STOP</td><td>AAG</td><td>K</td><td>Vitvi12g02051.CDS3</td><td>LOC104881</td><td></td></tr><tr><td>chr 13</td><td>23,444,751 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi13g01459.t01</td><td>LOC100263567</td><td>lowering time control protein FY [Vitis vinifera]</td></tr><tr><td>chr 14</td><td>8,620,581 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi14g00503</td><td>LOC100250122</td><td>phosphoribosylaminoimidazole carboxylase, chloroplastic</td></tr><tr><td>chr 15</td><td>13,526,354 L2</td><td>CDS</td><td></td><td>forward</td><td>AAC</td><td>N</td><td>AGC</td><td>S</td><td>Vitvi15g00455.CDS2</td><td>LOC100260889</td><td>homeobox-leucine zipper protein ANTHOCY-ANINLESS 2</td></tr><tr><td>chr 15</td><td>19,297,087 L1</td><td>CDS</td><td></td><td>forward</td><td>TCC</td><td>S</td><td>TTC</td><td>F</td><td>Vitvi15g00839.CDS59</td><td>LOC100245390</td><td>rho GTPase-activating protein 5 [Vitis vinifera]</td></tr><tr><td>chr 15</td><td>20,162,161 L1</td><td>CDS</td><td></td><td>reverse</td><td>ACC</td><td>T</td><td>ACT</td><td>T</td><td>Vitvi15g00907.CDS5</td><td>LOC100245390</td><td>autophagy-related protein 11</td></tr><tr><td>chr 15</td><td>20,298,163 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi00g01964.t01</td><td>LOC100267388</td><td>autophagy-related protein 11</td></tr><tr><td>chr 17</td><td>3,559,600 L2</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi17g00302.t01</td><td>LOC100251937</td><td>zinc finger protein 346</td></tr><tr><td>chr 17</td><td>7,610,259 L2</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi17g01488.t01</td><td>LOC100250198</td><td>serine/threonine-protein kinase ATM</td></tr><tr><td>chr 17</td><td>9,516,464 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi17g00813.t01</td><td>LOC100255285</td><td>callose synthase 10</td></tr><tr><td>chr 17</td><td>12,223,344 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi17g01006.t01</td><td>LOC100247288</td><td>acetyl-CoA carboxylase 1</td></tr><tr><td>chr 18</td><td>4,105,052 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi18g00368.t01</td><td>LOC100255806</td><td>WD and tetraricopeptide repeats protein 1</td></tr><tr><td>chr 18</td><td>5,367,698 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi18g00496</td><td>LOC100255806</td><td>purple acid phosphatase 2</td></tr><tr><td>chr 18</td><td>11,379,480 L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Vitvi18g02778.t01</td><td>LOC100247127</td><td>probable leucine-rich repeat receptor-like</td></tr><tr><td>chr 19</td><td>468,364 L2</td><td>CDS</td><td></td><td>reverse</td><td>CCG</td><td>A</td><td>GAG</td><td>E</td><td>Vitvi19g01801.CDS6</td><td>LOC100255520</td><td>serine/threonine-protein kinase At3g14840</td></tr><tr><td>chr 19</td><td>6,420,354 L1</td><td>CDS</td><td></td><td>reverse</td><td>GTG</td><td>V</td><td>GTA</td><td>V</td><td>Vitvi19g00471.CDS1</td><td></td><td></td></tr></tbody></table>

All chimeras included in a gene are shown in this table. Their exact location on Merlot-Root-Hap-CF and the cell layer are given. If the chimeras is located in a coding region it is marked “CDS” otherwise it is “gene”. When the position is in a coding region, the gene orientation either reversed or forward is given with the reference codon and amino acid, the nucleotide carrying the mutation is written in bold. The alternative sections show the impact of the mutation on the codon and the amino acid. Three last columns allow gene identification through gene name, LOC code and name.
## Conclusion

Through this study, whole genome DNA sequence was obtained using the latest genomic technologies and bioinformatics tools. Hifi long read sequencing, trio-binning, long read assembler, have all together allowed to obtain high quality, haplotype resolved pseudo-molecules. In addition, repeat masker tools, mapping and deep variant calling opened new possibilities in chimera detection. By comparing root and leaf samples and through severe selection it has been possible to identify hundred chimeras based on SNVs on both haplotypes. MIP validation has confirmed the presence of these chimeras. Other types of chimeras could be present, but we were not able to identify them. A functional interpretation was done through transferred annotation. Actual genomic tools open new doors in chimera detection, representing opportunities for perennial plant breeding. In addition, this high quality ‘Merlot’ genome, could also open new perspectives such as structural variants identifications, but could also serve as a basis for a study of intra-varietal variability for this cultivar.

## Materials and methods

### DNA sequencing

'Magdeleine Noire des Charentes' leaves were harvested from INRAE Vassal-Montpellier grape collection (Marseillan, FRANCE), while 'Merlot' clone 343 leafs and roots as well as 'Cabernet Franc' were harvested from IFV collection, Domaine de l'Espiguette (Grau du Roi, FRANCE). Two young leaves about 10 cm wide were collected, carefully rolled over and placed in a 13 ml tube. Secondary lateral roots from the same 'Merlot' clone 343 plant were also collected on the same day. This plant was not grafted and was destined to be pulled out which made it possible to collect its roots. All samples were conserved in a -80°C freezer until the DNA extraction process. DNA was extracted following the Tip 100 Qiagen Genomic kit with slight modifications. Lysis was performed 3 hours at 50°C on 0.5 g of ground plant material with 9.5 ml of G2 buffer supplemented with 1% PVP-40, 19 µl of RNase A and 500 µl of proteinase K. After tip filtrations, DNA was precipitated with isopropanol, centrifuged 15 min à 5000 g, washed with Ethanol 70° and re-suspended in 50 µl of TE buffer. DNA quality and high molecular weight were controlled. DO 260/280 ratio between 1.8 and 2.0 and DO 260/230 ratio between 2.0 and 2.2 were confirmed and an Agilent Genomic DNA Screen Tape was performed. Fifteen µg of high quality DNA were then used to carry out the sequencing. Samples were sequenced using Single Molecule Real Time PacBio SEQUEL II hifi long reads at INRAE Clermont-Ferrand GENTYANE platform (France).

### Assembly and building pseudo-molecules

DNA consensus call sequences obtained under BAM format were converted to fastq using bam2fastq tool from SMRTLink v9.0.0 PacBio library. The HIFI sequencing DNA quality was verified using FastQC version 0.11.7.

Figure 5 illustrates the whole bioinformatics workflow to build pseudo-molecules and transfer annotation. Paternal and maternal kmers were identified using the parental reads with yak-0.1 software. The outputs were then used in hifiasm-0.13 with default parameters to bin 'Merlot' long reads and assemble both haplotypes. This was done on both organs (leaf and root).

For each haplotype, contigs were aligned on PN40024. v4 using minimap2 version 2.17 [65]. Best contigs alignments were used to build an AGP file and from there reconstruct each pseudo-molecule. In order to refine the pseudo-molecules, we then reexecuted the same process starting with an alignment of each haplotype on the other previously reconstructed.

The embryophyta_odb10 lineage from BUSCO 5.3.1 software was carried out in genome mode to estimate the completeness of all assemblies [59]. BUSCO was also performed on protein sequences using "prot" option, protein sequences were obtained from the pseudomolecules by using gffread tool version 0.12.6 with default parameters.

Liftoff 1.6.1 tool with default parameters was used to transfer the annotation of PN40024 Vcost.v3 reference genome to the pseudomolecules [60].

### Chimera detection

Reads were mapped on Merlot-root-hap-CF and Merlot-root-hap-MG pseudomolecules with Minimap2 version 2.17 [65] with the option -x map-hifi and variant calling was performed with DeepVariant software version 1.1.0 [66] using PacBio model and default parameters. Finally filtering variants was done with vcftools 0.1.16 version [67].

Chimera detection was processed by filtering vcf from Merlot-leaf on Merlot-root pseudomolecule. We only conserved variants with more than 10 depth coverage, "PASS" quality flag and genotype quality (GQ) over 20. Non homozygous positions on all other sequences were excluded. Repeated sequences were identified by building a specific 'Merlot' library with repeatmodeler/2.0.2a-bin [68] and then using repeatmasker/4.1.1 software [69], and all chimeras in repeated sequences were excluded. Both repeatmodeler and repeatmasker were used with default parameters. Only single nucleotide variants were kept. Finally, Tables 7, 8 and 9 were manually checked site per site by visualisation in Integrative Genome Viewer (IGV 2.12.3) that allows a larger overview of the region on several samples [70]. These sites were crossed with the annotation file with intersect Bed function of BEDtools/2.30.0 [71] and Table 10 and 11. were completed.
<table><caption>**Table 12** Description of the periclinal chimeras on Merlot-Hap-MG located in a gene region</caption><thead><tr><th>Chr</th><th>Position on<br>Merlot Root<br>Hap MG</th><th>Cell<br>layer</th><th>Coding<br>region</th><th>Gene<br>orientation</th><th>REF<br>Codon</th><th>REF<br>Amino<br>Acid</th><th>ALT<br>Codon</th><th>ALT<br>Amino<br>Acid</th><th>Gene name</th><th>LOC</th><th>NAME</th></tr></thead><tbody><tr><td>chr 03</td><td>3,419,274</td><td>L1</td><td>CDS</td><td>reverse</td><td>ATC</td><td>I</td><td>TTC</td><td>F</td><td>Vitvi03g01484CDS2</td><td>LOC100243450</td><td>isoflavone reductase-like protein</td></tr><tr><td>chr 03</td><td>17,695,770</td><td>L2</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi03g01777</td><td></td><td></td></tr><tr><td>chr 04</td><td>27,055,780</td><td>L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi04g01695.t01</td><td>LOC100259141</td><td>GDP-mannose transporter GONST1</td></tr><tr><td>chr 05</td><td>5,349,008</td><td>L2</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi05g00464</td><td>LOC100261408</td><td>chloride channel CIC6</td></tr><tr><td>chr 05</td><td>8,794,703</td><td>L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi05g00751.t01</td><td></td><td></td></tr><tr><td>chr 05</td><td>8,797,586</td><td>L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi05g00751.t01</td><td></td><td></td></tr><tr><td>chr 05</td><td>12,819,049</td><td>L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi05g00975.t01</td><td>LOC100242520</td><td>dihydroxy-acid dehydratase, chloroplastic</td></tr><tr><td>chr 05</td><td>18,874,107</td><td>L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi05g01197.t01</td><td></td><td></td></tr><tr><td>chr 08</td><td>5,971,148</td><td>L2</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi08g00297</td><td>LOC100245232</td><td>uncharacterized</td></tr><tr><td>chr 08</td><td>13,911,900</td><td>L1</td><td>CDS</td><td>reverse</td><td>CGG</td><td>R</td><td>CTG</td><td>L</td><td>Vitvi08g00910CDS1</td><td>LOC100249911</td><td>uncharacterized LOC100249911</td></tr><tr><td>chr 10</td><td>2,278,377</td><td>L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi00g01047.t01</td><td>LOC100853350</td><td>sister chromatid cohesion protein PDS5 homolog B</td></tr><tr><td>chr 10</td><td>5,337,832</td><td>L1</td><td>CDS</td><td>forward</td><td>AGA</td><td>R</td><td>AGG</td><td>R</td><td>Vitvi10g00473CDS1</td><td>LOC100258824</td><td>putative nuclear RNA export factor SDE5</td></tr><tr><td>chr 10</td><td>10,968,103</td><td>L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi10g00847.t01</td><td>LOC104880520</td><td>uncharacterized LOC104880520</td></tr><tr><td>chr 10</td><td>22,333,673</td><td>L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi10g01425.t01</td><td></td><td></td></tr><tr><td>chr 11</td><td>13,732,865</td><td>L1</td><td>CDS</td><td>forward</td><td>AGT</td><td>S</td><td>AGC</td><td>S</td><td>Vitvi11g00979CDS8</td><td>LOC100246466</td><td>DEAD-box ATP-dependent RNA helicase 10</td></tr><tr><td>chr 11</td><td>13,734,413</td><td>L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi11g00979.t01</td><td>LOC100246466</td><td>DEAD-box ATP-dependent RNA helicase 10</td></tr><tr><td>chr 12</td><td>7,859,048</td><td>L2</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi12g00621</td><td>LOC100263032</td><td>5' exonuclease Apollo</td></tr><tr><td>chr 13</td><td>18,415,015</td><td>L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi13g01203.t01</td><td>LOC104881233</td><td>uncharacterized LOC104881233</td></tr><tr><td>chr 13</td><td>24,830,909</td><td>L1</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi13g02403.t01</td><td>LOC100261918</td><td>annexin D4-like</td></tr><tr><td>chr 13</td><td>25,521,690</td><td>L2</td><td>gene</td><td></td><td></td><td></td><td></td><td></td><td>Vitvi13g02429</td><td>LOC100244776</td><td>uncharacterized</td></tr><tr><td>chr 15</td><td>20,391,775</td><td>L2</td><td>CDS</td><td>forward</td><td>TCA</td><td>S</td><td>GCA</td><td>A</td><td>Vitvi15g00944exon15</td><td>LOC100255258</td><td>uncharacterized</td></tr><tr><td>chr 17</td><td>7,800,769</td><td>L1</td><td>CDS</td><td>forward</td><td>TGC</td><td>C</td><td>TAC</td><td>Y</td><td>Vitvi17g00653CDS1</td><td></td><td></td></tr><tr><td>chr 17</td><td>20,480,355</td><td>L1</td><td>CDS</td><td>reverse</td><td>GTA</td><td>V</td><td>ATA</td><td>I</td><td>Vitvi17g01265CDS1</td><td></td><td></td></tr></tbody></table>

All chimeras included in a gene are shown in this table. Their exact location on Merlot-Root-Hap-MG and the cell layer are given. If the chimaera is located in a coding region it is marked “CDS” otherwise it is “gene”. When the position is in a coding region, the gene orientation either reversed or forward is given with the reference codon and amino acid, the nucleotide carrying the mutation is written in bold. The alternative sections show the impact of the mutation on the codon and the amino acid. Three last columns allow gene identification through gene name, LOC code and name
<figure><img src="image_7.png" /><figcaption>Step 2. Trio-binning using parental sequences</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Step 3. Hifiasm assembly for each haplotype</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Step 4. Building pseudo-molecule on chromosome scale for each haplotype</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>Step 5. Transfered annotation from reference with Liftoff software</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 5** Bioinformatic workflow applied in this study. The workflow is described step by step. Step 1 is the kmers specific selection on each parent, these are variable size sequences that allow to specifically recognize reads from on parent. In step 2 these kmers are used to sort child reads out in two haplotypes that are each specific to one parent. If reads can't be attributed to one parent, they are considered to be in both. Step 3 is the assembly of several reads into contigs for each haplotype. Step 4 is the building of pseudomolecules using multiple alignments. Step 5 is the transfer of annotation from the reference genome to the pseudomolecules</figcaption></figure>
## Chimera validation

MIPGEN software [51] was used to design mips specific target regions previously identified as chimeras with following parameters: tag sizes 0.8 to introduce UMI (Unique Molecular Identifiers) to filter out duplicate reads and PCR errors, minimum ligature length 20, extension minimum length 16, arm length sums 36, 37, 38, 39, 40, minimum capture size 120, maximum capture size 150 and trf option was activated. DNA samples were adjusted in quantity using the previous DNA extraction and used in adapted MIP library protocol previously described [72] with some modifications. 100 ng of DNA template was added to a hybridization mix together with the oligo MIP pool (final concentration of 0.025 pM per probe) in 0.85x Ampligase buffer (Epicentre). Mix was incubated in a thermal cycler at 95 °C for 10 min, followed by a 60 °C cycle overnight. Products were mixed with dNTPs (Jena Bioscience, 15 pM), Betaine (Sigma-Aldrich, 375 mM), NAD+ (New England, Biolabs 1 mM), additional Ampligase buffer (0.75×), Ampligase (Epicentre, 1.25 U) and Klentaq (New England Biolabs, 0.16U). Mixture was incubated at 56 °C for 60 min followed by 72 °C for 20 min. Enzymatic digestion of linear probes was performed at 37 °C for 2 h, followed by 80 °C for 20 min by adding Exonuclease I (New England Biolabs, 8 U) and Exonuclease III (New England Biolabs, 50 U). Final product was amplified using Q5 Hot start High-Fidelity DNA Polymerase (New England Biolabs, 8 U) with different index combinations. PCR cycling conditions were an initial denaturation step for 2 min at 98 °C, followed by 20 cycles of 30 s at 98 °C, 20 s at 60 °C, and 20 s at 72 °C. PCR Samples were pooled and clean up using AMPureXP beads (Beckman Coulter) at 0.8× ratio. Samples were sequenced in 2×150 bp paired end mode using a MiSeq (Illumina) platforms with custom sequencing primers. UMI were extracted from obtained reads using umi_tools version 1.1.4 extract [73] with --extract-method=string and --bc-pattern=NNNNNNNNNNNN. Adapters were trimmed using cutadapt version 3.5 [74] with following parameters: -q 30 -m 100 -e 0.10 -a ACACTACCGTCCGGATCGTGCGTGT -A CTTCAGCTTCCCGATTACGGATCTCGTATG. SNP calling was done using process_reseq from VCFhunter version 2.2.0 with -s acefg option [75]. Finally, variant calling file was filtered when depth was below 10 for at least one sample.

### Abbreviations

RT PCR Real-Time Polymerase Chain Reaction

SSR Single Sequence Repeat

SNV Single nucleotide Variant

## Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s12864-023-09453-8.

**Additional file 1:** Additional_file_1.xls. Comparative data between whole genome sequences already published and the pseudomolecules built in this study.

**Additional file 2:** Additional_file_2.xls. PacBio and MIP data for each chimera found on Merlot-hap-MG.

**Additional file 3:** Additional_file_3.xls. PacBio and MIP data for each chimera found on Merlot-hap-CF.

## Acknowledgements

We are grateful to the members of the AGAP-DAAV team and especially Charles Romieu, for their help and their contribution in discussions. We acknowledge Laurent Torregrosa for his help in understanding grapevine cell layer physiology. We sincerely thank the Petrus team: Jean-Claude Berrouet, Olivier Berrouet, Fabienne Caillon, Lionel Caillon, Gilles Rabeyroux, Michael Paiva, Pierre-Jean Dalesme, Johann Ventre, Catherine Gaillard and Emilie Verral for their support during this work. We would like to thank Cécile Marchal and Sandrine Dedet for supplying samples from INRAE grapevine collection Vassal-Montpellier. We highly acknowledge the work of Carine Satgé from Toulouse CNRGV for adapting the DNA extraction protocol to grapevine.

## Authors' contributions

L.T., T.P., L.L., conceptualised the project. L.L., collected some samples. L.V., R.C., performed DNA extraction and sample preparation, S.V., S.G., G.N., B.P.F., built and executed bioinformatic work. M.P. and R.M. performed the MIP validation from samples preparation to Miseq sequencing. All authors contributed to writing the manuscript. The final manuscript was read and approved by the authors.

## Funding

This work is a collaboration between SC Petrus wine producer in Pomerol (France), National Research Institute for Agriculture, Food and Environment (INRAE) and the French institute for grapevine and wine (IFV). Fees and Phd student salary were covered by Petrus while INRAE, CIRAD and IFV contributed through researcher's salaries.

## Data Availability

Raw reads of PacBio sequencing are available on the European Nucleotide Archive repository, under the project named PRJEB59893: https://www.ebi.ac.uk/ena/browser/view/PRJEB59893.

Contigs, AGP file, chromosome scale assembly and annotations transferred from PN40024 VCost.v3 with liftoff are available on Recherche Data Gouv: https://doi.org/10.57745/OJ07SN.

PN40024 sequence and annotation used in this study are available on INTEGRAPE platform, https://integrape.eu/resources/genes-genomes/genome-accessions/.

Merlot clone 343 is available at INRAE, domaine de Vassal, under the code 0Mtp2399.

## Declarations

### Conflict of interest

The authors declare that they have no competing interest.

### Ethics approval and consent to participate

Not applicable.

### Consent for publication

Not applicable.

Received: 4 September 2022 / Accepted: 9 June 2023

Published online: 14 July 2023

## References

1. Bauhin C. Illustrated exposition of plants (the Pinax theatri botanici): sumptibus [et] typis Ludovici Regis; 1598.
2. Nati P. Petri Nati... Florentina phytologica obseruatio de malo limonia citrata-aurantia florentiae vulgo la bizzaria. typis Hippolyti de Naue; 1644.

3. Darwin C, Wallace A. On the variation of organic beings in a state of nature. J Proc Linn Soc Lond (Zoology). 1858;3:45–52.

4. Marcotrigiano M. Chimeras and variegation: patterns of Deceit. HortScience. 1997;32(5):773–84.

5. Winkler H. About grafted hybrids and plant chimeras. Ber dtsch bot Ges. 1907;25:568–76.

6. Baur E. Untersuchungen über die Erblichkeitsverhältnisse einer nur in Bастardform lebensfähigen Sippe von Antirrhinum majus. Bornträger; 1907.

7. Blakeslee AF, Avery AG, Bergner AD, Satina SA, Sinnott EW. Induction of periclinal chimeras in Datura stramonium by colchicine treatment. Science. 1939;89(2314):402.

8. Satina S, Blakeslee AF. Periclinal chimeras in datura stramonium in relation to development of leaf and flower. Am J Bot. 1941;28(10):862–71.

9. Satina S. Periclinal chimeras in Datura in relation to the development and structure of the ovule. Am J Bot 1945:72–81.

10. Szymkowiak EJ, Sussex IM. What chimeras can tell us about plant development. Annu Rev Plant Biol. 1996;47(1):351–76.

11. Satina S, Blakeslee AF, Avery AG. Demonstration of the three germ layers in the shoot apex of datura by means of induced polyploidy in periclinal chimeras. Am J Bot. 1940;27(10):895–905.

12. Thompson MM, Olmo H. Cytohistological studies of cytochimeric and tetraploid grapes. Am J Bot. 1963;50(9):901–6.

13. Frank MH, Chitwood DH. Plant chimeras: the good, the bad, and the 'Bizzaria'. Dev Biol. 2016;419(1):41–53.

14. Franks T, Botta R, Thomas MR, Franks J. Chimerism in grapevines: implications for cultivar identity, ancestry and genetic improvement. Theor Appl Genet. 2002;104(2–3):192–9.

15. Marcotrigiano M. Genetic mosaics and the analysis of Leaf Development. Int J Plant Sci. 2001;162:513–25.

16. Torregrosa L, Fernandez L, Bouquet A, Boursiquot J-M, Pelsy F, Martínez-Zapater JM. Origins and consequences of somatic variation in grapevine. Genet genomics Breed grapes. 2011;68:92.

17. Kazemian M, Mohajel Kazemi E, Kolahi M, Omran V. Floral ontogeny and molecular evaluation of anthocyanin biosynthesis pathway in pinwheel phenotype of Saintpaulia inthona Wendl. Periclinal chimera. Sci Hort. 2020;263:109142.

18. D'Amato F. Role of somatic mutations in the evolution of higher plants. Caryologia. 1997;50(1):1–15.

19. Carbonell-Bejerano P, Royo C, Mauri N, Ibáñez J. Miguel Martínez Zapater J. Somatic Variation and Cultivar Innovation in Grapevine. In: IntechOpen; 2019.

20. Vezzulli S, Leonardelli L, Malossini U, Stefanini M, Velasco R, Moser C. Pinot blanc and Pinot gris arose as independent somatic mutations of Pinot noir. J Exp Bot. 2012;63(18):6359–69.

21. Pelsy F, Dumas V, Bévilacqua L, Hocquigny S, Merdinoglu D. Chromosome replacement and deletion lead to Clonal Polymorphism of Berry Color in Grapevine. PLoS Genet. 2015;11(4):e1005081.

22. Röckel F, Moock C, Braun U, Schwander F, Cousins P, Maul E, Töpfer R, Hausmann L. Color Intensity of the Red-Fleshed Berry phenotype of Vitis vinifera Teinturier grapes varies due to a 408 bp duplication in the promoter of VvmybA1. Genes. 2020;11(8):891.

23. Boss PK, Thomas MR. Association of dwarfism and floral induction with a grape 'green revolution' mutation. Nature. 2002;416(6883):847–50.

24. Pratt C, Einset J, Zahur M. Radiation damage in Apple shoot Apices. Am J Bot. 1959;46(7):537–44.

25. Jailon. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007;449(7161):463–7.

26. Canaguier A, Grimplet J, Di Gaspero G, Scalabrin S, Duchêne E, Choisne N, Mohellibi N, Guichard C, Rombauts S, Le Clainche I, et al. A new version of the grapevine reference genome assembly (12X.v2) and of its annotation (VCost. v3). Genomics Data. 2017;14:56–62.

27. Di Genova A, Almeida AM, Muñoz-Espinoza C, Vizoso P, Travisany D, Moraga C, Pinto M, Hinrichsen P, Orellana A, Maass A. Whole genome comparison between table and wine grapes reveals a comprehensive catalog of structural variants. BMC Plant Biol. 2014;14(1):7.

28. Chin C-S, Peluso P, Sedlazeck FJ, Nattestad M, Concepcion GT, Clum A, Dunn C, O'Malley R, Figueroa-Balderas R, Morales-Cruz A, et al. Phased diploid genome assembly with single-molecule real-time sequencing. Nat Methods. 2016;13(12):1050–4.

29. Roach MJ, Johnson DL, Bohlmann J, van Vuuren HJJ, Jones SJM, Pretorius IS, Schmidt SA, Borneman AR. Population sequencing reveals clonal diversity

and ancestral inbreeding in the grapevine cultivar chardonnay. PLoS Genet. 2018;14(11):e1007807.

30. Minio A, Massonnet M, Figueroa-Balderas R, Castro A, Cantu D. Diploid genome Assembly of the wine grape Carménère. G3 Genes|Genomes|Genetics. 2019;9(5):1331–7.

31. Girollet N, Rubio B, Lopez-Roques C, Valière S, Ollat N, Bert P-F. Author correction: De novo phased assembly of the Vitis riparia grape genome. Sci Data 2019, 6(1).

32. Massonnet M, Cochetel N, Minio A, Vondras AM, Lin J, Muyle A, Garcia JF, Zhou Y, Delledonne M, Riaz S et al. The genetic basis of sex determination in grapes. Nat Commun 2020, 11(1).

33. Zou C, Massonnet M, Minio A, Patel S, Llaca V, Karn A, Gouker F, Cadle-Davidson L, Reisch B, Fennell A et al. Multiple independent recombinations led to hermaphroditism in grapevine. Proceedings of the National Academy of Sciences 2021, 118(15):e2023548118.

34. Cheng H, Concepcion GT, Feng X, Zhang H, Li H. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nat Methods. 2021;18(2):170–5.

35. Marx V. Long road to long-read assembly. Nat Methods. 2021;18(2):125–9.

36. Sugawara K, Oowada A, Moriguchi T, Omura M. Identification of Citrus chimeras by RAPD markers. HortScience. 1995;30(6):1276–8.

37. Riaz S, Garrison KE, Dangl GS, Boursiquot J-M, Meredith CP. Genetic divergence and chimerism within ancient asexually propagated Winegrape Cultivars. J Am Soc Hortic Sci. 2002;127(4):508–14.

38. Noh J-H, Park K-S, Yun H, Do G-R, Hur Y, Seung Hui K, Lee H-C, Ryou M-S, Park S-J, Jung SM. Determination of Chimera types and Ploidy Level of Sports from 'Campbell Early' grape (Vitis labruscana). Korean J Hortic Sci Technol 2010, 28.

39. Faize M, Faize L, Burgos L. Using quantitative real-time PCR to detect chimeras in transgenic tobacco and apricot and to monitor their dissociation. BMC Biotechnol. 2010;10(1):53.

40. Hocquigny S, Pelsy F, Dumas V, Kindt S, Heloir MC, Merdinoglu D. Diversification within grapevine cultivars goes through chimeric states. Genome. 2004;47(3):579–89.

41. Gambino G, Dal Molin A, Boccacci P, Minio A, Chitarra W, Avanzato CG, Tononi P, Perrone I, Raimondi S, Schneider A et al. Whole-genome sequencing and SNV genotyping of 'Nebbiolo' (Vitis vinifera L.) clones. Sci Rep 2017, 7(1).

42. Hou B-H, Tsai Y-H, Chiang M-H, Tsao S-M, Huang S-H, Chao C-P, Chen H-M. Cultivar-specific markers, mutations, and chimerism of Cavendish banana somaclonal variants resistant to Fusarium oxysporum f. sp. cubense tropical race 4. BMC Genomics 2022, 23(1).

43. Rohlin A, Wernersson J, Engwall Y, Wiklund L, Björk J, Nordling M. Parallel sequencing used in detection of mosaic mutations: comparison with four diagnostic DNA screening techniques. Hum Mutat. 2009;30(6):1012–20.

44. Hardenbol P, Banér J, Jain M, Nilsson M, Namsaraev EA, Karlin-Neumann GA, Fakhrai-Rad H, Ronaghi M, Willis TD, Landegren U, et al. Multiplexed genotyping with sequence-tagged molecular inversion probes. Nat Biotechnol. 2003;21(6):673–8.

45. Biezuner T, Brilon Y, Arye AB, Oron B, Kadam A, Danin A, Furer N, Minden Mark D, Hwan Kim DD, Shapira S et al. An improved molecular inversion probe based targeted sequencing approach for low variant allele frequency. NAR Genomics and Bioinformatics 2022, 4(1).

46. Wang Y, Moorhead M, Karlin-Neumann G, Falkowski M, Chen C, Siddiqui F, Davis RW, Willis TD, Faham M. Allele quantification using molecular inversion probes (MIP). Nucleic Acids Res. 2005;33(21):e183–3.

47. Waalkes A, Smith N, Penewit K, Hempelmann J, Konnick EQ, Hause RJ, Pritchard CC, Salipante SJ. Accurate Pan-Cancer Molecular diagnosis of micro-satellite instability by single-molecule Molecular Inversion Probe capture and high-throughput sequencing. Clin Chem. 2018;64(6):950–8.

48. Andersen EF, Paxton CN, O'Malley DP, Louissaint A Jr, Hornick JL, Griffin GK, Fedoriw Y, Kim YS, Weiss LM, Perkins SL, et al. Genomic analysis of follicular dendritic cell sarcoma by molecular inversion probe array reveals tumor suppressor-driven biology. Mod Pathol. 2017;30(9):1321–34.

49. Lau HY, Palanisamy R, Trau M, Botella JR. Molecular Inversion probe: a New Tool for highly specific detection of Plant Pathogens. PLoS ONE. 2014;9(10):e111182.

50. Wang H, Campbell B, Happ M, McConaughy S, Lorenz A, Amundsen K, Song Q, Pantalone V, Hyten D. Development of molecular inversion probes for soybean progeny genomic selection genotyping. The Plant Genome 2023, 16(1).

51. Boyle EA, O'Roak BJ, Martin BK, Kumar A, Shendure J. MIPgen: optimized modeling and design of molecular inversion probes for targeted resequencing. Bioinformatics. 2014;30(18):2670–2.
52. Almomani R, Marchi M, Sopacua M, Lindsey P, Salvi E, Koning BD, Santoro S, Magri S, Smeets HJM, Martinelli Boneschi F, et al. Evaluation of molecular inversion probe versus TruSeq® custom methods for targeted next-generation sequencing. PLoS ONE. 2020;15(9):e0238467.

53. Wu L, Chu X, Zheng J, Xiao C, Zhang Z, Huang G, Li D, Zhan J, Huang D, Hu P, et al. Targeted capture and sequencing of 1245 SNPs for forensic applications. Forensic Sci International: Genet. 2019;42:227–34.

54. Boursiquot J-M, Lacombe t. Laucou v, Julliard s, Perrin f-x, Lanier n, Legrand d, Meredith c, this p: parentage of Merlot and related winegrape cultivars of southwestern France: discovery of the missing link. Aust J Grape Wine Res. 2009;15(2):144–55.

55. : OIV, Release P. http://www.oiv.int/public/medias/5681/en-communiquédepresse-octobre-2017.pdf, (Accessed 19 September 2018). 2017.

56. IFV INRAE, Montpellier IA. Pl@ntGrape, Catalogue of Vines Cultivated in France, IFV - INRAE - Institut Agro Montpellier, 2009–2022. 2022.

57. Robinson J, Harding J, Vouillamoz J. Wine grapes: a complete guide to 1,368 vine varieties, including their origins and flavours. Penguin UK; 2013.

58. van Leeuwen C, Destrac-Irvine A, Dubernet M, Duchêne E, Gowdy M, Marguerit E, Pieri P, Parker A, de Rességuier L, Ollat N. An update on the impact of climate change in viticulture and potential adaptations. Agronomy. 2019;9(9):S14.

59. Mosé Manni MRB, Mathieu Seppey, Felipe A, Simão, Evgeny M, Zdobnov. BUSCO Update: Novel and Streamlined Workflows along with broader and deeper phylogenetic Coverage for Scoring of Eukaryotic, Prokaryotic, and viral genomes. Mol Biol Evol. 2021;38(10):4647–54.

60. Shumate A, Salzberg SL. Liftoff: accurate mapping of gene annotations. Bioinformatics. 2021;37(12):1639–43.

61. Ryan Poplin P-CC, Alexander D, Schwartz S, Colthurst T, Ku A, Newburger D, Dijamco J, Nguyen N, Afshar PT, Gross SS, Dorfman L, McLean CY. DePristo: DeepVariant: a universal SNP and small-indel variant caller using deep neural networks. Nat Biotechnol. 2018;36:983–7.

62. Velasco R, Zharkikh A, Troggio M, Cartwright DA, Cestaro A, Pruss D, Pindo M, Fitzgerald LM, Vezzulli S, Reid J, et al. A high Quality Draft Consensus sequence of the genome of a heterozygous Grapevine Variety. PLoS ONE. 2007;2(12):e1326.

63. Maestri S, Gambino G, Lopatriello G, Minio A, Perrone I, Cosentino E, Giovan none B, Marcolungo L, Alfano M, Rombauts S et al. 'Nebbiolo' genome assembly allows surveying the occurrence and functional implications of genomic structural variations in grapevines (Vitis vinifera L.). BMC Genomics 2022, 23(1).

64. Burian A. Does shoot apical Meristem function as the germline in safeguard- ing against excess of mutations? Front Plant Sci 2021, 12.

65. Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018;34:3094–100.

66. Poplin R, Chang P-C, Alexander D, Schwartz S, Colthurst T, Ku A, Newburger D, Dijamco J, Nguyen N, Afshar PT, et al. A universal SNP and small-indel variant caller using deep neural networks. Nat Biotechnol. 2018;36(10):983–7.

67. Danecek P, Auton A, Abecasis G, Albers CA, Banks E, Depristo MA, Handsaker RE, Lunter G, Marth GT, Sherry ST, et al. The variant call format and VCFtools. Bioinformatics. 2011;27(15):2156–8.

68. Flynn JM, Hubley R, Goubert C, Rosen J, Clark AG, Feschotte C, Smit AF. RepeatModeler2 for automated genomic discovery of transposable element families. *Proceedings of the National Academy of Sciences* 2020, 117(17):9451–9457.

69. RepeatMasker. http://repeatmasker.org.

70. Robinson JT, Thorvaldsdóttir H, Turner D, Mesirov JP. igv. js: an embeddable JavaScript implementation of the Integrative Genomics Viewer (IGV). bioRxiv 2020.

71. Quinlan ARHI. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2.

72. Hiatt JB, Pritchard CC, Salipante SJ, O’Roak BJ, Shendure J. Single molecule molecular inversion probes for targeted, high-accuracy detection of low-frequency variation. Genome Res. 2013;23(5):843–54.

73. Smith T, Heger A, Sudbery I. UMI-tools: modeling sequencing errors in Unique Molecular Identifiers to improve quantification accuracy. Genome Res. 2017;27(3):491–9.

74. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal. 2011;17(1):10.

75. Garsmeur O, Droc G, Antonise R, Grimwood J, Potier B, Aitken K, Jenkins J, Martin G, Charron C, Hervouet C et al. A mosaic monoploid reference sequence for the highly complex genome of sugarcane. Nat Commun 2018, 9(1).

## Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.