Loss and retention of resistance genes in five species of the Brassicaceae family

The nuclear-localized RNA helicase 13 is essential for chloroplast development in Arabidopsis thaliana

The chloroplast is a semi-autonomous organelle with a double membrane structure, and its structural stability is a prerequisite for its correct function. Chloroplast development is regulated by known nuclear-encoded chloroplast proteins or proteins encoded within the chloroplast itself. However, the mechanism of chloroplast development regulated by other organelles remains largely unknown. Here, we report that the nuclear-localized DEAD-box RNA helicase 13 (RH13) is essential for chloroplast development in Arabidopsis thaliana. RH13 is widely expressed in tissues and localized to the nucleolus. A homozygous rh13 mutant shows abnormal chloroplast structure and leaf morphogenesis. Proteomic analysis showed that the expression levels of photosynthesis-related proteins in chloroplasts were reduced due to loss of RH13. Furthermore, RNA-sequencing and proteomics data revealed decreases in the expression levels of these chloroplast-related genes, which undergo alternative splicing events in the rh13 mutant. Taken together, we propose that nucleolus-localized RH13 is critical for chloroplast development in Arabidopsis.

The Arabidopsis effector-triggered immunity landscape is conserved in oilseed crops

The bacterial phytopathogen Pseudomonas syringae causes disease on a wide array of plants, including the model plant Arabidopsis thaliana and its agronomically important relatives in the Brassicaceae family. To cause disease, P. syringae delivers effector proteins into plant cells through a type III secretion system. In response, plant nucleotide-binding leucine-rich repeat proteins recognize specific effectors and mount effector-triggered immunity (ETI). While ETI is pervasive across A. thaliana, with at least 19 families of P. syringae effectors recognized in this model species, the ETI landscapes of crop species have yet to be systematically studied. Here, we investigated the conservation of the A. thaliana ETI landscape in two closely related oilseed crops, Brassica napus (canola) and Camelina sativa (false flax). We show that the level of immune conservation is inversely related to the degree of evolutionary divergence from A. thaliana, with the more closely related C. sativa losing ETI responses to only one of the 19 P. syringae effectors tested, while the more distantly related B. napus loses ETI responses to four effectors. In contrast to the qualitative conservation of immune response, the quantitative rank order is not as well-maintained across the three species and diverges increasingly with evolutionary distance from A. thaliana. Overall, our results indicate that the A. thaliana ETI profile is qualitatively conserved in oilseed crops, but quantitatively distinct.

Comparison and evolutionary analysis of Brassica nucleotide binding site leucine rich repeat (NLR) genes and importance for disease resistance breeding

The Brassica genus contains many agriculturally significant oilseed and vegetable crops, however the crop yield is threatened by a range of fungal and bacterial pathogens. Nucleotide Binding Site Leucine Rich Repeat (NLR) genes play important roles in plant innate immunity. The evolution of NLR genes is influenced by genomic processes and pathogen selection. At the whole genome level, whole genome duplications (WGDs) generate abundant gene copies, most of which are lost during genome fractionation. At sub-genomic levels, some retained copies undergo duplication forming clusters which facilitate rapid evolution through recombination. The number, distribution and genetic variations of the NLR genes vary among Brassica species and within populations suggesting differential selection pressure exerted by pathogen populations throughout the evolutionary history. A study of the evolution of disease resistance genes in agriculturally important plants such as Brassicas helps gain insights into their function and inform the identification of resistance genes for breeding of resistant lines.

Genome sequence analysis provides insights on genomic variation and late blight resistance genes in potato somatic hybrid (parents and progeny)

