Systematic Analysis and Comparison of Nucleotide-Binding Site Disease Resistance Genes in a Diploid Cotton Gossypium raimondii

Plant disease resistance genes are a key component of defending plants from a range of pathogens. The majority of these resistance genes belong to the super-family that harbors a Nucleotide-binding site (NBS). A number of studies have focused on NBS-encoding genes in disease resistant breeding programs for diverse plants. However, little information has been reported with an emphasis on systematic analysis and comparison of NBS-encoding genes in cotton. To fill this gap of knowledge, in this study, we identified and investigated the NBS-encoding resistance genes in cotton using the whole genome sequence information of Gossypium raimondii. Totally, 355 NBS-encoding resistance genes were identified. Analyses of the conserved motifs and structural diversity showed that the most two distinct features for these genes are the high proportion of non-regular NBS genes and the high diversity of N-termini domains. Analyses of the physical locations and duplications of NBS-encoding genes showed that gene duplication of disease resistance genes could play an important role in cotton by leading to an increase in the functional diversity of the cotton NBS-encoding genes. Analyses of phylogenetic comparisons indicated that, in cotton, the NBS-encoding genes with TIR domain not only have their own evolution pattern different from those of genes without TIR domain, but also have their own species-specific pattern that differs from those of TIR genes in other plants. Analyses of the correlation between disease resistance QTL and NBS-encoding resistance genes showed that there could be more than half of the disease resistance QTL associated to the NBS-encoding genes in cotton, which agrees with previous studies establishing that more than half of plant resistance genes are NBS-encoding genes.     

Genome-wide identification of NBS genes in<Emphasis Type="Italic"> japonica</Emphasis> rice reveals significant expansion of divergent non-TIR NBS-LRR genes

A complete set of candidate disease resistance ( R) genes encoding nucleotide-binding sites (NBSs) was identified in the genome sequence of japonica rice ( Oryza sativa L. var. Nipponbare). These putative R genes were characterized with respect to structural diversity, phylogenetic relationships and chromosomal distribution, and compared with those in Arabidopsis thaliana. We found 535 NBS-coding sequences, including 480 non-TIR (Toll/IL-1 receptor) NBS-LRR (Leucine Rich Repeat) genes. TIR NBS-LRR genes, which are common in A. thaliana, have not been identified in the rice genome. The number of non-TIR NBS-LRR genes in rice is 8.7 times higher than that in A. thaliana, and they account for about 1% of all of predicted ORFs in the rice genome. Some 76% of the NBS genes were located in 44 gene clusters or in 57 tandem arrays, and 16 apparent gene duplications were detected in these regions. Phylogenetic analyses based both NBS and N-terminal regions classified the genes into about 200 groups, but no deep clades were detected, in contrast to the two distinct clusters found in A. thaliana. The structural and genetic diversity that exists among NBS-LRR proteins in rice is remarkable, and suggests that diversifying selection has played an important role in the evolution of R genes in this agronomically important species. (Supplemental material is available online at .)Communicated by R. HagemannThe first three authors contributed equally to this work     

玉米与其他六种植物NBS类抗病基因的进化比较分析

具有核苷酸结合位点(nucleotide binding site,NBS)的抗病基因在植物抵抗各种病原菌侵染中起关键作用。对玉米全基因组中具有NBS结构的基因进行鉴定和分析,并结合水稻、高粱、拟南芥、百脉根、苜蓿和杨树的NBS类基因比较其在数量、复制、染色体定位和亲缘关系上的进化差异。发现玉米NBS类基因数量、复制数和成簇基因数均明显少于其他植物。低复制频率可能导致玉米NBS类基因较少,并推测可能导致其功能具有多样性。在基因染色体定位上,除高梁外,玉米与其他五种植物相似,呈不均衡分布。此外,进化树分析表明玉米NBS类基因与高粱的亲缘关系最近,与拟南芥的最远,在物种间表现出较高的保守性。结果对掲示玉米NBS基因的进化特点与发掘有益的NBS类抗病基因提供了重要的理论依据。     

A large scale analysis of resistance gene homologues in<Emphasis Type="Italic"> Arachis</Emphasis>

