RIL Population Research Articles

An efficient genomic prediction method without the direct inverse of the genomic relationship matrix.

GBLUP, the most widely used genomic prediction (GP) method, consumes large and increasing amounts of computational resources as the training population size increases due to the inverse of the genomic relationship matrix (GRM). Therefore, in this study, we developed a new genomic prediction method (RHEPCG) that avoids the direct inverse of the GRM by combining randomized Haseman-Elston (HE) regression (RHE-reg) and a preconditioned conjugate gradient (PCG). The simulation results demonstrate that RHEPCG, in most cases, not only achieves similar predictive accuracy with GBLUP but also significantly reduces computational time. As for the real data, RHEPCG shows similar or better predictive accuracy for seven traits of the Arabidopsis thaliana F2 population and four traits of the Sorghum bicolor RIL population compared with GBLUP. This indicates that RHEPCG is a practical alternative to GBLUP and has better computational efficiency.

Leaf rolling in bread wheat (Triticum aestivum L.) is controlled by the upregulation of a pair of closely linked/duplicate zinc finger homeodomain class transcription factors during moisture stress conditions.

Zinc finger-homeodomain (ZF-HDs) class IV transcriptional factors (TFs) is a plant-specific transcription factor and play a key role in stress responses, plant growth, development, and hormonal signaling. In this study, two new leaf rolling TFs genes, namely TaZHD1 and TaZHD10, were identified in wheat using comparative genomic analysis of the target region that carried a major QTL for leaf rolling identified through multi-environment phenotyping and high throughput genotyping of a RIL population. Structural and functional annotation of the candidate ZHD genes with its closest rice orthologs reflects the species-specific evolution and, undoubtedly, validates the notions of remote-distance homology concept. Meanwhile, the morphological analysis resulted in contrasting difference for leaf rolling in extreme RILs between parental lines HD2012 and NI5439 at booting and heading stages. Transcriptome-wide expression profiling revealed that TaZHD10 transcripts showed significantly higher expression levels than TaZHD1 in all leaf tissues upon drought stress. The relative expression of these genes was further validated by qRT-PCR analysis, which also showed consistent results across the studied genotypes at the booting and anthesis stage. The contrasting modulation of these genes under drought conditions and the available evidenced for its epigenetic behavior that might involve the regulation of metabolic and gene regulatory networks. Prediction of miRNAs resulted in five Tae-miRs that could be associated with RNAi mediated control of TaZHD1 and TaZHD10 putatively involved in the metabolic pathway controlling rolled leaf phenotype. Gene interaction network analysis indicated that TaZHD1 and TaZHD10 showed pleiotropic effects and might also involve other functions in wheat in addition to leaf rolling. Overall, the results increase our understanding of TaZHD genes and provide valuable information as robust candidate genes for future functional genomics research aiming for the breeding of wheat varieties tolerant to leaf rolling.

Identification of drought tolerant entries based on stress tolerant indices and physiological traits in RIL population of cotton (Gossypium hirsutum)

Identification of drought tolerant entries based on stress tolerant indices and physiological traits in RIL population of cotton (Gossypium hirsutum)

Bulked segregant transcriptome analysis in pea identifies key expression markers for resistance to Peyronellaea pinodes

Peyronellaea pinodes is a devastating pathogen of pea crop. Quantitative trait loci (QTL) associated with resistance have been identified, as well as genes differentially expressed between resistant and susceptible pea lines. The key question is which of these many genes located into these QTLs, or differentially expressed, are the key genes that distinguish resistant from susceptible plants and could be used as markers. To identify these key genes, in the present study we applied MACE (Massive Analysis of cDNA Ends) -Seq to a whole Recombinant Inbred Line population segregating for resistance to this disease and their parental lines and identified those genes which expression was more correlated with the level of resistance. We also compared gene expression profiles between the most resistant and the most susceptible families of the RIL population. A total of 6780 transcripts were differentially expressed between the parental lines after inoculation. Of them, 803 showed the same expression pattern in the bulks formed by the most resistant and most susceptible RIL families. These genes, showing a consistent expression pattern, could be used as expression markers to distinguish resistant from susceptible plants. The analysis of these genes also discovered the crucial mechanisms acting against P. pinodes.

