Genome wide re-sequencing of newly developed Rice Lines from common wild rice (Oryza rufipogon Griff.) for the identification of NBS-LRR genes

The protective roles of Oryza glumaepatula and phytohormone in enhancing rice tolerance to cadmium stress by regulating gene expression, morphological, physiological, and antioxidant defense system

The heavy metal cadmium (Cd) is highly poisonous and has received significant attention from environmental scientists due to its harmful impacts on plants. Oryza glumaepatula is a wild rice that contains useful genes against biotic and abiotic stresses. Therefore, the current study used SG007, a single-segment substitution line (SSSL), generated by crossing O. glumaepatula with an elite rice cultivar (HJX74), to evaluate the resistance potential against Cd. Moreover, we assessed the efficacy of strigolactone GR24 (1 μM) against Cd toxicity (100 μM) by investigating physiological, biochemical, and molecular mechanisms in both cultivars (i.e., SG007 and HJX74). The findings of this study revealed that Cd toxicity declined the chlorophyll a, chlorophyll b, and carotenoids by 50%, 20%, and 44% in SG007, and 58%, 39%, and 59% in HJX74 by enhancing electrolyte leakage (EL), malondialdehyde (MDA), and hydrogen peroxide (H2O2) by 113%, 184%, and 119% in SG007 and 248%, 273% and 195% in HJX74, respectively. GR24 improved growth under Cd stress in both cultivars, and SG007 exhibited better plant growth parameters, antioxidant enzymatic activities, nitric oxide synthase (NOS), and nitric oxide (NO) levels than HJX74 under Cd toxicity. GR24 with SG007 regulated expressions of Cd transporters and reduced the cytological disruptions in cell organelles. The combined utilization of SG007 and GR24 reduced Cd accumulation and oxidative stress and improved plant growth parameters and enzymatic activities. In conclusion, our study highlights the potential of utilizing SG007 in conjunction with GR24 as a practical strategy to mitigate Cd pollution in rice. The results not only underscore the beneficial effects of strigolactone GR24 in alleviating Cd-induced stress but also emphasize the valuable genetic traits of O. glumaepatula in developing rice lines with enhanced tolerance to heavy metals, offering broader implications for sustainable agriculture and crop improvement in contaminated environments.

Genetic dissection of the soybean dwarf mutant dm with integrated genomic, transcriptomic and methylomic analyses

Plant height affects crop production and breeding practices, while genetic control of dwarfism draws a broad interest of researchers. Dwarfism in soybean (Glycine max) is mainly unexplored. Here, we characterized a dwarf mutant dm screened from ethyl methanesulfonate (EMS) mutated seeds of the soybean cultivar Zhongpin 661(ZP). Phenotypically, dm showed shorter and thinner stems, smaller leaves, and more nodes than ZP under greenhouse conditions. Genetically, whole-genome sequencing and comparison revealed that 210K variants of SNPs and InDel in ZP relative to the soybean reference genome Williams82, and EMS mutagenesis affected 636 genes with variants predicted to have a large impact on protein function in dm. Whole-genome methylation sequencing found 704 differentially methylated regions in dm. Further whole-genome RNA-Seq based transcriptomic comparison between ZP and dm leaves revealed 687 differentially expressed genes (DEGs), including 263 up-regulated and 424 down-regulated genes. Integrated omics analyses revealed 11 genes with both differential expressions and DNA variants, one gene with differential expression and differential methylation, and three genes with differential methylation and sequence variation, worthy of future investigation. Genes in cellulose, fatty acids, and energy-associated processes could be the key candidate genes for the dwarf phenotype. This study provides genetic clues for further understanding of the genetic control of dwarfism in soybean. The genetic resources could help to inbreed new cultivars with a desirable dwarf characteristic.

Identification of QTL under Brassinosteroid-Combined Cold Treatment at Seedling Stage in Rice Using Genotyping-by-Sequencing (GBS)

