The plasticity of NBS resistance genes in sorghum is driven by multiple evolutionary processes

A Sorghum BAK1/SERK4 Homolog Functions in Pathogen-Associated Molecular Patterns-Triggered Immunity and Cell Death in Response to Colletotrichum sublineola Infection

Sorghum bicolor is the fifth most important cereal crop and expected to gain prominence due to its versatility, low input requirements, and tolerance to hot and dry conditions. In warm and humid environments, the productivity of sorghum is severely limited by the hemibiotrophic fungal pathogen Colletotrichum sublineola, the causal agent of anthracnose. Cultivating anthracnose-resistant accessions is the most effective and environmentally benign way to safeguard yield. A previous genome-wide association study for anthracnose resistance in the Sorghum Association Panel uncovered single-nucleotide polymorphisms on chromosome 5 associated with resistance to anthracnose, including one located within the coding region of gene Sobic.005G182400. In this study, we investigated the molecular function of Sobic.005G182400 in response to C. sublineola infection. Conserved domain, phylogenetic, and structural analyses revealed that the protein encoded by Sobic.005G182400 shares significant structural similarity with the Arabidopsis BRASSINOSTEROID INSENSITIVE1-ASSOCIATED RECEPTOR KINASE1 (BAK1)/SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE4 (SERK4). Although sequence analysis of four sorghum accessions showed no substantial variation in the coding region, accession SC1330, which carries the resistance allele, exhibited significantly higher expression of Sobic.005G182400 during early infection (≤24 h). Co-expression network analysis identified that the module associated with Sobic.005G182400 was enriched in genes involved in endocytosis, autophagy, and vesicle transport. Gene regulatory network analysis further suggested that Sobic.005G182400 regulates genes required for BAK1/SERK4-mediated cell death via protein glycosylation. Together, these findings indicate that Sobic.005G182400 encodes a protein with similarity to Arabidopsis BAK1/SERK4 that enables pathogen-associated molecular patterns-triggered immunity and regulates cell death.

Haplotypes at the sorghum ARG4 and ARG5 NLR loci confer resistance to anthracnose

Sorghum anthracnose caused by the fungus Colletotrichum sublineola (Cs) is a damaging disease of the crop. Here, we describe the identification of ANTHRACNOSE RESISTANCE GENES (ARG4 and ARG5) encoding canonical nucleotide-binding leucine-rich repeat (NLR) receptors. ARG4 and ARG5 are dominant resistance genes identified in the sorghum lines SAP135 and P9830, respectively, that show broad-spectrum resistance to Cs. Independent genetic studies using populations generated by crossing SAP135 and P9830 with TAM428, fine mapping using molecular markers, comparative genomics and gene expression studies determined that ARG4 and ARG5 are resistance genes against Cs strains. Interestingly, ARG4 and ARG5 are both located within clusters of duplicate NLR genes at linked loci separated by ~1 Mb genomic region. SAP135 and P9830 each carry only one of the ARG genes while having the recessive allele at the second locus. Only two copies of the ARG5 candidate genes were present in the resistant P9830 line while five non-functional copies were identified in the susceptible line. The resistant parents and their recombinant inbred lines carrying either ARG4 or ARG5 are resistant to strains Csgl1 and Csgrg suggesting that these genes have overlapping specificities. The role of ARG4 and ARG5 in resistance was validated through sorghum lines carrying independent recessive alleles that show increased susceptibility. ARG4 and ARG5 are located within complex loci displaying interesting haplotype structures and copy number variation that may have resulted from duplication. Overall, the identification of anthracnose resistance genes with unique haplotype stucture provides a foundation for genetic studies and resistance breeding.

Multi-genome comprehensive identification of SSR/SV and development of molecular markers database to serve Sorghum bicolor (L.) breeding

BACKGROUND

As an important food and cash crop, identification of DNA molecular markers is of great significance for molecular marker-assisted breeding of Sorghum (Sorghum bicolor (L.) moench). Although some sorghum-related mutation databases have been published, the special SSR and SV databases still need to be constructed and updated.

RESULTS