Identification of QTLs and joint QTL segments of leaflet traits at different canopy layers in an interspecific RIL population of soybean.

A leaflet trait on different canopy layers may have different QTLs; leaflet trait QTLs may cluster to form joint QTL segments; all canopy layer QTLs form a complete QTL system for a leaflet trait. As the main part of the plant canopy structure, leaf/leaflet size and shape affect the plant architecture and yield. To explore the leaflet trait QTL system, a population composed of 199 recombinant inbred lines derived from Changling (annual wild, narrow leaflet) and Yiqianli (landrace, broad leaflet) with their parents was tested for leaflet length (LL), width (LW) and length to width (LLW). The population was genotyped with specific-locus amplified fragment sequencing (SLAF-seq) and applied for linkage mapping of the leaflet traits. The results showed that the leaflet traits varied greatly even within a plant, which supported a stratified leaflet sampling strategy to evaluate these traits at top, middle and bottom canopy layers. Altogether, 13 LL, 10 LW and 9 LLW in a total of 32 plus 3 duplicated QTLs were identified, in which, 17 QTLs were new ones, and 48.6%, 28.6% and 22.8% of QTLs were from the top, middle and bottom layers, respectively, indicating the genetic importance of the top layer leaves. Since a leaflet trait may have layer-specific QTLs, all layer QTLs form a complete QTL system. Five QTL clusters each with their QTL supporting intervals overlapped were designated as joint QTL segments (JQSs). In JQS-16, with its linkage map further validated using PCR markers, two QTLs, qLW-16-1 and qLLW-16-1 of the top layer leaflet, were identified six QTL·times. Six candidate genes were predicted, with Glyma.16G127900 as the most potential one for LW and LLW. Three PCR markers, Gm16PAV0653, BARCSOYSSR_16_0796 and YC-16-3, were suggested for marker-assisted selection for LW and LLW in JQS-16.

Discovery and Validation of Grain Shape Loci in U.S. Rice Germplasm Through Haplotype Characterization.

Rice grain shape is a major determinant of rice market value and the end-use. We mapped quantitative trait loci (QTL) for grain shape traits in a bi-parental recombinant inbred line population (Trenasse/Jupiter) and discovered two major grain length QTLs—qGL3.1 and qGL7.1. Previously, a major grain shape gene GS3 was reported in the qGL3.1 region and grain length gene GL7 was reported to be encompassing qGL7.1 locus. The re-sequencing SNP data on the International Rice Research Institute (IRRI) 3K Rice Genome Project (RGP) panel were obtained from the IRRI SNP-Seek database for both genes and haplotype diversity was characterized for each gene in this diverse panel. United States rice germplasm was not well represented in the IRRI 3K RGP database. Therefore, a minimum SNP set was identified for each gene that could differentiate all the characterized haplotypes. These haplotypes in the 3K RGP panel were screened across 323 elite U.S. genotypes using the minimum SNP set. The screening of haplotypes and phenotype association confirmed the role of GS3 under qGL3.1. However, screening of the GL7 haplotypes in the U.S. germplasm panel showed that GL7 did not play a role in qGL7.1, and in addition, GL7.1 did not segregate in the Trenasse/Jupiter RIL population. This concluded that qGL7.1 is a novel QTL discovered on chr7 for grain shape in the Trenasse/Jupiter RIL population. A high-throughput KASP-based SNP marker for each locus (GS3 and qGL7.1) was identified and validated in elite U.S. rice germplasm to be used in an applied rice breeding program.

Functional analysis of PHYB polymorphisms in Arabidopsis thaliana collected in Patagonia.

