QTL Mapping and Diurnal Transcriptome Analysis Identify Candidate Genes Regulating Brassica napus Flowering Time

Timely flowering is important for seed formation and maximization of rapeseed (Brassica napus) yield. Here, we performed flowering-time quantitative trait loci (QTL) mapping using a double haploid (DH) population grown in three environments to study the genetic architecture. Brassica 60 K Illumina Infinium™ single nucleotide polymorphism (SNP) array and simple sequence repeat (SSR) markers were used for genotyping of the DH population, and a high-density genetic linkage map was constructed. QTL analysis of flowering time from the three environments revealed five consensus QTLs, including two major QTLs. A major QTL located on chromosome A03 was detected specifically in the semi-winter rapeseed growing region, and the one on chromosome C08 was detected in all environments. Ribonucleic acid sequencing (RNA-seq) was performed on the parents’ leaves at seven time-points in a day to determine differentially expressed genes (DEGs). The biological processes and pathways with significant enrichment of DEGs were obtained. The DEGs in the QTL intervals were analyzed, and four flowering time-related candidate genes were found. These results lay a foundation for the genetic regulation of rapeseed flowering time and create a rapeseed gene expression library for seven time-points in a day.     

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.     

Fine Mapping and Functional Analysis of Major Regulatory Genes of Soluble Solids Content in Wax Gourd (Benincasa hispida)

Soluble solids content (SSC) is an important quality trait of wax gourd, but reports about its regulatory genes are scarce. In this study, the SSC regulatory gene BhSSC2.1 in wax gourd was mined via quantitative trait locus (QTL) mapping based on high-density genetic mapping containing 12 linkage groups (LG) and bulked segregant analysis (BSA)-seq. QTL mapping and BSA-seq revealed for the first time that the SSC QTL (107.658–108.176 cM) of wax gourd was on Chr2 (LG2). The interpretable phenotypic variation rate and maximum LOD were 16.033% and 6.454, respectively. The QTL interval contained 13 genes. Real-time fluorescence quantitative expression analysis, functional annotation, and sequence analysis suggested that Bch02G016960, named BhSSC2.1, was a candidate regulatory gene of the SSC in wax gourd. Functional annotation of this gene showed that it codes for a NADP-dependent malic enzyme. According to BhSSC2.1 sequence variation, an InDel marker was developed for molecular marker-assisted breeding of wax gourd. This study will lay the foundation for future studies regarding breeding and understanding genetic mechanisms of wax gourd.     

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.     

Analysis of QTLs and Candidate Genes for Tassel Symptoms in Maize Infected with Sporisorium reilianum

Heat smut is a fungal soil-borne disease caused by Sporisorium reilianum, and affects the development of male and female tassels. Our previous research found that the tassel symptoms in maize infected with Sporisorium reilianum significantly differed in inbred lines with Sipingtou blood, and exhibited stable heredity over time at multiple locations. In this study, cytological analysis demonstrated that the cellular organization structures of three typical inbred lines (Huangzao4, Jing7, and Chang7-2) showed significant discrepancies at the VT stage. QTLs that control the different symptoms of maize tassels infected with Sporisorium reilianum were located in two F₂ populations, which were constructed using three typical inbred lines. The BSA (bulked segregation analysis) method was used to construct mixed gene pools based on typical tassel symptoms. The QTLs of different symptoms of maize tassels infected with Sporisorium reilianum were detected with 869 SSR markers covering the whole maize genome. The mixed gene pools were screened with polymorphic markers between the parents. Additional SSR markers were added near the above marker to detect genotypes in partially single plants in F₂ populations. The QTL controlling tassel symptoms in the Huangzao4 and Jing7 lines was located on the bin 1.06 region, between the markers of umc1590 and bnlg1598, and explained 21.12% of the phenotypic variation with an additive effect of 0.6524. The QTL controlling the tassel symptoms of the Jing7 and Chang7-2 lines was located on the bin 2.07 region, between the markers of umc1042 and bnlg1335, and explained 11.26% phenotypic variation with an additive effect of 0.4355. Two candidate genes (ZmABP2 and Zm00001D006403) were identified by a conjoint analysis of label-free quantification proteome sequencings.     