Cold stress is a major threat to the sustainability of rice yield. Brassinosteroids (BR) application can enhance cold tolerance in rice. However, the regulatory mechanism related to cold tolerance and the BR signaling pathway in rice has not been clarified. In the current study, the seedling shoot length (SSL), seedling root length (SRL), seedling dry weight (SDW), and seedling wet weight (SWW) were used as the indices for identifying cold tolerance under cold stress and BR-combined cold treatment in a backcross recombinant inbred lines (BRIL) population. According to the phenotypic characterization for cold tolerance and a high-resolution SNP genetic map obtained from the GBS technique, a total of 114 QTLs were identified, of which 27 QTLs were detected under cold stress and 87 QTLs under BR-combined cold treatment. Among them, the intervals of many QTLs were coincident under different treatments, as well as different traits. A total of 13 candidate genes associated with cold tolerance or BR pathway, such as BRASSINAZOLE RESISTANT1 (OsBZR1), OsWRKY77, AP2 domain-containing protein, zinc finger proteins, basic helix-loop-helix (bHLH) protein, and auxin-induced protein, were predicted. Among these, the expression levels of 10 candidate genes were identified under different treatments in the parents and representative BRIL individuals. These results were helpful in understanding the regulation relationship between cold tolerance and BR pathway in rice.

Whole-genome resequencing reveals genomic footprints of Italian sweet and hot pepper heirlooms giving insight into genes underlying key agronomic and qualitative traits

BACKGROUND

Pepper is a major crop species of the Solanaceae family, largely appreciated for its high nutritional and healthy contribution to human diets. In the Mediterranean basin, the favorable pedoclimatic conditions enhanced the selection of several diversified landraces cultivated pepper (Capsicum annuum), for whom Italy can be considered a main pole of diversification. Hence, a survey of traditional C. annuum genetic resources is essential for deep understanding of such diversity and for applications in genomics assisted breeding. Here, we report whole-genome resequencing analyses of two sweet and two pungent genotypes highly diffused in South Italy and representative of the variability for shape, colour and nutritional properties.

RESULTS

The four genomes were reconstructed at a chromosomal scale using a reference-guided approach, based on a dataset of 2.6 billion paired-end reads, corresponding to 20× genome coverage and a mapping rate above 99% for a final genomes size of approximately 3 Gb. After five iterations of variant calling, a total of 29,258,818 single nucleotide polymorphisms (SNPs) and 1,879,112 InDels, were identified. Substantial differences were observed among the four genomes based on geographical origin, with chromosomes 9 and 11 showing more polymorphisms in the accessions with higher fruit weight and absence of pungency. Among the identified variants, a small private indel (T - > TA) shared between sweet and big fruits accessions induces a frameshift with the generation of a new stop codon in a gene annotated as extensin, whereas two private SNPs within hot types were identified in 1-aminocyclopropane-1-carboxylate oxidase (ACO), a key gene involved in fruit ripening. The estimation of repetitive elements highlights a preponderant presence of Long Terminal Repeats (LTRs), the majority of which belonged to Gypsy superfamily. By comparing the four genomes with publicly available references including 'CM334' and Zunla-1 highlight the presence of 49,475 shared gene families.

CONCLUSIONS

The new genomic sequences aim to enrich the whole genome information of pepper local varieties, providing a valuable tool for precision gene mapping, marker discovery, comparative studies. Such knowledge widens the frontiers to understand the selection history of Italian pepper landraces toward the recognition of specificity local agri-food products marks.

Omics: a tool for resilient rice genetic improvement strategies

