Deciphering the Genetic Basis of Root and Biomass Traits in Rapeseed (Brassica napus L.) through the Integration of GWAS and RNA-Seq under Nitrogen Stress

An excellent root system is responsible for crops with high nitrogen-use efficiency (NUE). The current study evaluated the natural variations in 13 root- and biomass-related traits under a low nitrogen (LN) treatment in a rapeseed association panel. The studied traits exhibited significant phenotypic differences with heritabilities ranging from 0.53 to 0.66, and most of the traits showed significant correlations with each other. The genome-wide association study (GWAS) found 51 significant and 30 suggestive trait–SNP associations that integrated into 14 valid quantitative trait loci (QTL) clusters and explained 5.7–21.2% phenotypic variance. In addition, RNA sequencing was performed at two time points to examine the differential expression of genes (DEGs) between high and low NUE lines. In total, 245, 540, and 399 DEGs were identified as LN stress-specific, high nitrogen (HN) condition-specific, and HNLN common DEGs, respectively. An integrated analysis of GWAS, weighted gene co-expression network, and DEGs revealed 16 genes involved in rapeseed root development under LN stress. Previous studies have reported that the homologs of seven out of sixteen potential genes control root growth and NUE. These findings revealed the genetic basis underlying nitrogen stress and provided worthwhile SNPs/genes information for the genetic improvement of NUE in rapeseed.     

The Combination of Conventional QTL Analysis,Bulked-Segregant Analysis,and RNA-Sequencing Provide New Genetic Insights into Maize Mesocotyl Elongation under Multiple Deep-Seeding Environments

Mesocotyl length (MES) is an important trait that affects the emergence of maize seedlings after deep-seeding and is closely associated with abiotic stress. The elucidation of constitutive-QTLs (cQTLs) and candidate genes for MES and tightly molecular markers are thus of great importance in marker-assisted selection (MAS) breeding. Therefore, the objective of this study was to perform detailed genetic analysis of maize MES across 346 F_2:3 families, 30/30 extreme bulks of an F₂ population, and two parents by conventional QTL analysis, bulked-segregation analysis (BSA), and RNA-sequencing when maize was sown at the depths of 3, 15, and 20 cm, respectively. QTL analysis identified four major QTLs in Bin 1.09, Bin 3.04, Bin 4.06–4.07, and Bin 6.01 under two or more environments, which explained 2.89–13.97% of the phenotypic variance within a single environment. BSA results revealed the presence of seven significantly linked SNP/InDel regions on chromosomes 1 and 4, and six SNP/InDel regions and the major QTL of qMES4-1 overlapped and formed a cQTL, cQMES4, within the 160.98–176.22 Mb region. In total, 18,001 differentially expressed genes (DEGs) were identified across two parents by RNA-sequencing, and 24 of these genes were conserved core DEGs. Finally, we validated 15 candidate genes in cQMES4 to involve in cell wall structure, lignin biosyntheis, phytohormones (auxin, abscisic acid, brassinosteroid) signal transduction, circadian clock, and plant organ formation and development. Our findings provide a basis for MAS breeding and enhance our understanding of the deep-seeding tolerance of maize.     

Combined QTL Mapping across Multiple Environments and Co-Expression Network Analysis Identified Key Genes for Embryogenic Callus Induction from Immature Maize Embryos

The ability of immature embryos to induce embryogenic callus (EC) is crucial for genetic transformation in maize, which is highly genotype-dependent. To dissect the genetic basis of maize EC induction, we conducted QTL mapping for four EC induction-related traits, the rate of embryogenic callus induction (REC), rate of shoot formation (RSF), length of shoot (LS), and diameter of callus (DC) under three environments by using an IBM Syn10 DH population derived from a cross of B73 and Mo17. These EC induction traits showed high broad-sense heritability (>80%), and significantly negative correlations were observed between REC and each of the other traits across multiple environments. A total of 41 QTLs for EC induction were identified, among which 13, 12, 10, and 6 QTLs were responsible for DC, RSF, LS, and REC, respectively. Among them, three major QTLs accounted for >10% of the phenotypic variation, including qLS1-1 (11.54%), qLS1-3 (10.68%), and qREC4-1 (11.45%). Based on the expression data of the 215 candidate genes located in these QTL intervals, we performed a weighted gene co-expression network analysis (WGCNA). A combined use of KEGG pathway enrichment and eigengene-based connectivity (KME) values identified the EC induction-associated module and four hub genes (Zm00001d028477, Zm00001d047896, Zm00001d034388, and Zm00001d022542). Gene-based association analyses validated that the variations in Zm00001d028477 and Zm00001d034388, which were involved in tryptophan biosynthesis and metabolism, respectively, significantly affected EC induction ability among different inbred lines. Our study brings novel insights into the genetic and molecular mechanisms of EC induction and helps to promote marker-assisted selection of high-REC varieties in maize.     