Arachis hypogaea L., commonly known as the peanut or groundnut, is an important and widespread food legume. Because the crop has a narrow genetic base, genetic diversity in A. hypogaea is low and it lacks sources of resistance to many pests and diseases. In contrast, wild diploid Arachis species are genetically diverse and are rich sources of disease resistance genes. The majority of known plant disease resistance genes encode proteins with a nucleotide binding site domain (NBS). In this study, degenerate PCR primers designed to bind to DNA regions encoding conserved motifs within this domain were used to amplify NBS-encoding regions from Arachis spp. The Arachis spp. used were A. hypogaea var. Tatu and wild species that are known to be sources of disease resistance: A. cardenasii, A. duranensis , A. stenosperma and A. simpsonii. A total of 78 complete NBS-encoding regions were isolated, of which 63 had uninterrupted ORFs. Phylogenetic analysis of the Arachis NBS sequences derived in this study and other NBS sequences from Arabidopsis thaliana, Medicago trunculata , Glycine max , Lotus japonicus and Phaseolus vulgaris that are available in public databases This analysis indicates that most Arachis NBS sequences fall within legume-specific clades, some of which appear to have undergone extensive copy number expansions in the legumes. In addition, NBS motifs from A. thaliana and legumes were characterized. Differences in the TIR and non-TIR motifs were identified. The likely effect of these differences on the amplification of NBS-encoding sequences by PCR is discussed.Electronic Supplementary Material Supplementary material is available for this article if you access the article at . A link in the frame on the left on that page takes you directly to the supplementary material.Communicated by M.-A. Grandbastien     

Plant disease resistance genes encode members of an ancient and diverse protein family within the nucleotide-binding superfamily

The nucleotide binding site (NBS) is a characteristic domain of many plant resistance gene products. An increasing number of NBS-encoding sequences are being identified through gene cloning, PCR amplification with degenerate primers, and genome sequencing projects. The NBS domain was analyzed from 14 known plant resistance genes and more than 400 homologs, representing 26 genera of monocotyledonous, dicotyle-donous and one coniferous species. Two distinct groups of diverse sequences were identified, indicating divergence during evolution and an ancient origin for these sequences. One group was comprised of sequences encoding an N-terminal domain with Toll/Interleukin-1 receptor homology (TIR), including the known resistance genes, N, M, L6, RPP1 and RPP5. Surprisingly, this group was entirely absent from monocot species in searches of both random genomic sequences and large collections of ESTs. A second group contained monocot and dicot sequences, including the known resistance genes, RPS2, RPM1, I2, Mi, Dm3, Pi-B, Xa1, RPP8, RPS5 and Prf. Amino acid signatures in the conserved motifs comprising the NBS domain clearly distinguished these two groups. The Arabidopsis genome is estimated to contain approximately 200 genes that encode related NBS motifs; TIR sequences were more abundant and outnumber non-TIR sequences threefold. The Arabidopsis NBS sequences currently in the databases are located in approximately 21 genomic clusters and 14 isolated loci. NBS-encoding sequences may be more prevalent in rice. The wide distribution of these sequences in the plant kingdom and their prevalence in the Arabidopsis and rice genomes indicate that they are ancient, diverse and common in plants. Sequence inferences suggest that these genes encode a novel class of nucleotide-binding proteins.     

Systematic analysis and comparison of nucleotide-binding site disease resistance genes in maize

Nucleotide-binding site (NBS) disease resistance genes play an integral role in defending plants from a range of pathogens and insect pests. Consequently, a number of recent studies have focused on NBS-encoding genes in molecular disease resistance breeding programmes for several important plant species. Little information, however, has been reported with an emphasis on systematic analysis and a comparison of NBS-encoding genes in maize. In the present study, 109 NBS-encoding genes were identified based on the complete genome sequence of maize (Zea mays cv. B73), classified as four different subgroups, and then characterized according to chromosomal locations, gene duplications, structural diversity and conserved protein motifs. Subsequent phylogenetic comparisons indicated that several maize NBS-encoding genes possessed high similarity to function-known NBS-encoding genes, and revealed the evolutionary relationships of NBS-encoding genes in maize comparede to those in other model plants. Analyses of the physical locations and duplications of NBS-encoding genes showed that gene duplication events of disease resistance genes were lower in maize than in other model plants, which may have led to an increase in the functional diversity of the maize NBS-encoding genes. Various expression patterns of maize NBS-encoding genes in different tissues were observed using an expressed-sequence tags database and, alternatively, after southern leaf blight infection or the application of exogenous salicylic acid. The results reported in the present study contribute to an improved understanding of the NBS-encoding gene family in maize.     