In this study, the quality of 18 different sorghum genomes was evaluated, and two genomes were assembled at chromosome level. Through the identification and comparative analysis of SSR loci in these genomes, the distribution characteristics of SSR in the above sorghum genomes were initially revealed. At the same time, five representative reference genomes were selected to identify the structural variation of sorghum. Finally, a convenient SSR/SV database of sorghum was constructed by integrating the above results ( http://www.sorghum.top:8079/ ; http://43.154.129.150:8079/ ; http://47.106.184.91:8079/ ). Users can query the information of related sites and primer pairs.

CONCLUSIONS

Anyway, our research provides convenience for sorghum researchers and will play an active role in sorghum molecular marker-assisted breeding.

PearMODB: a multiomics database for pear (Pyrus) genomics, genetics and breeding study

Pear (Pyrus ssp.) belongs to Rosaceae and is an important fruit tree widely cultivated around the world. Currently, challenges to cope with the burgeoning sets of multiomics data are rapidly increasing. Here, we constructed the Pear Multiomics Database (PearMODB) by integrating genome, transcriptome, epigenome and population variation data, and aimed to provide a portal for accessing and analyzing pear multiomics data. A variety of online tools were built including gene search, BLAST, JBrowse, expression heatmap, synteny analysis and primer design. The information of DNA methylation sites and single-nucleotide polymorphisms can be retrieved through the custom JBrowse, providing an opportunity to explore the genetic polymorphisms linked to phenotype variation. Moreover, different gene families involving transcription factors, transcription regulators and disease resistance (nucleotide-binding site leucine-rich repeat) were identified and compiled for quick search. In particular, biosynthetic gene clusters (BGCs) were identified in pear genomes, and specialized webpages were set up to show detailed information of BGCs, laying a foundation for studying metabolic diversity among different pear varieties. Overall, PearMODB provides an important platform for pear genomics, genetics and breeding studies. Database URL http://pearomics.njau.edu.cn.

The role of gene presence-absence variations on genetic incompatibility in Asian rice

Genetic incompatibilities are widespread between species. However, it remains unclear whether they all originated after population divergence as suggested by the Bateson-Dobzhansky-Muller model, and if not, what is their prevalence and distribution within populations. The gene presence-absence variations (PAVs) provide an opportunity for investigating gene-gene incompatibility. Here, we searched for the repulsion of coexistence between gene PAVs to identify the negative interaction of gene functions separately in two Oryza sativa subspecies. Many PAVs are involved in subspecies-specific negative epistasis and segregate at low-to-intermediate frequencies in focal subspecies but at low or high frequencies in the other subspecies. Incompatible PAVs are enriched in two functional groups, defense response and protein phosphorylation, which are associated with plant immunity and consistent with autoimmunity being a known mechanism of hybrid incompatibility in plants. Genes in the two enriched functional groups are older and seldom directly interact with each other. Instead, they interact with other younger gene PAVs with diverse functions. Our results illustrate the landscape of genetic incompatibility at gene PAVs in rice, where many incompatible pairs have already segregated as polymorphisms within subspecies, and many are novel negative interactions between older defense-related genes and younger genes with diverse functions.

Bioinformatics analysis and function prediction of NBS-LRR gene family in Broussonetia papyrifera

Most of the currently available disease resistance (R) genes have NBS (nucleotide-binding site) and LRR (leucine-rich-repeat) domain which belongs to the NBS-LRR gene family. The whole genome sequencing of Broussonetia papyrifera provides an important bioinformatics database for the study of the NBS-LRR gene family. In this study, 328 NBS-LRR family genes were identified and classified in B. papyrifera according to different classification schemes, where there are 92 N types, 47 CN type, 54 CNL type, 29 NL types, 55 TN type, and 51 TNL type. Subsequently, we conducted bioinformatics analysis of the NBS-LRR gene family. Classification, motif analysis of protein sequences, and phylogenetic tree studies of the NBS-LRR genes in B. papyrifera provide important basis for the functional study of NBS-LRR family genes. Additionally, we performed structural analysis of the chromosomal location, physicochemical properties, and sequences identified by genetic characterization. In addition, through the analysis of GO enrichment, it was found that NBS-LRR genes were involved in defense responses and were significantly enriched in biological stimulation, immune response, and abiotic stress. In addition, we found that Bp06g0955 was the most sensitive to low temperature and encoded the RPM1 protein by analyzing the low temperature transcriptome data of B. papyrifera. Quantitative results of gene expression after 48 h of Fusarium infection showed that Bp01g3293 increased 14 times after infection, which encodes RPM1 protein. The potential of NBS-LRR gene responsive to biotic and abiotic stresses can be exploited to improve the resistance of B. papyrifera.

Genome-wide identification, characterization, and evolutionary analysis of NBS genes and their association with disease resistance in Musa spp

Banana is an important food crop that is susceptible to a wide range of pests and diseases that can reduce yield and quality. The primary objective of banana breeding programs is to increase disease resistance, which requires the identification of resistance (R) genes. Despite the fact that resistant sources have been identified in bananas, the genes, particularly the nucleotide-binding site (NBS) family, which play an important role in protecting plants against pathogens, have received little attention. As a result, this study included a thorough examination of the NBS disease resistance gene family's classification, phylogenetic analysis, genome organization, evolution, cis-elements, differential expression, regulation by microRNAs, and protein-protein interaction. A total of 116 and 43 putative NBS genes from M. acuminata and M. balbisiana, respectively, were identified and characterized, and were classified into seven sub-families. Structural analysis of NBS genes revealed the presence of signal peptides, their sub-cellular localization, molecular weight and pI. Eight commonly conserved motifs were found, and NBS genes were unevenly distributed across multiple chromosomes, with the majority of NBS genes being located in chr3 and chr1 of the A and B genomes, respectively. Tandem duplication occurrences have helped bananas' NBS genes spread throughout evolution. Transcriptome analysis of NBS genes revealed significant differences in expression between resistant and susceptible cultivars of fusarium wilt, eumusae leaf spot, root lesion nematode, and drought, implying that they can be used as candidate resistant genes. Ninety miRNAs were discovered to have targets in 104 NBS genes from the A genome, providing important insights into NBS gene expression regulation. Overall, this study offers a valuable genomic resource and understanding of the function and evolution of NBS genes in relation to rapidly evolving pathogens, as well as providing breeders with selection targets for fast-tracking breeding of banana varieties with more durable resistance to pathogens.

Multilevel evolution shapes the function of NB-LRR encoding genes in plant innate immunity

A sophisticated innate immune system based on diverse pathogen receptor genes (PRGs) evolved in the history of plant life. To reconstruct the direction and magnitude of evolutionary trajectories of a given gene family, it is critical to detect the ancestral signatures. The rearrangement of functional domains made up the diversification found in PRG repertoires. Structural rearrangement of ancient domains mediated the NB-LRR evolutionary path from an initial set of modular proteins. Events such as domain acquisition, sequence modification and temporary or stable associations are prominent among rapidly evolving innate immune receptors. Over time PRGs are continuously shaped by different forces to find their optimal arrangement along the genome. The immune system is controlled by a robust regulatory system that works at different scales. It is important to understand how the PRG interaction network can be adjusted to meet specific needs. The high plasticity of the innate immune system is based on a sophisticated functional architecture and multi-level control. Due to the complexity of interacting with diverse pathogens, multiple defense lines have been organized into interconnected groups. Genomic architecture, gene expression regulation and functional arrangement of PRGs allow the deployment of an appropriate innate immunity response.

Chromosome-scale assembly and analysis of biomass crop Miscanthus lutarioriparius genome

Miscanthus, a rhizomatous perennial plant, has great potential for bioenergy production for its high biomass and stress tolerance. We report a chromosome-scale assembly of Miscanthus lutarioriparius genome by combining Oxford Nanopore sequencing and Hi-C technologies. The 2.07-Gb assembly covers 96.64% of the genome, with contig N50 of 1.71 Mb. The centromere and telomere sequences are assembled for all 19 chromosomes and chromosome 10, respectively. Allotetraploid origin of the M. lutarioriparius is confirmed using centromeric satellite repeats. The tetraploid genome structure and several chromosomal rearrangements relative to sorghum are clearly demonstrated. Tandem duplicate genes of M. lutarioriparius are functional enriched not only in terms related to stress response, but cell wall biosynthesis. Gene families related to disease resistance, cell wall biosynthesis and metal ion transport are greatly expanded and evolved. The expansion of these families may be an important genomic basis for the enhancement of remarkable traits of M. lutarioriparius.

Understanding the Sorghum-Colletotrichum sublineola Interactions for Enhanced Host Resistance

Improving sorghum resistance is a sustainable method to reduce yield losses due to anthracnose, a devastating disease caused by Colletotrichum sublineola. Elucidating the molecular mechanisms of sorghum-C. sublineola interactions would help identify biomarkers for rapid and efficient identification of novel sources for host-plant resistance improvement, understanding the pathogen virulence, and facilitating resistance breeding. Despite concerted efforts to identify resistance sources, the knowledge about sorghum-anthracnose interactions remains scanty. Hence, in this review, we presented an overview of the current knowledge on the mechanisms of sorghum-C. sublineola molecular interactions, sources of resistance for sorghum breeding, quantitative trait loci (QTL), and major (R-) resistance gene sequences as well as defense-related genes associated with anthracnose resistance. We summarized current knowledge about C. sublineola populations and its virulence. Illustration of the sorghum-C. sublineola interaction model based on the current understanding is also provided. We highlighted the importance of genomic resources of both organisms for integrated omics research to unravel the key molecular components underpinning compatible and incompatible sorghum-anthracnose interactions. Furthermore, sorghum-breeding strategy employing rapid sorghum germplasm screening, systems biology, and molecular tools is presented.

Elucidating Anthracnose Resistance Mechanisms in Sorghum-A Review

Sorghum (Sorghum bicolor) is the fifth most cultivated cereal crop in the world, traditionally providing food, feed, and fodder, but more recently also fermentable sugars for the production of renewable fuels and chemicals. The hemibiotrophic fungal pathogen Colletotrichum sublineola, the causal agent of anthracnose disease in sorghum, is prevalent in the warm and humid climates where much of the sorghum is cultivated and poses a serious threat to sorghum production. The use of anthracnose-resistant sorghum germplasm is the most environmentally and economically sustainable way to protect sorghum against this pathogen. Even though multiple anthracnose resistance loci have been mapped in diverse sorghum germplasm in recent years, the diversity in C. sublineola pathotypes at the local and regional levels means that these resistance genes are not equally effective in different areas of cultivation. This review summarizes the genetic and cytological data underlying sorghum's defense response and describes recent developments that will enable a better understanding of the interactions between sorghum and C. sublineola at the molecular level. This includes releases of the sorghum genome and the draft genome of C. sublineola, the use of next-generation sequencing technologies to identify gene expression networks activated in response to infection, and improvements in methodologies to validate resistance genes, notably virus-induced and transgenic gene silencing approaches.

Contrasting genetic variation and positive selection followed the divergence of NBS-encoding genes in Asian and European pears

BACKGROUND

The NBS disease-related gene family coordinates the inherent immune system in plants in response to pathogen infections. Previous studies have identified NBS-encoding genes in Pyrus bretschneideri ('Dangshansuli', an Asian pear) and Pyrus communis ('Bartlett', a European pear) genomes, but the patterns of genetic variation and selection pressure on these genes during pear domestication have remained unsolved.

RESULTS

In this study, 338 and 412 NBS-encoding genes were identified from Asian and European pear genomes. This difference between the two pear species was the result of proximal duplications. About 15.79% orthologous gene pairs had Ka/Ks ratio more than one, indicating two pear species undergo strong positive selection after the divergence of Asian and European pear. We identified 21 and 15 NBS-encoding genes under fire blight and black spot disease-related QTL, respectively, suggesting their importance in disease resistance. Domestication caused decreased nucleotide diversity across NBS genes in Asian cultivars (cultivated 6.23E-03; wild 6.47E-03), but opposite trend (cultivated 6.48E-03; wild 5.91E-03) appeared in European pears. Many NBS-encoding coding regions showed Ka/Ks ratio of greater than 1, indicating the role of positive selection in shaping diversity of NBS-encoding genes in pear. Furthermore, we detected 295 and 122 significantly different SNPs between wild and domesticated accessions in Asian and European pear populations. Two NBS genes (Pbr025269.1 and Pbr019876.1) with significantly different SNPs showed >5x upregulation between wild and cultivated pear accessions, and > 2x upregulation in Pyrus calleryana after inoculation with Alternaria alternata. We propose that positively selected and significantly different SNPs of an NBS-encoding gene (Pbr025269.1) regulate gene expression differences in the wild and cultivated groups, which may affect resistance in pear against A. alternata.

CONCLUSION

Proximal duplication mainly led to the different number of NBS-encoding genes in P. bretschneideri and P. communis genomes. The patterns of genetic diversity and positive selection pressure differed between Asian and European pear populations, most likely due to their independent domestication events. This analysis helps us understand the evolution, diversity, and selection pressure in the NBS-encoding gene family in Asian and European populations, and provides opportunities to study mechanisms of disease resistance in pear.

Identification and Distribution of NBS-Encoding Resistance Genes of Dactylis glomerata L. and Its Expression Under Abiotic and Biotic Stress

Orchardgrass (Dactylis glomerata L.) is drought resistant and tolerant to barren landscapes, making it one of the most important forages for animal husbandry, as well as ecological restoration of rocky landscapes that are undergoing desertification. However, orchardgrass is susceptible to rust, which can significantly reduce its yield and quality. Therefore, understanding the genes that underlie resistance against rust in orchardgrass is critical. The evolution, cloning of plant disease resistance genes, and the analysis of pathogenic bacteria induced expression patterns are important contents in the study of interaction between microorganisms and plants. Genes with nucleotide binding site (NBS) structure are disease-resistant genes ubiquitous in plants and play an important role in plant attacks against various pathogens. Using sequence analysis and re-annotation, we identified 413 NBS resistance genes in orchardgrass. Similar to previous studies, NBS resistance genes containing TIR (toll/interleukin-1 receptor) domain were not found in orchardgrass. The NBS resistance genes can be divided into four types: NBS (up to 264 homologous genes, accounting for 64% of the total number of NBS genes in orchardgrass), NBS-LRR, CC-NBS, and CC-NBS-LRR (minimum of 26 homologous genes, only 6% of the total number of NBS genes in orchardgrass). These 413 NBS resistance genes were unevenly distributed across seven chromosomes where chromosome 5 had up to 99 NBS resistance genes. There were 224 (54%) NBS resistance genes expressed in different tissues (roots, stems, leaves, flowers, and spikes), and we did not detect expression for 45 genes (11%). The remaining 145 (35%) were expressed in some tissues. And we found that 11 NBS resistance genes were differentially expressed under waterlogging stress, 5 NBS resistance genes were differentially expressed under waterlogging and drought stress, and 1 NBS resistance was is differentially expressed under waterlogging and heat stress. Most importantly, we found that 65 NBS resistance genes were significantly expressed in different control groups. On the 7th day of inoculation, 23 NBS resistance genes were differentially expressed in high resistance materials alone, of which 7 NBS resistance genes regulate the "plant-pathogen interaction" pathway by encoding RPM1. At the same time, 2 NBS resistance genes that were differentially expressed in the high resistance material after inoculation were also differentially expressed in abiotic stress. In summary, the NBS resistance gene plays a crucial role in the resistance of orchardgrass to rust.

Genome-Wide Association Mapping of Anthracnose (Colletotrichum sublineolum) Resistance in NPGS Ethiopian Sorghum Germplasm

The National Plant Germplasm System (NPGS) Ethiopian sorghum [Sorghum bicolor (L.) Moench] collection of the United States is an important genetic resource for sorghum improvement. Anthracnose (Colletotrichum sublineolum) is one of the most harmful fungal diseases in humid sorghum production regions. Although multiple resistance sources have been identified in temperate-adapted germplasm in the Sorghum Association Panel (SAP), these resistance loci explain a limited portion of the total variation, and sources of resistance from tropical germplasm are not available for breeding programs at temperate regions. Using a core set of 335 previously genotyped NPGS Ethiopian accessions, we identified 169 accessions resistant to anthracnose. To identify resistance loci, we merged the genotypic and anthracnose response data for both NPGS Ethiopian germplasm and the SAP and performed genome-wide association scans using 219,037 single nucleotide polymorphisms and 617 accessions. The integrated data set enabled the detection of a locus on chromosome 9 present in the SAP at a low frequency. The locus explains a limited portion of the observed phenotypic variation (r² = 0.31), suggesting the presence of other resistance loci. The locus in chromosome 9 was constituted by three R genes clustered within a 47-kb region. The presence of multiple sources of resistance in NPGS Ethiopian germplasm and SAP requires the inclusion of other resistance response evaluation that could revealed others low frequency resistance alleles in the panel.

A new reference genome for Sorghum bicolor reveals high levels of sequence similarity between sweet and grain genotypes: implications for the genetics of sugar metabolism

BACKGROUND

The process of crop domestication often consists of two stages: initial domestication, where the wild species is first cultivated by humans, followed by diversification, when the domesticated species are subsequently adapted to more environments and specialized uses. Selective pressure to increase sugar accumulation in certain varieties of the cereal crop Sorghum bicolor is an excellent example of the latter; this has resulted in pronounced phenotypic divergence between sweet and grain-type sorghums, but the genetic mechanisms underlying these differences remain poorly understood.

RESULTS

Here we present a new reference genome based on an archetypal sweet sorghum line and compare it to the current grain sorghum reference, revealing a high rate of nonsynonymous and potential loss of function mutations, but few changes in gene content or overall genome structure. We also use comparative transcriptomics to highlight changes in gene expression correlated with high stalk sugar content and show that changes in the activity and possibly localization of transporters, along with the timing of sugar metabolism play a critical role in the sweet phenotype.

CONCLUSIONS

The high level of genomic similarity between sweet and grain sorghum reflects their historical relatedness, rather than their current phenotypic differences, but we find key changes in signaling molecules and transcriptional regulators that represent new candidates for understanding and improving sugar metabolism in this important crop.

Genome wide association analysis of sorghum mini core lines regarding anthracnose, downy mildew, and head smut

In previous studies, a sorghum mini core collection was scored over several years for response to Colletotrichum sublineola, Peronosclerospora sorghi, and Sporisorium reilianum, the causal agents of the disease anthracnose, downy mildew, and head smut, respectively. The screening results were combined with over 290,000 Single nucleotide polymorphic (SNP) loci from an updated version of a publicly available genotype by sequencing (GBS) dataset available for the mini core collection. GAPIT (Genome Association and Prediction Integrated Tool) R package was used to identify chromosomal locations that differ in disease response. When the top scoring SNPs were mapped to the most recent version of the published sorghum genome, in each case, a nearby and most often the closest annotated gene has precedence for a role in host defense.

The Sorghum QTL Atlas: a powerful tool for trait dissection, comparative genomics and crop improvement

We describe the development and application of the Sorghum QTL Atlas, a high-resolution, open-access research platform to facilitate candidate gene identification across three cereal species, sorghum, maize and rice. The mechanisms governing the genetic control of many quantitative traits are only poorly understood and have yet to be fully exploited. Over the last two decades, over a thousand QTL and GWAS studies have been published in the major cereal crops including sorghum, maize and rice. A large body of information has been generated on the genetic basis of quantitative traits, their genomic location, allelic effects and epistatic interactions. However, such QTL information has not been widely applied by cereal improvement programs and genetic researchers worldwide. In part this is due to the heterogeneous nature of QTL studies which leads QTL reliability variation from study to study. Using approaches to adjust the QTL confidence interval, this platform provides access to the most updated sorghum QTL information than any database available, spanning 23 years of research since 1995. The QTL database provides information on the predicted gene models underlying the QTL CI, across all sorghum genome assembly gene sets and maize and rice genome assemblies and also provides information on the diversity of the underlying genes and information on signatures of selection in sorghum. The resulting high-resolution, open-access research platform facilitates candidate gene identification across 3 cereal species, sorghum, maize and rice. Using a number of trait examples, we demonstrate the power and resolution of the resource to facilitate comparative genomics approaches to provide a bridge between genomics and applied breeding.

Connecting genome structural variation with complex traits in crop plants

Structural genome variation is a major determinant of useful trait diversity. We describe how genome analysis methods are enabling discovery of trait-associated structural variants and their potential impact on breeding. As our understanding of complex crop genomes continues to grow, there is growing evidence that structural genome variation plays a major role in determining traits important for breeding and agriculture. Identifying the extent and impact of structural variants in crop genomes is becoming increasingly feasible with ongoing advances in the sophistication of genome sequencing technologies, particularly as it becomes easier to generate accurate long sequence reads on a genome-wide scale. In this article, we discuss the origins of structural genome variation in crops from ancient and recent genome duplication and polyploidization events and review high-throughput methods to assay such variants in crop populations in order to find associations with phenotypic traits. There is increasing evidence from such studies that gene presence-absence and copy number variation resulting from segmental chromosome exchanges may be at the heart of adaptive variation of crops to counter abiotic and biotic stress factors. We present examples from major crops that demonstrate the potential of pangenomic diversity as a key resource for future plant breeding for resilience and sustainability.

Segmental and Tandem Duplications Driving the Recent NBS-LRR Gene Expansion in the Asparagus Genome

Garden asparagus is an important horticultural plant worldwide. It is, however, susceptible to a variety of diseases, which can affect the potential yield, spear quality, and lifespan of production fields. Screening studies have identified resistant germplasm. The genetic resistance is usually complex, and the genes underlying that resistance are still unknown. Most often, disease resistance is determined by resistance genes (R). The most predominant R-genes contain nucleotide binding site and leucine-rich repeat (NBS-LRR) domains. Using bioinformatics and data mining approaches, we identified and characterized 68 NBS predicted proteins encoded by 49 different loci in the asparagus genome. The NBS-encoding genes were grouped into seven distinct classes based on their domain architecture. The NBS genes are unevenly distributed through the genome and nearly 50% of the genes are present in clusters. Chromosome 6 is significantly NBS-enriched and one single cluster hosts 10% of the genes. Phylogenetic analysis points to their diversification into three families during their evolution. Recent duplications are likely to have dominated the NBS expansion with both tandem genes and duplication events across multiple chromosomes. Transcriptome sequencing data provided evidence for their transcription and tissue-specific expression. The total number of cis-regulatory elements as well as their relative positions within the NBS promoters suggests a complex transcriptional network regulating defense responses. Our study provides a strong groundwork for the isolation of candidate R-genes in garden asparagus.

Evaluation of Genetic Variation among Sorghum Varieties from Southwest China via Genome Resequencing

Little is known regarding genomic variation among glutinous sorghum [ (L.) Moench] varieties grown in southwest China, which are primarily used to brew the popular Jiang-flavor liquor. This study evaluated genomic variation among six representative sorghum accessions via whole-genome resequencing. The evaluation revealed 2365,363 single-nucleotide polymorphisms (SNPs), 394,365 insertions and deletions, and 47,567 copy number variations among the six genomes. Chromosomes 5 and 10 showed relatively high SNP densities, whereas whole-genome diversity in this population was low. In addition, some chromosomal loci exhibited obvious selection during the breeding process. Sorghum accessions from southwest China formed an elite germplasm population compared with the findings of other geographic populations, and the elite variety 'Hongyingzi' contained 79 unique genes primarily involved in basic metabolism. The six sorghum lines contained a large number of high-confidence genes, with Hongyingzi in particular possessing 104 unique genes. These findings advance our understanding of domestication of the sorghum genome, and Chinese sorghum accessions will be valuable resources for further research and breeding improvements.

Genome-wide scans of selection highlight the impact of biotic and abiotic constraints in natural populations of the model grass Brachypodium distachyon

Grasses are essential plants for ecosystem functioning. Quantifying the selective pressures that act on natural variation in grass species is therefore essential regarding biodiversity maintenance. In this study, we investigate the selection pressures that act on two distinct populations of the grass model Brachypodium distachyon without prior knowledge about the traits under selection. We took advantage of whole-genome sequencing data produced for 44 natural accessions of B. distachyon and used complementary genome-wide selection scans (GWSS) methods to detect genomic regions under balancing and positive selection. We show that selection is shaping genetic diversity at multiple temporal and spatial scales in this species, and affects different genomic regions across the two populations. Gene ontology annotation of candidate genes reveals that pathogens may constitute important factors of positive and balancing selection in B. distachyon. We eventually cross-validated our results with quantitative trait locus data available for leaf-rust resistance in this species and demonstrate that, when paired with classical trait mapping, GWSS can help pinpointing candidate genes for further molecular validation. Thanks to a near base-perfect reference genome and the large collection of freely available natural accessions collected across its natural range, B. distachyon appears as a prime system for studies in ecology, population genomics and evolutionary biology.

Rapid identification of candidate genes for resistance to tomato late blight disease using next-generation sequencing technologies

Tomato late blight caused by Phytophthora infestans (Mont.) de Bary, also known as the Irish famine pathogen, is one of the most destructive plant diseases. Wild relatives of tomato possess useful resistance genes against this disease, and could therefore be used in breeding to improve cultivated varieties. In the genome of a wild relative of tomato, Solanum habrochaites accession LA1777, we identified a new quantitative trait locus for resistance against blight caused by an aggressive Egyptian isolate of P. infestans. Using double-digest restriction site-associated DNA sequencing (ddRAD-Seq) technology, we determined 6,514 genome-wide SNP genotypes of an F2 population derived from an interspecific cross. Subsequent association analysis of genotypes and phenotypes of the mapping population revealed that a 6.8 Mb genome region on chromosome 6 was a candidate locus for disease resistance. Whole-genome resequencing analysis revealed that 298 genes in this region potentially had functional differences between the parental lines. Among of them, two genes with missense mutations, Solyc06g071810.1 and Solyc06g083640.3, were considered to be potential candidates for disease resistance. SNP and SSR markers linking to this region can be used in marker-assisted selection in future breeding programs for late blight disease, including introgression of new genetic loci from wild species. In addition, the approach developed in this study provides a model for identification of other genes for attractive agronomical traits.

Copy number variation and disease resistance in plants

Plant genome diversity varies from single nucleotide polymorphisms to large-scale deletions, insertions, duplications, or re-arrangements. These re-arrangements of sequences resulting from duplication, gains or losses of DNA segments are termed copy number variations (CNVs). During the last decade, numerous studies have emphasized the importance of CNVs as a factor affecting human phenotype; in particular, CNVs have been associated with risks for several severe diseases. In plants, the exploration of the extent and role of CNVs in resistance against pathogens and pests is just beginning. Since CNVs are likely to be associated with disease resistance in plants, an understanding of the distribution of CNVs could assist in the identification of novel plant disease-resistance genes. In this paper, we review existing information about CNVs; their importance, role and function, as well as their association with disease resistance in plants.

Genome-Wide Association Study Identifies NBS-LRR-Encoding Genes Related with Anthracnose and Common Bacterial Blight in the Common Bean

Nucleotide-binding site and leucine-rich repeat (NBS-LRR) genes represent the largest and most important disease resistance genes in plants. The genome sequence of the common bean (Phaseolus vulgaris L.) provides valuable data for determining the genomic organization of NBS-LRR genes. However, data on the NBS-LRR genes in the common bean are limited. In total, 178 NBS-LRR-type genes and 145 partial genes (with or without a NBS) located on 11 common bean chromosomes were identified from genome sequences database. Furthermore, 30 NBS-LRR genes were classified into Toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) types, and 148 NBS-LRR genes were classified into coiled-coil (CC)-NBS-LRR (CNL) types. Moreover, the phylogenetic tree supported the division of these PvNBS genes into two obvious groups, TNL types and CNL types. We also built expression profiles of NBS genes in response to anthracnose and common bacterial blight using qRT-PCR. Finally, we detected nine disease resistance loci for anthracnose (ANT) and seven for common bacterial blight (CBB) using the developed NBS-SSR markers. Among these loci, NSSR24, NSSR73, and NSSR265 may be located at new regions for ANT resistance, while NSSR65 and NSSR260 may be located at new regions for CBB resistance. Furthermore, we validated NSSR24, NSSR65, NSSR73, NSSR260, and NSSR265 using a new natural population. Our results provide useful information regarding the function of the NBS-LRR proteins and will accelerate the functional genomics and evolutionary studies of NBS-LRR genes in food legumes. NBS-SSR markers represent a wide-reaching resource for molecular breeding in the common bean and other food legumes. Collectively, our results should be of broad interest to bean scientists and breeders.