A Quantitative Genetic Study of Sclerotinia Head Rot Resistance Introgressed from the Wild Perennial Helianthus maximiliani into Cultivated Sunflower (Helianthus annuus L.)

Sclerotinia head rot (HR), caused by Sclerotinia sclerotiorum, is an economically important disease of sunflower with known detrimental effects on yield and quality in humid climates worldwide. The objective of this study was to gain insight into the genetic architecture of HR resistance from a sunflower line HR21 harboring HR resistance introgressed from the wild perennial Helianthus maximiliani. An F₂ population derived from the cross of HA 234 (susceptible-line)/HR21 (resistant-line) was evaluated for HR resistance at two locations during 2019–2020. Highly significant genetic variations (p < 0.001) were observed for HR disease incidence (DI) and disease severity (DS) in both individual and combined analyses. Broad sense heritability (H²) estimates across environments for DI and DS were 0.51 and 0.62, respectively. A high-density genetic map of 1420.287 cM was constructed with 6315 SNP/InDel markers developed using genotype-by-sequencing technology. A total of 16 genomic regions on eight sunflower chromosomes, 1, 2, 10, 12, 13, 14, 16 and 17 were associated with HR resistance, each explaining between 3.97 to 16.67% of the phenotypic variance for HR resistance. Eleven of these QTL had resistance alleles from the HR21 parent. Molecular markers flanking the QTL will facilitate marker-assisted selection breeding for HR resistance in sunflower.     

Identifying Soybean Pod Borer (Leguminivora glycinivorella) Resistance QTLs and the Mechanism of Induced Defense Using Linkage Mapping and RNA-Seq Analysis

The soybean pod borer (Leguminivora glycinivorella) (SPB) is a major cause of soybean (Glycine max L.) yield losses in northeast Asia, thus it is desirable to elucidate the resistance mechanisms involved in soybean response to the SPB. However, few studies have mapped SPB-resistant quantitative trait loci (QTLs) and deciphered the response mechanism in soybean. Here, we selected two soybean varieties, JY93 (SPB-resistant) and K6 (SPB-sensitive), to construct F₂ and F_2:3 populations for QTL mapping and collected pod shells before and after SPB larvae chewed on the two parents to perform RNA-Seq, which can identify stable QTLs and explore the response mechanism of soybean to the SPB. The results show that four QTLs underlying SPB damage to seeds were detected on chromosomes 4, 9, 13, and 15. Among them, qESP-9-1 was scanned in all environments, hence it can be considered a stable QTL. All QTLs explained 0.79 to 6.09% of the phenotypic variation. Meanwhile, 2298 and 3509 DEGs were identified for JY93 and K6, respectively, after the SPB attack, and most of these genes were upregulated. Gene Ontology enrichment results indicated that the SPB-induced and differently expressed genes in both parents are involved in biological processes such as wound response, signal transduction, immune response, and phytohormone pathways. Interestingly, secondary metabolic processes such as flavonoid synthesis were only significantly enriched in the upregulated genes of JY93 after SPB chewing compared with K6. Finally, we identified 18 candidate genes related to soybean pod borer resistance through the integration of QTL mapping and RNA-Seq analysis. Seven of these genes had similar expression patterns to the mapping parents in four additional soybean germplasm after feeding by the SPB. These results provide additional knowledge of the early response and induced defense mechanisms against the SPB in soybean, which could help in breeding SPB-resistant soybean accessions.     

Antistress Action of Melanocortin Derivatives Associated with Correction of Gene Expression Patterns in the Hippocampus of Male Rats Following Acute Stress

Natural melanocortins (MCs) have been used in the successful development of drugs with neuroprotective properties. Here, we studied the behavioral effects and molecular genetic mechanisms of two synthetic MC derivatives-ACTH(4–7)PGP (Semax) and ACTH(6–9)PGP under normal and acute restraint stress (ARS) conditions. Administration of Semax or ACTH(6–9)PGP (100 μg/kg) to rats 30 min before ARS attenuated ARS-induced behavioral alterations. Using high-throughput RNA sequencing (RNA-Seq), we identified 1359 differentially expressed genes (DEGs) in the hippocampus of vehicle-treated rats subjected to ARS, using a cutoff of >1.5 fold change and adjusted p-value (Padj) < 0.05, in samples collected 4.5 h after the ARS. Semax administration produced > 1500 DEGs, whereas ACTH(6–9)PGP administration led to <400 DEGs at 4.5 h after ARS. Nevertheless, ~250 overlapping DEGs were identified, and expression of these DEGs was changed unidirectionally by both peptides under ARS conditions. Modulation of the expression of genes associated with biogenesis, translation of RNA, DNA replication, and immune and nervous system function was produced by both peptides. Furthermore, both peptides upregulated the expression levels of many genes that displayed decreased expression after ARS, and vice versa, the MC peptides downregulated the expression levels of genes that were upregulated by ARS. Consequently, the antistress action of MC peptides may be associated with a correction of gene expression patterns that are disrupted during ARS.     