Large-scale analysis of NBS domain-encoding resistance gene analogs in Triticeae

Proteins containing nucleotide binding sites (NBS) encoded by plant resistance genes play an important role in the response of plants to a wide array of pathogens. In this paper, an in silico search was conducted in order to identify and characterize members of NBS-encoding gene family in the tribe of Triticeae. A final dataset of 199 sequences was obtained by four search methods. Motif analysis confirmed the general structural organization of the NBS domain in cereals, characterized by the presence of the six commonly conserved motifs: P-loop, RNBS-A, Kinase-2, Kinase-3a, RNBS-C and GLPL. We revealed the existence of 11 distinct distribution patterns of these motifs along the NBS domain. Four additional conserved motifs were shown to be significantly present in all 199 sequences. Phylogenetic analyses, based on genetic distance and parsimony, revealed a significant overlap between Triticeae sequences and Coiled coil-Nucleotide binding site-Leucine rich repeat (CNL)-type functional genes from monocotyledons. Furthermore, several Triticeae sequences belonged to clades containing functional homologs from non Triticeae species, which has allowed for these sequences to be functionally assigned. The findings reported, in this study, will provide a strong groundwork for the isolation of candidate R-genes in Triticeae crops and the understanding of their evolution.     

Genetic diversity of NBS–LRR class disease-resistance gene analogs in cultivated and wild eggplants

The genes encoding the nucleotide-binding site (NBS) and leucine-rich repeat (LRR) motifs constitute a large gene family in plants and have attracted much interest, because most of the plant disease-resistance genes that have been cloned are from this gene family. In this study, degenerate oligonucleotide primers, designed on the basis of conserved regions of the NBS domains from known plant resistance genes, were used to isolate resistance gene analogs (RGAs) from cultivated and wild eggplants, i.e., S. melongena, S. aethiopicum gr. Gilo, S. linnaeanum, S. integrifolium, S. sisymbriifolium, and S. khasianum. Sequence analysis indicated that the cloned eggplant RGAs belong to the non-TIR–NBS–LRR type, which are very similar to the R genes or the RGAs identified in other plant species, especially Solanaceae plants, suggesting the existence of common ancestors. Wide genetic diversity of eggplant RGAs was observed both in interspecific and intraspecific sequences, and eight distinct families of eggplant RGAs were identified. Further studies revealed a high average ratio of synonymous to non-synonymous substitution and a low level of recombination. These results suggest that NBS-encoding sequences of RGAs in cultivated and wild eggplants are subject to gradual accumulation of mutations leading to purifying selection. This is the first report of NBS–LRR class RGAs in eggplants.     

Genome-Wide Analysis in Wild and Cultivated <Emphasis Type="Italic">Oryza</Emphasis> Species Reveals Abundance of NBS Genes in Progenitors of Cultivated Rice

NBS-encoding genes play a critical role in the plant defense system. Wild relatives of crop plants are rich reservoirs of plant defense genes. Here, we performed a stringent genome-wide identification of NBS-encoding genes in three cultivated and eight wild Oryza species, representing three different genomes (AA, BB, and FF) from four continents. A total of 2688 NBS-encoding genes were identified from 11 Oryza genomes. All the three progenitor species of cultivated rice, namely O. barthii, O. rufipogon, and O. nivara, were the richest reservoir of NBS-encoding genes (214, 313, and 307 respectively). Interestingly, the two Asian cultivated species showed a contrasting pattern in the number of NBS-encoding genes. While indica subspecies maintained nearly equal number of NBS genes as its progenitor (309 and 313), the japonica subspecies had retained only two third in the course of evolution (213 and 307). Other major sources for NBS-encoding genes could be (i) O. longistaminata since it had the highest proportion of NBS-encoding genes and (ii) O. glumaepatula as it clustered distinctly away from the rest of the AA genome species. The present study thus revealed that NBS-encoding genes can be exploited from the primary gene pool for disease resistance breeding in rice.     

Loss/retention and evolution of NBS-encoding genes upon whole genome triplication of Brassica rapa

