Construction of Genetic Map and QTL Mapping for Seed Size and Quality Traits in Soybean (Glycine max L.)

Soybean (Glycine max L.) is the main source of vegetable protein and edible oil for humans, with an average content of about 40% crude protein and 20% crude fat. Soybean yield and quality traits are mostly quantitative traits controlled by multiple genes. The quantitative trait loci (QTL) mapping for yield and quality traits, as well as for the identification of mining-related candidate genes, is of great significance for the molecular breeding and understanding the genetic mechanism. In this study, 186 individual plants of the F2 generation derived from crosses between Changjiangchun 2 and Yushuxian 2 were selected as the mapping population to construct a molecular genetic linkage map. A genetic map containing 445 SSR markers with an average distance of 5.3 cM and a total length of 2375.6 cM was obtained. Based on constructed genetic map, 11 traits including hundred-seed weight (HSW), seed length (SL), seed width (SW), seed length-to-width ratio (SLW), oil content (OIL), protein content (PRO), oleic acid (OA), linoleic acid (LA), linolenic acid (LNA), palmitic acid (PA), stearic acid (SA) of yield and quality were detected by the multiple- d size traits and 113 QTLs related to quality were detected by the multiple QTL model (MQM) mapping method across generations F2, F2:3, F2:4, and F2:5. A total of 71 QTLs related to seed size traits and 113 QTLs related to quality traits were obtained in four generations. With those QTLs, 19 clusters for seed size traits and 20 QTL clusters for quality traits were summarized. Two promising clusters, one related to seed size traits and the other to quality traits, have been identified. The cluster associated with seed size traits spans from position 27876712 to 29009783 on Chromosome 16, while the cluster linked to quality traits spans from position 12575403 to 13875138 on Chromosome 6. Within these intervals, a reference genome of William82 was used for gene searching. A total of 36 candidate genes that may be involved in the regulation of soybean seed size and quality were screened by gene functional annotation and GO enrichment analysis. The results will lay the theoretical and technical foundation for molecularly assisted breeding in soybean.

Construction of a genetic map and QTL mapping of seed size traits in soybean

Soybean seed size and seed shape traits are closely related to plant yield and appearance quality. In this study, 186 individual plants of the F2 generation derived from crosses between Changjiang Chun 2 and JiYu 166 were selected as the mapping population to construct a molecular genetic linkage map, and the phenotypic data of hundred-grain weight, seed length, seed width, and seed length-to-width ratio of soybean under three generations of F2 single plants and F2:3 and F2:4 lines were combined to detect the QTL (quantitative trait loci) for the corresponding traits by ICIM mapping. A soybean genetic map containing 455 markers with an average distance of 6.15 cM and a total length of 2799.2 cM was obtained. Forty-nine QTLs related to the hundred-grain weight, seed length, seed width, and seed length-to-width ratio of soybean were obtained under three environmental conditions. A total of 10 QTLs were detected in more than two environments with a phenotypic variation of over 10%. Twelve QTL clusters were identified on chromosomes 1, 2, 5, 6, 8, 13, 18, and 19, with the majority of the overlapping intervals for hundred-grain weight and seed width. These results will lay the theoretical and technical foundation for molecularly assisted breeding in soybean seed weight and seed shape. Eighteen candidate genes that may be involved in the regulation of soybean seed size were screened by gene functional annotation and GO enrichment analysis.

A QTL approach in faba bean highlights the conservation of genetic control of frost tolerance among legume species

Frost is a major abiotic stress of winter type faba beans (Vica faba L.) and has adverse effects on crop yield. Climate change, far from reducing the incidence of frost events, is making these phenomena more and more common, severe, and prolonged. Despite the important interaction that the environment has in the tolerance of faba bean to frost, this trait seems to have good levels of heritability. Several QTLs for frost tolerance have already been reported, however, a more robust identification is needed to more precisely identify the genomic regions involved in faba bean tolerance to sub-zero temperatures. Several pea (Pisum sativum L.) and barrel medic (Medicago truncatula L.) frost tolerance QTLs appear to be conserved between these two species, furthering the hypothesis that the genetic control of frost tolerance in legume species might be more generally conserved. In this work, the QTL mapping in two faba bean recombinant inbred line (RIL) populations connected by a common winter-type parent has led to the identification of five genomic regions involved in the control of frost tolerance on linkage groups I, III, IV, and V. Among them, a major and robust QTL of great interest for marker-assisted selection was identified on the lower part of the long-arm of LGI. The synteny between the faba bean frost tolerance QTLs and those previously identified in other legume species such as barrel medic, pea or soybean highlighted at least partial conservation of the genetic control of frost tolerance among different faba bean genetic pools and legume species. Four novel RILs showing high and stable levels of tolerance and the ability to recover from freezing temperatures by accumulating frost tolerance QTLs are now available for breeding programs.