Whole-Genome Sequencing and RNA-Seq Reveal Differences in Genetic Mechanism for Flowering Response between Weedy Rice and Cultivated Rice

Flowering is a key agronomic trait that influences adaptation and productivity. Previous studies have indicated the genetic complexity associated with the flowering response in a photoinsensitive weedy rice accession PSRR-1 despite the presence of a photosensitive allele of a key flowering gene Hd1. In this study, we used whole-genome and RNA sequencing data from both cultivated and weedy rice to add further insights. The de novo assembly of unaligned sequences predicted 225 genes, in which 45 were specific to PSRR-1, including two genes associated with flowering. Comparison of the variants in PSRR-1 with the 3K rice genome (RG) dataset identified unique variants within the heading date QTLs. Analyses of the RNA-Seq result under both short-day (SD) and long-day (LD) conditions revealed that many differentially expressed genes (DEGs) colocalized with the flowering QTLs, and some DEGs such as Hd1, OsMADS56, Hd3a, and RFT1 had unique variants in PSRR-1. Ehd1, Hd1, OsMADS15, and OsMADS56 showed different alternate splicing (AS) events between genotypes and day length conditions. OsMADS56 was expressed in PSRR-1 but not in Cypress under both LD and SD conditions. Based on variations in both sequence and expression, the unique flowering response in PSRR-1 may be due to the high-impact variants of flowering genes, and OsMADS56 is proposed as a key regulator for its day-neutral flowering response.     

Quantitative Trait Locus Mapping of Salt Tolerance in Wild Rice Oryza longistaminata

Salt stress is one of the most severe adverse environments in rice production; increasing salinization is seriously endangering rice production around the world. In this study, a rice backcross inbred line (BIL) population derived from the cross of 9311 and wild rice Oryza longistaminata was employed to identify the favorable genetic loci of O. longistaminata for salt tolerance. A total of 27 quantitative trait loci (QTLs) related to salt tolerance were identified in 140 rice BILs, and 17 QTLs formed seven QTL clusters on different chromosomes, of which 18 QTLs were derived from O. longistaminata, and a QTL for salt injury score (SIS), water content of seedlings (WCS) under salt treatment, and relative water content of seedlings (RWCS) was repeatedly detected and colocalized at the same site on chromosome 2, and a cytochrome P450 86B1 (MH02t0466900) was suggested as the potential candidate gene responsible for the salt tolerance based on sequence and expression analysis. These findings laid the foundation for further improving rice salt tolerance through molecular breeding in the future.     

Linkage Mapping Reveals QTL for Flowering Time-Related Traits under Multiple Abiotic Stress Conditions in Maize

Variation in flowering plays a major role in maize photoperiod adaptation during long-term domestication. It is of high value to investigate the genetic basis of maize flowering under a wide range of environmental conditions in order to overcome photoperiod sensitivity or enhance stress tolerance. A recombinant inbred line (RIL) population derived from a cross between Huangzaosi and Mo17, composed of 121 lines and genotyped by 8329 specifically developed markers, was field evaluated in two consecutive years under two planting densities (67,500 and 120,000 plants ha⁻¹) and two water treatments (normal irrigation and drought stress at the flowering stage). The days to silking (DTS), days to anthesis (DTA), and anthesis to silking interval (ASI) were all evaluated. Within the RIL population, DTS and DTA expanded as planting density and water deficit increased. For DTA, DTS, ASI, and ASI-delay, a total of 22, 17, 21, and 11 QTLs were identified, respectively. More than two significant QTLs were identified in each of the nine chromosomal intervals. Under diverse conditions and locations, six QTLs (quantitative trait locus) for DTS and DTA were discovered in Chr. 8: 118.13–125.31 Mb. Three chromosome regions, Chr. 3: 196.14–199.89 Mb, Chr. 8: 169.02–172.46 Mb, and Chr. 9: 128.12–137.26 Mb, all had QTLs for ASI-delay under normal and stress conditions, suggesting their possible roles in stress tolerance enhancement. These QTL hotspots will promote early-maturing or multiple abiotic stress-tolerant maize breeding, as well as shed light on the development of maize varieties with a broad range of adaptations.     