Rice is pivotal pyramid of about half of the world population. Bearing small genome size and worldwide utmost food crop rice has been known as ideal cereal crop for genome research. Currently, decreasing water table and soil fatigue are big challenges and intense consequences in changing climate. Whole sequenced genome of rice sized 389 Mb of which 95% is covered with excellent mapping order. Sequenced rice genome helps in molecular biology and transcriptomics of cereals as it provides whole genome sequence of indica and japonica sub species. Through rice genome sequencing and functional genomics, QTLs or genes, genetic variability and halophyte blocks for agronomic characters were identified which have proved much more useful in molecular breeding and direct selection. There are different numbers of genes or QTLs identified for yield related traits i.e., 6 QTLs/genes for plant architecture, 6 for panicle characteristics, 4 for grain number, 1 gene/QTL for tiller, HGW, grain filling and shattering. QTLS/genes for grain quality, biotic stresses and for abiotic stresses are 7, 23 and 13 respectively. Low yield, inferior quality and susceptibility to biotic and abiotic stresses of a crop is due to narrow genetic background of new evolving rice verities. Wild rice provides genetic resources for improvement of these characters, molecular and genomics tool at different stages can overcome these stresses and improve yield and quality of rice crop.

Genome-Wide Association Mapping of Crown and Brown Rust Resistance in Perennial Ryegrass

A population of 239 perennial ryegrass (Lolium perenne L.) genotypes was analyzed to identify marker-trait associations for crown rust (Puccinia coronata f. sp. lolii) and brown rust (Puccinia graminis f. sp. loliina) resistance. Phenotypic data from field trials showed a low correlation (r = 0.17) between the two traits. Genotypes were resequenced, and a total of 14,538,978 SNPs were used to analyze population structure, linkage disequilibrium (LD), and for genome-wide association study. The SNP heritability (h²_SNP) was 0.4 and 0.8 for crown and brown rust resistance, respectively. The high-density SNP dataset allowed us to estimate LD decay with the highest possible precision to date for perennial ryegrass. Results showed a low LD extension with a rapid decay of r² value below 0.2 after 520 bp on average. Additionally, QTL regions for both traits were detected, as well as candidate genes by applying Genome Complex Trait Analysis and Multi-marker Analysis of GenoMic Annotation. Moreover, two significant genes, LpPc6 and LpPl6, were identified for crown and brown rust resistance, respectively, when SNPs were aggregated to the gene level. The two candidate genes encode proteins with phosphatase activity, which putatively can be induced by the host to perceive, amplify and transfer signals to downstream components, thus activating a plant defense response.

Genomic resources in plant breeding for sustainable agriculture

Climate change during the last 40 years has had a serious impact on agriculture and threatens global food and nutritional security. From over half a million plant species, cereals and legumes are the most important for food and nutritional security. Although systematic plant breeding has a relatively short history, conventional breeding coupled with advances in technology and crop management strategies has increased crop yields by 56 % globally between 1965-85, referred to as the Green Revolution. Nevertheless, increased demand for food, feed, fiber, and fuel necessitates the need to break existing yield barriers in many crop plants. In the first decade of the 21st century we witnessed rapid discovery, transformative technological development and declining costs of genomics technologies. In the second decade, the field turned towards making sense of the vast amount of genomic information and subsequently moved towards accurately predicting gene-to-phenotype associations and tailoring plants for climate resilience and global food security. In this review we focus on genomic resources, genome and germplasm sequencing, sequencing-based trait mapping, and genomics-assisted breeding approaches aimed at developing biotic stress resistant, abiotic stress tolerant and high nutrition varieties in six major cereals (rice, maize, wheat, barley, sorghum and pearl millet), and six major legumes (soybean, groundnut, cowpea, common bean, chickpea and pigeonpea). We further provide a perspective and way forward to use genomic breeding approaches including marker-assisted selection, marker-assisted backcrossing, haplotype based breeding and genomic prediction approaches coupled with machine learning and artificial intelligence, to speed breeding approaches. The overall goal is to accelerate genetic gains and deliver climate resilient and high nutrition crop varieties for sustainable agriculture.

Isolation of Antimicrobial Genes from Oryza rufipogon Griff by Using a Bacillus subtilis Expression System with Potential Antimicrobial Activities