A genome triplication took place in the ancestor of Brassiceae species after the split of the Arabidopsis lineage. The postfragmentation and shuffling of the genome turned the ancestral hexaploid back to diploids and caused the radiation of Brassiceae species. The course of speciation was accompanied by the loss of duplicate genes and also influenced the evolution of retained genes. Of all the genes, those encoding NBS domains are typical R genes that confer resistance to invading pathogens. In this study, using the genome of Arabidopsis thaliana as a reference, we examined the loss/retention of orthologous NBS-encoding loci in the tripled Brassica rapa genome and discovered differential loss/retention frequencies. Further analysis indicated that loci of different retention ratios showed different evolutionary patterns. The loci of classesII and III (maintaining two and three syntenic loci, respectively, multi-loci) show sharper expansions by tandem duplications, have faster evolutionary rates and have more potential to be associated with novel gene functions. On the other hand, the loci that are retained at the minimal rate (keeping only one locus, class I, single locus) showed opposite patterns. Phylogenetic analysis indicated that recombination and translocation events were common among multi-loci in B. rapa, and differential evolutionary patterns between multi- and single-loci are likely the consequence of recombination. Investigations towards other gene families demonstrated different evolutionary characteristics between different gene families. The evolution of genes is more likely determined by the property of each gene family, and the whole genome triplication provided only a specific condition.     

Genome-wide comparative analysis of NBS-encoding genes between Brassica species and Arabidopsis thaliana

Background

Plant disease resistance (R) genes with the nucleotide binding site (NBS) play an important role in offering resistance to pathogens. The availability of complete genome sequences of Brassica oleracea and Brassica rapa provides an important opportunity for researchers to identify and characterize NBS-encoding R genes in Brassica species and to compare with analogues in Arabidopsis thaliana based on a comparative genomics approach. However, little is known about the evolutionary fate of NBS-encoding genes in the Brassica lineage after split from A. thaliana.

Results

Here we present genome-wide analysis of NBS-encoding genes in B. oleracea, B. rapa and A. thaliana. Through the employment of HMM search and manual curation, we identified 157, 206 and 167 NBS-encoding genes in B. oleracea, B. rapa and A. thaliana genomes, respectively. Phylogenetic analysis among 3 species classified NBS-encoding genes into 6 subgroups. Tandem duplication and whole genome triplication (WGT) analyses revealed that after WGT of the Brassica ancestor, NBS-encoding homologous gene pairs on triplicated regions in Brassica ancestor were deleted or lost quickly, but NBS-encoding genes in Brassica species experienced species-specific gene amplification by tandem duplication after divergence of B. rapa and B. oleracea. Expression profiling of NBS-encoding orthologous gene pairs indicated the differential expression pattern of retained orthologous gene copies in B. oleracea and B. rapa. Furthermore, evolutionary analysis of CNL type NBS-encoding orthologous gene pairs among 3 species suggested that orthologous genes in B. rapa species have undergone stronger negative selection than those in B .oleracea species. But for TNL type, there are no significant differences in the orthologous gene pairs between the two species.

Conclusion

This study is first identification and characterization of NBS-encoding genes in B. rapa and B. oleracea based on whole genome sequences. Through tandem duplication and whole genome triplication analysis in B. oleracea, B. rapa and A. thaliana genomes, our study provides insight into the evolutionary history of NBS-encoding genes after divergence of A. thaliana and the Brassica lineage. These results together with expression pattern analysis of NBS-encoding orthologous genes provide useful resource for functional characterization of these genes and genetic improvement of relevant crops.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-3) contains supplementary material, which is available to authorized users.     

Phylogenetic and evolutionary analysis of NBS-encoding genes in Rutaceae fruit crops

The nucleotide-binding site leucine-rich repeat (NBS-LRR) genes are the largest class of disease resistance genes in plants. However, our understanding of the evolution of NBS-LRR genes in Rutaceae fruit crops is rather limited. We report an evolutionary study of 103 NBS-encoding genes isolated from Poncirus trifoliata (trifoliate orange), Citrus reticulata (tangerine) and their F₁ progeny. In all, 58 of the sequences contained a continuous open reading frame. Phylogenetic analysis classified the 58 NBS genes into nine clades, eight of which were genus specific. This was taken to imply that most of the ancestors of these NBS genes evolved after the genus split. The motif pattern of the 58 NBS-encoding genes was consistent with their phylogenetic profile. An extended phylogenetic analysis, incorporating citrus NBS genes from the public database, classified 95 citrus NBS genes into six clades, half of which were genus specific. RFLP analysis showed that citrus NBS-encoding genes have been evolving rapidly, and that they are unstable when passed through an intergeneric cross. Of 32 NBS-encoding genes tracked by gene-specific PCR, 24 showed segregation distortion among a set of 94 F₁ individuals. This study provides new insight into the evolution of Rutaceae NBS genes and their behaviour following an intergeneric cross.     