Genome Analysis Identified Novel Candidate Genes for Ascochyta Blight Resistance in Chickpea Using Whole Genome Re-sequencing Data

Ascochyta blight (AB) is a fungal disease that can significantly reduce chickpea production in Australia and other regions of the world. In this study, 69 chickpea genotypes were sequenced using whole genome re-sequencing (WGRS) methods. They included 48 Australian varieties differing in their resistance ranking to AB, 16 advanced breeding lines from the Australian chickpea breeding program, four landraces, and one accession representing the wild chickpea species Cicer reticulatum. More than 800,000 single nucleotide polymorphisms (SNPs) were identified. Population structure analysis revealed relatively narrow genetic diversity amongst recently released Australian varieties and two groups of varieties separated by the level of AB resistance. Several regions of the chickpea genome were under positive selection based on Tajima's D test. Both Fst genome- scan and genome-wide association studies (GWAS) identified a 100 kb region (AB4.1) on chromosome 4 that was significantly associated with AB resistance. The AB4.1 region co-located to a large QTL interval of 7 Mb∼30 Mb identified previously in three different mapping populations which were genotyped at relatively low density with SSR or SNP markers. The AB4.1 region was validated by GWAS in an additional collection of 132 advanced breeding lines from the Australian chickpea breeding program, genotyped with approximately 144,000 SNPs. The reduced level of nucleotide diversity and long extent of linkage disequilibrium also suggested the AB4.1 region may have gone through selective sweeps probably caused by selection of the AB resistance trait in breeding. In total, 12 predicted genes were located in the AB4.1 QTL region, including those annotated as: NBS-LRR receptor-like kinase, wall-associated kinase, zinc finger protein, and serine/threonine protein kinases. One significant SNP located in the conserved catalytic domain of a NBS-LRR receptor-like kinase led to amino acid substitution. Transcriptional analysis using qPCR showed that some predicted genes were significantly induced in resistant lines after inoculation compared to non-inoculated plants. This study demonstrates the power of combining WGRS data with relatively simple traits to rapidly develop "functional makers" for marker-assisted selection and genomic selection.