Arabidopsis thaliana shows a wide range of natural genetic variation in light responses. Shade avoidance syndrome is a strategy of major adaptive significance that includes seed germination, elongation of vegetative structures, leaf hyponasty, and acceleration of flowering. Previously, we found that the southernmost Arabidopsis accession, collected in the south of Patagonia (Pat), is hyposensitive to light and displays a reduced response to shade light. This work aimed to explore the genetic basis of the shade avoidance response (SAR) for hypocotyl growth by QTL mapping in a recently developed 162 RIL population between Col-0 and Pat. We mapped four QTL for seedling hypocotyl growth: WL1 and WL2 QTL in white light, SHADE1 QTL in shade light, and SAR1 QTL for the SAR. PHYB is the strongest candidate gene for SAR1 QTL. Here we studied the function of two polymorphic indels in the promoter region, a GGGR deletion, and three non-synonymous polymorphisms on the PHYB coding region compared with the Col-0 reference genome. To decipher the contribution and relevance of each PHYB-Pat polymorphism, we constructed transgenic lines with single or double polymorphisms by using Col-0 as a reference genome. We found that single polymorphisms in the coding region of PHYB have discrete functions in seed germination, seedling development, and shade avoidance response. These results suggest distinct functions for each PHYB polymorphism to the adjustment of plant development to variable light conditions.

Identification of Quantitative Trait Loci Related to Salt Tolerance of Indica Rice RIL Population in Different Growth Stages

Identification of Quantitative Trait Loci Related to Salt Tolerance of Indica Rice RIL Population in Different Growth Stages

Identification of Novel Genes Associated with Partial Resistance to Aphanomyces Root Rot in Field Pea by BSR-Seq Analysis.

Aphanomyces root rot, caused by Aphanomyces euteiches, causes severe yield loss in field pea (Pisum sativum). The identification of a pea germplasm resistant to this disease is an important breeding objective. Polygenetic resistance has been reported in the field pea cultivar ‘00-2067’. To facilitate marker-assisted selection (MAS), bulked segregant RNA-seq (BSR-seq) analysis was conducted using an F8 RIL population derived from the cross of ‘Carman’ × ‘00-2067’. Root rot development was assessed under controlled conditions in replicated experiments. Resistant (R) and susceptible (S) bulks were constructed based on the root rot severity in a greenhouse study. The BSR-seq analysis of the R bulks generated 44,595,510~51,658,688 reads, of which the aligned sequences were linked to 44,757 genes in a reference genome. In total, 2356 differentially expressed genes were identified, of which 44 were used for gene annotation, including defense-related pathways (jasmonate, ethylene and salicylate) and the GO biological process. A total of 344.1 K SNPs were identified between the R and S bulks, of which 395 variants were located in 31 candidate genes. The identification of novel genes associated with partial resistance to Aphanomyces root rot in field pea by BSR-seq may facilitate efforts to improve management of this important disease.

SSR Linkage Maps and Identification of QTL Controlling Morpho-Phenological Traits in Two Iranian Wheat RIL Populations.

Wheat is one of the essential grains grown in large areas. Identifying the genetic structure of agronomic and morphological traits of wheat can help to discover the genetic mechanisms of grain yield. In order to map the morpho-phenological traits, an experiment was conducted in the two cropping years of 2020 and 2021 on the university farm of the Faculty of Agriculture, GonbadKavous University. This study used two F8 populations, including 120 lines resulting from Gonbad × Zagros and Gonbad × Kuhdasht. The number of days to physiological maturity, number of days to flowering, number of germinated grains, number of tillers, number of tillers per plant, grain filling periods, plant height, peduncle length, spike length, awn length, spike weight, peduncle diameter, flag leaf length and weight, number of spikelets per spike, number of grains per spike, grain length, grain width, 1000-grain weight, biomass, grain yield, harvest index, straw-weight, and number of fertile spikelets per spike were measured. A total of 21 and 13 QTLs were identified for 11 and 13 traits in 2020 and 2021, respectively. In 2020, qGL-3D and qHI-1A were identified for grain length and harvest index on chromosomes 3D and 1A, explaining over 20% phenotypic variation, respectively. qNT-5B, qNTS-2D, and qSL-1D were identified on chromosomes 5B, 2D, and 1D with the LOD scores of 4.5, 4.13, and 3.89 in 2021, respectively.

Identification and mapping of genetic locus conferring resistance to multiple plant viruses in soybean.

