Microsatellite marker development, mapping and applications in rice genetics and breeding

Quantitative trait loci mapping for cadmium and microelements accumulation in rice using recombinant inbred lines

Cadmium (Cd) contamination in Oryza sativa L. rice grain poses a significant food safety risk. Cd uptake in rice occurs via transmembrane proteins that also transport divalent metal microelements such as iron (Fe), manganese (Mn), copper (Cu), and zinc (Zn). Identifying stable quantitative trait loci (QTLs) for Cd across different environments, along with co-localized QTLs for Cd and other elements, is crucial for understanding the genetic mechanisms of low Cd accumulation and facilitating fine mapping and gene cloning. A recombinant inbred line (RIL) population derived from D62B (a Cd-safe rice line) × Wujin4B (a high-Cd-accumulating rice line) was used to identify QTLs associated with Cd, Fe, Mn, Cu, and Zn accumulation. QTL mapping was conducted in three distinct environments to assess locus stability. A total of 61 QTLs were identified, including 23 QTLs related to Cd and 38 QTLs associated with other microelements. Two novel and stable QTLs for Cd (qLCA1 and qLCA12) were consistently detected across different environments, significantly reducing Cd concentration in brown rice and Cd translocation within the plant. Additionally, a previously unreported co-localized QTL (qCdA-4/qCuA-4) was identified, which significantly reduced both Cd and Cu accumulation in whole plants. This study identified two novel and stable QTLs for Cd (qLCA1 and qLCA12) that consistently reduced Cd concentration in brown rice and Cd translocation across different environments. Additionally, a newly discovered co-localized QTL (qCdA-4/qCuA-4) was found to regulate both Cd and Cu accumulation.

Genetic dissection for seedling root-related traits using multiple-methods in bread wheat (Triticum aestivum L.)

KEY MESSAGE

Several quantitative trait loci (QTL) and structural chromosome variations (SCVs) related to seedling root traits were identified using multiple methods, which provided valuable insights to assist breeding efforts in wheat. The root system of wheat affects water and fertilizer use efficiency, stress tolerance, and agronomic traits. Using association analysis and linkage mapping, QTL associated with 11 seedling-stage root traits were identified with single nucleotide polymorphisms (SNPs) and SCVs under both hydroponic nutrient solution culture experiment (NCE) and vermiculite culture experiment (VCE). Except for maximum root length (MRL), the root traits of seedlings under NCE and VCE differed significantly. Root fresh weight (RFW) and root dry weight (RDW) were significantly correlated with most agronomic traits and grain yield. Identification of RFW and RDW by NCE might provide a reference basis for VCE. Co-localization analysis revealed that NCE and VCE simultaneously detected SNP-loci viz. QRdw.sxau-6A, QRd.sxau-1B.2, and QDw.sxau-6A (5.56-8.76% of R2). The SCV-loci Mr1B-3, Mr3A-3 and Mr3A-4 were detected in both NCE and VCE (4.74-9.07% of R2). Furthermore, QRdw.sxau-6A, QSfw.sxau-6A and QRd.sxau-4A were detected using the mixed linear model (MLM), 3 Variance-component multi-locus random-SNP-effect Mixed Linear Mode (3VmrMLM) and rrBLUP. In the association panel, SNPs and SCVs co-localized to 14 MTAs, of which Mr5A-6 and QRd.sxau-5A were significantly associated with root diameter (RD). The association panel and doubled haploid (DH) population co-located 10 QTL, of which QDw.sxau-1D was stably detected. Finally, QDw.sxau-6A and Mg6A-9 overlapped in same genomic location containing candidate genes TraesCS6A02G372300, TraesCS6A02G382900 and TraesCS6A02G365100. The present study contributes novel insights into the genetics of root architecture in wheat.

Dynamic QTL mapping reveals the genetic architecture of stem diameter across developmental stages in foxtail millet

MAIN CONCLUSION

Dynamic QTL analysis identified the key QTL qSD-9-1 and candidate genes SiGPI8, SiCesA5 and SiNPC1 associated with the developmental stages of foxtail millet stem diameter. Abstact. Stem diameter (SD) is a critical agronomic trait influencing lodging and yield in foxtail millet, yet its potential in selective breeding remains under-explored. Additionally, limited research has examined the dynamics of SD across various developmental stages. To address this gap, this study utilized two RIL populations (F8), RYRIL and JYRIL, comprising 215 and 224 lines, respectively. Dynamic QTL analysis of SD traits was conducted in two distinct environments, encompassing five developmental stages, to comprehensively investigate the genetic architecture of SD. Results revealed that the hB2 values for SD ranged from 60.19% to 82.51%, peaking during the maturity stage in both populations. A total of 26 unconditional and 16 conditional QTLs were identified, explaining 1.09-17.99% and 3.10-16.55% of the phenotypic variations, respectively. Notably, qSD-9-1 was consistently detected in both unconditional and conditional QTL mappings across the two populations, accounting for 1.66-17.99% of phenotypic variation. Phenotype-genotype association analysis within the bin marker interval of qSD-9-1 revealed significant differences in SD among RIL lines carrying parental genotypes. Furthermore, by predicting candidate genes within the physical interval of qSD-9-1 and integrating phenotype-haplotype association analysis, SiGPI8, SiCesA5 and SiNPC1 were identified as key candidate genes potentially regulating SD. This study advances our understanding of the genetic basis of SD in foxtail millet and provides a theoretical foundation for marker-assisted selection (MAS) in breeding lodging-resistant varieties.

Identification of superior haplotypes and candidate gene for seed size-related traits in soybean (Glycine max L.)

Seed size is an economically important trait that directly determines the seed yield in soybean. In the current investigation, we used an integrated strategy of linkage mapping, association mapping, haplotype analysis and candidate gene analysis to determine the genetic makeup of four seed size-related traits viz., 100-seed weight (HSW), seed area (SA), seed length (SL), and seed width (SW) in soybean. Linkage mapping identified a total of 23 quantitative trait loci (QTL) associated with four seed size-related traits in the F2 population; among them, 17 were detected as novel QTLs, whereas the remaining six viz., qHSW3-1, qHSW4-1, qHSW18-1, qHSW19-1, qSL4-1 and qSW6-1 have been previously identified. Six out of 23 QTLs were major possessing phenotypic variation explained (PVE) ≥ 10%. Besides, the four QTL Clusters/QTL Hotspots harboring multiple QTLs for different seed size-related traits were identified on Chr.04, Chr.16, Chr.19 and Chr.20. Genome-wide association study (GWAS) identified a total of 62 SNPs significantly associated with the four seed size-related traits. Interestingly, the QTL viz., qHSW18-1 was identified by both linkage mapping and GWAS, and was regarded as the most stable loci regulating HSW in soybean. In-silico, sequencing and qRT-PCR analysis identified the Glyma.18G242400 as the most potential candidate gene underlying the qHSW18-1 for regulating HSW. Moreover, three haplotype blocks viz., Hap2, Hap6A and Hap6B were identified for the SW trait, and one haplotype was identified within the Glyma.18G242400 for the HSW. These four haplotypes harbor three to seven haplotype alleles across the association mapping panel of 350 soybean accessions, regulating the seed size from lowest to highest through intermediate phenotypes. Hence, the outcome of the current investigation can be utilized as a potential genetic and genomic resource for breeding the improved seed size in soybean.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01525-1.

Unraveling the genetic diversity, population structure, and stability for yield-related traits of rice genotypes in mid-hills of northeastern India

The cultivation of nearly 10,000 indigenous rice landraces in the North-Eastern Hill (NEH) region by various ethnic groups creates opportunities for the utilization of unique landraces through systematic evaluation of genetic variability. In the present study, a set of 102 rice landraces were assessed based on morphological and SSR markers, and five checks in augmented design vis-à-vis high-yielding rice genotypes with stable performance were identified. The presence of high estimates of heritability, genotypic coefficient of variation, and genetic advance over mean indicated the predominance of additive gene action, which necessitated the effectiveness of selection in augmenting productivity. A total of 83.73% of the total variation was accounted by the first five principal components. A total of 132 alleles were detected, with an average of 3 alleles per locus. The PIC values ranged from 0.01 to 0.70, with an average of 0.40. Based on FST value (5.1%), significant differences between the genotypes of Arunachal Pradesh and Sikkim were observed. The percentage of variation among the population, among individuals within the population, and within individuals was 5.14, 75.66, and 19.2%, respectively. Both Nei's genetic distance and model-based clustering have differentiated the genotypes into five distinct clusters. Principal coordinate analysis illustrated that the genotypes of Manipur were scattered in all quadrants, showing that they are highly diverse, while the genotypes of Nagaland, Sikkim, and Meghalaya were found together, which represent the chance of mixing of the population at a certain point in time. Markers, namely RM 474, OSR 13, RM 413, and RM 259, were found to be associated with key traits for increasing yielding ability of plant. In a stability evaluation based on AMMI analysis and multi-trait genotype-ideoptype distance matrix (MGIDI), genotypes, namely Jyotrirmayie, RCPL 1-411, Tsamum firri, Ching Phouren, Rato Bhan Joha, MN-47, and Tara bali, were selected with higher yield potential.