Domestication and the storage starch biosynthesis pathway: signatures of selection from a whole sorghum genome sequencing strategy

Next-generation sequencing of complete genomes has given researchers unprecedented levels of information to study the multifaceted evolutionary changes that have shaped elite plant germplasm. In conjunction with population genetic analytical techniques and detailed online databases, we can more accurately capture the effects of domestication on entire biological pathways of agronomic importance. In this study, we explore the genetic diversity and signatures of selection in all predicted gene models of the storage starch synthesis pathway of Sorghum bicolor, utilizing a diversity panel containing lines categorized as either 'Landraces' or 'Wild and Weedy' genotypes. Amongst a total of 114 genes involved in starch synthesis, 71 had at least a single signal of purifying selection and 62 a signal of balancing selection and others a mix of both. This included key genes such as STARCH PHOSPHORYLASE 2 (SbPHO2, under balancing selection), PULLULANASE (SbPUL, under balancing selection) and ADP-glucose pyrophosphorylases (SHRUNKEN2, SbSH2 under purifying selection). Effectively, many genes within the primary starch synthesis pathway had a clear reduction in nucleotide diversity between the Landraces and wild and weedy lines indicating that the ancestral effects of domestication are still clearly identifiable. There was evidence of the positional rate variation within the well-characterized primary starch synthesis pathway of sorghum, particularly in the Landraces, whereby low evolutionary rates upstream and high rates downstream in the metabolic pathway were expected. This observation did not extend to the wild and weedy lines or the minor starch synthesis pathways.

