Sequence composition of BAC clones and SSR markers mapped to Upland cotton chromosomes 11 and 21 targeting resistance to soil-borne pathogens

Small cardamom genome: development and utilization of microsatellite markers from a draft genome sequence of Elettaria cardamomum Maton

Small cardamom (Elettaria cardamomum Maton), the queen of spices, is the third most expensive spice in the world after saffron and vanilla, valued highly for its aroma and taste. This perennial herbaceous plant is a native of coastal parts of Southern India and displays a significant amount of morphological diversity. Its genetic potential has not been exploited due to lack of genomic resources limiting our understanding of the genome and important metabolic pathways which give it the economic advantage in the spice industry. Here, we report upon the de novo assembled, draft whole genome sequence of cardamom variety, Njallani Green Gold. We used a hybrid assembly strategy using the reads from the Oxford Nanopore, Illumina and 10x Genomics GemCode sequencing chemistries. The assembled genome length was 1.06 Gb (gigabases) which is close to the estimated genome size of cardamom. More than 75% of the genome was captured in 8000 scaffolds with a N50 of 0.15 Mb. The genome appears to have a high repeat content and 68055 gene models were predicted. The genome is close to Musa species and displays an expansion and contraction in different gene families. The draft assembly was used for in silico mining of simple sequence repeats (SSRs). A total of 2,50,571 SSRs were identified of which 2,18,270 were perfect SSRs and 32,301 were compound SSRs. Among the perfect SSRs, trinucleotides were most abundant (1,25,329) and hexanucleotide repeats appear least (2,380). From the 2,50,571 SSRs mined, 2,27,808 primer pairs were designed based on flanking sequence information. Wet lab validation was performed for 246 SSR loci and based on their amplification profiles, 60 SSR markers were used for diversity analysis of a set of 60 diverse cardamom accessions. The average number of alleles detected per locus were 14.57 with a minimum of 4 and maximum of 30 alleles. Population structure analysis revealed the presence of high degree of admixtures which could primarily be due to cross-pollination prevalent in this species. The SSR markers identified would help in the development of gene or trait-linked markers which can be subsequently used for marker-assisted breeding for crop improvement in cardamom. The information on utilization of the SSR loci for generation of markers has been developed into a public database, 'cardamomSSRdb' that is freely available for use by the cardamom community.

Mapping of Candidate Genes in Response to Low Nitrogen in Rice Seedlings

Nitrogen is not only a macronutrient essential for crop growth and development, but also one of the most critical nutrients in farmland ecosystem. Insufficient nitrogen supply will lead to crop yield reduction, while excessive application of nitrogen fertilizer will cause agricultural and eco-environment damage. Therefore, mining low-nitrogen tolerant rice genes and improving nitrogen use efficiency are of great significance to the sustainable development of agriculture. This study was conducted by Genome-wide association study on a basis of two root morphological traits (root length and root diameter) and 788,396 SNPs of a natural population of 295 rice varieties. The transcriptome of low-nitrogen tolerant variety (Longjing 31) and low-nitrogen sensitive variety (Songjing 10) were sequenced between low and high nitrogen treatments. A total of 35 QTLs containing 493 genes were mapped. 3085 differential expressed genes were identified. Among these 493 genes, 174 genes showed different haplotype patterns. There were significant phenotype differences among different haplotypes of 58 genes with haplotype differences. These 58 genes were hypothesized as candidate genes for low nitrogen tolerance related to root morphology. Finally, six genes (Os07g0471300, Os11g0230400, Os11g0229300, Os11g0229400, Os11g0618300 and Os11g0229333) which expressed differentially in Longjing 31 were defined as more valuable candidate genes for low-nitrogen tolerance. The results revealed the response characteristics of rice to low-nitrogen, and provided insights into regulatory mechanisms of rice to nitrogen deficiency.

Mechanism of cotton resistance to abiotic stress, and recent research advances in the osmoregulation related genes

Abiotic stress is an important factor affecting the normal growth and development of plants and crop yield. To reduce the impact of abiotic adversity on cotton growth and development, the material basis of cotton resistance and its physiological functions are analyzed at the molecular level. At the same time, the use of genetic engineering methods to recombine resistance genes has become a hot spot in cotton resistance research. This paper provides an overviews of the resistance mechanism of cotton against the threat of non-biological adversity, as well as the research progress of osmoregulation-related genes, protein-acting genes, and transcription regulatory factor genes in recent years, and outlines the explored gene resources in cotton resistance genetic engineering, with the aim to provide ideas and reference bases for future research on cotton resistance.

Transgressive resistance to Heterodera glycines in chromosome segment substitution lines derived from susceptible soybean parents