Antimicrobial genes are distributed in all forms of life and provide a primary defensive shield due to their unique broad-spectrum resistance activities. To better isolate these genes, we used the Bacillus subtilis expression system as the host cells to build Oryza rufipogon Griff cDNA libraries and screen potential candidate genes from the library at higher flux using built-in indicator bacteria. We observed that the antimicrobial peptides OrR214 and OrR935 have strong antimicrobial activity against a variety of Gram-positive and Gram-negative bacteria, as well as several fungal pathogens. Owing to their high thermal and enzymatic stabilities, these two peptides can also be used as field biocontrol agents. Furthermore, we also found that the peptide OrR214 (MIC 7.7-10.7 μM) can strongly inhibit bacterial growth compared to polymyxin B (MIC 5-25 μM) and OrR935 (MIC 33-44 μM). The cell flow analysis, reactive oxygen burst, and electron microscopy (scanning and transmission electron microscopy) observations showed that the cell membranes were targeted by peptides OrR214 and OrR935, which revealed the mode of action of bacteriostasis. Moreover, the hemolytic activity, toxicity, and salt sensitivity experiments demonstrated that these two peptides might have the potential to be used for clinical applications. Overall, OrR214 and OrR935 antimicrobial peptides have a high-throughput bacteriostatic activity that acts as a new form of antimicrobial agent and can be used as a raw material in the field of drug development.

Evaluation of Whole-Genome Sequence, Genetic Diversity, and Agronomic Traits of Basmati Rice (Oryza sativa L.)

Basmati is considered a unique varietal group of rice (Oryza sativa L.) because of its aroma and superior grain quality. Previous genetic analyses of rice showed that most of the Basmati varieties are classified into the aromatic group. Despite various efforts, genomic relationship of Basmati rice with other varietal groups and genomic variation in Basmati rice are yet to be understood. In the present study, we resequenced the whole genome of three traditional Basmati varieties at a coverage of more than 25X using Illumina HiSeq2500 and mapped the obtained sequences to the reference genome sequences of Nipponbare (japonica rice), Kasalath (aus rice), and Zhenshan 97 (indica rice). Comparison of these sequences revealed common single nucleotide polymorphisms (SNPs) in the genic regions of three Basmati varieties. Analysis of these SNPs revealed that Basmati varieties showed fewer sequence variations compared with the aus group than with the japonica and indica groups. Gene ontology (GO) enrichment analysis indicated that SNPs were present in genes with various biological, molecular, and cellular functions. Additionally, functional annotation of the Basmati mutated gene cluster shared by Nipponbare, Kasalath, and Zhenshan 97 was found to be associated with the metabolic process involved in the cellular aromatic compound, suggesting that aroma is an important specific genomic feature of Basmati varieties. Furthermore, 30 traditional Basmati varieties were classified into three different groups, aromatic (22 varieties), aus (four varieties), and indica (four varieties), based on genome-wide SNPs. All 22 aromatic Basmati varieties harbored the fragrant-inducing Badh2 allele. We also performed comparative analysis of 13 key agronomic and grain quality traits of Basmati rice and other rice varieties. Three traits including length-to-width ratio of grain (L/W ratio), panicle length (PL), and amylose content (AC) showed significant (P < 0.05 and P < 0.01) differences between the aromatic and indica/aus groups. Comparative analysis of genome structure, based on genome sequence variation and GO analysis, revealed that the Basmati genome was derived mostly from the aus and japonica groups. Overall, whole-genome sequence data and genetic diversity information obtained in this study will serve as an important resource for molecular breeding and genetic analysis of Basmati varieties.

Whole-genome resequencing analysis of 20 Micro-pigs

BACKGROUND

Miniature pigs have been increasingly used as mammalian model animals for biomedical research because of their similarity to human beings in terms of their metabolic features and proportional organ sizes. However, despite their importance, there is a severe lack of genome-wide studies on miniature pigs.

OBJECTIVE

In this study, we performed whole-genome sequencing analysis of 20 Micro-pigs obtained from Medi Kinetics to elucidate their genomic characteristics.

RESULTS