Whole-Genome Resequencing of a Cucumber Chromosome Segment Substitution Line and Its Recurrent Parent to Identify Candidate Genes Governing Powdery Mildew Resistance

Cucumber is an economically important vegetable crop worldwide. Powdery mildew (PM) is one of the most severe diseases that can affect cucumber crops. There have been several research efforts to isolate PM resistance genes for breeding PM-resistant cucumber. In the present study, we used a chromosome segment substitution line, SSL508-28, which carried PM resistance genes from the donor parent, JIN5-508, through twelve generations of backcrossing with a PM-susceptible inbred line, D8. We performed whole-genome resequencing of SSL508-28 and D8 to identify single nucleotide polymorphisms (SNPs), and insertions and deletions (indels). When compared against the reference genome of the inbred cucumber line 9930, a total of 468,616 SNPs and 67,259 indels were identified in SSL508-28, and 537,352 SNPs and 91,698 indels were identified in D8. Of these, 3,014 non-synonymous SNPs and 226 frameshift indels in SSL508-28, and 3,104 non-synonymous SNPs and 251 frameshift indels in D8, were identified. Bioinformatics analysis of these variations revealed a total of 15,682 SNPs and 6,262 indels between SSL508-28 and D8, among which 120 non-synonymous SNPs and 30 frameshift indels in 94 genes were detected between SSL508-28 and D8. Finally, out of these 94 genes, five resistance genes with nucleotide-binding sites and leucine-rich repeat domains were selected for qRT-PCR analysis. This revealed an upregulation of two transcripts, Csa2M435460.1 and Csa5M579560.1, in SSL508-28. Furthermore, the results of qRT-PCR analysis of these two genes in ten PM resistant and ten PM susceptible cucumber lines showed that when exposed to PM, Csa2M435460.1 and Csa5M579560.1 exhibited a higher expression level of resistant lines than susceptible lines. This indicates that Csa2M435460.1 and Csa5M579560.1 are candidate genes for PM resistance in cucumber. In addition, the non-synonymous SNPs in Csa2M435460.1 and Csa5M579560.1, identified in SSL508-28 and D8, might be the key to high PM-resistance in SSL508-28.