A reliable locus confers broad-spectrum resistance to multiple plant viruses in soybean under field conditions. Soybean mosaic disease (SMD) can be caused by a variety of viruses, most of which have been largely overlooked in breeding programs. Effective mitigation of the adverse of SMD might result from breeding cultivars with broad-spectrum resistance. However, reports on broad-spectrum resistance to multiple virus have been limited. To catalog viral community members behind SMD, virus samples were collected from symptomatic field plots, and pathogenicity of component strains was assessed. Preliminary ELISA and PCR detection revealed that 39.58% and 66.67% of samples contained two or more virus strains, respectively. Only three soybean accessions were completely asymptomatic, while 42% exhibited moderate or severe susceptibility, indicating that co-infection of multiple virus remains a significant threat in current soybean production systems. Further, a RIL population consisting of 150 F7:9 strains derived from two soybean genotypes with contrasting reactions to virus infection was constructed and explored for significant markers and resistance genes. QTL analysis returned a reliable locus, named GmRmv, on chromosome 13. Significance of GmRmv in imparting resistance to SMD was further confirmed in NIL lines and delimited into a 157-kb interval that contains 17 annotated genes. Among these genes, three, Glyma.13G190000, Glyma.13G190300 and Glyma.13G190400, each contained LRR domains, as well as significant variation in coding sequences between resistant and susceptible parents. Hence, these three genes are considered strong candidate genes for explaining GmRmv significance. In summary, this research opens a new avenue for formulating strategies to breed soybean varieties with broad-spectrum resistance to multiple virus associated with SMD.

Alternative Splicing (AS) Dynamics in Dwarf Soybean Derived from Cross of Glycine max and Glycine soja

Short crop height is the preferred breeding trait since there is a positive correlation between lodging resistance and a crop yield increase. Alternative splicing can alter transcriptome diversity and contribute to plant adaptation to environmental stress. We characterized the transcriptomes obtained from dwarf and normal soybean lines derived from a cross of Glycine max var. Peking (G. max) and G. soja var. IT182936 in an F7 RIL population to study the differences between the isoforms. Full-length mRNA derived from leaf tissues was sequenced using the PacBio RSII platform, generating 904,474 circular consensus sequence (CCS) reads. Using the Structural and Quality Annotation of Novel Transcript Isoforms (SQANTI) process, 42,582 and 44,762 high-quality isoforms, and 91 and 179 polished low-quality isoforms were obtained in dwarf and normal cells, respectively. As a result, 832 and 36,772 nonredundant transcripts were generated. Approximately 30% of the identified genes were estimated to produce two or more isoforms. We detected an average of 166,171 splice junctions (SJs), of which 93.8% were canonical SJs. We identified that novel isoforms accounted for 19% of all isoforms, among which 12% fell within coding regions. The dwarf soybean demonstrated a greater number of isoforms in most of the annotated genes, particularly in genes related to growth hormones and defense responses. Our study provides comprehensive isoform and gene information that may accelerate transcriptome research in G. max and provide a basis to further study the impact of these isoforms on plant growth.

Epistatic interaction effect between chromosome 1BL (Yr29) and a novel locus on 2AL facilitating resistance to stripe rust in Chinese wheat Changwu 357-9.

Four stable QTL for adult plant resistance were identified in wheat line Changwu 357-9, including a new QTL on 2AL showing significant interaction with Yr29 to reduce stripe rust severity. Stripe rust (yellow rust) is a serious disease of bread wheat (Triticum aestivum L.) worldwide. Genetic resistance is considered the most economical, effective and environmentally friendly method to control the disease and to minimize the use of fungicides. The current study focused on characterizing the components of stripe rust resistance and understanding the interactions in Changwu 357-9 (CW357-9)/Avocet S RIL population. A genetic linkage map constructed using a new GenoBaits Wheat 16K Panel and the 660K SNP array had 5104 polymorphic SNP markers spanning 3533.11cM. Four stable QTL, consistently identified across five environments, were detected on chromosome arms 1BL, 2AL, 3DS, and 6BS in Changwu357-9. The most effective QTL QYrCW357-1BL was Yr29. The 6BS QTL was identified as Yr78, which has been combined with the 1BL QTL in many wheat cultivars and breeding lines. The novel QTL on 2AL with moderate effect showed a stable and significant epistatic interaction with Yr29. The QTL on 3DL should be same as QYrsn.nwafu-3DL and enriches the overall stripe rust resistance gene pool for breeding. Polymorphisms of flanking AQP markers AX-110020417 (for QYrCW357-1BL), AX-110974948 (for QYrCW357-2AL), AX-109466386 (for QYrCW357-3DL), and AX-109995005 (for QYrCW357-6BS) were evaluated in a diversity panel including 225 wheat cultivars and breeding lines. These results suggested that these high-throughput markers could be used to introduce QYrCW357-1BL, QYrCW357-2AL, QYrCW357-3DL, and QYrCW357-6BS into commercial wheat cultivars. Combinations of these genes with other APR QTL should lead to higher levels of stripe rust resistance along with the beneficial effects of multi-disease resistance gene Yr29 on improving resistance to other diseases.