Wild Solanum species are the important resources for potato improvement. With the availability of potato genome and sequencing progress, knowledge about genomic resources is essential for novel genes discovery. Hence, the aim of this study was to decipher draft genome sequences of unique potato genotypes i.e. somatic hybrid P8 (J1), wild species S. pinnatisectum (J2), progeny MSH/14-112 (P8 × cv. Kufri Jyoti) (J3), and S. tuberosum dihaploid C-13 (J4). Draft genome sequencing using Illumina platform and reference-based assemblies with the potato genome yielded genome assembly size of 725.01 Mb (J1), 724.95 Mb (J2), 725.01 Mb (J3), and 809.59 Mb (J4). Further, 39,260 (J1), 25,711 (J2), 39,730 (J3) and 30,241 (J4) genes were identified and 17,411 genes were found common in the genotypes particularly late blight resistance genes (R3a, RGA2, RGA3, R1B-16, Rpi-blb2, Rpi and Rpi-vnt1). Gene ontology (GO) analysis showed that molecular function was predominant and signal transduction was major KEGG pathways. Further, gene enrichment analysis revealed dominance of metabolic process (GO: 0008152) in all the samples. Phylogeny analysis showed relatedness with potato and other plant species. Heterozygous single nucleotide polymorphism (SNP) was more than homozygous, and SNP in genic region was more than inter-genic region. Copy number variation (CNV) analysis indicated greater number of deletions than duplications. Sequence diversity and conserved motifs analysis revealed variation for late blight resistance genes. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis showed differential expression of late blight resistance genes. Our study provides insights on genome sequence, structural variation and late blight resistance genes in potato somatic hybrid (parents and progeny) for future research.

Recent Findings Unravel Genes and Genetic Factors Underlying Leptosphaeria maculans Resistance in Brassica napus and Its Relatives

Among the Brassica oilseeds, canola (Brassica napus) is the most economically significant globally. However, its production can be limited by blackleg disease, caused by the fungal pathogen Lepstosphaeria maculans. The deployment of resistance genes has been implemented as one of the key strategies to manage the disease. Genetic resistance against blackleg comes in two forms: qualitative resistance, controlled by a single, major resistance gene (R gene), and quantitative resistance (QR), controlled by numerous, small effect loci. R-gene-mediated blackleg resistance has been extensively studied, wherein several genomic regions harbouring R genes against L. maculans have been identified and three of these genes were cloned. These studies advance our understanding of the mechanism of R gene and pathogen avirulence (Avr) gene interaction. Notably, these studies revealed a more complex interaction than originally thought. Advances in genomics help unravel these complexities, providing insights into the genes and genetic factors towards improving blackleg resistance. Here, we aim to discuss the existing R-gene-mediated resistance, make a summary of candidate R genes against the disease, and emphasise the role of players involved in the pathogenicity and resistance. The comprehensive result will allow breeders to improve resistance to L. maculans, thereby increasing yield.

Immunodiversity of the Arabidopsis ZAR1 NLR Is Conveyed by Receptor-Like Cytoplasmic Kinase Sensors

The Arabidopsis nucleotide-binding leucine-rich repeat protein ZAR1 can recognize at least six distinct families of pathogenic effector proteins to mount an effector-triggered immune response. This remarkable immunodiversity appears to be conveyed by receptor-like cytoplasmic kinase (RLCK) complexes, which associate with ZAR1 to sense several effector-induced kinase perturbations. Here we show that the recently identified ZAR1-mediated immune responses against the HopX1, HopO1, and HopBA1 effector families of Pseudomonas syringae rely on an expanded diversity of RLCK sensors. We show that individual sensors can recognize distinct effector families, thereby contributing to the expanded surveillance potential of ZAR1 and supporting its role as a guardian of the plant kinome.

Evolutionary balance between LRR domain loss and young NBS-LRR genes production governs disease resistance in Arachis hypogaea cv. Tifrunner

BACKGROUND

Cultivated peanut (Arachis hypogaea L.) is an important oil and protein crop, but it has low disease resistance; therefore, it is important to reveal the number, sequence features, function, and evolution of genes that confer resistance. Nucleotide-binding site-leucine-rich repeats (NBS-LRRs) are resistance genes that are involved in response to various pathogens.

RESULTS