Whole genome sequencing of a banana wild relative Musa itinerans provides insights into lineage-specific diversification of the Musa genus

Crop wild relatives are valuable resources for future genetic improvement. Here, we report the de novo genome assembly of Musa itinerans, a disease-resistant wild banana relative in subtropical China. The assembled genome size was 462.1 Mb, covering 75.2% of the genome (615.2Mb) and containing 32, 456 predicted protein-coding genes. Since the approximate divergence around 5.8 million years ago, the genomes of Musa itinerans and Musa acuminata have shown conserved collinearity. Gene family expansions and contractions enrichment analysis revealed that some pathways were associated with phenotypic or physiological innovations. These include a transition from wood to herbaceous in the ancestral Musaceae, intensification of cold and drought tolerances, and reduced diseases resistance genes for subtropical marginally distributed Musa species. Prevalent purifying selection and transposed duplications were found to facilitate the diversification of NBS-encoding gene families for two Musa species. The population genome history analysis of M. itinerans revealed that the fluctuated population sizes were caused by the Pleistocene climate oscillations, and that the formation of Qiongzhou Strait might facilitate the population downsizing on the isolated Hainan Island about 10.3 Kya. The qualified assembly of the M. itinerans genome provides deep insights into the lineage-specific diversification and also valuable resources for future banana breeding.