Construction of a high-density genetic map and localization of grazing-tolerant QTLs in Medicago falcata L

BACKGROUND

Using genomic DNA from 79 F1 plants resulted from a crossing between parents with strong and weak grazing tolerance in Medicago falcata L., we generated an EcoRI restriction site-associated DNA (RAD) sequencing library. After sequencing and assembly, a high-density genetic map with high-quality SNP markers was constructed, with a total length of 1312.238 cM and an average density of 0.844 SNP/cM.

METHODS

The phenotypic traits of 79 F1 families were observed and the QTLS of 6 traits were analyzed by interval mapping.

RESULTS

Sixty three QTLs were identified for seven traits with LOD values from 3 to 6 and the contribution rates from 15% to 30%. Among the 63 QTLs, 17 were for natural shoot height, 12 for rhizome Length, 10 for Shoot canopy diameter, 9 for Basal plant diameter, 6 for stem number, 5 for absolute shoot height, and 4 for rhizome width. These QTLs were concentrated on LG2, LG4, LG5, LG7, and LG8. LG6 had only 6 QTLs. According to the results of QTL mapping, comparison of reference genomes, and functional annotation, 10 candidate genes that may be related to grazing tolerance were screened. qRT-PCR analysis showed that two candidate genes (LOC11412291 and LOC11440209) may be the key genes related to grazing tolerance of M. falcata.

CONCLUSION

The identified trait-associated QTLs and candidate genes in this study will provide a solid foundation for future molecular breeding for enhanced grazing-tolerance in M. falcata.

Genetic Linkage and Physical Mapping for an Oyster Mushroom (Pleurotus cornucopiae) and Quantitative Trait Locus Analysis for Cap Color

Oyster mushrooms are grown commercially worldwide, especially in many developing countries, for their easy cultivation and high biological efficiency. Pleurotus cornucopiae is one of the main oyster mushroom species because of its gastronomic value and nutraceutical properties. Cap color is an important trait, since consumers prefer dark mushrooms, which are now represented by only a small portion of the commercial varieties. Breeding efforts are required to improve quality-related traits to satisfy various demands of consumers. Here, we present a saturated genetic linkage map of P. cornucopiae constructed by using a segregating population of 122 monokaryons and 3,449 single nucleotide polymorphism (SNP) markers generated by the 2b-RAD approach. The map contains 11 linkage groups covering 961.6 centimorgans (cM), with an average marker spacing of 0.27 cM. The genome of P. cornucopiae was de novo sequenced, resulting in 425 scaffolds (>1,000 bp) with a total genome size of 35.1 Mb. The scaffolds were assembled to the pseudochromosome level with the assistance of the genetic linkage map. A total of 97% SNP markers (3,357) were physically localized on 140 scaffolds that were assigned to 11 pseudochromosomes, with a total of 32.5 Mb, representing 92.5% of the whole genome. Six quantitative trait loci (QTL) controlling cap color of P. cornucopiae were detected, accounting for a total phenotypic variation of 65.6%, with the highest value for the QTL on pseudochromosome 5 (18%). The results of our study provide a solid base for marker-assisted breeding for agronomic traits and especially for studies on biological mechanisms controlling cap color in oyster mushrooms. IMPORTANCE Oyster mushrooms are produced and consumed all over the world. Pleurotus cornucopiae is one of the main oyster mushroom species. Dark-cap oyster mushrooms are becoming more and more popular with consumers, but dark varieties are rare on the market. Prerequisites for efficient breeding programs are the availability of high-quality whole genomes and genetic linkage maps. Genetic studies to fulfill some of these prerequisites have hardly been done for P. cornucopiae. In this study, we de novo sequenced the genome and constructed a saturated genetic linkage map for P. cornucopiae. The genetic linkage map was effectively used to assist the genome assembly and identify QTL that genetically control the trait cap color. As well, the genome characteristics of P. cornucopiae were compared to the closely related species Pleurotus ostreatus. The results provided a basis for understanding the genetic background and marker-assisted breeding of this economically important mushroom species.

High-Density Genetic Linkage Mapping of Lepidium Based on Genotyping-by-Sequencing SNPs and Segregating Contig Tag Haplotypes