QTL mapping and candidate gene mining of seed size and seed weight in castor plant (Ricinus communis L.)

BACKGROUND

Castor (Ricinus communis L., 2n = 2x = 20) is an important industrial crop, due to its oil is very important to the global special chemical industry. Seed size and seed weight are fundamentally important in determining castor yield, while little is known about it. In this study, QTL analysis and candidate gene mining of castor seed size and seed weight were conducted with composite interval mapping (CIM), inclusive composite interval mapping (ICIM) and marker enrichment strategy in 4 populations, i.e., populations F2, BC1, S1-1 and S1-2, derived from 2 accessions with significant phenotypic differences.

RESULTS

In the QTL primary mapping, 2 novel QTL clusters were detected in marker intervals RCM520-RCM76 and RCM915-RCM950. In order to verify their accuracy and to narrow their intervals, QTL remapping was carried out in populations F2 and BC1. Among them, 44 and 30 QTLs underlying seed size and seed weight were detected in F2 population using methods CIM and ICIM-ADD respectively, including 4-9 and 3-5 ones conferring each trait were identified with a phenotypic variation explained ranged from 37.92 to 115.81% and 32.86-45.98% respectively. The remapping results in BC1 population were consistent with those in F2 population. Importantly, 3 QTL clusters (i.e. QTL-cluster1, QTL-cluster2 and QTL-cluster3) were found in marker intervals RCM74-RCM76 (37.1 kb), RCM930-RCM950 (259.8 kb) and RCM918-RCM920 (172.9 kb) respectively; in addition, all of them were detected again, the former one was found in the S1-2 population, and the latter two were found simultaneously in the populations S1-1 and S1-2. Finally, 6 candidate genes (i.e. LOC8266555, LOC8281168, LOC8281151, LOC8259066, LOC8258591 and LOC8270077) were screened in the above QTL clusters, they were differentially expressed in multiple seed tissues of both parents, signifying the potential role in regulating seed size and seed weight.

CONCLUSION

The above results not only provide new insights into the genetic structure of seed size and seed weight in castor, but also lay the foundation for the functional identification of these candidate genes.

Pomological and Molecular Characterization of Apple Cultivars in the German Fruit Genebank

Traditional varieties are a valuable tool in modern apple breeding. However, the use of synonyms and missing source documentation hinder an effective identification and conservation of relevant cultivars. During several projects, the authenticity and diversity of the apple cultivar collection of the German Fruit Genebank (GFG) was evaluated extensively. The trueness-to-type of 7890 apple trees was assessed on a pomological and molecular level. Pomological evaluations were performed by at least two experienced experts to identify the original cultivar names. On the molecular level, a set of 17 SSR markers was used to determine a unique genetic profile for each apple cultivar. The pomological and molecular characterization was expressed in terms of a comprehensive trueness-to-type criterion and the results were previously published as a well-curated dataset. In this study, the published dataset was analyzed to evaluate the quality and diversity of the apple collection of the GFG and highlight new findings based on phylogenetic and parentage analysis. The dataset contains 1404 unique genetic profiles corresponding to unambiguous cultivar names. Of these 1404 cultivars, 74% were assessed as true-to-type. The collection of diploid apple cultivars showed a high degree of expected heterozygosity (Hexp = 0.84). Genetic diversity in terms of year and location of origin was investigated with a STRUCTURE analysis. It was hypothesized that genetic diversity might decline overtime due to restrictive breeding programs. The results showed a shift dynamic between older and newer cultivars in one specific cluster, but no significant decrease in genetic diversity was observed in this study. Lastly, a parentage analysis was performed to check parental relationships based on historical research. Out of 128 parent-child trios, 110 trios resulted in significant relationships and reconfirmed the information from the literature. In some cases, the information from the literature was disproven. This analysis also allowed for readjusting the trueness-to-type criteria for previously undetermined cultivars. Overall, the importance of authenticity evaluations for gene bank cultivars was highlighted. Furthermore, the direct use of the dataset was shown by relevant investigations on the genetic diversity and structure of the apple cultivar collections of the GFG.

Quantitative Trait Loci Mapping of Heading Date in Wheat under Phosphorus Stress Conditions

Wheat (Triticum aestivum L.) is a crucial cereal crop, contributing around 20% of global caloric intake. However, challenges such as diminishing arable land, water shortages, and climate change threaten wheat production, making yield enhancement crucial for global food security. The heading date (HD) is a critical factor influencing wheat's growth cycle, harvest timing, climate adaptability, and yield. Understanding the genetic determinants of HD is essential for developing high-yield and stable wheat varieties. This study used a doubled haploid (DH) population from a cross between Jinmai 47 and Jinmai 84. QTL analysis of HD was performed under three phosphorus (P) treatments (low, medium, and normal) across six environments, using Wheat15K high-density SNP technology. The study identified 39 QTLs for HD, distributed across ten chromosomes, accounting for 2.39% to 29.52% of the phenotypic variance. Notably, five stable and major QTLs (Qhd.saw-3A.7, Qhd.saw-3A.8, Qhd.saw-3A.9, Qhd.saw-4A.4, and Qhd.saw-4D.3) were consistently detected across varying P conditions. The additive effects of these major QTLs showed that favorable alleles significantly delayed HD. There was a clear trend of increasing HD delay as the number of favorable alleles increased. Among them, Qhd.saw-3A.8, Qhd.saw-3A.9, and Qhd.saw-4D.3 were identified as novel QTLs with no prior reports of HD QTLs/genes in their respective intervals. Candidate gene analysis highlighted seven highly expressed genes related to Ca2+ transport, hormone signaling, glycosylation, and zinc finger proteins, likely involved in HD regulation. This research elucidates the genetic basis of wheat HD under P stress, providing critical insights for breeding high-yield, stable wheat varieties suited to low-P environments.

Identification of QTLs and candidate genes for water-soluble protein content in soybean seeds

Soybean represents a vital source of premium plant-based proteins for human nutrition. Importantly, the level of water-soluble protein (WSP) is crucial for determining the overall quality and nutritional value of such crops. Enhancing WSP levels in soybean plants is a high-priority goal in crop improvement. This study aimed to elucidate the genetic basis of WSP content in soybean seeds by identifying quantitative trait loci (QTLs) and set the foundation for subsequent gene cloning and functional analysis. Using 180 F10 recombinant inbred lines generated by crossing the high-protein soybean cultivar JiDou 12 with the wild variety Ye 9, our researcher team mapped the QTLs influencing protein levels, integrating Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and gene expression profiling to identify candidate genes. During the 2020 and 2022 growing seasons, a standard bell-shaped distribution of protein content trait data was observed in these soybean lines. Eight QTLs affecting protein content were found across eight chromosomes, with LOD scores ranging from 2.59 to 7.30, explaining 4.15-11.74% of the phenotypic variance. Notably, two QTLs were newly discovered, one with a elite allele at qWSPC-15 from Ye 9. The major QTL, qWSPC-19, on chromosome 19 was stable across conditions and contained genes involved in nitrogen metabolism, amino acid biosynthesis, and signaling. Two genes from this QTL, Glyma.19G185700 and Glyma.19G186000, exhibited distinct expression patterns at maturity, highlighting the influence of these genes on protein content. This research revealed eight QTLs for WSP content in soybean seeds and proposed a gene for the key QTL qWSPC-19, laying groundwork for gene isolation and enhanced soybean breeding through the use of molecular markers. These insights are instrumental for developing protein-rich soybean cultivars.

Mapping of dynamic quantitative trait loci for plant height in a RIL population of foxtail millet (Setaria italica L.)