SorghumFDB: sorghum functional genomics database with multidimensional network analysis

Sorghum (Sorghum bicolor [L.] Moench) has excellent agronomic traits and biological properties, such as heat and drought-tolerance. It is a C₄ grass and potential bioenergy-producing plant, which makes it an important crop worldwide. With the sorghum genome sequence released, it is essential to establish a sorghum functional genomics data mining platform. We collected genomic data and some functional annotations to construct a sorghum functional genomics database (SorghumFDB). SorghumFDB integrated knowledge of sorghum gene family classifications (transcription regulators/factors, carbohydrate-active enzymes, protein kinases, ubiquitins, cytochrome P450, monolignol biosynthesis related enzymes, R-genes and organelle-genes), detailed gene annotations, miRNA and target gene information, orthologous pairs in the model plants Arabidopsis, rice and maize, gene loci conversions and a genome browser. We further constructed a dynamic network of multidimensional biological relationships, comprised of the co-expression data, protein–protein interactions and miRNA-target pairs. We took effective measures to combine the network, gene set enrichment and motif analyses to determine the key regulators that participate in related metabolic pathways, such as the lignin pathway, which is a major biological process in bioenergy-producing plants.

Database URL:http://structuralbiology.cau.edu.cn/sorghum/index.html.

Directional Selection from Host Plants Is a Major Force Driving Host Specificity in Magnaporthe Species