Lepidium campestre has been targeted for domestication as future oilseed and catch crop. Three hundred eighty plants comprising genotypes of L. campestre, Lepidium heterophyllum, and their interspecific F2 mapping population were genotyped using genotyping by sequencing (GBS), and the generated polymorphic markers were used for the construction of high-density genetic linkage map. TASSEL-GBS, a reference genome-based pipeline, was used for this analysis using a draft L. campestre whole genome sequence. The analysis resulted in 120,438 biallelic single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) above 0.01. The construction of genetic linkage map was conducted using MSTMap based on phased SNPs segregating in 1:2:1 ratio for the F2 individuals, followed by genetic mapping of segregating contig tag haplotypes as dominant markers against the linkage map. The final linkage map consisted of eight linkage groups (LGs) containing 2,330 SNP markers and spanned 881 Kosambi cM. Contigs (10,302) were genetically mapped to the eight LGs, which were assembled into pseudomolecules that covered a total of ∼120.6 Mbp. The final size of the pseudomolecules ranged from 9.4 Mbp (LG-4) to 20.4 Mpb (LG-7). The following major correspondence between the eight Lepidium LGs (LG-1 to LG-8) and the five Arabidopsis thaliana (At) chromosomes (Atx-1-Atx-5) was revealed through comparative genomics analysis: LG-1&2_Atx-1, LG-3_Atx-2&3, LG-4_Atx-2, LG-5_Atx-2&Atx-3, LG-6_Atx-4&5, LG-7_Atx-4, and LG-8_Atx-5. This analysis revealed that at least 66% of the sequences of the LGs showed high collinearity with At chromosomes. The sequence identity between the corresponding regions of the LGs and At chromosomes ranged from 80.6% (LG-6) to 86.4% (LG-8) with overall mean of 82.9%. The map positions on Lepidium LGs of the homologs of 24 genes that regulate various traits in A. thaliana were also identified. The eight LGs revealed in this study confirm the previously reported (1) haploid chromosome number of eight in L. campestre and L. heterophyllum and (2) chromosomal fusion, translocation, and inversion events during the evolution of n = 8 karyotype in ancestral species shared by Lepidium and Arabidopsis to n = 5 karyotype in A. thaliana. This study generated highly useful genomic tools and resources for Lepidium that can be used to accelerate its domestication.

Genetic mapping in Diversity Outbred mice identifies a Trpa1 variant influencing late-phase formalin response

Identification of genetic variants that influence susceptibility to pain is key to identifying molecular mechanisms and targets for effective and safe therapeutic alternatives to opioids. To identify genes and variants associated with persistent pain, we measured late-phase response to formalin injection in 275 male and female Diversity Outbred mice genotyped for over 70,000 single nucleotide polymorphisms. One quantitative trait locus reached genome-wide significance on chromosome 1 with a support interval of 3.1 Mb. This locus, Nociq4 (nociceptive sensitivity quantitative trait locus 4; MGI: 5661503), harbors the well-known pain gene Trpa1 (transient receptor potential cation channel, subfamily A, member 1). Trpa1 is a cation channel known to play an important role in acute and chronic pain in both humans and mice. Analysis of Diversity Outbred founder strain allele effects revealed a significant effect of the CAST/EiJ allele at Trpa1, with CAST/EiJ carrier mice showing an early, but not late, response to formalin relative to carriers of the 7 other inbred founder alleles (A/J, C57BL/6J, 129S1/SvImJ, NOD/ShiLtJ, NZO/HlLtJ, PWK/PhJ, and WSB/EiJ). We characterized possible functional consequences of sequence variants in Trpa1 by assessing channel conductance, TRPA1-TRPV1 interactions, and isoform expression. The phenotypic differences observed in CAST/EiJ relative to C57BL/6J carriers were best explained by Trpa1 isoform expression differences, implicating a splice junction variant as the causal functional variant. This study demonstrates the utility of advanced, high-precision genetic mapping populations in resolving specific molecular mechanisms of variation in pain sensitivity.

Introgression of Aegilops speltoides segments in Triticum aestivum and the effect of the gametocidal genes

BACKGROUND AND AIMS

Bread wheat (Triticum aestivum) has been through a severe genetic bottleneck as a result of its evolution and domestication. It is therefore essential that new sources of genetic variation are generated and utilized. This study aimed to generate genome-wide introgressed segments from Aegilops speltoides. Introgressions generated from this research will be made available for phenotypic analysis.

METHODS