Plant height (PH) is a crucial trait for strengthening lodging resistance and boosting yield in foxtail millet. To identify quantitative trait loci (QTL) and candidate genes associated with PH, we first developed a genetic map using a recombinant inbred line (RIL) population derived from a cross between Aininghuang and Jingu 21. Then, PH phenotyping data and four variations of best linear unbiased prediction (BLUP) were collected from nine environments and three development stages. Next, QTL mapping was conducted using both unconditional and conditional QTL methods. Subsequently, candidate genes were predicted via transcriptome analysis of parental samples at three developmental stages. The results revealed that the genetic map, based on re-sequencing, consisted of 4,360 bin markers spanning 1,016.06 cM with an average genetic distance of 0.23 cM. A total of 19 unconditional QTL, accounting for 5.23%-35.36% of the phenotypic variation explained (PVE), which included 7 major and 4 stable QTL, were identified. Meanwhile, 13 conditional QTL, explaining 5.88%-40.35% of PVE, including 5 major and 3 stable QTL, were discovered. Furthermore, four consistent and stable QTL were identified. Finally, eight candidate genes were predicted through RNA-seq and weighted gene co-expression network analysis (WGCNA). Those findings provide a crucial foundation for understanding the genetic mechanisms underlying PH development and facilitate molecular marker-assisted breeding of ideal plant types in foxtail millet.

Development of a molecular marker for the Pi1 gene based on the association of the SNAP protocol with the touch-up gradient amplification method

Genetic resistance is the most effective and eco-friendly approach to combat rice blast. The application of resistance genes may be facilitated by the availability of molecular markers that allow marker-assisted selection during the breeding process. The Pi1 gene, considered to be a broad-spectrum resistance gene, might contribute to enhancing resistance to rice blast, but it lacks a suitable marker that can be used. In this study, we investigated nucleotide polymorphism in the Pik locus and combined the SNAP protocol with the touch-up gradient amplification method to develop a SNAP marker. The Pi1 SNAP marker could distinguish Pi1 from Pik alleles, and when used for screening a germplasm bank and an F2 population, it consistently identified germplasms carrying the Pi1 gene. The P1 SNAP marker offers as advantages to involve only the presence/absence analysis of PCR amplicons resolved on an agarose gel.

Exploring Pan-Genomes: An Overview of Resources and Tools for Unraveling Structure, Function, and Evolution of Crop Genes and Genomes

The availability of multiple sequenced genomes from a single species made it possible to explore intra- and inter-specific genomic comparisons at higher resolution and build clade-specific pan-genomes of several crops. The pan-genomes of crops constructed from various cultivars, accessions, landraces, and wild ancestral species represent a compendium of genes and structural variations and allow researchers to search for the novel genes and alleles that were inadvertently lost in domesticated crops during the historical process of crop domestication or in the process of extensive plant breeding. Fortunately, many valuable genes and alleles associated with desirable traits like disease resistance, abiotic stress tolerance, plant architecture, and nutrition qualities exist in landraces, ancestral species, and crop wild relatives. The novel genes from the wild ancestors and landraces can be introduced back to high-yielding varieties of modern crops by implementing classical plant breeding, genomic selection, and transgenic/gene editing approaches. Thus, pan-genomic represents a great leap in plant research and offers new avenues for targeted breeding to mitigate the impact of global climate change. Here, we summarize the tools used for pan-genome assembly and annotations, web-portals hosting plant pan-genomes, etc. Furthermore, we highlight a few discoveries made in crops using the pan-genomic approach and future potential of this emerging field of study.

Combining transcriptomics and metabolomics to identify key response genes for aluminum toxicity in the root system of Brassica napus L. seedlings

By integrating QTL mapping, transcriptomics and metabolomics, 138 hub genes were identified in rapeseed root response to aluminum stress and mainly involved in metabolism of lipids, carbohydrates and secondary metabolites. Aluminum (Al) toxicity has become one of the important abiotic stress factors in areas with acid soil, which hinders the absorption of water and nutrients by roots, and consequently retards the growth of crops. A deeper understanding of the stress-response mechanism of Brassica napus may allow us to identify the tolerance gene(s) and use this information in breeding-resistant crop varieties. In this study, a population of 138 recombinant inbred lines (RILs) was subjected to aluminum stress, and QTL (quantitative trait locus) mapping was used to preliminarily locate quantitative trait loci related to aluminum stress. Root tissues from seedlings of an aluminum-resistant (R) line and an aluminum-sensitive (S) line from the RIL population were harvested for transcriptome sequencing and metabolome determination. By combining the data on quantitative trait genes (QTGs), differentially expressed genes (DEGs), and differentially accumulated metabolites (DAMs), key candidate genes related to aluminum tolerance in rapeseed were determined. The results showed that there were 3186 QTGs in the RIL population, 14,232 DEGs and 457 DAMs in the comparison between R and S lines. Lastly, 138 hub genes were selected to have a strong positive or negative correlation with 30 important metabolites (|R|≥ 0.95). These genes were mainly involved in the metabolism of lipids, carbohydrates and secondary metabolites in response to Al toxicity stress. In summary, this study provides an effective method for screening key genes by combining QTLs, transcriptome sequencing and metabolomic analysis, but also lists key genes for exploring the molecular mechanism of Al tolerance in rapeseed seedling roots.

Identification of major quantitative trait loci and candidate genes for seed weight in soybean

Four major quantitative trait loci for 100-seed weight were identified in a soybean RIL population under five environments, and the most likely candidate genes underlying these loci were identified. Seed weight is an important target of soybean breeding. However, the genes underlying the major quantitative trait loci (QTL) controlling seed weight remain largely unknown. In this study, a soybean population of 300 recombinant inbred lines (RILs) derived from a cross between PI595843 (PI) and WH was used to map the QTL and identify candidate genes for seed weight. The RIL population was genotyped through whole genome resequencing, and phenotyped for 100-seed weight under five environments. A total of 38 QTL were detected, and four major QTL, each explained at least 10% of the variation in 100-seed weight, were identified. Six candidate genes within these four major QTL regions were identified by analyses of their tissue expression patterns, gene annotations, and differential gene expression levels in soybean seeds during four developmental stages between two parental lines. Further sequence variation analyses revealed a C to T substitution in the first exon of the Glyma.19G143300, resulting in an amino acid change between PI and WH, and thus leading to a different predicted kinase domain, which might affect its protein function. Glyma.19G143300 is highly expressed in soybean seeds and encodes a leucine-rich repeat receptor-like protein kinase (LRR-RLK). Its predicted protein has typical domains of LRR-RLK family, and phylogenetic analyses reveled its similarity with the known LRR-RLK protein XIAO (LOC_Os04g48760), which is involved in controlling seed size. The major QTL and candidate genes identified in this study provide useful information for molecular breeding of new soybean cultivars with desirable seed weight.

Identification and validation of quantitative trait loci for chlorophyll content of flag leaf in wheat under different phosphorus treatments

In wheat, the leaf chlorophyll content in flag leaves is closely related to the degree of phosphorus stress. Identifying major genes/loci associated with chlorophyll content in flag leaves under different phosphorus conditions is critical for breeding wheat varieties resistant to low phosphorus (P). Under normal, medium, and low phosphorus conditions, the chlorophyll content of flag leaves was investigated by a double haploid (DH) population derived from a cross between two popular wheat varieties Jinmai 47 and Jinmai 84, at different grain filling stages. Chlorophyll content of the DH population and parents decreased gradually during the S1 to the S3 stages and rapidly at the S4 stage. At the S4 stage, the chlorophyll content of the DH population under low phosphorus conditions was significantly lower than under normal phosphate conditions. Using a wheat 15K single-nucleotide polymorphism (SNP) panel, a total of 157 QTLs were found to be associated with chlorophyll content in flag leaf and were identified under three phosphorus conditions. The phenotypic variation explained (PVE) ranged from 3.07 to 31.66%. Under three different phosphorus conditions, 36, 30, and 48 QTLs for chlorophyll content were identified, respectively. Six major QTLs Qchl.saw-2B.1, Qchl.saw-3B.1, Qchl.saw-4D.1, Qchl.saw-4D.2, Qchl.saw-5A.9 and Qchl.saw-6A.4 could be detected under multiple phosphorus conditions in which Qchl.saw-4D.1, Qchl.saw-4D.2, and Qchl.saw-6A.4 were revealed to be novel major QTLs. Moreover, the closely linked SNP markers of Qchl.saw-4D.1 and Qchl.saw-4D.2 were validated as KASP markers in a DH population sharing the common parent Jinmai 84, showed extreme significance (P <0.01) in more than three environments under different phosphorus conditions, which has the potential to be utilized in molecular marker-assisted breeding for low phosphorus tolerance in wheat.

Assessment of the Genetic Distinctiveness and Uniformity of Pre-Basic Seed Stocks of Italian Ryegrass Varieties

