- Research article
- Open Access
Variability among Cucurbitaceae species (melon, cucumber and watermelon) in a genomic region containing a cluster of NBS-LRR genes
BMC Genomics volume 18, Article number: 138 (2017)
Cucurbitaceae species contain a significantly lower number of genes coding for proteins with similarity to plant resistance genes belonging to the NBS-LRR family than other plant species of similar genome size. A large proportion of these genes are organized in clusters that appear to be hotspots of variability. The genomes of the Cucurbitaceae species measured until now are intermediate in size (between 350 and 450 Mb) and they apparently have not undergone any genome duplications beside those at the origin of eudicots. The cluster containing the largest number of NBS-LRR genes has previously been analyzed in melon and related species and showed a high degree of interspecific and intraspecific variability. It was of interest to study whether similar behavior occurred in other cluster of the same family of genes.
The cluster of NBS-LRR genes located in melon chromosome 9 was analyzed and compared with the syntenic regions in other cucurbit genomes. This is the second cluster in number within this species and it contains nine sequences with a NBS-LRR annotation including two genes, Fom1 and Prv, providing resistance against Fusarium and Ppapaya ring-spot virus (PRSV). The variability within the melon species appears to consist essentially of single nucleotide polymorphisms. Clusters of similar genes are present in the syntenic regions of the two species of Cucurbitaceae that were sequenced, cucumber and watermelon. Most of the genes in the syntenic clusters can be aligned between species and a hypothesis of generation of the cluster is proposed. The number of genes in the watermelon cluster is similar to that in melon while a higher number of genes (12) is present in cucumber, a species with a smaller genome than melon. After comparing genome resequencing data of 115 cucumber varieties, deletion of a group of genes is observed in a group of varieties of Indian origin.
Clusters of genes coding for NBS-LRR proteins in cucurbits appear to have specific variability in different regions of the genome and between different species. This observation is in favour of considering that the adaptation of plant species to changing environments is based upon the variability that may occur at any location in the genome and that has been produced by specific mechanisms of sequence variation acting on plant genomes. This information could be useful both to understand the evolution of species and for plant breeding.
The evolution of genes related to pathogen resistance in plants has been the object of intense research for different reasons. These genes are probably related to plant adaptation to different environments where they are confronted with evolving pathogens. These are genes that have also been under selective pressure on domestication, as pathogen resistance is one of the most important traits for plant breeders in many crops. Among these genes, those coding for the NBS-LRR class of proteins have attracted a significant amount of research because they have been associated with effector-triggered immunity, an important component of plant resistance to pathogens, as has recently been reviewed . They are also interesting examples of the evolution of plant gene sequences: the mechanisms that produce genome variability, such as single nucleotide polymorphisms and copy number variation that may occur even without large genome duplication events, appear to act upon these sequences in those species that have been studied.
With the completion of genome sequences of major plant species, it appeared that the structure of gene families with similarity to resistance genes were very variable when comparing different species. In a recent study in angiosperms, the evolution of the lineages of different classes of NBS-LRR genes indicates the importance of whole genome duplication events . The number of genes coding for these proteins varies considerably between species and it has been shown that, in some species, these genes occur in clusters that may be the result of amplification of the gene families. This is the case, for instance, in well studied species such as rice .
The case of Cucurbitaceae has attracted attention in this respect as they have been shown to contain a reduced number of sequences belonging to the family of NBS-LRR proteins compared to other plants. This has been observed in the main Cucurbit species sequenced so far, such as cucumber , melon  and watermelon , where the numbers (104, 89 and 54, respectively), have been found to be significantly lower when compared to species having a similar genome size, such as those belonging to the Prunus family . In peach, for instance, with a smaller genome than the three cucurbit species studied, more than 400 NBS-LRR sequences have been identified . The report of a correlation between a high number of NBS-LRR genes and the existence of miRNAs controlling their expression  is a further indication of the importance of controlling the activity of these genes for the fitness of the species. These observations indicate how variability in the structure and expression of NBS-LRR genes may be involved in the adaptation of plants to different environments. It has also been shown that, at least in melon, a group of sequences coding for these proteins is a hotspot of variability . As an example, a detailed study of a large cluster of NBS-LRR genes in chromosome 5 of melon has shown a high interspecific and intraspecific variability in the number of genes at this location . It was of interest to study whether this observation could be generalized to homologous gene families and to other Cucurbit species. In melon, a second large cluster of NBS-LRR has been observed in chromosome 9, with two major resistance genes (Fom1 and Prv) that have been characterized . The analysis of the variability of this cluster in melon and in the related Cucurbit species was the object of the present study.
Results and discussion
It has been observed in melon, as in other cases [12–14], that many of the genes coding for NBS-LRR sequences are located in clusters within the genome . These clusters in melon are among those with the highest variability in the genome in terms of presence/absence of genes . As a possible mechanism for adaptation of plants, this has prompted our interest in analysing these clusters in detail in melon and in the syntenic regions of the two other Cucurbit species whose genome sequences are available (cucumber and watermelon).
The largest of the NBS-LRR gene clusters in the melon genome is in chromosome 5 and includes the gene Vat, responsible for aphid resistance . The genome of this region in melon has been sequenced and the variability of the genes in this cluster has been compared with other cucurbits . Another cluster of sequences coding for NBS-LRR proteins have been reported in chromosome 9 . There is particular interest in this cluster as it contains two genes, Fom1 and Prv, that have been shown to be responsible for resistance to Fusarium and Papaya ring-spot virus, respectively. These genes have been identified from a molecular and genetic point of view and have been found to have an unusual head-to-tail structure . It was therefore of interest to study whether the level of sequence variability in this cluster was similar to that observed in chromosome 5, both between cucurbit species and within these species. In order to answer this question, we carried out a bioinformatic comparative analysis of this region in the different genomic sequences available in melon. This analysis was extended to the syntenic regions of cucumber and watermelon.
Figure 1 shows the structure of the cluster of chromosome 9 and it is compared to the syntenic regions in cucumber and watermelon. Genes with similarity to NBS-LRR proteins (using RGH nomenclature and the numbering in the melon genome) and the other annotated genes are schematically indicated. This sequence corresponds to that reported for part of the cluster previously described  but extended with data from the published melon genome sequences . The two sequences correspond to Piel de Sapo melon varieties as whole genome sequencing was with the DHL92 double haploid line of a cross between a Piel de Sapo and a Korean variety, and some regions correspond to each of the parental lines. The region analyzed corresponds to a Piel de Sapo fragment . It appears that the cluster has nine TIR-NBS-LRR genes, which is less than the cluster in chromosome 5 previously studied, that contains more than 20 genes, depending on the melon variety. A number of transposable elements have been observed between the genes and within an intron of the Prv gene . In Fig. 1, the melon sequence is aligned with the corresponding syntenic loci of cucumber and watermelon.
When comparing the syntenic areas in melon and cucumber, the length of the genome region of cucumber is shorter, a fact that may correlate with the difference in genome size of the two species: 450 Mb in the case of melon and 350 in the case of cucumber. Nevertheless, the number of genes in the cluster, surprisingly, is much greater in cucumber according to the annotation in the reference genomes published. In the syntenic region located in chromosome 5 of this species, at least 12 of the sequences annotated as NBS-LRR in the reference cucumber genome  can be observed. In the case of watermelon, the number of sequences homologous to the melon genes is similar to that in melon but lower than in cucumber. The density of genes in this watermelon region appears to be lower than in melon. However that may be an effect of the quality of the available watermelon sequence, as, for instance, a large insertion between two NBS-LRR genes could be an artefact of the sequencing methods employed.
The NBS-LRR sequences of the cluster in chromosome 9 from melon were compared with the corresponding sequences in cucumber and watermelon and a phylogenetic analysis of the sequences between the three species and within each species was carried out. The tree shown in Fig. 2 indicates a possible pathway of how the sequences could be generated in the three species as a result of this analysis. The alignment was compared with that published on the PhylomeDB  and Plaza websites . The alignments essentially agree: the exception is in the Plaza system of genes homologous to Fom1 and Prv in watermelon, shown in a different colour in the figure. However the two genes in melon and cucumber are both in the opposite orientation, indicating that they were generated before divergence of the two species. In the other cases, complex patterns of relations, including genes being expressed in opposite directions, were observed. It should be noted that a complete transposon of the CUMULE family  was present near Fom1, and this may indicate an independent origin of this sequence under the influence of the transposon. Amplification of the gene clusters clearly occurs independently in the different species, in particular in cucumber, as will be discussed later.
A number of genome sequences are currently available for melon varieties, and they have been used for analysis of intraspecific variability within the species. They have been particularly useful for studying variations in the number of genes present in clusters . In a similar analysis, the reference genome sequence was compared. Single nucleotide polymorphisms, and copy number were recorded along chromosome 9 (Fig. 3). These polymorphisms were unevenly distributed along the chromosome, with an average value around 200 polymorphisms per 10 kb. The value of the genome interval containing the NBS-LRR recorded was very similar. No change in the structure of the cluster was observed in the melon varieties analysed, in contrast to the significant number of presence/absence of genes observed in the cluster of NBS-LRR genes in chromosome 5 . A specific difference was observed between melon varieties in the gene MGH9 that corresponded to Fom1. The cluster of NBS-LRR sequences in chromosome 9 of melon contained two gene sequences that have been shown to be responsible for two genes providing resistance to pathogens, Fom1 and Prv and it has been reported that they have an unusual structure . When analysing the sequence of these two genes (shown in Fig. 1), in the published genome sequence of the DHL92 double haploid line, a unique insertion was found that corresponds to the cultivated Piel de Sapo variety. The insertion was in one of the introns of the gene. Searching in the REPET annotation , it was found to contain a complex pattern of several incomplete transposable elements.
A significant difference was observed in the syntenic locus in chromosome 5 of the genome of cucumber. Between the sequences with similarity to those in melon, a number of new sequences belonging to the same family of genes were observed. In this region, 12 sequences with similarity to NBS-LRR genes were found, while only nine were present in melon. Upon further examination, it appeared that the new sequences also coded for NBS-LRR proteins, but were more similar to each other than to the melon genes. These sequences were interspersed with the previous gene sequences (see Fig. 1). A possible hypothesis on how these sequences were generated is based on phylogenetic analysis (Fig. 2). Specific amplification of the cluster in this locus appears to have occurred in cucumber, which has a smaller genome than melon. The generation of new NBS-LRR sequences could have occurred in the different species at specific loci after the separation of the cucurbits during the speciation process. It was of interest to study whether there was variability in the sequences in this cluster in chromosome 5 of cucumber, by comparing the sequences of the two genomes published so far. This result is given in Additional file 1: Figure S1. The two genomic sequences were compared to each other and to the structure of the locus proposed by Lin et al., . The fragmented structure of the genes from the published cucumber genomes made it difficult to accurately compare the two sequences, but no major rearrangement were observed between them.
Nevertheless, the study was extended to the other cucumber genomic sequences that have been published . This is a collection of 115 sequences and they have been compared with the reference cucumber genome. It appears that, in a subset of accessions corresponding to Indian harwickii cucumber varieties, a deletion can be observed indicating that presence-absence mechanisms are also active in cucumber, as has been observed in melon. Varieties where PAV polymorphisms are observed are presented in Fig. 4 and the full set of varieties are listed in Additional file 2: File S3. It is interesting to note that four of these PAV polymorphisms include tandem Fom1 and Prv orthologous genes, thus indicating the relevance of genomic variability of these clusters in establishing resistances to pathogens.
In the Cucurbitaceae family, the number of genes coding for the NBS-LRR group of proteins, which are generally related to disease resistance, is consistently lower than in other plant families, and a large proportion of these genes are present in clusters. It has previously been shown that the largest cluster in melon on chromosome 5 is one of the most variable regions in the genome of this species, in terms of presence and absence of genes . Examination of a smaller gene cluster in chromosome 9 showed that the variability in terms of single nucleotide polymorphisms was within the average range observed in the whole chromosome. Apparently there is no major change in the number of genes in the cluster or in the structure of the cluster itself, at least in those melon varieties whose genome has been resequenced to date .
In the syntenic region of the two other cucurbit species, cucumber and watermelon, whose genomes have already been sequenced, a cluster of NBS-LRR gene sequences have also been found. The most interesting observation in this cluster was that, while there were no changes in the number of genes in melon, the number of NBS-LRR genes in cucumber was found to be amplified when the locus was compared with the syntenic one. In thesespecies, the cluster was larger than in melon, in spite of the smaller size of the genome. In cucumber, a number of genes conferring resistance to pathogens have been found in chromosome 5, although not in the region of the cluster except a QTL of pathogen resistance to powdery mildew located in this region . Therefore it is probable that there are genes active in providing pathogen resistance in this cluster. Unfortunately the quality of the published sequences of the cucumber genome does not allow the variability in this region to be studied. However when a large number of cucumber genomes is analyzed the absence of a number of genes in the cluster is observed in some varieties from Indian origin. These varieties are included in the hardwickii subgroup of cucumber varieties that clearly separates from other varieties after the analysis of their genomes (see Supplementary Figure 6 in ref. 36). It has to be taken into account that the genetic bases of cultivated cucumber is relatively low and therefore the level of variability is expected to be lower than in other species, including melon.
It may be concluded that clusters of genes coding for proteins having the NBS-LRR features in cucurbits have specific strategies of variability in different regions and between species. It has to be remembered that these are species that have not undergone recent whole genome duplications and that contain a low number of sequences with similarity to the family of resistance genes when compared with other plants. Amplification of NBS-LRR gene numbers have occurred in one locus in chromosome 9 in the case of melon and in another one in the syntenic chromosome 5 in the case of cucumber. This leads to considering that the adaptation of plant species to changing environments may be based upon the variability that may occur at any location in the genome and that has been produced by any of the mechanisms of sequence variation acting on plant genomes.
The protein sequences (MELO3C022143 - MELO3C022157), described as a cluster of resistance genes in chromosome 9 of melon , were manually inspected and re-annotated with Augustus  and FGENESH , with Arabidopsis as the model organism. For each of the 9 resistance genes in the cluster, a protein sequence was decided based on the output of the re-annotation, the protein sequences described in BAC sequencing studies of the region , and taking into account the analysis and criteria (exon structure, domain composition, sequence integrity, etc.) of a study of resistance genes in several Cucurbitacea . No additional resistance genes in the cluster were found after scanning the DRAGO database  and running a chromosome-wide prediction of NBS, LRR and TIR domains with hmmscan  against Pfam profiles database .
For watermelon and cucumber (Gy14 and Chinese Long 9930) protein fasta sequences and annotation gff3 were collected from ICUGI  and Phytozome . All were downloaded on January 15th, 2015. Cucumber scaffolds where assembled in linkage groups based on a previously published, high resolution genetic map .
The phylogeny of the resistance genes for each species was analysed with the NBS domain sequence (Pfam ID: PF00931) predicted with hmmscan and the “One Click” option of the phylogeny.fr web server, based on MUSCLE alignments, PhyML phylogenic analysis and TreeDyn tree rendering .
The synthenic regions between each pair of species were determined with MCScanX , as described elsewhere . The orthologous relationships between each pair of species was calculated with InParanoid4.1  with default parameters, bootstraping and no outgroup analysis. Further parsing and filtering was done with in-house scripts. Blastp was used for double-check prediction for watermelon orthologous pairs. Melon-cucumber and melon-watermelon orthologous pairs were verified with the PhylomeDB 4  and Plaza 3.0 database , respectively.
For annotation of transposon elements in melon, the REPET package v2.2  was used. The SNPs and indels for seven resequences of melon with respect the reference genome  have been previously published . The insertion at the first intron in MELO3C022145 gene was confirmed to be unique with respect to previously published sequences [11, 16].
In order to check for presence/absence variation in cucumber, we retrieved 115 published resequences  of several varieties, publicly available at ENA  (Study: PRJNA171718). The fastq files were trimmed and filtered with skewer  and aligned to the Gy14 reference genome with bwa mem . Large deletions were detected with Delly . PAVs were retained when ten or more paired end reads with a minimum mapping quality of one supported the deletion.
Jones JDG, Vance RE, Dangl JL. Intracellular innate immune surveillance devices in plants and animals. Sci. 2016;354(6316):1117.
Shao ZQ, Xue JY, Wu P, Zhang YM, Wu Y, Hang YY, Wang B, Chen JQ. Large-scale analyses of angiosperm nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes reveal three anciently diverged classes with distinct evolutionary patterns. Plant Physiol. 2016;170(4):2095–109.
Zhou T, Wang Y, Chen JQ, Araki H, Jing Z, Jiang K, Shen J, Tian D. Genome-wide identification of NBS genes in japonica rice reveals significant expansion of divergent non-TIR NBS-LRR genes. Mol Genet Genomics. 2004;271(4):402–15.
Huang S, Li R, Zhang Z, Li L, Gu X, Fan W, Lucas WJ, Wang X, Xie B, Ni P, Ren Y, et al. The genome of the cucumber, Cucumis sativus L. Nat Genet. 2009;41:1275–81.
Garcia-Mas J, Benjak A, Sanseverino W, Bourgeois M, Mir G, González VM, Hénaff E, Câmara F, Cozzuto L, Lowy E, Alioto T, Capella-Guitérrez S, Blanca J, Cañizares J, Ziarsolo P, Gonzalez-Ibeas D, Rodríguez-Moreno L, Droege M, Du L, Alvarez-Tejado M, Lorente-Galdós B, Melé M, et al. The genome of melon (Cucumis melo L.). PNAS. 2012;109(29):11872–7.
Guo S, Zhang J, Sun H, Salse J, Lucas WJ, Zhang H, Zheng Y, Mao L, Ren Y, Wang Z, et al. The draft genome of watermelon (Citrullus lanatus) and resequencing of 20 diverse accessions. Nat Genet. 2013;45:51–8.
Jia X, Yuan Y, Zhang Y, Yang S, Zhang X. Extreme expansion of NBS-encoding genes in Rosaceae. BMC Genet. 2015;16:48.
González VM, Müller S, Baulcombe D, Puigdomenech P. Evolution of NBS-LRR gene copies among Dicot plants and its regulation by members of the miR482/2119 superfamily of miRNAs. Mol Plant. 2015;8:329–31.
González VM, Aventín N, Centeno E, Puigdomènech P. High presence/absence gene variability in defense-related gene clusters of Cucumis melo. BMC Genomics. 2013;14:782–95.
González VM, Aventín N, Centeno E, Puigdomènech P. Interspecific and intraspecific gene variability in a 1-Mb region containing the highest density of NBS-LRR genes found in the melon genome. BMC Genomics. 2014;15:1131.
Brotman Y, Normantovich M, Goldenberg Z, Zvirin Z, Kovalski I, Stovbun N, Doniger T, Bolger AM, Troadec C, Bendahmane A, Cohen R, Katzir N, Pitrat M, Dogimont C, Perl-Treves R. Dual resistance of melon to Fusarium oxysporum races 0 and 2 and to Papaya ring-spot virus is controlled by a pair of head-to-head-oriented NB-LRR genes of unusual architecture. Mol Plant. 2013;6(1):235–8.
Guo YL, Fitz J, Schneeberger K, Ossowski S, Cao J, Weigel D. Genome-wide comparison of nucleotide-binding site-leucine-rich repeat-encoding genes in Arabidopsis. Plant Physiol. 2011;157(2):757–69.
Jupe F, Pritchard L, Etherington GJ, Mackenzie K, Cock PJ, Wright F, Sharma SK, Bolser D, Bryan GJ, Jones JD, Hein I. Identification and localisation of the NB-LRR gene family within the potato genome. BMC Genomics. 2012;13:75.
Christopoulou M, Wo SR, Kozik A, McHale LK, Truco MJ, Wroblewski T, Michelmore RW. Genome-wide architecture of disease resistance genes in lettuce. G3 (Bethesda). 2015;5(12):2655–69.
Dogimont C, Chovelon V, Pauquet J, Boualem A, Bendahmane A. The Vat locus encodes for a CC-NBS-LRR protein that confers resistance to Aphis gossypii infestation and A. gossypii-mediated virus resistance. Plant J. 2014;80(6):993–1004.
van Leeuwen H, Garcia-Mas J, Coca M, Puigdoménech P, Monfort A. Analysis of the melon genome in regions encompassing TIR-NBS-LRR resistance genes. Mol Genet Genomics. 2005;273(3):240–51.
Huerta-Cepas J, Capella-Gutiérrez S, Pryszcz LP, Marcet-Houben M, Gabaldón T. PhylomeDB v4: zooming into the plurality of evolutionary histories of a genome. Nucleic Acids Res. 2014;42:D897–902.
Proost S, Van Bel M, Vaneecchoutte D, Van de Peer Y, Inze D, Mueller-Roeber B, Vandepoele K. PLAZA 3.0: an access point for plant comparative genomics. Nucleic Acids Res. 2015;43:D974–81.
van Leeuwen H, Monfort A, Puigdomenech P. Mutator-like elements identified in melon, Arabidopsis and rice contain ULP1 protease domains. Mol Genet Genomics. 2007;277(4):357–64.
Flutre T, Duprat E, Feuillet C, Quesneville H. Considering transposable element diversification in de novo annotation approaches. PLoS One. 2011;6(1):e16526.
Lin X, Zhang Y, Kuang H, Chen J. Frequent loss of lineages and deficient duplications accounted for low copy number of disease resistance genes in Cucurbitaceae. BMC Genomics. 2013;14:335.
European Nucleotide Archive. http://www.ebi.ac.uk/ena/.
Sanseverino W, Hénaff E, Vives C, Pinosio S, Burgos-Paz W, Morgante M, Ramos-Onsins SE, Garcia-Mas J, Casacuberta JM. Transposon insertions, structural variations, and SNPs contribute to the evolution of the melon genome. Mol Biol Evol. 2015;32(10):2760–74.
He X, Li X, Pandey S, Yandell BS, Pathak M, Weng Y. QTL mapping of powdery mildew resistance in WI 2757 cucumber (Cucumis sativus L.). Theor Appl Genet. 2013;126:2149–61.
Stanke M, Morgenstern B. AUGUSTUS: a web server for gene prediction in eukaryotes that allows user-defined constraints. Nucl Acid Res. 2005;33:W465–7.
Solovyev V, Kosarev P, Seledsov I, Vorobyev D. Automatic annotation of eukaryotic genes, pseudogenes and promoters. Genome Biol. 2006;7(1):10. 1–12.
Sanseverino W, Roma G, De Simone M, Faino L, Melito S, Stupka E, Frusciante L, Ercolano MR. PRGdb: a bioinformatics platform for plant resistance gene analysis. Nucleic Acids Res. 2010;38:D814–21.
Finn RD, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching. Nucl Acid Res. 2011;39:W29–37.
Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, Sonnhammer ELL, Tate J, Punta M. The Pfam protein families database. Nucleic Acids Res. 2014;42:D222–30.
Cucurbit Genomes Database. http://icugi.org. Accessed 15 Jan 2015.
Phytozome 11, The Plant Genomics Resource. http://phytozome.jgi.doe.gov. Accessed 15 Jan 2015.
Luming Y, Dawei L, Yuhong L, Xingfang G, Sanwen H, Jordi G-M, Yiqun W. A 1,681-locus consensus genetic map of cultivated cucumber including 67 NB-LRR resistance gene homolog and ten gene loci. BMC Plant Biol. 2013;13:53.
Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, Dufayard JF, Guindon S, Lefort V, Lescot M, Claverie JM, Gascuel O. Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucl Acid Res. 2008;36:W465–9.
Wang Y, Tang H, DeBarry JD, Tan X, Li J, Wang X, Lee TH, Jin H, Marler B, Guo H, Kissinger JC, Paterson AH. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucl Acid Res. 2012;40(7):e49.
O’Brien KP, Remm M, Sonnhammer EL. InParanoid: a comprehensive database of eukaryotic orthologs. Nucleic Acids Res. 2005;33:D476–80.
Qi J, Liu X, Shen D, Miao H, Xie B, Li X, Zeng P, Wang S, Shang Y, Gu X, Du Y, Li Y, Lin T, Yuan J, Yang X, Chen J, Chen H, Xiong X, Huang K, Fei Z, Mao L, Tian L, Städler T, Renner SS, Kamoun S, Lucas WJ, Zhang Z, Huang S. A genomic variation map provides insights into the genetic basis of cucumber domestication and diversity. Nat Genet. 2013;45(12):1510–5. Accessed 4 Nov 2016.
Jiang H, Lei R, Ding SW, Zhu S. Skewer: a fast and accurate adapter trimmer for next-generation sequencing paired-end reads. BMC Bioinformatics. 2014;15:182.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler Transform. Bioinformatics. 2009;25:1754–60.
Rausch T, Zichner T, Schlattl A, Stuetz AM, Benes V, Korbel JO. Delly: structural variant discovery by integrated paired-end and split-read analysis. Bioinformatics. 2012;28:i333–9.
The authors thank doctors Jordi Garcia-Mas and Josep Casacuberta (CRAG) for critical reading of the manuscript.
We acknowledge financial support from the Spanish Ministry of Economy and Competitiveness, through the “Severo Ochoa Programme for Centres of Excellence in R&D” 2016–2019 (SEV-2015-0533) and project MICINN AGL2013-43244 to PP.
Availability of data and materials
Information on the resistance genes for melon, watermelon and cucumber used in this study can be found in the supplementary material (Additional file 3).
PP conceived the study, supervised the research and wrote the manuscript. JM analysed the bioinformatic data and participated in drafting the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Distribution of the resistance gene cluster in chromosome 5 in the published genomes of cucumber: 9930 (upper) and Gy14 (middle). Dashed lines represent identical regions. Differences are due basically to Ns stretches. Partial and pseudogene annotation (bottom) according to . (PDF 70 kb)
Detailed list of PAV polymorphisms in cucumber chromosome 5 cluster of NBS-LRR genes. Orange boxes indicate PAV polymorphisms. (XLS 45 kb)
List of the resistance genes in melon, cucumber and watermelon used in this study, along with the gene code, resistance gene class and position in the genome. Protein fasta sequences are available in Additional File 4. Melon gene nomenclature has three codes: gene code , resistance gene code [11, 16] and gene name . These genes are shown in Fig. 2. LG = Linkage Group. (XLS 11 kb)
About this article
Cite this article
Morata, J., Puigdomènech, P. Variability among Cucurbitaceae species (melon, cucumber and watermelon) in a genomic region containing a cluster of NBS-LRR genes. BMC Genomics 18, 138 (2017). https://doi.org/10.1186/s12864-017-3529-5
- Syntenic Region
- Cucumber Genome
- Cucurbit Species
- Melon Genome