Approximately 595 gigabase pairs (Gb) of sequence reads were generated to be mapped to the swine reference genome assembly (Sus scrofa 10.2); on average, the sequence reads covered 99.15% of the reference genome at an average of 9.6-fold coverage. We detected a total of 19,518,548 SNPs, of which 8.7% were found to be novel. With further annotation of all of the SNPs, we retrieved 144,507 nonsynonymous SNPs (nsSNPs); of these, 5968 were found in all 20 individuals used in this study. SIFT prediction for these SNPs identified that 812 nsSNPs in 402 genes were deleterious. Among these 402 genes, we identified some genes that could potentially affect traits of interest in Micro-pigs, such as RHEB and FRAS1. Furthermore, we performed runs of homozygosity analysis to locate potential selection signatures in the genome, detecting several loci that might be involved in phenotypic characteristics in Micro-pigs, such as MSTN, GDF5, and GDF11.

CONCLUSION

In this study, we identified numerous nsSNPs that could be used as candidate genetic markers with involvement in traits of interest. Furthermore, we detected putative selection footprints that might be associated with recent selection applied to miniature pigs.

Discovery of Functional SNPs via Genome-Wide Exploration of Malaysian Pigmented Rice Varieties

Recently, rice breeding program has shown increased interests on the pigmented rice varieties due to their benefits to human health. However, the genetic variation of pigmented rice varieties is still scarce and remains unexplored. Hence, we performed genome-wide SNP analysis from the genome resequencing of four Malaysian pigmented rice varieties, representing two black and two red rice varieties. The genome of four pigmented varieties was mapped against Nipponbare reference genome sequences, and 1.9 million SNPs were discovered. Of these, 622 SNPs with polymorphic sites were identified in 258 protein-coding genes related to metabolism, stress response, and transporter. Comparative analysis of 622 SNPs with polymorphic sites against six rice SNP datasets from the Ensembl Plants variation database was performed, and 70 SNPs were identified as novel SNPs. Analysis of SNPs in the flavonoid biosynthetic genes revealed 40 nonsynonymous SNPs, which has potential as molecular markers for rice seed colour identification. The highlighted SNPs in this study show effort in producing valuable genomic resources for application in the rice breeding program, towards the genetic improvement of new and improved pigmented rice varieties.

Comparative study of the mycorrhizal root transcriptomes of wild and cultivated rice in response to the pathogen Magnaporthe oryzae

BACKGROUND

Rice, which serves as a staple food for more than half of the world's population, is very susceptible to the pathogenic fungus, Magnaporthe oryzae. However, common wild rice (Oryza rufipogon), which is the ancestor of Asian cultivated rice (O. sativa), has significant potential as a genetic source of resistance to M. oryzae. Recent studies have shown that the domestication of rice has altered its relationship to symbiotic arbuscular mycorrhizae. A comparative response of wild and domestic rice inhabited by mycorrhizae to infection by M. oryzae has not been documented.

RESULTS

In the current study, roots of wild and cultivated rice colonized with the arbuscular mycorrhizal (AM) fungus (AMF) Rhizoglomus intraradices were used to compare the transcriptomic responses of the two species to infection by M. oryzae. Phenotypic analysis indicated that the colonization of wild and cultivated rice with R. intraradices improved the resistance of both genotypes to M. oryzae. Wild AM rice, however, was more resistant to M. oryzae than the cultivated AM rice, as well as nonmycorrhizal roots of wild rice. Transcriptome analysis indicated that the mechanisms regulating the responses of wild and cultivated AM rice to M. oryzae invasion were significantly different. The expression of a greater number of genes was changed in wild AM rice than in cultivated AM rice in response to the pathogen. Both wild and cultivated AM rice exhibited a shared response to M. oryzae which included genes related to the auxin and salicylic acid pathways; all of these play important roles in pathogenesis-related protein synthesis. In wild AM rice, secondary metabolic and biotic stress-related analyses indicated that the jasmonic acid synthesis-related α-linolenic acid pathway, the phenolic and terpenoid pathways, as well as the phenolic and terpenoid syntheses-related mevalonate (MVA) pathway were more affected by the pathogen. Genes related to these pathways were more significantly enriched in wild AM rice than in cultivated AM rice in response to M. oryzae. On the other hand, genes associated with the 'brassinosteroid biosynthesis' were more enriched in cultivated AM rice.