Lolium multiflorum Lam., commonly known as Italian ryegrass, is a forage grass mostly valued for its high palatability and digestibility, along with its high productivity. However, Italian ryegrass has an outbreeding nature and therefore has high genetic heterogeneity within each variety. Consequently, the exclusive use of morphological descriptors in the existing varietal identification and registration process based on the Distinctness, Uniformity, and Stability (DUS) test results in an inadequately precise assessment. The primary objective of this work was to effectively test whether the uniformity observed at the phenological level within each population of Italian ryegrass was confirmed at the genetic level through an SSR marker analysis. In this research, using 12 polymorphic SSR loci, we analyzed 672 samples belonging to 14 different Italian ryegrass commercial varieties to determine the pairwise genetic similarity (GS), verified the distribution of genetic diversity within and among varieties, and investigated the population structure. Although the fourteen commercial varieties did not show elevated genetic differentiation, with only 13% of the total variation attributable to among-cultivar genetic variation, when analyzed as a core, each variety constitutes a genetic cluster on its own, resulting in distinct characteristics from the others, except for two varieties. In this way, by combining a genetic tool with the traditional morphological approach, we were able to limit biases linked to the environmental effect of field trials, assessing the real source of diversity among varieties and concretely answering the key requisites of the Plant Variety Protection (PVP) system.

Dynamic Quantitative Trait Loci Mapping for Plant Height in Recombinant Inbred Line Population of Upland Cotton

Plant height (PH) is a key plant architecture trait for improving the biological productivity of cotton. Ideal PH of cotton is conducive to lodging resistance and mechanized harvesting. To detect quantitative trait loci (QTL) and candidate genes of PH in cotton, a genetic map was constructed with a recombinant inbred line (RIL) population of upland cotton. PH phenotype data under nine environments and three best linear unbiased predictions (BLUPs) were used for QTL analyses. Based on restriction-site-associated DNA sequence (RAD-seq), the genetic map contained 5,850 single-nucleotide polymorphism (SNP) markers, covering 2,747.12 cM with an average genetic distance of 0.47 cM. Thirty-seven unconditional QTL explaining 1.03-12.50% of phenotypic variance, including four major QTL and seven stable QTL, were identified. Twenty-eight conditional QTL explaining 3.27-28.87% of phenotypic variance, including 1 major QTL, were identified. Importantly, five QTL, including 4 stable QTL, were both unconditional and conditional QTL. Among the 60 PH QTL (including 39 newly identified), none of them were involved in the whole period of PH growth, indicating that QTL related to cotton PH development have dynamic expression characteristics. Based on the functional annotation of Arabidopsis homologous genes and transcriptome data of upland cotton TM-1, 14 candidate genes were predicted within 10 QTL. Our research provides valuable information for understanding the genetic mechanism of PH development, which also increases the economic production of cotton.

Identification of quantitative trait loci (QTLs) and candidate genes of seed Iron and zinc content in soybean [Glycine max (L.) Merr.]

BACKGROUND

Deciphering the hereditary mechanism of seed iron (Fe) and zinc (Zn) content in soybean is important and sustainable to address the "hidden hunger" that presently affects approximately 2 billion people worldwide. Therefore, in order to detect genomic regions related to soybean seed Fe and Zn content, a recombinant inbred line (RIL) population with 248 lines was assessed in four environments to detect Quantitative Trait Loci (QTLs) related to soybean seed Fe and Zn content.

RESULT

Wide variation was found in seed Fe and Zn content in four environments, and genotype, environment, and genotype × environment interactions had significant influences on both the seed Fe and Zn content. A positive correlation was observed between seed Fe content and seed Zn content, and broad-sense heritability (H2) of seed Fe and Zn content were 0.73 and 0.75, respectively. In this study, five QTLs for seed Fe content were detected with 4.57 - 32.71% of phenotypic variation explained (PVE) and logarithm of odds (LOD) scores ranging from 3.60 to 33.79. Five QTLs controlling the seed Zn content were detected, and they individually explained 3.35 to 26.48% of the phenotypic variation, with LOD scores ranging from 3.64 to 20.4. Meanwhile, 409,541 high-quality single-nucleotide variants (SNVs) and 85,102 InDels (except intergenic regions) between two bi-parental lines were identified by whole genome resequencing. A total of 12 candidate genes were reported in one major QTL for seed Fe content and two major QTLs for seed Zn content, with the help of RNA-Seq analysis, gene ontology (GO) enrichment, gene annotation, and bi-parental whole genome sequencing (WGS) data.

CONCLUSIONS

Limited studies were performed about microelement of soybean, so these results may play an important role in the biofortification of Fe and Zn and accelerate the development of marker-assisted selection (MAS) for breeding soybeans fortified with iron and zinc.

Genome survey sequencing and characterization of simple sequence repeat (SSR) markers in Platostoma palustre (Blume) A.J.Paton (Chinese mesona)

Platostoma palustre (Blume) A.J.Paton is an annual herbaceous persistent plant of the Labiatae family. However, there is a lack of genomic data for this plant, which severely restricts its genetic improvement. In this study, we performed genome survey sequencing of P. palustre and developed simple sequence repeat (SSR) markers based on the resulting sequence. K-mer analysis revealed that the assembled genome size was approximately 1.21 Gb. A total of 15,498 SSR motifs were identified and characterized in this study; among them, dinucleotide, and hexanucleotide repeats had the highest and lowest, respectively. Among the dinucleotide repeat motifs, AT/TA repeat motifs were the most abundant, and GC/CG repeat motifs were rather rare, accounting for 44.28% and 0.63%, respectively. Genetic similarity coefficient analysis by the UPMGA methods clustered 12 clones, of P. palustre and related species into two subgroups. These results provide helpful information for further research on P. palustre resources and variety improvements.

A novel miR167a-OsARF6-OsAUX3 module regulates grain length and weight in rice

Grain size is one of the most important factors that control rice yield, as it is associated with grain weight (GW). To date, dozens of rice genes that regulate grain size have been isolated; however, the regulatory mechanism underlying GW control is not fully understood. Here, the quantitative trait locus qGL5 for grain length (GL) and GW was identified in recombinant inbred lines of 9311 and Nipponbare (NPB) and fine mapped to a candidate gene, OsAUX3. Sequence variations between 9311 and NPB in the OsAUX3 promoter and loss of function of OsAUX3 led to higher GL and GW. RNA sequencing, gene expression quantification, dual-luciferase reporter assays, chromatin immunoprecipitation-quantitative PCR, and yeast one-hybrid assays demonstrated that OsARF6 is an upstream transcription factor regulating the expression of OsAUX3. OsARF6 binds directly to the auxin response elements of the OsAUX3 promoter, covering a single-nucleotide polymorphism site between 9311 and NPB/Dongjin/Hwayoung, and thereby controls GL by altering longitudinal expansion and auxin distribution/content in glume cells. Furthermore, we showed that miR167a positively regulate GL and GW by directing OsARF6 mRNA silencing. Taken together, our study reveals that a novel miR167a-OsARF6-OsAUX3 module regulates GL and GW in rice, providing a potential target for the improvement of rice yield.

GhUBX controlling helical growth results in production of stronger cotton fiber

Cotton fiber is an excellent model for studying plant cell elongation and cell wall biogenesis as well because they are highly polarized and use conserved polarized diffuse growth mechanism. Fiber strength is an important trait among cotton fiber qualities due to ongoing changes in spinning technology. However, the molecular mechanism of fiber strength forming is obscure. Through map-based cloning, we identified the fiber strength gene GhUBX. Increasing its expression, the fiber strength of the transgenic cotton was significantly enhanced compared to the receptor W0 and the helices number of the transgenic fiber was remarkably increased. Additionally, we proved that GhUBX regulates the fiber helical growth by degrading the GhSPL1 via the ubiquitin 26S–proteasome pathway. Taken together, we revealed the internal relationship between fiber helices and fiber stronger. It will be useful for improving the fiber quality in cotton breeding and illustrating the molecular mechanism for plant twisted growth.

•

Isolation of the first fiber strength gene GhUBX using map-based cloning strategy

•

Verification of the function of GhUBX experimentally in transgenic cotton

•

Link helices to the cotton fiber strength, that more helices make fiber stronger

•

An ubiquitin–proteasome system regulating the development of cotton fiber

Identification and QTL Analysis of Flavonoids and Carotenoids in Tetraploid Roses Based on an Ultra-High-Density Genetic Map