Identification and Fine Mapping of a Locus Related to Leaf Up-Curling Trait (Bnuc3) in Brassica napus

Leaf trait is an important target trait in crop breeding programs. Moderate leaf curling may be a help for improving crop yield by minimizing the shadowing by leaves. Mining locus for leaf curling trait is of significance for plant genetics and breeding researches. The present study identified a novel rapeseed accession with up-curling leaf, analyzed the up-curling leaf trait inheritance, and fine mapped the locus for up-curling leaf property (Bnuc3) in Brassica napus. Genetic analysis revealed that the up-curling leaf trait is controlled by a single dominant locus, named BnUC3. We performed an association study of BnUC3 with single nucleotide polymorphism (SNP) markers using a backcross population derived from the homozygous up-curling leaf line NJAU-M1295 and the canola variety ‘zhongshuang11’ with typical flat leaves, and mapped the BnUC3 locus in a 1.92 Mb interval of chromosome A02 of B. napus. To further map BnUC3, 232 simple sequence repeat (SSR) primers and four pairs of Insertion/Deletion (InDel) primers were developed for the mapping interval. Among them, five SSR markers and two InDel markers were polymorphic. By these markers, the mapping interval was narrowed to 92.0 kb using another F₂ population. This fine mapping interval has 11 annotated genes among which BnaA02T0157000ZS were inferred to be candidate casual genes for up-curling leaf based on the cloned sequence analysis, gene functionality, and gene expression analysis. The current study laid a foundational basis for further elucidating the mechanism of BnUC3 and breeding of variety with up-curling leaf.     

Integrated mRNA and miRNA Transcriptome Analysis Suggests a Regulatory Network for UV–B-Controlled Terpenoid Synthesis in Fragrant Woodfern (Dryopteris fragrans)

Fragrant woodfern (Dryopteris fragrans) is a medicinal plant rich in terpenoids. Ultraviolet-B (UV–B) light could increase concentration of terpenoids. The aim of this study was to analyze how UV–B regulates the terpenoid synthesis of the molecular regulatory mechanism in fragrant woodfern. In this study, compared with the control group, the content of the terpenes was significantly higher in fragrant woodfern leaves under UV–B treatment for 4 days (d). In order to identify how UV–B regulates the terpenoid metabolic mechanism in fragrant woodfern, we examined the mRNAs and small RNAs in fragrant woodfern leaves under UV–B treatment. mRNA and miRNA–seq identified 4533 DEGs and 17 DEMs in the control group compared with fragrant woodfern leaves under UV–B treatment for 4 d. mRNA–miRNA analysis identified miRNA target gene pairs consisting of 8 DEMs and 115 miRNAs. The target genes were subjected to GO and KEGG analyses. The results showed that the target genes were mainly enriched in diterpene biosynthesis, terpenoid backbone biosynthesis, plant hormone signal transduction, MEP pathway and MVA pathway, in which miR156 and miR160 regulate these pathways by targeting DfSPL and DfARF, respectively. The mRNA and miRNA datasets identified a subset of candidate genes. It provides the theoretical basis that UV–B regulates the terpenoid synthesis of the molecular regulatory mechanism in fragrant woodfern.     

Development and Validation of SNP and InDel Markers for Pod-Shattering Tolerance in Soybean