CONCLUSIONS

The AMF R. intraradices-colonized rice plants exhibited greater resistance to M. oryzae than non-AMF-colonized plants. The findings of the current study demonstrate the potential effects of crop domestication on the benefits received by the host via root colonization with AMF(s), and provide new information on the underlying molecular mechanisms. In addition, results of this study can also help develop guidelines for the applications of AMF(s) when planting rice.

Re-sequencing and transcriptome analysis reveal rich DNA variations and differential expressions of fertility-related genes in neo-tetraploid rice

Autotetraploid rice is a useful germplasm for polyploid rice breeding, however, low seed setting is the major barrier in commercial utilization of autotetraploid rice. Our research group has developed neo-tetraploid rice lines, which have the characteristics of high fertility and heterosis when crossed with autotetraploid rice. In the present study, re-sequencing and RNA-seq were employed to detect global DNA variations and differentially expressed genes (DEGs) during meiosis stage in three neo-tetraploid rice lines compared to their parents, respectively. Here, a total of 4109881 SNPs and 640592 InDels were detected in neo-tetraploid lines compared to the reference genome, and 1805 specific presence/absence variations (PAVs) were detected in three lines. Approximately 12% and 0.5% of the total SNPs and InDels identified in three lines were located in genic regions, respectively. A total of 28 genes, harboring at least one of the large-effect SNP and/or InDel which affect the integrity of the encoded protein, were identified in the three lines. Together, 324 specific mutation genes, including 52 meiosis-related genes and 8 epigenetics-related genes were detected in neo-tetraploid rice compared to their parents. Of these 324 genes, five meiosis-related and three epigenetics-related genes displayed differential expressions during meiosis stage. Notably, 498 specific transcripts, 48 differentially expressed transposons and 245 differentially expressed ncRNAs were also detected in neo-tetraploid rice. Our results suggested that genomic structural reprogramming, DNA variations and differential expressions of some important meiosis and epigenetics related genes might be associated with high fertility in neo-tetraploid rice.

A catalog of annotated high-confidence SNPs from exome capture and sequencing reveals highly polymorphic genes in Norway spruce (Picea abies)

BACKGROUND

Norway spruce [Picea abies (L.) Karst.] is ecologically and economically one of the most important conifer worldwide. Our main goal was to develop a large catalog of annotated high confidence gene SNPs that should sustain the development of genomic tools for the conservation of natural and domesticated genetic diversity resources, and hasten tree breeding efforts in this species.

RESULTS

Targeted sequencing was achieved by capturing P. abies exome with probes previously designed from the sequenced transcriptome of white spruce (Picea glauca (Moench) Voss). Capture efficiency was high (74.5%) given a high level of exome conservation between the two species. Using stringent criteria, we delimited a set of 61,771 high-confidence SNPs across 13,543 genes. To validate SNPs, a high-throughput genotyping array was developed for a subset of 5571 predicted SNPs representing as many different gene loci, and was used to genotype over 1000 trees. The estimated true positive rate of the resource was 84.2%, which was comparable with the genotyping success rate obtained for P. abies control SNPs recycled from previous genotyping efforts. We also analyzed SNP abundance across various gene functional categories. Several GO terms and gene families involved in stress response were found over-represented in highly polymorphic genes.

CONCLUSION

The annotated high-confidence SNP catalog developed herein represents a valuable genomic resource, being representative of over 13 K genes distributed across the P. abies genome. This resource should serve a variety of population genomics and breeding applications in Norway spruce.

Genome-wide analysis of DNA polymorphisms, the methylome and transcriptome revealed that multiple factors are associated with low pollen fertility in autotetraploid rice