Aegilops speltoides was crossed as the male parent to T. aestivum 'Paragon'. The interspecific hybrids were then backcrossed to Paragon. Introgressions were detected and characterized using the Affymetrix Axiom Array and genomic in situ hybridization (GISH).

KEY RESULTS

Recombination in the gametes of the F1 hybrids was at a level where it was possible to generate a genetic linkage map of Ae. speltoides. This was used to identify 294 wheat/Ae. speltoides introgressions. Introgressions from all seven linkage groups of Ae. speltoides were found, including both large and small segments. Comparative analysis showed that overall macro-synteny is conserved between Ae. speltoides and T. aestivum, but that Ae. speltoides does not contain the 4A/5A/7B translocations present in wheat. Aegilops speltoides has been reported to carry gametocidal genes, i.e. genes that ensure their transmission through the gametes to the next generation. Transmission rates of the seven Ae. speltoides linkage groups introgressed into wheat varied. A 100 % transmission rate of linkage group 2 demonstrates the presence of the gametocidal genes on this chromosome.

CONCLUSIONS

A high level of recombination occurs between the chromosomes of wheat and Ae. speltoides, leading to the generation of large numbers of introgressions with the potential for exploitation in breeding programmes. Due to the gametocidal genes, all germplasm developed will always contain a segment from Ae. speltoides linkage group 2S, in addition to an introgression from any other linkage group.

Multigenic Control of Pod Shattering Resistance in Chinese Rapeseed Germplasm Revealed by Genome-Wide Association and Linkage Analyses

The majority of rapeseed cultivars shatter seeds upon maturity especially under hot-dry and windy conditions, reducing yield and gross margin return to growers. Here, we identified quantitative trait loci (QTL) for resistance to pod shatter in an unstructured diverse panel of 143 rapeseed accessions, and two structured populations derived from bi-parental doubled haploid (DH) and inter-mated (IF2) crosses derived from R1 (resistant to pod shattering) and R2 (prone to pod shattering) accessions. Genome-wide association analysis identified six significant QTL for resistance to pod shatter located on chromosomes A01, A06, A07, A09, C02, and C05. Two of the QTL, qSRI.A09 delimited with the SNP marker Bn-A09-p30171993 (A09) and qSRI.A06 delimited with the SNP marker Bn-A06-p115948 (A06) could be repeatedly detected across environments in a diversity panel, DH and IF2 populations, suggesting that at least two loci on chromosomes A06 and A09 were the main contributors to pod shatter resistance in Chinese germplasm. Significant SNP markers identified in this study especially those that appeared repeatedly across environments provide a cost-effective and an efficient method for introgression and pyramiding of favorable alleles for pod shatter resistance via marker-assisted selection in rapeseed improvement programs.

Combining powers of linkage and association mapping for precise dissection of QTL controlling resistance to gray leaf spot disease in maize (Zea mays L.)

BACKGROUND

Gray Leaf Spot (GLS causal agents Cercospora zeae-maydis and Cercospora zeina) is one of the most important foliar diseases of maize in all areas where the crop is being cultivated. Although in the USA the situation with GLS severity is not as critical as in sub-Saharan Africa or Brazil, the evidence of climate change, increasing corn monoculture as well as the narrow genetic base of North American resistant germplasm can turn the disease into a serious threat to US corn production. The development of GLS resistant cultivars is one way to control the disease. In this study we combined the high QTL detection power of genetic linkage mapping with the high resolution power of genome-wide association study (GWAS) to precisely dissect QTL controlling GLS resistance and identify closely linked molecular markers for robust marker-assisted selection and trait introgression.

RESULTS

Using genetic linkage analysis with a small bi-parental mapping population, we identified four GLS resistance QTL on chromosomes 1, 6, 7, and 8, which were validated by GWAS. GWAS enabled us to dramatically increase the resolution within the confidence intervals of the above-mentioned QTL. Particularly, GWAS revealed that QTLGLSchr8, detected by genetic linkage mapping as a locus with major effect, was likely represented by two QTL with smaller effects. Conducted in parallel, GWAS of days-to-silking demonstrated the co-localization of flowering time QTL with GLS resistance QTL on chromosome 7 indicating that either QTLGLSchr7 is a flowering time QTL or it is a GLS resistance QTL that co-segregates with the latter. As a result, this genetic linkage - GWAS hybrid mapping system enabled us to identify one novel GLS resistance QTL (QTLGLSchr8a) and confirm with more refined positions four more previously mapped QTL (QTLGLSchr1, QTLGLSchr6, QTLGLSchr7, and QTLGLSchr8b). Through the novel Single Donor vs. Elite Panel method we were able to identify within QTL confidence intervals SNP markers that would be suitable for marker-assisted selection of gray leaf spot resistant genotypes containing the above-mentioned GLS resistance QTL.