Pod-shattering causes a significant yield loss in many soybean cultivars. Shattering-tolerant cultivars provide the most effective approach to minimizing this loss. We developed molecular markers for pod-shattering and validated them in soybeans with diverse genetic backgrounds. The genes Glyma.16g141200, Glyma.16g141500, and Glyma.16g076600, identified in our previous study by quantitative trait locus (QTL) mapping and whole-genome resequencing, were selected for marker development. The whole-genome resequencing of three parental lines (one shattering-tolerant and two shattering-susceptible) identified single nucleotide polymorphism (SNP) and/or insertion/deletion (InDel) regions within or near the selected genes. Two SNPs and one InDel were converted to Kompetitive Allele-Specific PCR (KASP) and InDel markers, respectively. The accuracy of the markers was examined in the two recombinant inbred line populations used for the QTL mapping, as well as the 120 varieties and elite lines, through allelic discrimination and phenotyping by the oven-drying method. Both types of markers successfully discriminated the pod shattering-tolerant and shattering-susceptible genotypes. The prediction accuracy, which was as high as 90.9% for the RILs and was 100% for the varieties and elite lines, also supported the accuracy and usefulness of these markers. Thus, the markers can be used effectively for genetic and genomic studies and the marker-assisted selection for pod-shattering tolerance in soybean.     

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 a Candidate restorer-of-fertility Gene Rf3 Encoding a Pentatricopeptide Repeat Protein for the Cytoplasmic Male Sterility in Soybean

The cytoplasmic male sterility/restorer-of-fertility (CMS/Rf) system plays a vital role in high-efficiency hybrid seed production in crops, including soybean (Glycine max (L.) Merr.). The markers linked to fertility restoration and the restorer-of-fertility (Rf) genes are essential because they can facilitate the breeding of new CMS lines and production of commercial hybrid soybean seeds. To date, several soybean Rf genes have been mapped to various genetic loci in diverse genetic populations. However, the mapping range of restorer genes remains narrow, with relatively limited practical applicability. Therefore, in the present study, F₂ and F₃ segregating populations derived from the CMS line JLCMS5A crossed with the restorer line JLR2 were developed and used for Rf3 gene fine mapping. Genetic investigation indicated that the restorer line JLR2 was controlled by a single dominant gene, Rf3. By integrating bulk-segregant analysis and next-generation sequencing, a 4 Mb region on chromosome 9 was identified, which was most likely the target region harboring the candidate gene responsible for fertility restoration. This region was further narrowed down to 86.44 Kb via fine mapping in F₂ and F₃ populations using SSR, InDel, and dCAPS markers. This region contained 10 putative genes (Glyma.09G171100–Glyma.09G172000). Finally, Glyma.09G171200, which encodes a mitochondria-targeted pentatricopeptide repeat protein, was proposed as the potential candidate for Rf3 using sequence alignment and expression analysis in restorer and CMS lines. Based on single-nucleotide polymorphisms in Glyma.09G171200, a CAPS marker co-segregated with Rf3 named CAPS1712 was developed. Our results will be fundamental in the assisted selection and creation of potent lines for the production and rapid selection of novel restorer lines.     

Yield-Related QTL Clusters and the Potential Candidate Genes in Two Wheat DH Populations

In the present study, four large-scale field trials using two doubled haploid wheat populations were conducted in different environments for two years. Grain protein content (GPC) and 21 other yield-related traits were investigated. A total of 227 QTL were mapped on 18 chromosomes, which formed 35 QTL clusters. The potential candidate genes underlying the QTL clusters were suggested. Furthermore, adding to the significant correlations between yield and its related traits, correlation variations were clearly shown within the QTL clusters. The QTL clusters with consistently positive correlations were suggested to be directly utilized in wheat breeding, including 1B.2, 2A.2, 2B (4.9–16.5 Mb), 2B.3, 3B (68.9–214.5 Mb), 4A.2, 4B.2, 4D, 5A.1, 5A.2, 5B.1, and 5D. The QTL clusters with negative alignments between traits may also have potential value for yield or GPC improvement in specific environments, including 1A.1, 2B.1, 1B.3, 5A.3, 5B.2 (612.1–613.6 Mb), 7A.1, 7A.2, 7B.1, and 7B.2. One GPC QTL (5B.2: 671.3–672.9 Mb) contributed by cultivar Spitfire was positively associated with nitrogen use efficiency or grain protein yield and is highly recommended for breeding use. Another GPC QTL without negatively pleiotropic effects on 2A (50.0–56.3 Mb), 2D, 4D, and 6B is suggested for quality wheat breeding.