Roses are highly valuable within the flower industry. The metabolites of anthocyanins, flavonols, and carotenoids in rose petals are not only responsible for the various visible petal colors but also important bioactive compounds that are important for human health. In this study, we performed a QTL analysis on pigment contents to locate major loci that determine the flower color traits. An F1 population of tetraploid roses segregating for flower color was used to construct an ultra-high-density genetic linkage map using whole-genome resequencing technology to detect genome-wide SNPs. Previously developed SSR and SNP markers were also utilized to increase the marker density. Thus, a total of 9,259 markers were mapped onto seven linkage groups (LGs). The final length of the integrated map was 1285.11 cM, with an average distance of 0.14 cM between adjacent markers. The contents of anthocyanins, flavonols and carotenoids of the population were assayed to enable QTL analysis. Across the 33 components, 46 QTLs were detected, explaining 11.85-47.72% of the phenotypic variation. The mapped QTLs were physically clustered and primarily distributed on four linkage groups, namely LG2, LG4, LG6, and LG7. These results improve the basis for flower color marker-assisted breeding of tetraploid roses and guide the development of rose products.

Identify of Fast-Growing Related Genes Especially in Height Growth by Combining QTL Analysis and Transcriptome in Salix matsudana (Koidz)

The study on the fast-growing traits of trees, mainly valued by tree height (TH) and diameter at breast height (DBH), is of great significance to promote the development of the forest industry. Quantitative trait locus (QTL) mapping based on high-density genetic maps is an efficient approach to identify genetic regions for fast-growing traits. In our study, a high-density genetic map for the F₁ population was constructed. The genetic map had a total size of 5,484.07 centimorgan (cM), containing 5,956 single nucleotide polymorphisms (SNPs) based on Specific Length Amplified Fragment sequencing. Six fast-growing related stable QTL were identified on six chromosomes, and five stable QTL were identified by a principal component analysis (PCA). By combining the RNA-seq analysis for the two parents and two progenies with the qRT-PCR analysis, four candidate genes, annotated as DnaJ, 1-aminocyclopropane-1-carboxylate oxidase 1 (ACO1), Caffeic acid 3-O-methyltransferase 1 (COMT1), and Dirigent protein 6 (DIR6), that may regulate height growth were identified. Several lignin biosynthesis-related genes that may take part in height growth were detected. In addition, 21 hotspots in this population were found. The results of this study will provide an important foundation for further studies on the molecular and genetic regulation of TH and DBH.

Complete chloroplast genome of a rare and endangered plant species Osteomeles subrotunda: genomic features and phylogenetic relationships with other Rosaceae plants

Osteomeles subrotunda is a rare and endangered plant species with extremely small populations. In our study, we sequenced the complete chloroplast (CP) genome of O. subrotunda and described its structural organization, and performed comparative genomic analyses with other Rosaceae CP genomes. The plastome of O. subrotunda was 159,902 bp in length with 36.6% GC content and contained a pair of inverted repeats of 26,367 bp which separated a large single-copy region of 87,933 bp and a small single-copy region of 19,235 bp. The CP genome included 130 genes, of which 85 were protein-coding genes, 37 were transfer RNAs, and eight were ribosomal RNAs. Two genes, rps19 and ycf1, which are located at the borders of IRB/SSC and IRB/LSC, were presumed to be pseudogenes. A total of 61 SSRs were detected, of which, 59 loci were mono-nucleotide repeats, and two were di-nucleotide repeats. The phylogenic analysis indicated that the 14 Rosaceae species were divided into three groups, among which O. subrotunda grouped with P. rupicola, E. japonica, P. pashia, C. japonica, S. torminalis, and M. florentina, and it was found to be a sister clade to C. japonica. Our newly sequenced CP genome of O. subrotunda will provide essential data for further studies on population genetics and biodiversity.

Identification of quantitative trait nucleotides and candidate genes for soybean seed weight by multiple models of genome-wide association study

BACKGROUND

Seed weight is a complex yield-related trait with a lot of quantitative trait loci (QTL) reported through linkage mapping studies. Integration of QTL from linkage mapping into breeding program is challenging due to numerous limitations, therefore, Genome-wide association study (GWAS) provides more precise location of QTL due to higher resolution and diverse genetic diversity in un-related individuals.

RESULTS

The present study utilized 573 breeding lines population with 61,166 single nucleotide polymorphisms (SNPs) to identify quantitative trait nucleotides (QTNs) and candidate genes for seed weight in Chinese summer-sowing soybean. GWAS was conducted with two single-locus models (SLMs) and six multi-locus models (MLMs). Thirty-nine SNPs were detected by the two SLMs while 209 SNPs were detected by the six MLMs. In all, two hundred and thirty-one QTNs were found to be associated with seed weight in YHSBLP with various effects. Out of these, seventy SNPs were concurrently detected by both SLMs and MLMs on 8 chromosomes. Ninety-four QTNs co-localized with previously reported QTL/QTN by linkage/association mapping studies. A total of 36 candidate genes were predicted. Out of these candidate genes, four hub genes (Glyma06g44510, Glyma08g06420, Glyma12g33280 and Glyma19g28070) were identified by the integration of co-expression network. Among them, three were relatively expressed higher in the high HSW genotypes at R5 stage compared with low HSW genotypes except Glyma12g33280. Our results show that using more models especially MLMs are effective to find important QTNs, and the identified HSW QTNs/genes could be utilized in molecular breeding work for soybean seed weight and yield.

CONCLUSION

Application of two single-locus plus six multi-locus models of GWAS identified 231 QTNs. Four hub genes (Glyma06g44510, Glyma08g06420, Glyma12g33280 & Glyma19g28070) detected via integration of co-expression network among the predicted candidate genes.

Functional Markers for Precision Plant Breeding

Advances in molecular biology including genomics, high-throughput sequencing, and genome editing enable increasingly faster and more precise cultivar development. Identifying genes and functional markers (FMs) that are highly associated with plant phenotypic variation is a grand challenge. Functional genomics approaches such as transcriptomics, targeting induced local lesions in genomes (TILLING), homologous recombinant (HR), association mapping, and allele mining are all strategies to identify FMs for breeding goals, such as agronomic traits and biotic and abiotic stress resistance. The advantage of FMs over other markers used in plant breeding is the close genomic association of an FM with a phenotype. Thereby, FMs may facilitate the direct selection of genes associated with phenotypic traits, which serves to increase selection efficiencies to develop varieties. Herein, we review the latest methods in FM development and how FMs are being used in precision breeding for agronomic and quality traits as well as in breeding for biotic and abiotic stress resistance using marker assisted selection (MAS) methods. In summary, this article describes the use of FMs in breeding for development of elite crop cultivars to enhance global food security goals.

Genetic Analysis for Cooking and Eating Quality of Super Rice and Fine Mapping of a Novel Locus qGC10 for Gel Consistency

Rice (Oryza sativa L.) is an important cereal that provides food for more than half of the world's population. Besides grain yield, improving grain quality is also essential to rice breeders. Amylose content (AC), gelatinization temperature (GT) and gel consistency (GC) are considered to be three indicators for cooking and eating quality in rice. Using a genetic map of RILs derived from the super rice Liang-You-Pei-Jiu with high-density SNPs, we detected 3 QTLs for AC, 3 QTLs for GT, and 8 QTLs for GC on chromosomes 3, 4, 5, 6, 10, and 12. Wx locus, an important determinator for AC and GC, resided in one QTL cluster for AC and GC, qAC6 and qGC6 here. And a novel major QTL qGC10 on chromosome 10 was identified in both Lingshui and Hangzhou. With the BC₄F₂ population derived from a CSSL harboring the segment for qGC10 from 93-11 in PA64s background, it was fine mapped between two molecular markers within 181 kb region with 27 annotated genes. Quantitative real-time PCR results showed that eight genes were differentially expressed in endosperm of two parents. After DNA sequencing, only LOC_Os10g04900, which encodes a F-box domain containing protein, has 2 bp deletion in the exon of PA64s, resulting in a premature stop codon. Therefore, LOC_Os10g04900 is considered to be the most likely candidate gene for qGC10 associated with gel consistency. Identification of qGC10 provides a new genetic resource for improvement of rice quality.

Genetic Mapping Combined with a Transcriptome Analysis to Screen for Candidate Genes Responsive to Abscisic Acid Treatment in Brassica napus Embryos During Seed Germination