CONCLUSION

The application of a genetic linkage - GWAS hybrid mapping system enabled us to dramatically increase the resolution within the confidence interval of GLS resistance QTL by-passing labor- and time-intensive fine mapping. This method appears to have a great potential to accelerate the pace of QTL mapping projects. It is universal and can be used in the QTL mapping projects in any crops.

SNP markers-based map construction and genome-wide linkage analysis in Brassica napus

An Illumina Infinium array comprising 5306 single nucleotide polymorphism (SNP) markers was used to genotype 175 individuals of a doubled haploid population derived from a cross between Skipton and Ag-Spectrum, two Australian cultivars of rapeseed (Brassica napus L.). A genetic linkage map based on 613 SNP and 228 non-SNP (DArT, SSR, SRAP and candidate gene markers) covering 2514.8 cM was constructed and further utilized to identify loci associated with flowering time and resistance to blackleg, a disease caused by the fungus Leptosphaeria maculans. Comparison between genetic map positions of SNP markers and the sequenced Brassica rapa (A) and Brassica oleracea (C) genome scaffolds showed several genomic rearrangements in the B. napus genome. A major locus controlling resistance to L. maculans was identified at both seedling and adult plant stages on chromosome A07. QTL analyses revealed that up to 40.2% of genetic variation for flowering time was accounted for by loci having quantitative effects. Comparative mapping showed Arabidopsis and Brassica flowering genes such as Phytochrome A/D, Flowering Locus C and agamous-Like MADS box gene AGL1 map within marker intervals associated with flowering time in a DH population from Skipton/Ag-Spectrum. Genomic regions associated with flowering time and resistance to L. maculans had several SNP markers mapped within 10 cM. Our results suggest that SNP markers will be suitable for various applications such as trait introgression, comparative mapping and high-resolution mapping of loci in B. napus.

Construction of integrated linkage map of a recombinant inbred line population of white lupin (Lupinus albus L.)

We report the development of a Diversity Arrays Technology (DArT) marker panel and its utilisation in the development of an integrated genetic linkage map of white lupin (Lupinus albus L.) using an F8 recombinant inbred line population derived from Kiev Mutant/P27174. One hundred and thirty-six DArT markers were merged into the first genetic linkage map composed of 220 amplified fragment length polymorphisms (AFLPs) and 105 genic markers. The integrated map consists of 38 linkage groups of 441 markers and spans a total length of 2,169 cM, with an average interval size of 4.6 cM. The DArT markers exhibited good genome coverage and were associated with previously identified genic and AFLP markers linked with quantitative trait loci for anthracnose resistance, flowering time and alkaloid content. The improved genetic linkage map of white lupin will aid in the identification of markers for traits of interest and future syntenic studies.

Molecular Mapping of D₁, D₂ and ms5 Revealed Linkage between the Cotyledon Color Locus D₂ and the Male-Sterile Locus ms5 in Soybean

In soybean, genic male sterility can be utilized as a tool to develop hybrid seed. Several male-sterile, female-fertile mutants have been identified in soybean. The male-sterile, female-fertile ms5 mutant was selected after fast neutron irradiation. Male-sterility due to ms5 was associated with the “stay-green” cotyledon color mutation. The cotyledon color trait in soybean is controlled by two loci, D₁ and D₂. Association between cotyledon color and male-sterility can be instrumental in early phenotypic selection of sterility for hybrid seed production. The use of such selection methods saves time, money, and space, as fewer seeds need to be planted and screened for sterility. The objectives of this study were to compare anther development between male-fertile and male-sterile plants, to investigate the possible linkages among the Ms5, D₁ and D₂ loci, and to determine if any of the d₁ or d₂ mutations can be applied in hybrid seed production. The cytological analysis during anther development displayed optically clear, disintegrating microspores and enlarged, engorged pollen in the male-sterile, female-fertile ms5ms5 plants, a common characteristic of male-sterile mutants. The D₁ locus was mapped to molecular linkage group (MLG) D1a and was flanked by Satt408 and BARCSOYSSR_01_1622. The ms5 and D₂ loci were mapped to MLG B1 with a genetic distance ~12.8 cM between them. These results suggest that use of the d₂ mutant in the selection of male-sterile line may attenuate the cost hybrid seed production in soybean.