We identified 713 full-length NBS-LRRs in A. hypogaea cv. Tifrunner. Genetic exchange events occurred on NBS-LRRs in A. hypogaea cv. Tifrunner, which were detected in the same subgenomes and also found in different subgenomes. Relaxed selection acted on NBS-LRR proteins and LRR domains in A. hypogaea cv. Tifrunner. Using quantitative trait loci (QTL), we found that NBS-LRRs were involved in response to late leaf spot, tomato spotted wilt virus, and bacterial wilt in A. duranensis (2 NBS-LRRs), A. ipaensis (39 NBS-LRRs), and A. hypogaea cv. Tifrunner (113 NBS-LRRs). In A. hypogaea cv. Tifrunner, 113 NBS-LRRs were classified as 75 young and 38 old NBS-LRRs, indicating that young NBS-LRRs were involved in response to disease after tetraploidization. However, compared to A. duranensis and A. ipaensis, fewer LRR domains were found in A. hypogaea cv. Tifrunner NBS-LRR proteins, partly explaining the lower disease resistance of the cultivated peanut.

CONCLUSIONS

Although relaxed selection acted on NBS-LRR proteins and LRR domains, LRR domains were preferentially lost in A. hypogaea cv. Tifrunner compared to A. duranensis and A. ipaensis. The QTL results suggested that young NBS-LRRs were important for resistance against diseases in A. hypogaea cv. Tifrunner. Our results provid insight into the greater susceptibility of A. hypogaea cv. Tifrunner to disease compared to A. duranensis and A. ipaensis.

An EDS1-SAG101 Complex Is Essential for TNL-Mediated Immunity in Nicotiana benthamiana

Heterodimeric complexes containing the lipase-like protein ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1) are regarded as central regulators of plant innate immunity. In this context, a complex of EDS1 with PHYTOALEXIN DEFICIENT4 (PAD4) is required for basal resistance and signaling downstream of immune receptors containing an N-terminal Toll-interleukin-1 receptor-like domain (TNLs) in Arabidopsis (Arabidopsis thaliana). Here we analyze EDS1 functions in the model Solanaceous plant Nicotiana benthamiana (Nb). Stable Nb mutants deficient in EDS1 complexes are not impaired in basal resistance, a finding which contradicts a general role for EDS1 in immunity. In Nb, PAD4 demonstrated no detectable immune functions, but TNL-mediated resistance responses required EDS1 complexes incorporating a SENESCENCE ASSOCIATED GENE101 (SAG101) isoform. Intriguingly, SAG101 is restricted to those genomes also encoding TNL receptors, and we propose it may be required for TNL-mediated immune signaling in most plants, except the Brassicaceae. Transient complementation in Nb was used for accelerated mutational analyses while avoiding complex biotic interactions. We identify a large surface essential for EDS1-SAG101 immune functions that extends from the N-terminal lipase domains to the C-terminal EDS1-PAD4 domains and might mediate interaction partner recruitment. Furthermore, this work demonstrates the value of genetic resources in Nb, which will facilitate elucidation of EDS1 functions.

A Species-Wide Inventory of NLR Genes and Alleles in Arabidopsis thaliana

Infectious disease is both a major force of selection in nature and a prime cause of yield loss in agriculture. In plants, disease resistance is often conferred by nucleotide-binding leucine-rich repeat (NLR) proteins, intracellular immune receptors that recognize pathogen proteins and their effects on the host. Consistent with extensive balancing and positive selection, NLRs are encoded by one of the most variable gene families in plants, but the true extent of intraspecific NLR diversity has been unclear. Here, we define a nearly complete species-wide pan-NLRome in Arabidopsis thaliana based on sequence enrichment and long-read sequencing. The pan-NLRome largely saturates with approximately 40 well-chosen wild strains, with half of the pan-NLRome being present in most accessions. We chart NLR architectural diversity, identify new architectures, and quantify selective forces that act on specific NLRs and NLR domains. Our study provides a blueprint for defining pan-NLRomes.