Brassica napus embryos contain precursor tissues for the leaves, stem, and root, as well as the cotyledons, and these precursor tissues play key roles in seed germination, seedling survival, and subsequent seedling growth. Abscisic acid (ABA) plays a prominent role in the inhibition of seed germination. The underlying molecular mechanisms of the embryo responses to ABA stress followed by inhibited seed germination have not been reported in B. napus to date. In this study, we conducted quantitative trait locus (QTL) analysis of B. napus seed in response to ABA stress using 170 recombinant inbred lines. Furthermore, we performed transcriptome sequencing (RNA-seq) analyses by using B. napus ZS11 embryos under sterile deionized water (control) and 10 mg/L (10A), 20 mg/L (20A), and 30 mg/L (30A) ABA treatment conditions. In total, 10 QTLs were screened for explaining 2.70-6.73% of the phenotypic variation under ABA stress. In addition, 1495, 3332, and 3868 differentially expressed genes (DEGs) were identified in the "control vs 10A," "control vs 20A," and "control vs 30A" comparisons, respectively. Gene Ontology (GO) enrichment analysis indicated that DEG functions are mainly related to response to stimuli, response to oxygen-containing compounds, response to lipids, and the transport and seed dormancy processes. These DEGs mainly participated in the response to plant hormone signal transduction, starch and sucrose metabolism, cutin, suberine, and wax biosynthesis, and phenylpropanoid biosynthesis processes pathways. Our results provide a foundation for further explorations of the molecular regulatory mechanisms of B. napus embryos in response to abiotic stress during the seed germination stage.

High-Resolution Mapping in Two RIL Populations Refines Major "QTL Hotspot" Regions for Seed Size and Shape in Soybean (Glycine max L.)

Seed size and shape are important traits determining yield and quality in soybean. However, the genetic mechanism and genes underlying these traits remain largely unexplored. In this regard, this study used two related recombinant inbred line (RIL) populations (ZY and K3N) evaluated in multiple environments to identify main and epistatic-effect quantitative trait loci (QTLs) for six seed size and shape traits in soybean. A total of 88 and 48 QTLs were detected through composite interval mapping (CIM) and mixed-model-based composite interval mapping (MCIM), respectively, and 15 QTLs were common among both methods; two of them were major (R2 > 10%) and novel QTLs (viz., qSW-1-1ZN and qSLT-20-1K3N). Additionally, 51 and 27 QTLs were identified for the first time through CIM and MCIM methods, respectively. Colocalization of QTLs occurred in four major QTL hotspots/clusters, viz., "QTL Hotspot A", "QTL Hotspot B", "QTL Hotspot C", and "QTL Hotspot D" located on Chr06, Chr10, Chr13, and Chr20, respectively. Based on gene annotation, gene ontology (GO) enrichment, and RNA-Seq analysis, 23 genes within four "QTL Hotspots" were predicted as possible candidates, regulating soybean seed size and shape. Network analyses demonstrated that 15 QTLs showed significant additive x environment (AE) effects, and 16 pairs of QTLs showing epistatic effects were also detected. However, except three epistatic QTLs, viz., qSL-13-3ZY, qSL-13-4ZY, and qSW-13-4ZY, all the remaining QTLs depicted no main effects. Hence, the present study is a detailed and comprehensive investigation uncovering the genetic basis of seed size and shape in soybeans. The use of a high-density map identified new genomic regions providing valuable information and could be the primary target for further fine mapping, candidate gene identification, and marker-assisted breeding (MAB).

Target sequencing reveals genetic diversity, population structure, core-SNP markers, and fruit shape-associated loci in pepper varieties

BACKGROUND

The widely cultivated pepper (Capsicum spp.) is one of the most diverse vegetables; however, little research has focused on characterizing the genetic diversity and relatedness of commercial varieties grown in China. In this study, a panel of 92 perfect single-nucleotide polymorphisms (SNPs) was identified using re-sequencing data from 35 different C. annuum lines. Based on this panel, a Target SNP-seq genotyping method was designed, which combined multiplex amplification of perfect SNPs with Illumina sequencing, to detect polymorphisms across 271 commercial pepper varieties.

RESULTS

The perfect SNPs panel had a high discriminating capacity due to the average value of polymorphism information content, observed heterozygosity, expected heterozygosity, and minor allele frequency, which were 0.31, 0.28, 0.4, and 0.31, respectively. Notably, the studied pepper varieties were morphologically categorized based on fruit shape as blocky-, long horn-, short horn-, and linear-fruited. The long horn-fruited population exhibited the most genetic diversity followed by the short horn-, linear-, and blocky-fruited populations. A set of 35 core SNPs were then used as kompetitive allele-specific PCR (KASPar) markers, another robust genotyping technique for variety identification. Analysis of genetic relatedness using principal component analysis and phylogenetic tree construction indicated that the four fruit shape populations clustered separately with limited overlaps. Based on STRUCTURE clustering, it was possible to divide the varieties into five subpopulations, which correlated with fruit shape. Further, the subpopulations were statistically different according to a randomization test and F_st statistics. Nine loci, located on chromosomes 1, 2, 3, 4, 6, and 12, were identified to be significantly associated with the fruit shape index (p < 0.0001).

CONCLUSIONS

Target SNP-seq developed in this study appears as an efficient power tool to detect the genetic diversity, population relatedness and molecular breeding in pepper. Moreover, this study demonstrates that the genetic structure of Chinese pepper varieties is significantly influenced by breeding programs focused on fruit shape.

The indica nitrate reductase gene OsNR2 allele enhances rice yield potential and nitrogen use efficiency

The indica and japonica rice (Oryza sativa) subspecies differ in nitrate (NO₃⁻) assimilation capacity and nitrogen (N) use efficiency (NUE). Here, we show that a major component of this difference is conferred by allelic variation at OsNR2, a gene encoding a NADH/NADPH-dependent NO₃⁻ reductase (NR). Selection-driven allelic divergence has resulted in variant indica and japonica OsNR2 alleles encoding structurally distinct OsNR2 proteins, with indica OsNR2 exhibiting greater NR activity. Indica OsNR2 also promotes NO₃⁻ uptake via feed-forward interaction with OsNRT1.1B, a gene encoding a NO₃⁻ uptake transporter. These properties enable indica OsNR2 to confer increased effective tiller number, grain yield and NUE on japonica rice, effects enhanced by interaction with an additionally introgressed indica OsNRT1.1B allele. In consequence, indica OsNR2 provides an important breeding resource for the sustainable increases in japonica rice yields necessary for future global food security.

Indica rice has higher nitrate assimilation and nitrogen use efficiency (NUE) than japonica rice, but the mechanism is unclear. Here, the authors reveal that the difference is partly due to allelic variation of a nitrate reductase encoding gene and this indica allele can increase yield potential and NUE.

Variability Assessment for Root and Drought Tolerance Traits and Genetic Diversity Analysis of Rice Germplasm using SSR Markers

The studies on genetic variation, diversity and population structure of rice germplasm of North East India could be an important step for improvements of abiotic and biotic stress tolerance in rice. Genetic diversity and genetic relatedness among 114 rice genotypes of North East India were assessed using genotypic data of 65 SSR markers and phenotypic data. The phenotypic diversity analysis showed the considerable variation across genotypes for root, shoot and drought tolerance traits. The principal component analysis (PCA) revealed the fresh shoot weight, root volume, dry shoot weight, fresh root weight and drought score as a major contributor to diversity. Genotyping of 114 rice genotypes using 65 SSR markers detected 147 alleles with the average polymorphic information content (PIC) value of 0.51. Population structure analysis using the Bayesian clustering model approach, distance-based neighbor-joining cluster and principal coordinate analysis using genotypic data grouped the accession into three sub-populations. Population structure analysis revealed that rice accession was moderately structured based on FST value estimates. Analysis of molecular variance (AMOVA) and pairwise FST values showed significant differentiation among all the pairs of sub-population ranging from 0.152 to 0.222 suggesting that all the three subpopulations were significantly different from each other. AMOVA revealed that most of the variation in rice accession mainly occurred among individuals. The present study suggests that diverse germplasm of NE India could be used for the improvement of root and drought tolerance in rice breeding programmes.

Effect of marker segregation distortion on high density linkage map construction and QTL mapping in Soybean (Glycine max L.)

Marker segregation distortion is a natural phenomenon. Severely distorted markers are usually excluded in the construction of linkage maps. We investigated the effect of marker segregation distortion on linkage map construction and quantitative trait locus (QTL) mapping. A total of 519 recombinant inbred lines of soybean from orthogonal and reciprocal crosses between LSZZH and NN493-1 were genotyped by specific length amplified fragment markers and seed linoleic acid content was measured in three environments. As a result, twenty linkage groups were constructed with 11,846 markers, including 1513 (12.77%) significantly distorted markers, on 20 chromosomes, and the map length was 2475.86 cM with an average marker-interval of 0.21 cM. The inclusion of distorted markers in the analysis was shown to not only improve the grouping of the markers from the same chromosomes, and the consistency of linkage maps with genome, but also increase genome coverage by markers. Combining genotypic data from both orthogonal and reciprocal crosses decreased the proportion of distorted markers and then improved the quality of linkage maps. Validation of the linkage maps was confirmed by the high collinearity between positions of markers in the soybean reference genome and in linkage maps and by the high consistency of 24 QTL regions in this study compared with the previously reported QTLs and lipid metabolism related genes. Additionally, linkage maps that include distorted markers could add more information to the outputs from QTL mapping. These results provide important information for linkage mapping, gene cloning and marker-assisted selection in soybean.