Chromosome segment substitution lines (CSSLs) are valuable genetic resources for quantitative trait loci (QTL) mapping of complex agronomic traits especially suitable for minor effect QTL. Here, 162 BC3 F7 -BC7 F3 CSSLs derived from crossing two susceptible parent lines, soybean [Glycine max (L.) Merr.] 'Suinong14' (recurrent parent) × wild soybean (G. soja Siebold & Zucc.) ZYD00006, were used for QTL mapping of soybean cyst nematode (SCN, Heterodera glycine Ichinohe) resistance based on female index (FI) and cysts per gram root (CGR) through phenotypic screening and whole-genome resequencing of CSSLs. Phenotypic results displayed a wide range of distribution and transgressive lines in both HG Type 2.5.7 FI and CGR and demonstrated a higher correlation between CGR and root weight (R2 = .5424) compared with than between FI and CGR (R2 = .0018). Using the single-marker analysis nonparametric mapping test, 33 significant QTL were detected on 18 chromosomes contributing resistance to FI and CGR. Fourteen QTL contributing 5.6-15.5% phenotypic variance (PVE) to FI were revealed on 11 chromosomes, and 16 QTL accounting for 6.1-36.2% PVE in CGR were detected on 14 chromosomes with strong additive effect by multiple-QTL model (MQM) mapping. Twenty-five and 13 out of all 38 QTL identified for FI and CGR on 20 chromosomes were from ZYD00006 and Suinong14, respectively. The CSSLs with the combination of positive alleles for FI, CGR, and root weight exhibited low nematode reproduction. For the first time, QTL associated with CGR have been detected, and both FI and CGR should be considered for breeding purposes in the absence of strong resistance genes such as rhg1 and Rhg4.

Genome-wide identification and expression analysis of Raffinose synthetase family in cotton

Background

The Raffinose synthetase (RAFS) genes superfamily is critical for the synthesis of raffinose, which accumulates in plant leaves under abiotic stress. However, it remains unclear whether RAFS contributes to resistance to abiotic stress in plants, specifically in the Gossypium species.

Results

In this study, we identified 74 RAFS genes from G. hirsutum, G. barbadense, G. arboreum and G. raimondii by using a series of bioinformatic methods. Phylogenetic analysis showed that the RAFS gene family in the four Gossypium species could be divided into four major clades; the relatively uniform distribution of the gene number in each species ranged from 12 to 25 based on species ploidy, most likely resulting from an ancient whole-genome polyploidization. Gene motif analysis showed that the RAFS gene structure was relatively conservative. Promoter analysis for cis-regulatory elements showed that some RAFS genes might be regulated by gibberellins and abscisic acid, which might influence their expression levels. Moreover, we further examined the functions of RAFS under cold, heat, salt and drought stress conditions, based on the expression profile and co-expression network of RAFS genes in Gossypium species. Transcriptome analysis suggested that RAFS genes in clade III are highly expressed in organs such as seed, root, cotyledon, ovule and fiber, and under abiotic stress in particular, indicating the involvement of genes belonging to clade III in resistance to abiotic stress. Gene co-expressed network analysis showed that GhRFS2A-GhRFS6A, GhRFS6D, GhRFS7D and GhRFS8A-GhRFS11A were key genes, with high expression levels under salt, drought, cold and heat stress.

Conclusion

The findings may provide insights into the evolutionary relationships and expression patterns of RAFS genes in Gossypium species and a theoretical basis for the identification of stress resistance materials in cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04276-4.

Association Mapping of Verticillium Wilt Disease in a Worldwide Collection of Cotton (Gossypium hirsutum L.)