Autotetraploid rice is a useful germplasm with high biomass production; however, low fertility is the main barrier in commercial utilization. In our previous study, differential expression of meiosis-related miRNAs was found to be involved in the pollen sterility of autotetraploid rice. However, genome-wide DNA variations and methylomes associated with low fertility of autotetraploid rice are still poorly understood. Here, we measured both global DNA variations and the methylome and compared them with the transcriptome during pollen development in autotetraploid rice by high-throughput sequencing. A total of 34416 SNPs, 6993 InDels, 1003 SVs and 25 CNVs were detected, and 11367 and 41117 differentially methylated regions showed hypermethylation and hypomethylation in 02428-4x. In total, 1110 genes displayed differentially expression in 02428-4x during meiosis, of these six harbored CNVs, including four upregulated genes with gain CNVs, such as LOC_Os11g38620. We identified 122 genes by comparing with the previous data that might be associated with low fertility during pollen development in 02428-4x. Of the 122 gens, 98 were displayed methylation and differential expression, including OsMADS98, CYP703A3 and OsABCG26. The downregulation of these three genes were confirmed by qPCR during meiosis of 02428-4x, which played pivotal roles in pollen fertility. These results indicate that the low fertility of autotetraploid rice is not only caused by the differential expression of genes involved in pollen development, but also by sequence variation and differential methylation, suggesting that the reason for pollen sterility in autotetraploid rice is complex and might be affected by multiple factors.

Genome-Wide Analysis of Genetic Variations and the Detection of Rich Variants of NBS-LRR Encoding Genes in Common Wild Rice Lines

Common wild rice (Oryza rufipogon Griff.) is invaluable genetic resource for rice resistance breeding. Whole-genome re-sequencing was conducted to systematically analyze the variations in two new inbred lines (Huaye 3 and Huaye 4) developed from a common wild rice. A total of 4,841,127 SNPs, 1,170,479 InDels, 24,080 structural variations (SVs), and 298 copy number variations (CNVs) were identified in three materials. Approximately 16.24 and 5.64% of the total SNPs and InDels of Huaye 3 and Huaye 4 were located in genic regions, respectively. Together, 12,486 and 15,925 large-effect SNPs, and 12,417 and 14,513 large-effect InDels, which affect the integrity of the encoded protein, were identified in Huaye 3 and Huaye 4, respectively. The distribution map of 194 and 245 NBS-LRR encoding homologs was constructed across 12 rice chromosomes. Further, GO enrichment analysis of the homologs with identical genotype variations in Huaye 3 and Huaye 4 revealed 67, 82, and 58 homologs involved in cell death, response to stress, and both terms, respectively. Comparative analysis displayed that 550 out of 652 SNPs and 129 out of 147 InDels were present in a widely used blast-susceptible rice variety (LTH). Protein-protein interaction analysis revealed a strong interaction between NBS-LRR candidates and several known R genes. One homolog of disease resistance protein (RPM1) was involved in the plant-pathogen interaction pathway. Artificial inoculation of disease/insect displayed resistance phenotypes against rice blast and brown planthopper in two lines. The results will provide allele-specific markers for rice molecular breeding.

Differential expression of NBS-LRR-encoding genes in the root transcriptomes of two Solanum phureja genotypes with contrasting resistance to Globodera rostochiensis

BACKGROUND

The characterization of major resistance genes (R genes) in the potato remains an important task for molecular breeding. However, R genes are rapidly evolving and frequently occur in genomes as clusters with complex structures, and their precise mapping and identification are complicated and time consuming.

RESULTS

Comparative analysis of root transcriptomes of Solanum phureja genotypes with contrasting resistance to Globodera rostochiensis revealed a number of differentially expressed genes. However, compiling a list of candidate R genes for further segregation analysis was hampered by their scarce annotation. Nevertheless, combination of transcriptomic analysis with data on predicted potato NBS-LRR-encoding genes considerably improved the quality of the results and provided a reasonable number of candidate genes that provide S. phureja with strong resistance to the potato golden cyst nematode.

CONCLUSION

Combination of comparative analyses of tissue-specific transcriptomes in resistant and susceptible genotypes may be used as an approach for the rapid identification of candidate potato R genes for co-segregation analysis and may be used in parallel with more sophisticated studies based on genome resequencing.