Utilization of Interspecific High-Density Genetic Map of RIL Population for the QTL Detection and Candidate Gene Mining for 100-Seed Weight in Soybean

Seed-weight is one of the most important traits determining soybean yield. Hence, it is prerequisite to have detailed understanding of the genetic basis regulating seed-weight for the development of improved cultivars. In this regard, the present study used high-density interspecific linkage map of NJIR4P recombinant inbred population evaluated in four different environments to detect stable Quantitative trait loci (QTLs) as well as mine candidate genes for 100-seed weight. In total, 19 QTLs distributed on 12 chromosomes were identified in all individual environments plus combined environment, out of which seven were novel and eight are stable identified in more than one environment. However, all the novel QTLs were minor (R 2 < 10%). The remaining 12 QTLs detected in this study were co-localized with the earlier reported QTLs with narrow genomic regions, and out of these only 2 QTLs were major (R 2 > 10%) viz., qSW-17-1 and qSW-17-4. Beneficial alleles of all identified QTLs were derived from cultivated soybean parent (Nannong493-1). Based on Protein ANalysis THrough Evolutionary Relationships, gene annotation information, and literature search, 29 genes within 5 stable QTLs were predicted to be possible candidate genes that might regulate seed-weight/size in soybean. However, it needs further validation to confirm their role in seed development. In conclusion, the present study provides better understanding of trait genetics and candidate gene information through the use high-density inter-specific bin map, and also revealed considerable scope for genetic improvement of 100-seed weight in soybean using marker-assisted breeding.

Whole-genome mapping identified novel "QTL hotspots regions" for seed storability in soybean (Glycine max L.)

BACKGROUND

Seed aging in soybean is a serious challenge for agronomic production and germplasm preservation. However, its genetic basis remains largely unclear in soybean. Unraveling the genetic mechanism involved in seed aging, and enhancing seed storability is an imperative goal for soybean breeding. The aim of this study is to identify quantitative trait loci (QTLs) using high-density genetic linkage maps of soybean for seed storability. In this regard, two recombinant inbred line (RIL) populations derived from Zhengyanghuangdou × Meng 8206 (ZM6) and Linhefenqingdou × Meng 8206 (LM6) crosses were evaluated for three seed-germination related traits viz., germination rate (GR), normal seedling length (SL) and normal seedling fresh weight (FW) under natural and artificial aging conditions to map QTLs for seed storability.

RESULTS

A total of 34 QTLs, including 13 QTLs for GR, 11 QTLs for SL and 10 QTLs for FW, were identified on 11 chromosomes with the phenotypic variation ranged from 7.30 to 23.16% under both aging conditions. All these QTLs were novel, and 21 of these QTLs were clustered in five QTL-rich regions on four different chromosomes viz., Chr3, Chr5, Chr17 &Chr18, among them the highest concentration of seven and six QTLs were found in "QTL hotspot A" (Chr17) and "QTL hotspot B" (Chr5), respectively. Furthermore, QTLs within all the five QTL clusters are linked to at least two studied traits, which is also supported by highly significant correlation between the three germination-related traits. QTLs for seed-germination related traits in "QTL hotspot B" were found in both RIL populations and aging conditions, and also QTLs underlying "QTL hotspot A" are identified in both RIL populations under artificial aging condition. These are the stable genomic regions governing the inheritance of seed storability in soybean, and will be the main focus for soybean breeders.

CONCLUSION

This study uncovers the genetic basis of seed storability in soybean. The newly identified QTLs provides valuable information, and will be main targets for fine mapping, candidate gene identification and marker-assisted breeding. Hence, the present study is the first report for the comprehensive and detailed investigation of genetic architecture of seed storability in soybean.

Genetic Diversity of Torch Ginger (Etlingera elatior) Germplasm Revealed by ISSR and SSR Markers

Fifty-seven accessions of torch ginger (Etlingera elatior) collected from seven states in Peninsular Malaysia were evaluated for their molecular characteristics using ISSR and SSR markers to assess the pattern of genetic diversity and association among the characteristics. Diversity study through molecular characterization showed that high variability existed among the 57 torch ginger accessions. ISSR and SSR molecular markers revealed the presence of high genetic variability among the torch ginger accessions. The combination of different molecular markers offered reliable and convincing information about the genetic diversity of torch ginger germplasm. This study found that SSR marker was more informative compared to ISSR marker in determination of gene diversity, polymorphic information content (PIC), and heterozygosity in this population. SSR also revealed high ability in evaluating diversity levels, genetic structure, and relationships of torch ginger due to their codominance and rich allelic diversity. High level of genetic diversity discovered by SSR markers showed the effectiveness of this marker to detect the polymorphism in this germplasm collection.

Development and Application of InDel Markers for Capsicum spp. Based on Whole-Genome Re-Sequencing

Genome-wide identification of Insertion/Deletion polymorphisms (InDels) in Capsicum spp. was performed through comparing whole-genome re-sequencing data from two Capsicum accessions, C. annuum cv. G29 and C. frutescens cv. PBC688, with the reference genome sequence of C. annuum cv. CM334. In total, we identified 1,664,770 InDels between CM334 and PBC688, 533,523 between CM334 and G29, and 1,651,856 between PBC688 and G29. From these InDels, 1605 markers of 3-49 bp in length difference between PBC688 and G29 were selected for experimental validation: 1262 (78.6%) showed polymorphisms, 90 (5.6%) failed to amplify, and 298 (18.6%) were monomorphic. For further validation of these InDels, 288 markers were screened across five accessions representing five domesticated species. Of these assayed markers, 194 (67.4%) were polymorphic, 87 (30.2%) monomorphic and 7 (2.4%) failed. We developed three interspecific InDels, which associated with three genes and showed specific amplification in five domesticated species and clearly differentiated the interspecific hybrids. Thus, our novel PCR-based InDel markers provide high application value in germplasm classification, genetic research and marker-assisted breeding in Capsicum species.

Population genetics and evolutionary history of the wild rice species Oryza rufipogon and O. nivara in Sri Lanka

Genetic diversity and population genetic structure of the wild rice species Oryza rufipogon and O. nivara in Sri Lanka were studied using 33 microsatellite markers. A total of 315 individuals of 11 natural populations collected from the wet, intermediate, and dry zones of the country were used in the study. We found a moderate to high level of genetic diversity at the population level, with the polymorphic loci (P) ranging from 60.6% to 100% (average 81.8%) and the expected heterozygosity (H E) varying from 0.294 to 0.481 (average 0.369). A significant genetic differentiation between species and strong genetic structure within species were also observed. Based on species distribution modeling, we detected the dynamics of the preferred habitats for the two species in Sri Lanka and demonstrated that both O. rufipogon and O. nivara populations have expanded substantially since the last internal glacial. In addition, we showed that the geographical distribution of the two species corresponded to the climate zones and identified a few of key environmental variables that contribute to the distribution of the two species, implying the potential mechanism for ecological adaptation of these two species in Sri Lanka. These studies provided important insights into the population genetics and evolution of these wild species in Sri Lanka and are of great significance to the in situ conservation and utilization of these wild resources in genetic improvement of rice.

Agro-morphological description, genetic diversity and population structure of sugarcane varieties from sub-tropical India