Development of candidate gene-based markers and map-based cloning of a dominant rust resistance gene in cultivated groundnut (Arachis hypogaea L.)

A dominant rust resistance gene, VG 9514-Rgene was isolated through map-based cloning. Sequence analysis revealed non-synonymous mutations in the TIR, NBS and LRR region of the R-protein. Candidate gene-based markers from these SNPs revealed complete co-segregation of the isolated VG 9514-Rgene with rust resistance in a RIL population and confirmed their map position in between FRS 72 and SSR_GO340445 markers in arahy03 chromosome. Blastp search of VG 9514-Rprotein detected Arahy.T6DCA5 with >80.0% identity that localized at 142,544,745.0.142,549,184 in arahy03 chromosome. Ka/Ks calculation revealed that VG 9514-Rgene had undergone positive selection compared to four homologous genes in the groundnut genome. Homology based structure modelling of this R-protein revealed a typical consensus three-dimensional folding of TIR-NBS-LRR protein. Non-synonymous mutations in susceptible version of R-protein were mapped and found E268Q mutation in hhGRExE motif, Y309F in RNBS-A motif and I579T in MHD motif of NB-ARC domain are probable candidates for loss of function.

Identification of Drought Tolerant Recombinant Inbred Lines and Estimation of Genetic Diversity by Clustering and Multivariate Analysis in Wheat RILs Populations

The current experiments were carried out in alpha lattice design at the Seed Breeding Farm, Department of Plant Breeding and Genetics, Jawaharlal Nehru Krishi Vishwa Vidyalaya, Jabalpur, Madhya Pradesh state, India, during 2019−20 and 2020−21 crop seasons, using 236 RILs population. The aim of current experiment is to assessment of genetic diversity and identification of promising recombinant inbred lines for irrigated and restricted irrigated conditions through PCA and cluster analysis. The drought selection indices viz. Relative drought index, yield stability index, drought resistance index, mean productivity, stress tolerance index and yield index were performed from two year grain yield pooled data. The principal component analysis and cluster analysis were performed through drought selection indices. Drought selection indices viz. relative drought index, yield stability index, yield index were confirm strong positively associated with grain yield under restricted irrigated conditions while mean productivity, stress tolerance index and yield index were strong positively associated with grain yield under irrigated conditions. Moreover, high cluster mean, for grain yield under restricted irrigated condition with associated selection indices was confirmed by 77 inbred lines from cluster I. similarly 128 superior inbred lines were found for irrigated condition. The highest inter cluster distance was observed between cluster I and cluster III therefore the inbred line occupy by cluster I and cluster III were considered as most diverse lines and could be used in farther breeding program to achieve more recombination for drought tolerance.

Comprehensive transcriptomic analysis of two RIL parents with contrasting salt responsiveness identifies polyadenylated and non-polyadenylated flower lncRNAs in chickpea.