Dissection of the Genetic Basis of Yield Traits in Line per se and Testcross Populations and Identification of Candidate Genes for Hybrid Performance in Maize

Dissecting the genetic basis of yield traits in hybrid populations and identifying the candidate genes are important for molecular crop breeding. In this study, a BC1F3:4 population, the line per se (LPS) population, was constructed by using elite inbred lines Zheng58 and PH4CV as the parental lines. The population was genotyped with 55,000 SNPs and testcrossed to Chang7-2 and PH6WC (two testers) to construct two testcross (TC) populations. The three populations were evaluated for hundred kernel weight (HKW) and yield per plant (YPP) in multiple environments. Marker–trait association analysis (MTA) identified 24 to 151 significant SNPs in the three populations. Comparison of the significant SNPs identified common and specific quantitative trait locus/loci (QTL) in the LPS and TC populations. Genetic feature analysis of these significant SNPs proved that these SNPs were associated with the tested traits and could be used to predict trait performance of both LPS and TC populations. RNA-seq analysis was performed using maize hybrid varieties and their parental lines, and differentially expressed genes (DEGs) between hybrid varieties and parental lines were identified. Comparison of the chromosome positions of DEGs with those of significant SNPs detected in the TC population identified potential candidate genes that might be related to hybrid performance. Combining RNA-seq analysis and MTA results identified candidate genes for hybrid performance, providing information that could be useful for maize hybrid breeding.     

Discovery of Genomic Regions and Candidate Genes Controlling Root Development Using a Recombinant Inbred Line Population in Rapeseed (Brassica napus L.)

Marker-assisted selection enables breeders to quickly select excellent root architectural variations, which play an essential role in plant productivity. Here, ten root-related and shoot biomass traits of a new F₆ recombinant inbred line (RIL) population were investigated under hydroponics and resulted in high heritabilities from 0.61 to 0.83. A high-density linkage map of the RIL population was constructed using a Brassica napus 50k Illumina single nucleotide polymorphism (SNP) array. A total of 86 quantitative trait loci (QTLs) explaining 4.16–14.1% of the phenotypic variances were detected and integrated into eight stable QTL clusters, which were repeatedly detected in different experiments. The codominant markers were developed to be tightly linked with three major QTL clusters, qcA09-2, qcC08-2, and qcC08-3, which controlled both root-related and shoot biomass traits and had phenotypic contributions greater than 10%. Among these, qcA09-2, renamed RT.A09, was further fine-mapped to a 129-kb interval with 19 annotated genes in the B. napus reference genome. By integrating the results of real-time PCR and comparative sequencing, five genes with expression differences and/or amino acid differences were identified as important candidate genes for RT.A09. Our findings laid the foundation for revealing the molecular mechanism of root development and developed valuable markers for root genetic improvement in rapeseed.     

Identification of Genomic Regions Associated with Vine Growth and Plant Height of Soybean

Vining growth (VG) and high plant height (PH) are the physiological traits of wild soybean that preclude their utilization for domesticated soybean breeding and improvement. To identify VG- and PH-related quantitative trait loci (QTLs) in different genetic resources, two populations of recombinant inbred lines (RILs) were developed by crossing a cultivated soybean, Zhonghuang39 (ZH39), with two wild soybean accessions, NY27-38 and NY36-87. Each line from the two crosses was evaluated for VG and PH. Three QTLs for VG and three for PH, detected in the ZH39 × NY27-38 population of the RILs, co-located on chromosomes 2, 17 and 19. The VG- and PH-related QTL in the ZH39 × NY36-87 population co-located on chromosome 19. A common QTL shared by the two populations was located on chromosome 19, suggesting that this major QTL was consistently selected for in different genetic backgrounds. The results suggest that different loci are involved in the domestication or adaptations of soybean of various genetic backgrounds. The molecular markers presented here would benefit the fine mapping and cloning of candidate genes underlying the VG and PH co-localized regions and thus facilitate the utilization of wild resources in breeding by avoiding undesirable traits.     

Characterization of the Resistance to Powdery Mildew and Leaf Rust Carried by the Bread Wheat Cultivar Victo