A Species-Wide Inventory of NLR Genes and Alleles in Arabidopsis thaliana

Infectious disease is both a major force of selection in nature and a prime cause of yield loss in agriculture. In plants, disease resistance is often conferred by nucleotide-binding leucine-rich repeat (NLR) proteins, intracellular immune receptors that recognize pathogen proteins and their effects on the host. Consistent with extensive balancing and positive selection, NLRs are encoded by one of the most variable gene families in plants, but the true extent of intraspecific NLR diversity has been unclear. Here, we define a nearly complete species-wide pan-NLRome in Arabidopsis thaliana based on sequence enrichment and long-read sequencing. The pan-NLRome largely saturates with approximately 40 well-chosen wild strains, with half of the pan-NLRome being present in most accessions. We chart NLR architectural diversity, identify new architectures, and quantify selective forces that act on specific NLRs and NLR domains. Our study provides a blueprint for defining pan-NLRomes.

Two Arabidopsis Receptor-like Cytoplasmic Kinases SZE1 and SZE2 Associate with the ZAR1-ZED1 Complex and Are Required for Effector-Triggered Immunity

Plants utilize intracellular nucleotide-binding leucine-rich repeat domain-containing receptors (NLRs) to recognize pathogen effectors and induce a robust defense response named effector-triggered immunity (ETI). The Arabidopsis NLR protein HOPZ-ACTIVATED RESISTANCE 1 (ZAR1) forms a precomplex with HOPZ-ETI-DEFICIENT 1 (ZED1), a receptor-like cytoplasmic kinase (RLCK) XII-2 subfamily member, to recognize the Pseudomonas syringae effector HopZ1a. We previously described a dominant mutant of Arabidopsis ZED1, zed1-D, which displays temperature-sensitive autoimmunity in a ZAR1-dependent manner. Here, we report that the RLCKs SUPPRESSOR OF ZED1-D1 (SZE1) and SZE2 associate with the ZAR1-ZED1 complex and are required for the ZED1-D-activated autoimmune response and HopZ1a-triggered immunity. We show that SZE1 but not SZE2 has autophosphorylation activity, and that the N-terminal myristoylation of both SZE1 and SZE2 is critical for their plasma membrane localization and ZED1-D-activated autoimmunity. Furthermore, we demonstrate that SZE1 and SZE2 both interact with ZAR1 to form a functional complex and are required for resistance against P. syringae pv. tomato DC3000 expressing HopZ1a. We also provide evidence that SZE1 and SZE2 interact with HopZ1a and function together with ZED1 to change the intramolecular interactions of ZAR1, leading to its activation. Taken together, our results reveal SZE1 and SZE2 as critical signaling components of HopZ1a-triggered immunity.

Identification and expression profiling analysis of NBS-LRR genes involved in Fusarium oxysporum f.sp. conglutinans resistance in cabbage

As one of the most important resistance (R) gene families in plants, the NBS-LRR genes, encoding proteins with nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains, play significant roles in resisting pathogens. The published genomic data for cabbage (Brassica oleracea L.) provide valuable data to identify and characterize the genomic organization of cabbage NBS-LRR genes. Ultimately, we identified 105 TIR (N-terminal Toll/interleukin-1 receptor)-NBS-LRR (TNL) genes and 33 CC (coiled-coil)-NBS-LRR (CNL) genes. Further research indicated that 50.7% of the 138 NBS-LRR genes exist in 27 clusters and there are large differences among the gene structures and protein characteristics. Conserved motif and phylogenetic analysis showed that the structures of TNLs and CNLs were similar, with some differences. These NBS-LRRs are evolved under negative selection and mostly arose from whole-genome duplication events during evolution. Tissue-expression profiling of NBS-LRR genes revealed that 37.1% of the TNL genes are highly or specifically expressed in roots, especially the genes on chromosome 7 (76.5%). Digital gene expression and reverse transcription PCR analyses revealed the expression patterns of the NBS-LRR genes upon challenge by Fusarium oxysporum f.sp. conglutinans: nine genes were upregulated, and five were downregulated. The major resistance gene Foc1 probably works together with the other four genes in the same cluster to resist F. oxysporum infection.