SummarySalinity severely affects the yield of chickpea. Understanding the role of lncRNAs can shed light on chickpea salt tolerance mechanisms. However, because lncRNAs are encoded by multiple sites within the genome, their classification to reveal functional versatility at the transcriptional and the post‐transcriptional levels is challenging. To address this, we deep sequenced 24 salt‐challenged flower transcriptomes from two parental genotypes of a RIL population that significantly differ in salt tolerance ability. The transcriptomes for the first time included 12 polyadenylated and 12 non‐polyadenylated RNA libraries to a sequencing depth of ~50 million reads. The ab initio transcriptome assembly comprised ~34 082 transcripts from three biological replicates of salt‐tolerant (JG11) and salt‐sensitive (ICCV2) flowers. A total of 9419 lncRNAs responding to salt stress were identified, 2345 of which were novel lncRNAs specific to chickpea. The expression of poly(A+) lncRNAs and naturally antisense transcribed RNAs suggest their role in post‐transcriptional modification and gene silencing. Notably, 178 differentially expressed lncRNAs were induced in the tolerant genotype but repressed in the sensitive genotype. Co‐expression network analysis revealed that the induced lncRNAs interacted with the FLOWERING LOCUS (FLC), chromatin remodelling and DNA methylation genes, thus inducing flowering during salt stress. Furthermore, 26 lncRNAs showed homology with reported lncRNAs such as COOLAIR, IPS1 and AT4, thus confirming the role of chickpea lncRNAs in controlling flowering time as a crucial salt tolerance mechanism in tolerant chickpea genotype. These robust set of differentially expressed lncRNAs provide a deeper insight into the regulatory mechanisms controlled by lncRNAs under salt stress.

Genome-wide identification of quantitative trait loci for morpho-agronomic and yield-related traits in foxtail millet (Setaria italica) across multi-environments

Foxtail millet (Setaria italica) is an ideal model of genetic system for functional genomics of the Panicoideae crop. Identification of QTL responsible for morpho-agronomic and yield-related traits facilitates dissection of genetic control and breeding in cereal crops. Here, based on a Yugu1 × Longgu7 RIL population and genome-wide resequencing data, an updated linkage map harboring 2297 bin and 74 SSR markers was constructed, spanning 1315.1 cM with an average distance of 0.56 cM between adjacent markers. A total of 221 QTL for 17 morpho-agronomic and yield-related traits explaining 5.5 ~ 36% of phenotypic variation were identified across multi-environments. Of these, 109 QTL were detected in two to nine environments, including the most stable qLMS6.1 harboring a promising candidate gene Seita.6G250500, of which 70 were repeatedly identified in different trials in the same geographic location, suggesting that foxtail millet has more identical genetic modules under the similar ecological environment. One hundred-thirty QTL with overlapping intervals formed 22 QTL clusters. Furthermore, six superior recombinant inbred lines, RIL35, RIL48, RIL77, RIL80, RIL115 and RIL125 with transgressive inheritance and enrichment of favorable alleles in plant height, tiller, panicle morphology and yield related-traits were screened by hierarchical cluster. These identified QTL, QTL clusters and superior lines lay ground for further gene-trait association studies and breeding practice in foxtail millet.

Residual effect of defeated stripe rust resistance genes/QTLs in bread wheat against prevalent pathotypes of Puccinia striiformis f. sp. tritici.

The periodic breakdowns of stripe rust resistance due to emergence of new virulent and more aggressive pathotypes of Puccinia striiformis f. sp. tritici have resulted in severe epidemics in India. This necessitates the search for new and more durable resistance sources against stripe rust. The three bread wheat cultivars PBW 343 (carries Yr9 and Yr27), PBW 621 (carries Yr17) and HD 2967 (gene not known) were highly popular among the farmers after their release in 2011. But presently all three cultivars are highly susceptible to stripe rust at seedling as well as at adult plant stages as their resistance has been broken down due to emergence of new pathotypes of the pathogen (110S119, 238S119). In previous study, the crosses of PBW 621 with PBW 343 and HD 2967 and evaluation of further generations (up to F4) against pathotype 78S84 resulted in resistant segregants. In the present study, the F5 and F6 RIL populations have been evaluated against new pathotypes of Pst. The RILs categorized based on the disease severity on the P (Penultimate leaf) and F (flag) leaf into three categories i.e., high, moderate and low level of APR (adult plant resistance) having 1-200, 201-400 and >400 values of AUDPC, respectively, upon infection with stripe rust. The various APR components (latent period, lesion growth rate, spore production and uredial density) were studied on each category, i.e., resistant, moderately resistant and susceptible. The values of APR parameters decreased as the level of resistance increased. Based on molecular analysis, the lines (representing different categories of cross PBW 621 X PBW 343) containing the genes Yr9 and Yr17 due to their interactive effect provide resistance. Based on BSA using 35k SNPs and KASP markers association with phenotypic data of the RIL population (PBW 621 X HD 2967) showed the presence of two QTLs (Q.Pst.pau-6B, Q.Pst.pau-5B) responsible for the residual resistance and two SNPs AX-94891670 and AX-94454107 were found to be associated with the trait of interest on chromosome 6B and 5B respectively. The present study concludes that in the population of both the crosses (PBW 621 X PBW 343 and PBW 621 X HD 2967) major defeated gene contributed towards residual resistance by interacting with minor gene/QTLs.