One major threat to global food security that requires immediate attention, is the increasing incidence of host shift and host expansion in growing number of pathogenic fungi and emergence of new pathogens. The threat is more alarming because, yield quality and quantity improvement efforts are encouraging the cultivation of uniform plants with low genetic diversity that are increasingly susceptible to emerging pathogens. However, the influence of host genome differentiation on pathogen genome differentiation and its contribution to emergence and adaptability is still obscure. Here, we compared genome sequence of 6 isolates of Magnaporthe species obtained from three different host plants. We demonstrated the evolutionary relationship between Magnaporthe species and the influence of host differentiation on pathogens. Phylogenetic analysis showed that evolution of pathogen directly corresponds with host divergence, suggesting that host-pathogen interaction has led to co-evolution. Furthermore, we identified an asymmetric selection pressure on Magnaporthe species. Oryza sativa-infecting isolates showed higher directional selection from host and subsequently tends to lower the genetic diversity in its genome. We concluded that, frequent gene loss or gain, new transposon acquisition and sequence divergence are host adaptability mechanisms for Magnaporthe species, and this coevolution processes is greatly driven by directional selection from host plants.

Comparative Analysis and Identification of miRNAs and Their Target Genes Responsive to Salt Stress in Diploid and Tetraploid Paulownia fortunei Seedlings

Salt stress is a global environmental problem that affects plant growth and development. Paulownia fortunei is an adaptable and fast-growing deciduous tree native to China that is environmentally and economically important. MicroRNAs (miRNAs) play important regulatory roles in growth, development, and stress responses in plants. MiRNAs that respond to biotic stresses have been identified; however, how miRNAs in P. fortunei respond to salt stress has not yet been reported. To identify salt-stress-responsive miRNAs and predict their target genes, four small RNA and four degradome libraries were constructed from NaCl-treated and NaCl-free leaves of P. fortunei seedlings. The results indicated that salt stress had different physiological effects on diploid and tetraploid P. fortunei. We detected 53 conserved miRNAs belonging to 17 miRNA families and 134 novel miRNAs in P. fortunei. Comparing their expression levels in diploid and tetraploid P. fortunei, we found 10 conserved and 10 novel miRNAs that were significantly differentially expressed under salt treatment, among them eight were identified as miRNAs probably associated with higher salt tolerance in tetraploid P. fortunei than in diploid P. fortunei. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed to predict the functions of the target genes of the conserved and novel miRNAs. The expressions of 10 differentially expressed miRNAs were validated by quantitative real-time polymerase chain reaction (qRT-PCR). This is the first report on P. fortunei miRNAs and their target genes under salt stress. The results provided information at the physiological and molecular levels for further research into the response mechanisms of P. fortunei to salt stress.

Genome-Wide Analysis of NBS-LRR Genes in Sorghum Genome Revealed Several Events Contributing to NBS-LRR Gene Evolution in Grass Species

The nucleotide-binding site (NBS)–leucine-rich repeat (LRR) gene family is crucially important for offering resistance to pathogens. To explore evolutionary conservation and variability of NBS-LRR genes across grass species, we identified 88, 107, 24, and 44 full-length NBS-LRR genes in sorghum, rice, maize, and Brachypodium, respectively. A comprehensive analysis was performed on classification, genome organization, evolution, expression, and regulation of these NBS-LRR genes using sorghum as a representative of grass species. In general, the full-length NBS-LRR genes are highly clustered and duplicated in sorghum genome mainly due to local duplications. NBS-LRR genes have basal expression levels and are highly potentially targeted by miRNA. The number of NBS-LRR genes in the four grass species is positively correlated with the gene clustering rate. The results provided a valuable genomic resource and insights for functional and evolutionary studies of NBS-LRR genes in grass species.

Disease Resistance Gene Analogs (RGAs) in Plants

Plants have developed effective mechanisms to recognize and respond to infections caused by pathogens. Plant resistance gene analogs (RGAs), as resistance (R) gene candidates, have conserved domains and motifs that play specific roles in pathogens' resistance. Well-known RGAs are nucleotide binding site leucine rich repeats, receptor like kinases, and receptor like proteins. Others include pentatricopeptide repeats and apoplastic peroxidases. RGAs can be detected using bioinformatics tools based on their conserved structural features. Thousands of RGAs have been identified from sequenced plant genomes. High-density genome-wide RGA genetic maps are useful for designing diagnostic markers and identifying quantitative trait loci (QTL) or markers associated with plant disease resistance. This review focuses on recent advances in structures and mechanisms of RGAs, and their identification from sequenced genomes using bioinformatics tools. Applications in enhancing fine mapping and cloning of plant disease resistance genes are also discussed.

Two distinct classes of QTL determine rust resistance in sorghum

BACKGROUND

Agriculture is facing enormous challenges to feed a growing population in the face of rapidly evolving pests and pathogens. The rusts, in particular, are a major pathogen of cereal crops with the potential to cause large reductions in yield. Improving stable disease resistance is an on-going major and challenging focus for many plant breeding programs, due to the rapidly evolving nature of the pathogen. Sorghum is a major summer cereal crop that is also a host for a rust pathogen Puccinia purpurea, which occurs in almost all sorghum growing areas of the world, causing direct and indirect yield losses in sorghum worldwide, however knowledge about its genetic control is still limited. In order to further investigate this issue, QTL and association mapping methods were implemented to study rust resistance in three bi-parental populations and an association mapping set of elite breeding lines in different environments.

RESULTS

In total, 64 significant or highly significant QTL and 21 suggestive rust resistance QTL were identified representing 55 unique genomic regions. Comparisons across populations within the current study and with rust QTL identified previously in both sorghum and maize revealed a high degree of correspondence in QTL location. Negative phenotypic correlations were observed between rust, maturity and height, indicating a trend for both early maturing and shorter genotypes to be more susceptible to rust.

CONCLUSIONS

The significant amount of QTL co-location across traits, in addition to the consistency in the direction of QTL allele effects, has provided evidence to support pleiotropic QTL action across rust, height, maturity and stay-green, supporting the role of carbon stress in susceptibility to rust. Classical rust resistance QTL regions that did not co-locate with height, maturity or stay-green QTL were found to be significantly enriched for the defence-related NBS-encoding gene family, in contrast to the lack of defence-related gene enrichment in multi-trait effect rust resistance QTL. The distinction of disease resistance QTL hot-spots, enriched with defence-related gene families from QTL which impact on development and partitioning, provides plant breeders with knowledge which will allow for fast-tracking varieties with both durable pathogen resistance and appropriate adaptive traits.