Genome-Wide Identification and Characterization of Nucleotide-Binding Site (NBS) Resistance Genes in Pineapple

Pineapple is a major tropical fruit and the most important crop processing CAM photosynthesis. It originated in southwest Brazil and northeast Paraguay and survived the harsh, semi-arid environment. Disease resistance genes have contributed to the survival and thriving of this species. The largest class of disease resistance (R) genes in plants consists of genes encoding nucleotide-binding site (NBS) domains. The sequenced genome of pineapple (Ananas comosus (L.) Merr.) provides a resource for analyzing the NBS-encoding genes in this species. A total of 177 NBS-encoding genes were identified using automated and manual analysis criteria, and these represent about 0.6 % of the total number of predicted pineapple genes. Five genes identified here contained the N-terminal Toll/Interleukin-l receptor (TIR) domain, and 46 genes carried the N-terminal Coiled-Coil (CC) motif. A majority of these NBS-encoding genes (84 %) contained a leucine-rich repeat (LRR) domain. A total of 130 of 177 (73 %) of these NBS-encoding genes were distributed across 20 pineapple linkage groups. The identification and characterization of NBS genes in pineapple yielded a valuable genomic resource and improved understanding of R genes in pineapple, which will facilitate the development of disease resistant pineapple cultivars.     

Genome-wide identification of NBS-encoding resistance genes in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Nucleotide-binding site (NBS)-encoding resistance genes are key plant disease-resistance genes and are abundant in plant genomes, comprising up to 2% of all genes. The availability of genome sequences from several plant models enables the identification and cloning of NBS-encoding genes from closely related species based on a comparative genomics approach. In this study, we used the genome sequence of Brassica rapa to identify NBS-encoding genes in the Brassica genome. We identified 92 non-redundant NBS-encoding genes [30 CC-NBS-LRR (CNL) and 62 TIR-NBS-LRR (TNL) genes] in approximately 100 Mbp of B. rapa euchromatic genome sequence. Despite the fact that B. rapa has a significantly larger genome than Arabidopsis thaliana due to a recent whole genome triplication event after speciation, B. rapa contains relatively small number of NBS-encoding genes compared to A. thaliana, presumably because of deletion of redundant genes related to genome diploidization. Phylogenetic and evolutionary analyses suggest that relatively higher relaxation of selective constraints on the TNL group after the old duplication event resulted in greater accumulation of TNLs than CNLs in both Arabidopsis and Brassica genomes. Recent tandem duplication and ectopic deletion are likely to have played a role in the generation of novel Brassica lineage-specific resistance genes.     

Isolation,Cloning and Characterization of Resistance Gene Analogues in Pearl Millet Based on Conserved Nucleotide‐binding Sites

Sequence analysis of plant disease resistance genes shows similarity among themselves, with the presence of conserved motifs common to the nucleotide‐binding site (NBS). Oligonucleotide degenerate primers designed from the conserved NBS motifs encoded by several plant disease resistance genes were used to amplify resistance gene analogues (RGAs) corresponding to the NBS sequences from the genomic DNA of various plant species. Using specific primers designed from the conserved NBS regions, 22 RGAs were cloned and sequenced from pearl millet (Pennisetum glaucum L. Br.). Phylogenetic analysis of the predicted amino acid sequences grouped the RGAs into nine distinct classes. GenBank database searches with the consensus protein sequences of each of the nine classes revealed their conserved NBS domains and similarity to other known R genes of various crop species. One RGA 213 was mapped onto LG1 and LG7 in the pearl millet linkage map. This is the first report of the isolation and characterization of RGAs from pearl millet, which will facilitate the improvement of marker‐assisted breeding strategies.     

Characterization of the rice blast resistance gene Pik cloned from Kanto51