Leaf rust and powdery mildew are two important foliar diseases in wheat. A recombinant inbred line (RIL) population, obtained by crossing two bread wheat cultivars (‘Victo’ and ‘Spada’), was evaluated for resistance to the two pathogens at seedling stage. Upon developing a genetic map of 8726 SNP loci, linkage analysis identified three resistance Quantitative Trait Loci (QTLs), with ‘Victo’ contributing the resistant alleles to all loci. One major QTL (QPm.gb-7A) was detected in response to Blumeria graminis on chromosome 7A, which explained 90% of phenotypic variation (PV). The co-positional relationship with known powdery mildew (Pm) resistance loci suggested that a new source of resistance was identified in T. aestivum. Two QTLs were detected in response to Puccinia triticina: a major gene on chromosome 5D (QLr.gb-5D), explaining a total PV of about 59%, and a minor QTL on chromosome 2B (QLr.gb-2B). A positional relationship was observed between the QLr.gb-5D with the known Lr1 gene, but polymorphisms were found between the cloned Lr1 and the corresponding ‘Victo’ allele, suggesting that QLr.gb-5D could represent a new functional Lr1 allele. Lastly, upon anchoring the QTL on the T. aestivum reference genome, candidate genes were hypothesized on the basis of gene annotation and in silico gene expression analysis.     

Combined BSA-Seq Based Mapping and RNA-Seq Profiling Reveal Candidate Genes Associated with Plant Architecture in Brassica napus

Plant architecture involves important agronomic traits affecting crop yield, resistance to lodging, and fitness for mechanical harvesting in Brassica napus. Breeding high-yield varieties with plant architecture suitable for mechanical harvesting is the main goal of rapeseed breeders. Here, we report an accession of B. napus (4942C-5), which has a dwarf and compact plant architecture in contrast to cultivated varieties. A BC₈ population was constructed by crossing a normal plant architecture line, 8008, with the recurrent parent 4942C-5. To investigate the molecular mechanisms underlying plant architecture, we performed phytohormone profiling, bulk segregant analysis sequencing (BSA-Seq), and RNA sequencing (RNA-Seq) in BC₈ plants with contrasting plant architecture. Genetic analysis indicated the plant architecture traits of 4942C-5 were recessive traits controlled by multiple genes. The content of auxin (IAA), gibberellin (GA), and abscisic acid (ABA) differed significantly between plants with contrasting plant architecture in the BC₈ population. Based on BSA-Seq analysis, we identified five candidate intervals on chromosome A01, namely those of 0 to 6.33 Mb, 6.45 to 6.48 Mb, 6.51 to 6.53 Mb, 6.77 to 6.79 Mb, and 7 to 7.01 Mb regions. The RNA-Seq analysis revealed a total of 4378 differentially expressed genes (DEGs), of which 2801 were up-regulated and 1577 were down-regulated. There, further analysis showed that genes involved in plant hormone biosynthesis and signal transduction, cell structure, and the phenylpropanoid pathway might play a pivotal role in the morphogenesis of plant architecture. Association analysis of BSA-Seq and RNA-Seq suggested that seven DEGs involved in plant hormone signal transduction and a WUSCHEL-related homeobox (WOX) gene (BnaA01g01910D) might be candidate genes responsible for the dwarf and compact phenotype in 4942C-5. These findings provide a foundation for elucidating the mechanisms underlying rapeseed plant architecture and should contribute to breed new varieties suitable for mechanization.     

GsRSS3L,a Candidate Gene Underlying Soybean Resistance to Seedcoat Mottling Derived from Wild Soybean (Glycine soja Sieb. and Zucc)