NBS-Encoding Genes in Brassica napus Evolved Rapidly After Allopolyploidization and Co-localize With Known Disease Resistance Loci

Genes containing nucleotide-binding sites (NBS) play an important role in pathogen resistance in plants. However, the evolutionary fate of NBS-encoding genes after formation of allotetraploid Brassica napus (A_nA_nC_nC_n, 2n = 38) is still unknown. We performed a genome-wide comparison of putatively functional NBS-encoding genes in B. napus and its progenitor species Brassica rapa (A_rA_r, 2n = 20) and Brassica oleracea (C_oC_o, 2n = 18), identifying 464, 202, and 146 putatively functional NBS-encoding genes respectively, with genes unevenly distributed in several clusters. The A_n-subgenome of B. napus possessed similar numbers of NBS-encoding genes (191 genes) to the A_r genome of B. rapa (202 genes) and similar clustering patterns. However, the C_n genome of B. napus had many more genes (273) than the B. oleracea C_o genome (146), with different clustering trends. Only 97 NBS-encoding genes (66.4%) in B. oleracea were homologous with NBS-encoding genes in B. napus, while 176 NBS-encoding genes (87.1%) were homologous between B. rapa and B. napus. These results suggest a greater diversification of NBS-encoding genes in the C genome may have occurred after formation of B. napus. Although most NBS-encoding genes in B. napus appeared to derive from the progenitors, the birth and death of several NBS-encoding genes was also putatively mediated by non-homologous recombination. The Ka/Ks values of most homologous pairs between B. napus and the progenitor species were less than 1, suggesting purifying selection during B. napus evolution. The majority of NBS-encoding genes (60% in all species) showed higher expression levels in root tissue (out of root, leaf, stem, seed and flower tissue types). Comparative analysis of NBS-encoding genes with mapped resistance QTL against three major diseases of B. napus (blackleg, clubroot and Sclerotinia stem rot) found 204 NBS-encoding genes in B. napus located within 71 resistance QTL intervals. The majority of NBS-encoding genes were co-located with resistance QTLs against a single disease, while 47 genes were co-located with QTLs against two diseases and 3 genes were co-located with QTLs against all three. Our results revealed significant variation as well as interesting evolutionary trajectories of NBS-encoding genes in the different Brassica subgenomes, while co-localization of NBS-encoding genes and resistance QTL may facilitate resistance breeding in oilseed rape.

Structural, Functional, and Genomic Diversity of Plant NLR Proteins: An Evolved Resource for Rational Engineering of Plant Immunity

Plants employ a diverse intracellular system of NLR (nucleotide binding-leucine-rich repeat) innate immune receptors to detect pathogens of all types. These receptors represent valuable agronomic traits that plant breeders rely on to maximize yield in the face of devastating pathogens. Despite their importance, the mechanistic underpinnings of NLR-based disease resistance remain obscure. The rapidly increasing numbers of plant genomes are revealing a diverse array of NLR-type immune receptors. In parallel, mechanistic studies are describing diverse functions for NLR immune receptors. In this review, we intend to broadly describe how the structural, functional, and genomic diversity of plant immune receptors can provide a valuable resource for rational engineering of plant immunity.

Evolutionary Divergence of TNL Disease-Resistant Proteins in Soybean (Glycine max) and Common Bean (Phaseolus vulgaris)