Cotton (Gossypium spp.) is the best plant fiber source in the world and provides the raw material for industry. Verticillium wilt caused by Verticillium dahliae Kleb. is accepted as a major disease of cotton production. The most practical way to deal with verticillium wilt is to develop resistant/tolerant varieties after cultural practices. One of the effective selections in plant breeding is the use of marker-assisted selection (MAS) via quantitative trait loci (QTL). Therefore, in this study, we aimed to discover the genetic markers associated with the disease. Through the association mapping analysis, common single nucleotide polymorphism (SNP) markers were obtained using 4730 SNP alleles. As a result, twenty-three markers were associated with defoliating (PYDV6 isolate) pathotype, twenty-one markers with non-defoliating (Vd11 isolate) pathotype, ten QTL with Disease Severity Index (DSI) of the leaves at the 50–60% boll opening period and eight markers were associated with DSI in the stem section. Some of the markers that show significant associations are located on protein coding genes such as protein Mpv17-like, 21 kDa protein-like, transcription factor MYB113-like, protein dehydration-induced 19 homolog 3-like, F-box protein CPR30-like, extracellular ribonuclease LE-like, putative E3 ubiquitin-protein ligase LIN, pentatricopeptide repeat-containing protein At3g62890-like, fructose-1,6-bisphosphatase, tubby-like F-box protein 8, endoglucanase 16-like, glucose-6-phosphate/phosphate translocator 2, metal tolerance protein 11-like, VAN3-binding protein-like, transformation/transcription domain-associated protein-like, pyruvate kinase isozyme A, ethylene-responsive transcription factor CRF2-like, molybdate transporter 2-like, IRK-interacting protein-like, glycosylphosphatidylinositol anchor attachment 1 protein, U3 small nucleolar RNA-associated protein 4-like, microtubule-associated protein futsch-like, transport and Golgi organization 2 homolog, splicing factor 3B subunit 3-like, mediator of RNA polymerase II transcription subunit 15a-like, putative ankyrin repeat protein, and protein networked 1D-like. It has been reported in previous studies that most of these genes are associated with biotic and abiotic stress factors. As a result, once validated, it would be possible to use the markers obtained in the study in Marker Assisted Selection (MAS) breeding.

Sequence Composition of Bacterial Chromosome Clones in a Transgressive Root-Knot Nematode Resistance Chromosome Region in Tetraploid Cotton

Plants evolve innate immunity including resistance genes to defend against pest and pathogen attack. Our previous studies in cotton (Gossypium spp.) revealed that one telomeric segment on chromosome (Chr) 11 in G. hirsutum cv. Acala NemX (rkn1 locus) contributed to transgressive resistance to the plant parasitic nematode Meloidogyne incognita, but the highly homologous segment on homoeologous Chr 21 had no resistance contribution. To better understand the resistance mechanism, a bacterial chromosome (BAC) library of Acala N901 (Acala NemX resistance source) was used to select, sequence, and analyze BAC clones associated with SSR markers in the complex rkn1 resistance region. Sequence alignment with the susceptible G. hirsutum cv. TM-1 genome indicated that 23 BACs mapped to TM-1-Chr11 and 18 BACs mapped to TM-1-Chr 21. Genetic and physical mapping confirmed less BAC sequence (53-84%) mapped with the TM-1 genome in the rkn1 region on Chr 11 than to the homologous region (>89%) on Chr 21. A 3.1-cM genetic distance between the rkn1 flanking markers CIR316 and CIR069 was mapped in a Pima S-7 × Acala NemX RIL population with a physical distance ∼1 Mbp in TM-1. NCBI Blast and Gene annotation indicated that both Chr 11 and Chr 21 harbor resistance gene-rich cluster regions, but more multiple homologous copies of Resistance (R) proteins and of adjacent transposable elements (TE) are present within Chr 11 than within Chr 21. (CC)-NB-LRR type R proteins were found in the rkn1 region close to CIR316, and (TIR)-NB-LRR type R proteins were identified in another resistance rich region 10 cM from CIR 316 (∼3.1 Mbp in the TM-1 genome). The identified unique insertion/deletion in NB-ARC domain, different copies of LRR domain, multiple copies or duplication of R proteins, adjacent protein kinases, or TE in the rkn1 region on Chr 11 might be major factors contributing to complex recombination and transgressive resistance.

Genome-wide association study of Gossypium arboreum resistance to reniform nematode

BACKGROUND

Reniform nematode (Rotylenchulus reniformis) has emerged as one of the most destructive root pathogens of upland cotton (Gossypium hirsutum) in the United States. Management of R. reniformis has been hindered by the lack of resistant G. hirsutum cultivars; however, resistance has been frequently identified in germplasm accessions from the G. arboreum collection. To determine the genetic basis of reniform nematode resistance, a genome-wide association study (GWAS) was performed using 246 G. arboreum germplasm accessions that were genotyped with 7220 single nucleotide polymorphic (SNP) sequence markers generated from genotyping-by-sequencing.

RESULTS

Fifteen SNPs representing 12 genomic loci distributed over eight chromosomes showed association with reniform nematode resistance. For 14 SNPs, major alleles were shown to be associated with resistance. From the 15 significantly associated SNPs, 146 genes containing or physically close to these loci were identified as putative reniform nematode resistance candidate genes. These genes are involved in a broad range of biological pathways, including plant innate immunity, transcriptional regulation, and redox reaction that may have a role in the expression of resistance. Eighteen of these genes corresponded to differentially expressed genes identified from G. hirsutum in response to reniform nematode infection.