Soybeans are a major crop that produce the best vegetable oil and protein for use in food and beverage products worldwide. However, one of the most well-known viral infections affecting soybeans is the Soybean Mosaic Virus (SMV), a member of the Potyviridae family. A crucial method for preventing SMV damage is the breeding of resistant soybean cultivars. Adult resistance and resistance of seedcoat mottling are two types of resistance to SMV. Most studies have focused on adult-plant resistance but not on the resistance to seedcoat mottling. In this study, chromosome segment-substituted lines derived from a cross between Suinong14 (cultivated soybean) and ZYD00006 (wild soybean) were used to identify the chromosome region and candidate genes underlying soybean resistance to seed coat mottling. Herein, two quantitative trait loci (QTLs) were found on chromosome 17, and eighteen genes were found in the QTL region. RNA-seq was used to evaluate the differentially expressed genes (DEGs) among the eighteen genes located in the QTLs. According to the obtained data, variations were observed in the expression of five genes following SMV infection. Furthermore, Nicotiana benthamiana was subjected to an Agrobacterium-mediated transient expression assay to investigate the role of the five candidate genes in SMV resistance. It has also been revealed that Glyma.17g238900 encoding a RICE SALT SENSITIVE 3-like protein (RSS3L) can inhibit the multiplication of SMV in N. benthamiana. Moreover, two nonsynonymous single-nucleotide polymorphisms (SNPs) were found in the coding sequence of Glyma.17g238900 derived from the wild soybean ZYD00006 (GsRSS3L), and the two amino acid mutants may be associated with SMV resistance. Hence, it has been suggested that GsRSS3L confers seedcoat mottling resistance, shedding light on the mechanism of soybean resistance to SMV.     

Mapping QTLs Controlling Beneficial Fatty Acids Based on the Embryo and Maternal Plant Genomes in Brassica napus L.

This study was mainly carried out to detect the quantitative trait loci (QTLs) located simultaneously in embryo and maternal plant nuclear genomes for beneficial fatty acid contents of rapeseed under different environments by using a mapping model and mapping populations from a cross of ‘Tapidor’ × ‘Ningyou7.’ Eight, three, six, seven and three QTLs respectively for palmitic, oleic, linoleic, linolenic and eicosenoic acid contents were identified and subsequently mapped on chromosomes A1, A3, A4, A6, A8, A9, C1, C3, C5, C7, C8 and C9. Among them, 12 QTLs were major ones that could respectively explain more than 10 % of phenotypic variation. Two specific genomic regions on A8 and C3 were associated with several QTLs relating to multiple fatty acid contents. Some of the markers such as HBr015 and JICB0633 showed tight linkage with the QTLs, which could be used for improving the fatty acid component(s) of rapeseed.     

Identification of a Major QTL and Validation of Related Genes for Tiller Angle in Rice Based on QTL Analysis

An ideal plant architecture is an important condition to achieve high crop yields. The tiller angle is an important and complex polygenic trait of rice (Oryza sativa L.) plant architecture. Therefore, the discovery and identification of tiller angle-related genes can aid in the improvement of crop architecture and yield. In the present study, 222 SSR markers were used to establish a high-density genetic map of rice doubled haploid population, and a total of 8 quantitative trait loci (QTLs) were detected based on the phenotypic data of the tiller angle and tiller crown width over 2 years. Among them, four QTLs (qTA9, qCW9, qTA9-1, and qCW9-1) were overlapped at marker interval RM6235–RM24288 on chromosome 9 with a large effect value regarded as a stable major QTL. The selected promising related genes were further identified by relative gene expression analysis, which gives us a basis for the future cloning of these genes. Finally, OsSAURq9, which belongs to the SMALL AUXIN UP RNA (SAUR), an auxin-responsive protein family, was selected as a target gene. Overall, this work will help broaden our knowledge of the genetic control of tiller angle and tiller crown width, and this study provides both a good theoretical basis and a new genetic resource for the breeding of ideal-type rice.     

Identification and Allele Combination Analysis of Rice Grain Shape-Related Genes by Genome-Wide Association Study

Grain shape is an important agronomic character of rice, which affects the appearance, processing, and the edible quality. Screening and identifying more new genes associated with grain shape is beneficial to further understanding the genetic basis of grain shape and provides more gene resources for genetic breeding. This study has a natural population containing 623 indica rice cultivars. Genome-wide association studies/GWAS of several traits related to grain shape (grain length/GL, grain width/GW, grain length to width ratio/GLWR, grain circumferences/GC, and grain size/grain area/GS) were conducted by combining phenotypic data from four environments and the second-generation resequencing data, which have identified 39 important Quantitative trait locus/QTLs. We analyzed the 39 QTLs using three methods: gene-based association analysis, haplotype analysis, and functional annotation and identified three cloned genes (GS3, GW5, OsDER1) and seven new candidate genes in the candidate interval. At the same time, to effectively utilize the genes in the grain shape-related gene bank, we have also analyzed the allelic combinations of the three cloned genes. Finally, the extreme allele combination corresponding to each trait was found through statistical analysis. This study’s novel candidate genes and allele combinations will provide a valuable reference for future breeding work.