To study similar, but distinct, plant disease resistance (R) specificities exhibited by allelic genes at the rice blast resistance locus Pik/Pikm, we cloned the Pik gene from rice cultivar Kanto51 and compared its molecular features with those of Pikm and of another Pik gene cloned from cv. Kusabue. Like Pikm, Pik is composed of two adjacent NBS-LRR (nucleotide-binding site, leucine-rich repeat) genes: the first gene, Pik1-KA, and the second gene, Pik2-KA. Pik from Kanto51 and Pik from Kusabue were not identical; although the predicted protein sequences of the second genes were identical, the sequences differed by three amino acids within the NBS domain of the first genes. The Pik proteins from Kanto51 and Kusabue differed from Pikm in eight and seven amino acids, respectively. Most of these substituted amino acids were within the coiled-coil (CC) and NBS domains encoded by the first gene. Of these substitutions, all within the CC domain were conserved between the two Pik proteins, whereas all within the NBS domain differed between them. Comparison of the two Pik proteins and Pikm suggests the importance of the CC domain in determining the resistance specificities of Pik and Pikm. This feature contrasts with that of most allelic or homologous NBS-LRR genes characterized to date, in which the major specificity determinant is believed to lie in the highly diverged LRR domain. In addition, our study revealed high evolutionary flexibility in the genome at the Pik locus, which may be relevant to the generation of new R specificities at this locus.     

Physcomitrella patens Has Kinase-LRR R Gene Homologs and Interacting Proteins

Plant disease resistance gene (R gene)-like sequences were screened from the Physcomitrella patens genome. We found 603 kinase-like, 475 Nucleotide Binding Site (NBS)-like and 8594 Leucine Rich Repeat (LRR)-like sequences by homology searching using the respective domains of PpC24 (Accession No. BAD38895), which is a candidate kinase-NBS-LRR (kinase-NL) type R-like gene, as a reference. The positions of these domains in the genome were compared and 17 kinase-NLs were predicted. We also found four TIR-NBS-LRR (TIR-NL) sequences with homology to Arabidopsis TIR-NL (NM_001125847), but three out of the four TIR-NLs had tetratricopeptide repeats or a zinc finger domain in their predicted C-terminus. We also searched for kinase-LRR (KLR) type sequences by homology with rice OsXa21 and Arabidopsis thaliana FLS2. As a result, 16 KLRs with similarity to OsXa21 were found. In phylogenetic analysis of these 16 KLRs, PpKLR36, PpKLR39, PpKLR40, and PpKLR43 formed a cluster with OsXa21. These four PpKLRs had deduced transmembrane domain sequences and expression of all four was confirmed. We also found 14 homologs of rice OsXB3, which is known to interact with OsXa21 and is involved in signal transduction. Protein–protein interaction was observed between the four PpKLRs and at least two of the XB3 homologs in Y2H analysis.     

The genomic architecture of disease resistance in lettuce

Genbank and The Compositae Genome Project database, containing over 42,000 lettuce unigenes from Lactuca sativa cv. Salinas and L. serriola accession UC96US23 were mined to identify 702 candidate genes involved in pathogen recognition (RGCs), resistance signal transduction, defense responses, and disease susceptibility. In addition, to identify sequences representing additional sub-families of nucleotide binding site (NBS)-leucine-rich repeat encoding genes; the major classes of resistance genes (R-genes), NBS-encoding sequences were amplified by PCR using degenerate oligonucleotides designed to NBS sub-families specific to the subclass Asteridae, which includes the Compositae family. These products were cloned and sequenced resulting in 18 novel NBS sequences from cv. Salinas and 15 novel NBS sequences from UC96US23. Using a variety of marker technologies, 294 of the 735 candidate disease resistance genes were mapped in our primary mapping population, which consisted of 119 F₇ recombinant inbred lines derived from an interspecific cross between cv. Salinas and UC96US23. Using markers shared across multiple genetic maps, 36 resistance phenotypic loci, including two new loci for resistance to downy mildew and two quantitative trait loci for resistance to anthracnose were positioned onto the reference map to provide a global view of the genomic architecture of disease resistance in lettuce and to identify candidate genes for resistance phenotypes. The majority but not all of the resistance phenotypes were genetically associated with RGCs. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. An erratum to this article can be found at