Unravelling the Genetic Architecture of Rust Resistance in the Common Bean (Phaseolus vulgaris L.) by Combining QTL-Seq and GWAS Analysis.

The common bean (Phaseolus vulgaris L.) is the most important legume crop directly used for human consumption worldwide. Bean rust, caused by Uromyces appendiculatus, is a devastating disease and usually causes severe loss of seed yield and pod quality. Deployment of resistant cultivars is the best strategy to combat this disease. However, despite being the largest snap bean-producing country, the genetic basis research of rust resistance has largely lagged in China. In this study, an RIL population and a diversity panel were evaluated for rust resistance against a purified rust isolate Cua-LS using a detached leaf assay. Deploying a QTL-Seq analysis in the RIL population, a 1.81 Mb interval on chromosome 4, a 2.73 Mb interval on chromosome 5 and a 1.26 Mb interval on chromosome 6 were identified as major QTLs for rust resistance, designated as Qur-1, Qur-2 and Qur-3, respectively. Through a GWAS diversity panel, 64 significant SNPs associated with rust resistance were detected, distributed in all 11 chromosomes and explaining 19–49% of the phenotypic variation. Synteny analysis showed that Qur-2 was validated in GWAS, but the rust QTL/SNPs detected in our study were different from the known genes, except Ur-11. A total of 114 candidate genes, including the typical NBS-LRR genes, protein kinase superfamily proteins and ABC transporter family proteins, were identified and proposed as the likely candidates. The identified 17 resistant accessions will enrich the resistant germplasm resources, and the detected QTLs/SNPs will facilitate the molecular breeding of rust resistance in the common bean.

Maize plant architecture trait QTL mapping and candidate gene identification based on multiple environments and double populations

BackgroundThe plant architecture traits of maize determine the yield. Plant height, ear position, leaf angle above the primary ear and internode length above the primary ear together determine the canopy structure and photosynthetic efficiency of maize and at the same time affect lodging and disease resistance. A flat and tall plant architecture confers an obvious advantage in the yield of a single plant but is not conducive to dense planting and results in high rates of lodging; thus, it has been gradually eliminated in production. Although using plants that are too compact, short and density tolerant can increase the yield per unit area to a certain extent, the photosynthetic efficiency of such plants is low, ultimately limiting yield increases. Genetic mapping is an effective method for the improvement of plant architecture to identify candidate genes for regulating plant architecture traits.ResultsTo find the best balance between the yield per plant and the yield per unit area of maize, in this study, the F2:3 pedigree population and a RIL population with the same male parent were used to identify QTL for plant height (PH), ear height (EH), leaf angle and internode length above the primary ear (LAE and ILE) in Changchun and Gongzhuling for 5 consecutive years (2016–2020). A total of 11, 13, 23 and 13 QTL were identified for PH, EH, LAE, and ILE, respectively. A pleiotropic consistent QTL for PH overlapped with that for EH on chromosome 3, with a phenotypic variation explanation rate from 6.809% to 21.96%. In addition, there were major consistent QTL for LAE and ILE, and the maximum phenotypic contribution rates were 24.226% and 30.748%, respectively. Three candidate genes were mined from the three consistent QTL regions and were involved in the gibberellin-activated signal pathway, brassinolide signal transduction pathway and auxin-activated signal pathway, respectively. Analysis of the expression levels of the three genes showed that they were actively expressed during the jointing stage of vigorous maize growth.ConclusionsIn this study, three consistent major QTL related to plant type traits were identified and three candidate genes were screened. These results lay a foundation for the cloning of related functional genes and marker-assisted breeding of related functional genes.