Disease-resistant genes (R genes) encode proteins that are involved in protecting plants from their pathogens and pests. Availability of complete genome sequences from soybean and common bean allowed us to perform a genome-wide identification and analysis of the Toll interleukin-1 receptor-like nucleotide-binding site leucine-rich repeat (TNL) proteins. Hidden Markov model (HMM) profiling of all protein sequences resulted in the identification of 117 and 77 regular TNL genes in soybean and common bean, respectively. We also identified TNL gene homologs with unique domains, and signal peptides as well as nuclear localization signals. The TNL genes in soybean formed 28 clusters located on 10 of the 20 chromosomes, with the majority found on chromosome 3, 6 and 16. Similarly, the TNL genes in common bean formed 14 clusters located on five of the 11 chromosomes, with the majority found on chromosome 10. Phylogenetic analyses of the TNL genes from Arabidopsis, soybean and common bean revealed less divergence within legumes relative to the divergence between legumes and Arabidopsis. Syntenic blocks were found between chromosomes Pv10 and Gm03, Pv07 and Gm10, as well as Pv01 and Gm14. The gene expression data revealed basal level expression and tissue specificity, while analysis of available microRNA data showed 37 predicted microRNA families involved in targeting the identified TNL genes in soybean and common bean.

Defended to the Nines: 25 Years of Resistance Gene Cloning Identifies Nine Mechanisms for R Protein Function

Plants have many, highly variable resistance (R) gene loci, which provide resistance to a variety of pathogens. The first R gene to be cloned, maize (Zea mays) Hm1, was published over 25 years ago, and since then, many different R genes have been identified and isolated. The encoded proteins have provided clues to the diverse molecular mechanisms underlying immunity. Here, we present a meta-analysis of 314 cloned R genes. The majority of R genes encode cell surface or intracellular receptors, and we distinguish nine molecular mechanisms by which R proteins can elevate or trigger disease resistance: direct (1) or indirect (2) perception of pathogen-derived molecules on the cell surface by receptor-like proteins and receptor-like kinases; direct (3) or indirect (4) intracellular detection of pathogen-derived molecules by nucleotide binding, leucine-rich repeat receptors, or detection through integrated domains (5); perception of transcription activator-like effectors through activation of executor genes (6); and active (7), passive (8), or host reprogramming-mediated (9) loss of susceptibility. Although the molecular mechanisms underlying the functions of R genes are only understood for a small proportion of known R genes, a clearer understanding of mechanisms is emerging and will be crucial for rational engineering and deployment of novel R genes.

Investigation of orthologous pathogen recognition gene-rich regions in solanaceous species

Pathogen receptor proteins such as receptor-like protein (RLP), receptor-like kinase (RLK), and nucleotide-binding leucine-rich repeat (NLR) play a leading role in plant immunity activation. The genome architecture of such genes has been extensively investigated in several plant species. However, we still know little about their elaborate reorganization that arose during the plant speciation process. Using recently released pepper and eggplant genome sequences, we were able to identify 1097 pathogen recognition genes (PRGs) in the cultivated pepper Zunla-1 and 775 in the eggplant line Nakate-Shinkuro. The retrieved genes were analysed for their tendency to cluster, using different methods to infer the means of grouping. Orthologous relationships among clustering loci were found, and interesting reshuffling within given loci was observed for each analysed species. The information obtained was integrated into a comparative map to highlight the evolutionary dynamics in which the PRG loci were involved. Diversification of 14 selected PRG-rich regions was also explored using a DNA target-enrichment approach. A large number of gene variants were found as well as rearrangements of sequences encoding single protein domain and changes in chromosome gene order among species. Gene duplication and transposition activity have clearly influenced plant genome R-gene architecture and diversification. Our findings contribute to addressing several biological questions concerning the parallel evolution that occurred between genomes of the family Solanaceae. Moreover, the integration of different methods proved a powerful approach to reconstruct the evolutionary history in plant families and to transfer important biology findings among plant genomes.

Comparative metabolite and genome analysis of tuber-bearing potato species