CONCLUSIONS

The identification of multiple genomic loci associated with reniform nematode resistance would indicate that the G. arboreum collection is a significant resource of novel resistance genes. The significantly associated markers identified from this GWAS can be used for the development of molecular tools for breeding improved reniform nematode resistant upland cotton with resistance introgressed from G. arboreum. Additionally, a greater understanding of the molecular mechanisms of reniform nematode resistance can be determined through genetic structure and functional analyses of candidate genes, which will aid in the pyramiding of multiple resistance genes.

Quantitative Trait Loci Mapping of Multiple Independent Loci for Resistance to Fusarium oxysporum f. sp. vasinfectum Races 1 and 4 in an Interspecific Cotton Population

Fusarium wilt, caused by the soilborne fungal pathogen Fusarium oxysporum f. sp. vasinfectum, is a vascular disease of cotton (Gossypium spp.). F. oxysporum f. sp. vasinfectum race 1 (FOV1) causes major plant injury and yield loss in G. hirsutum cultivars with coinfection with root-knot nematode (Meloidogyne incognita), while F. oxysporum f. sp. vasinfectum race 4 (FOV4) causes plant damage without nematode coinfection in G. hirsutum and in G. barbadense cultivars. Quantitative trait loci (QTL) analysis of the interspecific cross G. barbadense Pima S-7 × G. hirsutum Acala NemX revealed separate multiple loci determining resistance to FOV1 and FOV4, confirming that race specificity occurs in F. oxysporum f. sp. vasinfectum. Based on the area under the disease progress stairs, six major QTLs on chromosomes (Chrs) 1, 2, 12, 15 (2), and 21 contributing 7 to 15% to FOV1 resistance and two major QTLs on Chrs 14 and 17 contributing 12 to 33% to FOV4 resistance were identified. Minor-effect QTLs contributing to resistance to both FOV1 and FOV4 were also identified. These results define and establish a pathosystem of race-specific resistance under polygenic control. This research also validates the importance of previously reported markers and chromosome regions and adds new information for the location of F. oxysporum f. sp. vasinfectum resistance genes. Some F₈ recombinant inbred lines have resistance to both FOV1 and FOV4 and also to root-knot nematode, providing multiple resistance sources for breeding.

QTL Analysis of Transgressive Nematode Resistance in Tetraploid Cotton Reveals Complex Interactions in Chromosome 11 Regions

Transgressive segregation in cotton (Gossypium spp.) provides an important approach to enhance resistance to the major pest root-knot nematode (RKN) Meloidogyne incognita. Our previous studies reported transgressive RKN resistance in an intraspecific Gossypium hirsutum resistant NemX × susceptible SJ-2 recombinant inbred line (RIL) population and early generations of interspecific cross Gossypium barbadense (susceptible Pima S-7) × G. hirsutum (NemX). However, the underlying functional mechanisms for this phenomenon are not known. In this study, the region of RKN resistance gene rkn1 on chromosome (Chr) 11 and its homoeologous Chr 21 was fine mapped with G. raimondii D₅ genome reference sequence. Transgressive resistance was found in the later generation of a new RIL population F_2:7 (Pima S-7 × NemX) and one interspecific F₂ (susceptible Pima S-7 × susceptible SJ-2). QTL analysis revealed similar contributions to root-galling and egg-production resistance phenotypes associated with SSR marker CIR316 linked to resistance gene rkn1 in NemX on Chr 11 in all seven populations analyzed. In testcross NemX × F₁ (Pima S-7 × SJ-2) marker allele CIR069-271 from Pima S-7 linked to CIR316 contributed 63% of resistance to galling phenotype in the presence of rkn1. Similarly, in RIL population F_2:8 (NemX × SJ-2), SJ-2 markers closely linked to CIR316 contributed up to 82% of resistance to root-galling. These results were confirmed in BC₁F₁ SJ-2 × F₁ (NemX × SJ-2), F₂ (NemX × SJ-2), and F₂ (Pima S-7 × SJ-2) populations in which up to 44, 36, and 15% contribution in resistance to galling was found, respectively. Transgressive segregation for resistance was universal in all intra- and inter-specific populations, although stronger transgressive resistance occurred in later than in early generations in the intraspecific cross compared with the interspecific cross. Transgressive effects on progeny from susceptible parents are possibly provided in the rkn1 resistance region of chromosome 11 by tandemly arrayed allele (TAA) or gene (TAG) interactions contributing to transgressive resistance. Complex TAA and TAG recombination and interactions in the rkn1 resistance region provide three genes and a model to study disease and transgressive resistance in polyploid plants, and novel genotypes for plant breeding.