Genetic diversity in 92 sugarcane varieties of sub-tropical India was assessed using 30 morphological descriptors and 643 simple sequence repeat (SSR) marker loci. Out of the 30 morphological descriptors, 14 were found polymorphic, and significant variability was recorded for plant height, cane diameter and number of millable canes. Grouping traits like plant growth habit, leaf blade curvature and leaf sheath adherence were found to be predominantly monomorphic. There were a few pairs of varieties (e.g., CoP 9702 and CoP 9302, CoP 9301 and CoSe 01424, UP 05 and Co 1336, CoS 96258 and CoH 110) that showed similar DUS profiles except differing for a few descriptors. The STRUCTURE profile suggest that all the 92 sugarcane varieties had admixtures and no sub-group had a pure unblemished structure profile. An average Nei's genetic distance of 0.49 was found to be a better measure of diversity, whereas, the average band informativeness (Ib_av) value of all the 80 SSR primers was 0.434. Although, the mean Ib_av values for EST-SSR and genomic-SSR primers were same (0.43), the range of Ib_av of EST-SSR (0.04-0.85) was more compared to genomic-SSR (0.12-0.63) primers. The segregation of the varieties based on morphological traits was not in accordance with their geographical distribution or maturity groups, but principal component analysis was able to group the sugarcane varieties that had similar pedigree together. Results indicate that the SSRs have a potential use in the DNA fingerprinting of varieties to prevent any malpractice like unauthorised re-registration of a previously registered sugarcane variety under PPV&FR Act. The marker profiles could also be utilised for variety identification and release, since at present, it has been made mandatory to include it in addition to the morphological descriptors.

QTL mapping for resistance to and tolerance for the rice root-knot nematode, Meloidogyne graminicola

Background

The root-knot nematode Meloidogyne graminicola is an obligate biotrophic pathogen considered to be the most damaging nematode species that causes significant yield losses to upland and rainfed lowland rice production in South and Southeast Asia. Mapping and identification of quantitative trait loci (QTL) for resistance to and tolerance for M. graminicola may offer a safe and economic management option to farmers. In this study, resistance to and tolerance for M. graminicola in Asian rice (Oryza sativa L.) were studied in a mapping population consisting of 300 recombinant inbred lines (RILs) derived from IR78877–208-B-1-2, an aerobic rice genotype with improved resistance to and tolerance for M. graminicola, and IR64, a popular, high-yielding rice mega-variety susceptible to M. graminicola. RILs were phenotyped for resistance and tolerance in the dry seasons of 2012 and 2013. QTL analysis was performed using 131 single nucleotide polymorphism (SNP) and 33 simple sequence repeat (SSR) markers.

Results

Three QTLs with main effects on chromosomes 4 (qMGR_4.1), 7 (qMGR_7.1) and 9 (qMGR_9.1) and two epistatic interactions (qMGR_3.1/ qMGR_11.1 and qMGR_4.2/ qMGR_8.1) associated with nematode reproduction that were consistent in the two seasons were detected. A QTL affecting root galling was found on chromosomes 4 (qGR_4.1) and 8 (qGR_8.1), and QTLs for nematode tolerance were found on chromosomes 5 (qYR_5.1) and 11 (qYR_11.1). These QTLs were consistent in both seasons. A QTL for grain yield was found on chromosome 10 (qGYLD_10.1), a QTL affecting filled grains per panicle was detected on chromosome 11 (qFG_11.1) and a QTL for fresh root weight was found on chromosomes 2 (qFRWt_2.1), 8 (qFRWt_8.1) and 12 (qFRWt_12.1) in both seasons. The donor of the alleles for qMGR_4.1, qMGR_7.1, qMGR_9.1, qGR_4.1, qGR_8.1, qYR_5.1 and qFRWt_2.1 was IR78877–208-B-1-2, whereas for qYR_11.1, qGYLD_10.1 and qFG_11.1, qFRWt_8.1 and qFRWt_12.1 was IR64. Lines having favorable alleles for resistance, tolerance and yield provided better yield under nematode-infested conditions and could be a starting point of marker-assisted breeding (MAB) for the improvement of M. graminicola resistance and tolerance in Asian rice.

Conclusion

This study identified a total of 12 QTLs with main effects and two epistatic interactions in the 1st season and 2nd season related to M. graminicola resistance and tolerance, and other agronomic traits such as plant yield, percentage of filled grains, and fresh and dry root weight. Rice genotypes that have the favorable alleles for resistance (qMGR_4.1, qMGR_7.1, qMGR_9.1, qGR_4.1, qGR_8.1) and tolerance (qYR_5.1, and qYR_11.1,) QTLs, and which are either resistant or partially resistant and tolerant, were also selected. These selected genotypes and the identified QTLs are vital information in designing MAB for the improvement of high-yielding rice genotypes but are susceptible to M. graminicola infection.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0656-1) contains supplementary material, which is available to authorized users.

Genetic diversity and population structure of Urochloa grass accessions from Tanzania using simple sequence repeat (SSR) markers

Urochloa (syn.-Brachiaria s.s.) is one of the most important tropical forages that transformed livestock industries in Australia and South America. Farmers in Africa are increasingly interested in growing Urochloa to support the burgeoning livestock business, but the lack of cultivars adapted to African environments has been a major challenge. Therefore, this study examines genetic diversity of Tanzanian Urochloa accessions to provide essential information for establishing a Urochloa breeding program in Africa. A total of 36 historical Urochloa accessions initially collected from Tanzania in 1985 were analyzed for genetic variation using 24 SSR markers along with six South American commercial cultivars. These markers detected 407 alleles in the 36 Tanzania accessions and 6 commercial cultivars. Markers were highly informative with an average polymorphic information content of 0.79. The analysis of molecular variance revealed high genetic variation within individual accessions in a species (92%), fixation index of 0.05 and gene flow estimate of 4.77 showed a low genetic differentiation and a high level of gene flow among populations. An unweighted neighbor-joining tree grouped the 36 accessions and six commercial cultivars into three main clusters. The clustering of test accessions did not follow geographical origin. Similarly, population structure analysis grouped the 42 tested genotypes into three major gene pools. The results showed the Urochloa brizantha (A. Rich.) Stapf population has the highest genetic diversity (I = 0.94) with high utility in the Urochloa breeding and conservation program. As the Urochloa accessions analyzed in this study represented only 3 of 31 regions of Tanzania, further collection and characterization of materials from wider geographical areas are necessary to comprehend the whole Urochloa diversity in Tanzania.

Combined QTL mapping, physiological and transcriptomic analyses to identify candidate genes involved in Brassica napus seed aging

Seed aging is an inevitable problem in the germplasm conservation of oil crops. Thus, clarifying the genetic mechanism of seed aging is important for rapeseed breeding. In this study, Brassica napus seeds were exposed to an artificial aging environment (40 °C and 90% relative humidity). Using a population of 172 recombinant inbred lines, 13 QTLs were detected on 8 chromosomes, which explained ~ 9.05% of the total phenotypic variation. The QTLs q2015AGIA-C08 and q2016AGI-C08-2 identified in the two environments were considered the same QTL. After artificial aging, lower germination index, increased relative electrical conductivity, malondialdehyde and proline content, and reduced soluble sugar, protein content and antioxidant enzyme activities were detected. Furthermore, seeds of extreme lines that were either left untreated (R0 and S0) or subjected to 15 days of artificial aging (R15 and S15) were used for transcriptome sequencing. In total, 2843, 1084, 429 and 1055 differentially expressed genes were identified in R15 vs. R0, S15 vs. S0, R0 vs. S0 and R15 vs. S15, respectively. Through integrated QTL mapping and RNA-sequencing analyses, seven genes, such as BnaA03g37460D, encoding heat shock transcription factor C1, and BnaA03g40360D, encoding phosphofructokinase 4, were screened as candidate genes involved in seed aging. Further researches on these candidate genes could broaden our understanding of the regulatory mechanisms of seed aging.

Transcriptomic signature of drought response in pearl millet (Pennisetum glaucum (L.) and development of web-genomic resources

Pearl millet, (Pennisetum glaucum L.), an efficient (C4) crop of arid/semi-arid regions is known for hardiness. Crop is valuable for bio-fortification combating malnutrition and diabetes, higher caloric value and wider climatic resilience. Limited studies are done in pot-based experiments for drought response at gene-expression level, but field-based experiment mimicking drought by withdrawal of irrigation is still warranted. We report de novo assembly-based transcriptomic signature of drought response induced by irrigation withdrawal in pearl millet. We found 19983 differentially expressed genes, 7595 transcription factors, gene regulatory network having 45 hub genes controlling drought response. We report 34652 putative markers (4192 simple sequence repeats, 12111 SNPs and 6249 InDels). Study reveals role of purine and tryptophan metabolism in ABA accumulation mediating abiotic response in which MAPK acts as major intracellular signal sensing drought. Results were validated by qPCR of 13 randomly selected genes. We report the first web-based genomic resource ( http://webtom.cabgrid.res.in/pmdtdb/ ) which can be used for candidate genes-based SNP discovery programs and trait-based association studies. Looking at climatic change, nutritional and pharmaceutical importance of this crop, present investigation has immense value in understanding drought response in field condition. This is important in germplasm management and improvement in endeavour of pearl millet productivity.