The cultivated potato Solanum tuberosum is unrivalled among crop plants for its wild relatives, which potentially represent an important source of genetic diversity to improve the nutritional value of potato varieties and understand metabolism regulation. The main aim of this research was to profile human health-related metabolites in a number of clones from 13 Solanum species. Results from HPLC-DAD and LC-ESI-MS analyses highlighted a high interspecific variability in the level of metabolites analysed. Ascorbic acid was confirmed to be the most abundant antioxidant in potato and chlorogenic acid the primary polyphenol. Generally, metabolite-based hierarchical clustering (HCL) and correlation networks did not group clones of identical species in the same cluster. This might be due to various factors, including the outcrossing nature of potato species, gene expression level and metabolic profiling techniques. Access to the genome sequence of S. tuberosum and S. commersonii allowed comparison of the genes involved in ascorbic acid, aromatic amino acid, phenylpropanoid and glycoalkaloid biosynthesis and helped interpret their respective pathways.

Altools: a user friendly NGS data analyser

BACKGROUND

Genotyping by re-sequencing has become a standard approach to estimate single nucleotide polymorphism (SNP) diversity, haplotype structure and the biodiversity and has been defined as an efficient approach to address geographical population genomics of several model species. To access core SNPs and insertion/deletion polymorphisms (indels), and to infer the phyletic patterns of speciation, most such approaches map short reads to the reference genome. Variant calling is important to establish patterns of genome-wide association studies (GWAS) for quantitative trait loci (QTLs), and to determine the population and haplotype structure based on SNPs, thus allowing content-dependent trait and evolutionary analysis. Several tools have been developed to investigate such polymorphisms as well as more complex genomic rearrangements such as copy number variations, presence/absence variations and large deletions. The programs available for this purpose have different strengths (e.g. accuracy, sensitivity and specificity) and weaknesses (e.g. low computation speed, complex installation procedure and absence of a user-friendly interface). Here we introduce Altools, a software package that is easy to install and use, which allows the precise detection of polymorphisms and structural variations.

RESULTS

Altools uses the BWA/SAMtools/VarScan pipeline to call SNPs and indels, and the dnaCopy algorithm to achieve genome segmentation according to local coverage differences in order to identify copy number variations. It also uses insert size information from the alignment of paired-end reads and detects potential large deletions. A double mapping approach (BWA/BLASTn) identifies precise breakpoints while ensuring rapid elaboration. Finally, Altools implements several processes that yield deeper insight into the genes affected by the detected polymorphisms. Altools was used to analyse both simulated and real next-generation sequencing (NGS) data and performed satisfactorily in terms of positive predictive values, sensitivity, the identification of large deletion breakpoints and copy number detection.

CONCLUSIONS

Altools is fast, reliable and easy to use for the mining of NGS data. The software package also attempts to link identified polymorphisms and structural variants to their biological functions thus providing more valuable information than similar tools.

The Solanum commersonii Genome Sequence Provides Insights into Adaptation to Stress Conditions and Genome Evolution of Wild Potato Relatives

Here, we report the draft genome sequence of Solanum commersonii, which consists of ∼830 megabases with an N50 of 44,303 bp anchored to 12 chromosomes, using the potato (Solanum tuberosum) genome sequence as a reference. Compared with potato, S. commersonii shows a striking reduction in heterozygosity (1.5% versus 53 to 59%), and differences in genome sizes were mainly due to variations in intergenic sequence length. Gene annotation by ab initio prediction supported by RNA-seq data produced a catalog of 1703 predicted microRNAs, 18,882 long noncoding RNAs of which 20% are shown to target cold-responsive genes, and 39,290 protein-coding genes with a significant repertoire of nonredundant nucleotide binding site-encoding genes and 126 cold-related genes that are lacking in S. tuberosum. Phylogenetic analyses indicate that domesticated potato and S. commersonii lineages diverged ∼2.3 million years ago. Three duplication periods corresponding to genome enrichment for particular gene families related to response to salt stress, water transport, growth, and defense response were discovered. The draft genome sequence of S. commersonii substantially increases our understanding of the domesticated germplasm, facilitating translation of acquired knowledge into advances in crop stability in light of global climate and environmental changes.