Chromosomal regions associated with segregation distortion of molecular markers in F2, backcross, doubled haploid, and recombinant inbred populations in rice (Oryza sativa L.)

Mapping of quantitative trait loci and mining of candidate genes for seed viability in soybean [Glycine max (L.) Merr.]

Global oilseed crop soybean [Glycine max (L.) Merrill] contains 18%-20% oil, 40%-45% protein, and countless nutrients vital for human health. It is grown worldwide for food, feed, pharmaceutical, and industrial applications. However, inherent loss of seed viability during ambient storage poses serious bottleneck in the production and maintenance of quality seeds. Understanding inheritance and mapping of quantitative trait loci (QTLs) for seed viability would help in designing breeding program for developing varieties with higher viability of the seeds. In this study, attempt was made to map QTLs and identify candidate genes for seed viability in soybean. A high-viable genotype EC1023 (>90% germination after 1 year of storage) was hybridized with VLS61, a poor storing genotype (<70% germination after 1 year of storage), and the F1 seeds were advanced to the next generation. The F2:3 seeds were subjected to accelerated ageing (AA) by exposing it to 41°C at 100% RH for 72h followed by viability testing through germination test. After AA test, the germination of the parental genotypes EC1023 and VLS61 were 40% and 14%, respectively, and that of the F2:3 seeds ranged from 4.16% to 71.42% indicating wide variability in the viability of the seeds. Genetic polymorphism studied with 517 SSR markers indicated the polymorphism between the parental genotypes to be 20.35%; however, distribution of the polymorphism was not uniform across the chromosomes; Chr. 14 had 30.00% polymorphism as against 7.14% on Chrs.12. Through inclusive composite interval mapping approach, 8 QTL for seed viability, namely, qSv-6.1 and qSv-6.2, qSv-7.1, qSv-8.1, and qSv-8.2, qSv-10.1, qSv-13.1, and qSv-17.1 were mapped on Chrs. 6, 7, 8, 10, 13 and 17, respectively. The phenotypic variation explained (PVE) by the QTL were 1.97%-11.10%. Two QTL, namely, qSv-7.1 (PVE = 11.10%) and qSv13.1 (PVE = 11.08%) appeared to be major QTLs for seed viability and rest minor ones. All QTL except qSv8.2 appeared to be novel. The mapped QTLs were validated in 40 inter-specific RILs with varying level of seed viability. The SSR marker Satt538 linked to the QTL qSv8.2 could successfully (70%) separate the highly viable RILs from the poor-viable RILs. Similarly, SSR markers Sat_316 and Sat_173 were 80%-85% successful in separating the high and poor viable RILs. Based on Protein Analysis Through Evolutionary Relationships (PANTHER), gene annotation information, and literature search, more than 500 candidate genes for seed viability underlying the mapped QTL were identified. The mapped QTL and the identified candidate genes will pave the way for marker-assisted breeding of soybean to generate genotypes with improved seed viability.

Introgression and Mapping of a Novel Bacterial Blight Resistance Gene xa49(t) from Oryza rufipogon acc. CR100098A into O. sativa

Bacterial blight (BB) caused by Xanthomonas oryzae pv. oryzae is one of the epidemic diseases in rice. Rapid changes in the pathogenicity of the X. oryzae pv. oryzae pathogen demand the identification and characterization of novel BB resistance genes. Here, we report the transfer and mapping of a new BB resistance gene from Oryza rufipogon acc. CR100098A. Inheritance studies on the BC2F2 population, BC2F3 progenies, and backcross-derived recombinant inbred lines derived from a cross between Pusa44/O. rufipogon acc. CR100098A//2*PR114 showed that a single recessive gene confers resistance in O. rufipogon acc. CR100098A. Bulked segregant analysis using 203 simple sequence repeat (SSR) markers localized the BB resistance gene on chromosome 11 bracketed between two SSR markers, RM27235 and RM2136. Using PR114 and O. rufipogon acc. CR100098A genotyping by sequencing data, 86 KASP markers within the bracketed region were designed and tested for bulked segregant analysis. Only five KASP markers showed polymorphism between parents, and three were associated with the target gene. Seventy-seven new SSR markers were designed from the same interval. A total of 33 polymorphic markers were analyzed on the whole population and mapped the BB gene in an interval of 2.8 cM flanked by SSR markers PAU11_65 and PAU11_44 within a physical distance of 376.3 kb. The BB resistance gene mapped in this study is putatively new and designated as xa49(t). Fourteen putative candidate genes were identified within the xa49(t) region having a role in biotic stress resistance. The linked markers to the xa49(t) gene were validated in other rice cultivars for its successful deployment in BB resistance breeding.

Selfing Promotes Spread and Introgression of Segregation Distorters in Hermaphroditic Plants

Segregation distorters (SDs) are genetic elements that distort the Mendelian segregation ratio to favor their own transmission and are able to spread even when they incur fitness costs on organisms carrying them. Depending on the biology of the host organisms and the genetic architecture of the SDs, the population dynamics of SDs can be highly variable. Inbreeding is considered an effective mechanism for inhibiting the spread of SDs in populations, and can evolve as a defense mechanism against SDs in some systems. However, we show that inbreeding in the form of selfing in fact promotes the spread of SDs acting as pollen killers in a toxin-antidote system in hermaphroditic plants by two mechanisms: (i) By reducing the effective recombination rate between killer and antidote loci in the two-locus system and (ii) by increasing the proportion of SD alleles in individual flowers, rather than in the general gene-pool. We also show that in rice (Oryza sativa L.), a typical hermaphroditic plant, all molecularly characterized SDs associated with pollen killing were involved in population hybridization and have introgressed across different species. Paradoxically, these loci, which are associated with hybrid incompatibility and can be thought of as Bateson-Dobzhansky-Muller incompatibility loci are expected to reduce gene-flow between species, in fact cross species boundaries more frequently than random loci, and may act as important drivers of introgression.

Genetic linkage map construction and QTL analysis for plant height in proso millet (Panicum miliaceum L.)

KEY MESSAGE

A genetic linkage map representing proso millet genome was constructed with SSR markers, and a major QTL corresponding to plant height was mapped on chromosome 14 of this map. Proso millet (Panicum miliaceum L.) has the lowest water requirements of all cultivated cereal crops. However, the lack of a genetic map and the paucity of genomic resources for this species have limited the utility of proso millet for detailed genetic studies and hampered genetic improvement programs. In this study, 97,317 simple sequence repeat (SSR) markers were developed based on the genome sequence of the proso millet landrace Longmi 4. Using some of these markers in conjunction with previously identified SSRs, an SSR-based linkage map for proso millet was successfully constructed using a large mapping population (316 F2 offspring). In total, 186 SSR markers were assigned to 18 linkage groups corresponding to the haploid chromosomes. The constructed map had a total length of 3033.42 centimorgan (cM) covering 78.17% of the assembled reference genome. The length of the 18 linkage groups ranged from 88.89 cM (Chr. 15) to 274.82 cM (Chr. 16), with an average size of 168.17 cM. To our knowledge, this is the first genetic linkage map for proso millet based on SSR markers. Plant height is one of the most important traits in crop improvement. A major QTL was repeatedly detected in different environments, explaining 8.70-24.50% of the plant height variations. A candidate gene affecting auxin biosynthesis and transport, and ROS homeostasis regulation was predicted. Thus, the linkage map and QTL analysis provided herein will promote the development of gene mining and molecular breeding in proso millet.

Unveiling the Genetic Basis Underlying Rice Anther Culturability via Segregation Distortion Analysis in Doubled Haploid Population

Anther culture (AC) is a valuable technique in rice breeding. However, the genetic mechanisms underlying anther culturability remain elusive, which has hindered its widespread adoption in rice breeding programs. During AC, microspores carrying favorable alleles for AC are selectively regenerated, leading to segregation distortion (SD) of chromosomal regions linked to these alleles in the doubled haploid (DH) population. Using the AC method, a DH population was generated from the japonica hybrid rice Shenyou 26. A genetic map consisting of 470 SNPs was constructed using this DH population, and SD analysis was performed at both the single- and two-locus levels to dissect the genetic basis underlying anther culturability. Five segregation distortion loci (SDLs) potentially linked to anther culturability were identified. Among these, SDL5 exhibited an overrepresentation of alleles from the female parent, while SDL1.1, SDL1.2, SDL2, and SDL7 displayed an overrepresentation of alleles from the male parent. Furthermore, six pairs of epistatic interactions (EPIs) that influenced two-locus SDs in the DH population were discovered. A cluster of genetic loci, associated with EPI-1, EPI-3, EPI-4, and EPI-5, overlapped with SDL1.1, indicating that the SDL1.1 locus may play a role in regulating anther culturability via both additive and epistatic mechanisms. These findings provide valuable insights into the genetic control of anther culturability in rice and lay the foundation for future research focused on identifying the causal genes associated with anther culturability.

Identification of structural variations related to drought tolerance in wheat (Triticum aestivum L.)

Structural variations are common in plant genomes, affecting meiotic recombination and distorted segregation in wheat. And presence/absence variations can significantly affect drought tolerance in wheat. Drought is a major abiotic stress limiting wheat production. Common wheat has a complex genome with three sub-genomes, which host large numbers of structural variations (SVs). SVs play critical roles in understanding the genetic contributions of plant domestication and phenotypic plasticity, but little is known about their genomic characteristics and their effects on drought tolerance. In the present study, high-resolution karyotypes of 180 doubled haploids (DHs) were developed. Signal polymorphisms between the parents involved with 8 presence-absence variations (PAVs) of tandem repeats (TR) distributed on the 7 (2A, 4A, 5A, 7A, 3B, 7B, and 2D) of 21 chromosomes. Among them, PAV on chromosome 2D showed distorted segregation, others transmit normal conforming to a 1:1 segregation ration in the population; and a PAVs recombination occurred on chromosome 2A. Association analysis of PAV and phenotypic traits under different water regimes, we found PAVs on chromosomes 4A, 5A, and 7B showed negative effect on grain length (GL) and grain width (GW); PAV.7A had opposite effect on grain thickness (GT) and spike length (SL), with the effect on traits differing under different water regimes. PAVs on linkage group 2A, 4A, 7A, 2D, and 7B associated with the drought tolerance coefficients (DTCs), and significant negative effect on drought resistance values (D values) were detected in PAV.7B. Additionally, quantitative trait loci (QTL) associated with phenotypic traits using the 90 K SNP array showed QTL for DTCs and grain-related traits in chromosomes 4A, and 5A, 3B were co-localized in differential regions of PAVs. These PAVs can cause the differentiation of the target region of SNP and could be used for genetic improvement of agronomic traits under drought stress via marker-assisted selection (MAS) breeding.

Systematic trait dissection in oilseed rape provides a comprehensive view, further insight, and exact roadmap for yield determination

Background

Yield is the most important and complex trait that is influenced by numerous relevant traits with very complicated interrelations. While there are a large number of studies on the phenotypic relationship and genetic basis of yield traits, systematic studies with further dissection focusing on yield are limited. Therefore, there is still lack of a comprehensive and in-depth understanding of the determination of yield.

Results

In this study, yield was systematically dissected at the phenotypic, genetic to molecular levels in oilseed rape (Brassica napus L.). The analysis of correlation, network, and principal component for 21 traits in BnaZN-RIL population showed that yield was determined by a complex trait network with key contributors. The analysis of the constructed high-density single nucleotide polymorphism (SNP) linkage map revealed the concentrated distribution of distorted and heterozygous markers, likely due to selection on genes controlling the growth period and yield heterosis. A total of 134 consensus quantitative trait loci (QTL) were identified for 21 traits, of which all were incorporated into an interconnecting QTL network with dozens of hub-QTL. Four representative hub-QTL were further dissected to the target or candidate genes that governed the causal relationships between the relevant traits.

Conclusions

The highly consistent results at the phenotypic, genetic, and molecular dissecting demonstrated that yield was determined by a multilayer composite network that involved numerous traits and genes showing complex up/down-stream and positive/negative regulation. This provides a systematic view, further insight, and exact roadmap for yield determination, which represents a significant advance toward the understanding and dissection of complex traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02134-w.

Identification of QTLs for zinc deficiency tolerance in a recombinant inbred population of rice (Oryza sativa L.)

BACKGROUND

Deficiency of Zn is a major soil constraint in rice plant growth and yield. Edaphic factors such as Zn deficiency in soil in relation to plant performance are still poorly understood. Here, we report promising quantitative trait loci (QTL) conferring tolerance to Zn deficiency, which were identified through biparental mapping. The experiment was conducted using the 236 F₇ recombinant inbred line mapping population derived from the cross of Kinandang Patong (Zn deficiency sensitive) and A69-1 (Zn deficiency tolerant).

RESULTS

A total of six QTLs (qLB-2B, qLB-4B, qPM-4B, qPM-6B, qRZC-4B, qSZC-4B) on chromosomes 2, 4 and 6 were identified for environment 1, whereas five QTLs (qLB-2 N, qLB-4 N, qPM-4 N, qRZC-4 N, qSZC-4 N) on chromosomes 2 and 4 were detected for environment 2. Among these, five major (51.30, 48.70, 28.60, 56.00, 52.00 > 10 R² ) and one minor (5.40 < 10 R² ) QTLs for environment 1 and four major (51.48, 50.20, 53.00, 48.00 > 10 R² ) and one minor (4.44 < 10) QTLs for environment 2 for Zn deficiency tolerance with a logarithm of odd threshold value higher than 3 were identified. The QTLs (qLB-4B, qPM-4B, qRZC-4B, qSZC-4B, qLB-4 N, qPM-4 N, qRZC-4 N, qSZC-4 N) for leaf bronzing, plant mortality root zinc concentration and shoot zinc concentration identified on chromosome 4 were found to be the most promising and highly reproducible across the locations that explained phenotypic variation from 48.00% to 56.00% with the same marker interval RM6748-RM303.

CONCLUSION

The new QTLs and its linked markers identified in the present study can be utilized for Zn deficiency tolerance in elite cultivars using marker-assisted backcrossing. © 2022 Society of Chemical Industry.

Molecular mapping of QTLs for grain dimension traits in Basmati rice

Basmati rice is known for its extra-long slender grains, exceptional kernel dimensions after cooking, high volume expansion, and strong aroma. Developing high yielding Basmati rice varieties with good cooking quality is a gigantic task. Therefore, identifying the genomic regions governing the grain and cooked kernel dimension traits is of utmost importance for its use in marker-assisted breeding. Although several QTLs governing grain dimension traits have been reported, limited attempts have been made to map QTLs for grain and cooked kernel dimension traits of Basmati rice. In the current study, a population of recombinant inbred lines (RIL) was generated from a cross of Sonasal and Pusa Basmati 1121 (PB1121). In the RIL population, there was a significant positive correlation among the length (RRL: rough rice length, MRL: milled rice length, CKL: cooked kernel length) and breadth (RRB: rough rice breadth, MRB: milled rice breadth and CKB: cooked kernel breadth) of the related traits, while there was significant negative correlation between them. QTL mapping has led to the identification of four major genomic regions governing MRL and CKL. Two QTLs co-localize with the earlier reported major gene GS3 and a QTL qGRL7.1, while the remaining two QTLs viz., qCKL3.2 (qMRL3.2) and qCKL4.1 (qMRL4.1) were novel. The QTL qCKL3.2 has been bracketed to a genomic region of 0.78 Mb between the markers RM15247 and RM15281. Annotation of this region identified 18 gene models, of which the genes predicted to encode pentatricopeptides and brassinosteroid insensitive 1-associated receptor kinase 1 precursor may be the putative candidate genes. Furthermore, we identified a novel QTL qKER2.1 governing kernel elongation ratio (KER) in Basmati rice.

Pentaploidization Enriches the Genetic Diversity of Wheat by Enhancing the Recombination of AB Genomes

Allohexaploidization and continuous introgression play a key role in the origin and evolution of bread wheat. The genetic bottleneck of bread wheat resulting from limited germplasms involved in the origin and modern breeding may be compensated by gene flow from tetraploid wheat through introgressive hybridization. The inter-ploidy hybridization between hexaploid and tetraploid wheat generates pentaploid hybrids first, which absorbed genetic variations both from hexaploid and tetraploid wheat and have great potential for re-evolution and improvement in bread wheat. Therefore, understanding the effects of the pentaploid hybrid is of apparent significance in our understanding of the historic introgression and in informing breeding. In the current study, two sets of F2 populations of synthetic pentaploid wheat (SPW1 and SPW2) and synthetic hexaploid wheat (SHW1 and SHW2) were created to analyze differences in recombination frequency (RF) of AB genomes and distorted segregation of polymorphic SNP markers through SNP genotyping. Results suggested that (1) the recombination of AB genomes in the SPW populations was about 3- to 4-fold higher than that in the SHW populations, resulting from the significantly (P < 0.01) increased RF between adjacent and linked SNP loci, especially the variations that occurred in a pericentromeric region which would further enrich genetic diversity; (2) the crosses of hexaploid × tetraploid wheat could be an efficient way to produce pentaploid derivatives than the crosses of tetraploid × hexaploid wheat according to the higher germination rate found in the former crosses; (3) the high proportion of distorted segregation loci that skewed in favor of the female parent genotype/allele in the SPW populations might associate with the fitness and survival of the offspring. Based on the presented data, we propose that pentaploid hybrids should increasingly be used in wheat breeding. In addition, the contribution of gene flow from tetraploid wheat to bread wheat mediated by pentaploid introgressive hybridization also was discussed in the re-evolution of bread wheat.

Understanding the Nature of Hybrid Sterility and Divergence of Asian Cultivated Rice

Intraspecific hybrid sterility is a common form of postzygotic reproductive isolation in Asian cultivated rice, which is also the major obstacle to utilize the strong heterosis in the rice breeding program. Here, we review recent progress in classification and hybrid sterility in Asian cultivated rice. A genome-wide analysis of numerous wild relatives of rice and Asian cultivated rice has provided insights into the origin and differentiation of Asian cultivated rice, and divided Asian cultivated rice into five subgroups. More than 40 conserved and specific loci were identified to be responsible for the hybrid sterility between subgroup crosses by genetic mapping, which also contributed to the divergence of Asian cultivated rice. Most of the studies are focused on the sterile barriers between indica and japonica crosses, ignoring hybrid sterility among other subgroups, leading to neither a systematical understanding of the nature of hybrid sterility and subgroup divergence, nor effectively utilizing strong heterosis between the subgroups in Asian cultivated rice. Future studies will aim at identifying and characterizing genes for hybrid sterility and segregation distortion, comparing and understanding the molecular mechanism of hybrid sterility, and drawing a blueprint for intraspecific hybrid sterility loci derived from cross combinations among the five subgroups. These studies would provide scientific and accurate guidelines to overcome the intraspecific hybrid sterility according to the parent subgroup type identification, allowing the utilization of heterosis among subgroups, also helping us unlock the mysterious relationship between hybrid sterility and Asian cultivated rice divergence.

Mapping of dwarfing QTL of Ari1327, a semi-dwarf mutant of upland cotton

BACKGROUND

Upland Cotton (Gossypium hirsutum L.) has few cotton varieties suitable for mechanical harvesting. The plant height of the cultivar is one of the key features that need to modify. Hence, this study was planned to locate the QTL for plant height in a 60Co γ treated upland cotton semi-dwarf mutant Ari1327.

RESULTS

Interestingly, bulk segregant analysis (BSA) and genotyping by sequencing (GBS) methods exhibited that candidate QTL was co-located in the region of 5.80-9.66 Mb at D01 chromosome in two F2 populations. Using three InDel markers to genotype a population of 1241 individuals confirmed that the offspring's phenotype is consistent with the genotype. Comparative analysis of RNA-seq between the mutant and wild variety exhibited that Gh_D01G0592 was identified as the source of dwarfness from 200 genes. In addition, it was also revealed that the appropriate use of partial separation markers in QTL mapping can escalate linkage information.

CONCLUSIONS

Overwhelmingly, the results will provide the basis to reveal the function of candidate genes and the utilization of excellent dwarf genetic resources in the future.

Identification of Candidate Genes for Self-Compatibility in Perennial Ryegrass (Lolium perenne L.)

Self-incompatibility (SI) is a genetic mechanism preventing self-pollination in ~40% of plant species. Two multiallelic loci, called S and Z, control the gametophytic SI system of the grass family (Poaceae), which contains all major forage grasses. Loci independent from S and Z have been reported to disrupt SI and lead to self-compatibility (SC). A locus causing SC in perennial ryegrass (Lolium perenne L.) was previously mapped on linkage group (LG) 5 in an F2 population segregating for SC. Using a subset of the same population (n = 68), we first performed low-resolution quantitative trait locus (QTL) mapping to exclude the presence of additional, previously undetected contributors to SC. The previously reported QTL on LG 5 explained 38.4% of the phenotypic variation, and no significant contribution from other genomic regions was found. This was verified by the presence of significantly distorted markers in the region overlapping with the QTL. Second, we fine mapped the QTL to 0.26 centimorgan (cM) using additional 2,056 plants and 23 novel sequence-based markers. Using Italian ryegrass (Lolium multiflorum Lam.) genome assembly as a reference, the markers flanking SC were estimated to span a ~3 Mb region encoding for 57 predicted genes. Among these, seven genes were proposed as relevant candidate genes based on their annotation and function described in previous studies. Our study is a step forward to identify SC genes in forage grasses and provides diagnostic markers for marker-assisted introgression of SC into elite germplasm.

Identification of novel QTLs for grain fertility and associated traits to decipher poor grain filling of basal spikelets in dense panicle rice

High grain number is positively correlated with grain yield in rice, but it is compromised because of poor filling of basal spikelets in dense panicle bearing numerous spikelets. The phenomenon that turns the basal spikelets of compact panicle sterile in rice is largely unknown. In order to understand the factor(s) that possibly determines such spikelet sterility in compact panicle cultivars, QTLs and candidate genes were identified for spikelet fertility and associated traits like panicle compactness, and ethylene production that significantly influences the grain filling using recombinant inbred lines developed from a cross between indica rice cultivars, PDK Shriram (compact, high spikelet number) and Heera (lax, low spikelet number). Novel QTLs, qSFP1.1, qSFP3.1, and qSFP6.1 for spikelet fertility percentage; qIGS3.2 and qIGS4.1 for panicle compactness; and qETH1.2, qETH3.1, and qETH4.1 for ethylene production were consistently identified in both kharif seasons of 2017 and 2018. The comparative expression analysis of candidate genes like ERF3, AP2-like ethylene-responsive transcription factor, EREBP, GBSS1, E3 ubiquitin-protein ligase GW2, and LRR receptor-like serine/threonine-protein kinase ERL1 associated with identified QTLs revealed their role in poor grain filling of basal spikelets in a dense panicle. These candidate genes thus could be important for improving grain filling in compact-panicle rice cultivars through biotechnological interventions.

Multi-Environment Quantitative Trait Loci Mapping for Grain Iron and Zinc Content Using Bi-parental Recombinant Inbred Line Mapping Population in Pearl Millet

Pearl millet is a climate-resilient, nutritious crop with low input requirements that could provide economic returns in marginal agro-ecologies. In this study, we report quantitative trait loci (QTLs) for iron (Fe) and zinc (Zn) content from three distinct production environments. We generated a genetic linkage map using 210 F6 recombinant inbred line (RIL) population derived from the (PPMI 683 × PPMI 627) cross using genome-wide simple sequence repeats (SSRs). The molecular linkage map (seven linkage groups) of 151 loci was 3,273.1 cM length (Kosambi). The content of grain Fe in the RIL population ranged between 36 and 114 mg/Kg, and that of Zn from 20 to 106 mg/Kg across the 3 years (2014-2016) at over the three locations (Delhi, Dharwad, and Jodhpur). QTL analysis revealed a total of 22 QTLs for grain Fe and Zn, of which 14 were for Fe and eight were for Zn on three consecutive years at all locations. The observed phenotypic variance (R 2) explained by different QTLs for grain Fe and Zn content ranged from 2.85 (QGFe.E3.2014-2016_Q3) to 19.66% (QGFe.E1.2014-2016_Q3) and from 2.93 (QGZn.E3.2014-2016_Q3) to 25. 95% (QGZn.E1.2014-2016_Q1), respectively. Two constitutive expressing QTLs for both Fe and Zn co-mapped in this population, one on LG 2 and second one on LG 3. Inside the QTLs candidate genes such as Ferritin gene, Al3+ Transporter, K+ Transporters, Zn2+ transporters and Mg2+ transporters were identified using bioinformatics approaches. The identified QTLs and candidate genes could be useful in pearl millet population improvement programs, seed, restorer parents, and marker-assisted selection programs.

A new genetic locus for self-compatibility in the outcrossing grass species perennial ryegrass (Lolium perenne)

BACKGROUND

Self-incompatibility (SI) is a physiological mechanism that many flowering plants employ to prevent self-fertilization and maintain heterozygosity. In the grass family this is known to be controlled by a two locus (S-Z) system; however, the SI system is intrinsically leaky. Modifier genes of both the S and Z loci and a further locus, T, are known to override SI leading to self-fertilization and self-seed production. This has implications for the ecological and evolutionary success as well as the commercial breeding of grasses. Here we report a study where the genetic control of self-compatibility (SC) was determined from the results of self-pollinating an F2 population of perennial ryegrass from two independently derived inbred lines produced by single-seed descent.

METHODS

In vitro self-pollinations of 73 fertile plants were analysed. A genetic association analysis was made with a panel of 1863 single-nucleotide polymorphism (SNP) markers, generated through genotype-by-sequencing methodology. Markers were placed on a recombination map of seven linkage groups (LGs) created using Joinmap v.5. The seed set on self- and open-pollinated inflorescences was determined on 143 plants, including the 73 plants analysed for self-pollination response.

KEY RESULTS

Self-pollinations revealed a bimodal distribution of percentage SC with peaks at 50 and 100 %. A single quantitative trait locus (QTL) was identified with peak association for marker 6S14665z17875_11873 that mapped to LG 6. Peak position was associated with maximum marker segregation distortion. The self-compatible plants were equally fecund after self- and open pollination.

CONCLUSIONS

This is the first report in the Poaceae family of an SC locus located on LG 6. This new SC QTL discovery, as well as indicating the complex nature of the pollen-stigma recognition process and its evolutionary significance, provides an additional source of SC for breeding perennial ryegrass.

Substitution Mapping of a Locus Responsible for Hybrid Breakdown in Populations Derived From Interspecific Introgression Line

Hybrid breakdown, a form of postzygotic reproductive barrier, has been reported to hinder gene flow in many crosses between wild and cultivated rice. Here, the phenomenon of hybrid breakdown was observed as low-tillering (i.e., low tiller number) in some progeny of an interspecific cross produced in an attempt to introduce Oryza meridionalis Ng (W1625) chromosomal segments into Oryza sativa L. ssp. japonica "Taichung 65" (T65). Low-tillering lines were obtained in BC4-derived progeny from a cross between W1625 and "Taichung 65," but the locus for low-tillering could not be mapped in segregating populations. As a second approach to map the locus for low-tillering, we analyzed an F2 population derived from a cross between the low-tillering lines and a high-yielding indica cultivar, "Takanari." A major QTL for low-tillering, qLTN4, was detected between PCR-based markers MS10 and RM307 on the long arm of chromosome 4, with a LOD score of 15.6. The low-tillering phenotype was associated with weak growth and pale yellow phenotype; however, low-tillering plant had less reduction of grain fertility. In an F4 population (4896 plants), 563 recombinant plants were identified and the low-tillering locus was delimited to a 4.6-Mbp region between markers W1 and C5-indel3729. This region could not be further delimited because recombination is restricted in this region of qLTN4, which is near the centromere. Understanding the genetic basis of hybrid breakdown, including the low-tillering habit, will be important for improving varieties in rice breeding.

Identification and characterization of genes related to salt stress tolerance within segregation distortion regions of genetic map in F2 population of upland cotton

Segregation distortion (SD) is a genetic mechanism commonly found in segregating or stable populations. The principle behind this puzzles many researchers. The F₂ generation developed from wild Gossypium darwinii and G. hirsutum CCRI12 species was used to investigate the possible transcription factors within the segregation distortion regions (SDRs). The 384 out of 2763 markers were distorted in 29 SDRs on 18 chromosomes. Good collinearity was observed among genetic and physical maps of G. hirsutum and G. barbadense syntenic blocks. Total 568 genes were identified from SDRs of 18 chromosomes. Out of these genes, 128 belonged to three top-ranked salt-tolerant gene families. The DUF597 contained 8 uncharacterized genes linked to Pkinase (PF00069) gene family in the phylogenetic tree, while 15 uncharacterized genes clustered with the zinc finger gene family. Two hundred thirty four miRNAs targeted numerous genes, including ghr-miR156, ghr-miR399 and ghr-miR482, while others targeted top-ranked stress-responsive transcription factors. Moreover, these genes were involved in the regulation of numerous stress-responsive cis-regulatory elements. The RNA sequence data of fifteen upregulated genes were verified through the RT-qPCR. The expression profiles of two highly upregulated genes (Gh_D01G2015 and Gh_A01G1773) in salt-tolerant G. darwinii showed antagonistic expression in G. hirsutum. The results indicated that salt-tolerant genes have been possibly transferred from the wild G. darwinii species. A detailed functional analysis of these genes can be carried out which might be helpful in the future for gene cloning, transformation, gene editing and the development of salt-resistant cotton varieties.

Segregation distortion: high genetic load suggested by a Chinese shrimp family under high-intensity selection

Segregation distortion is a common phenomenon found in most genetic mapping studies and is an important resource to dissect the mechanism of action in gene loci that cause deviation. Marine animals possess high genetic diversity and genomic heterozygosity, they therefore are ideal model organisms to study segregation distortion induced by selection. In the present study, we constructed a full-sib family of Fenneropenaeus chinensis and exerted high-intensity selection on 10,000 incipient progenies. 2b-RAD method was employed in remaining 273 individuals to develop genome-wide SNPs for segregating analysis and 41,612 SNPs were developed. 50.77% of 32,229 high-quality representative markers deviated from the expected Mendelian ratio. Results showed that most of these distorted markers (91.57%) were influenced at zygotic level. Heterozygote excess (53.07%) and homozygous deletions (41.96%) may both play an important role, sum of which explained 95.03% of distortion after fertilization. However, further results identified highly probable linkage among deleterious alleles, which may account for a considerable portion of heterozygote excess rather than single locus with heterozygote advantage. Results of this study support a major role of deleterious alleles in genetic load, thus in favor of partial dominance hypothesis. It would also offer necessary recommendations for the formulation of breeding strategy in shrimps.

Mapping QTLs for enhancing early biomass derived from Aegilops tauschii in synthetic hexaploid wheat

Strong early vigour plays a crucial role in wheat yield improvement by enhancing resource utilization efficiency. Synthetic hexaploid wheat (SHW) combines the elite genes of tetraploid wheat with Aegilops tauschii and has been widely used in wheat genetic improvement for its abundant genetic variation. The two SHWs Syn79 and Syn80 were derived from the crossing of the same tetraploid wheat DOY1 with two different Ae. tauschii accessions, AT333 and AT428, respectively. The Syn80 possessed better early vigour traits than Syn79, theretically caused by their D genome from Ae. tauschii. To dissect their genetic basis in a hexaploid background, 203 recombinant inbred lines (RILs) derived from the cross of Syn79 x Syn80 were developed to detect quantitative trait loci (QTL) for four early biomass related traits: plant height (PH), tiller number (TN), shoot fresh weight (SFW) and shoot dry weight (SDW) per plant, under five different environmental conditions. Determined from the data of SNP markers, two genome regions on 1DS and 7D were stably associated with the four early biomass related traits showing pleiotropic effects. Four stable QTLs QPh.saas-1DS, QTn.saas-1DS, QSfw.saas-1DS and QSdw.saas-1DS explaining 7.92, 15.34, 9.64 and 10.15% of the phenotypic variation, respectively, were clustered in the region of 1DS from AX-94812958 to AX-110910133. Meanwhile, QPh.saas-7D, QTn.saas-7D, QSfw.saas-7D and QSdw.saas-7D were flanked by AX-109917900 and AX-110605376 on 7D, explaining 16.12, 24.35, 15.25 and 13.37% of the phenotypic variation on average, respectively. Moreover, these genomic QTLs on 1DS and 7D enhancing biomass in the parent Syn80 were from Ae. tauschii AT428. These findings suggest that these two QTLs from Ae. tauschii can be expressed stably in a hexaploid background at the jointing stage and be used for wheat improvement.

High-resolution genetic mapping of a novel bacterial blight resistance gene xa-45(t) identified from Oryza glaberrima and transferred to Oryza sativa

A novel recessive bacterial blight resistance locus designated as a xa-45(t) was identified from Oryza glaberrima accession IRGC 102600B, transferred to O. sativa and mapped to the long arm of chromosome 8 using ddRAD sequencing approach. The identified QTL spans 80 kb region on Nipponbare reference genome IRGSP-1.0 and contains 9 candidate genes. An STS marker developed from the locus LOC_Os08g42410 was found co-segregating with the trait and will be useful for marker-assisted transfer of this recessive resistance gene in breeding programs. Bacterial blight, caused by Xanthomonas oryzae pv. oryzae, is one of the major constraints of rice productivity in Southeast Asia. In spite of having 44 bacterial blight resistance genes from cultivated rice and wild species, the durability of resistance is always at stake due to the continually evolving nature of the pathogen and lack of suitable chemical control. Here, we report high-resolution genetic mapping of a novel bacterial blight resistance gene tentatively designated as a xa-45(t) from an introgression line derived from Oryza glaberrima accession IRGC 102600B. This introgression line was crossed with the susceptible rice indica cultivar cv. Pusa 44 to generate F₂ and F_2:3 populations for inheritance and mapping studies. The inheritance studies revealed the presence of single recessive locus controlling resistance to the Xanthomonas pathotype seven. A high-density linkage map was constructed using double-digest restriction-associated DNA sequencing of 96 F₂ populations along with the parents. The QTL mapping identified a major locus on the long arm of rice chromosome 8 with a LOD score of 33.22 between the SNP markers C8.26737175 and C8.26818765. The peak marker, C8.26810477, explains 49.8% of the total phenotypic variance and was positioned at 202.90 cM on the linkage map. This major locus spans 80 kb region on Nipponbare reference genome IRGSP-1.0 and contains 9 candidate genes. A co-segregating STS marker was developed from the LOC_Os08g42410 for efficient transfer of this novel gene to elite cultivars.

Genotyping-by-sequencing and multilocation evaluation of two interspecific backcross populations identify QTLs for yield-related traits in pigeonpea

This study has identified single-nucleotide polymorphism (SNP) markers associated with nine yield-related traits in pigeonpea by using two backcross populations (BP) developed through interspecific crosses and evaluating them at two locations and 3 years. In both the populations, markers have shown strong segregation distortion; therefore, a quantitative trait locus (QTL) mapping mixed model was used. A total of 86 QTLs explaining 12-21% phenotypic variation were detected in BP-1. On the other hand, 107 QTLs explaining 11-29% phenotypic variation were detected in BP-2. Although most QTLs were environment and trait specific, few stable and consistent QTLs were also detected. Interestingly, 11 QTLs in BP-2 were associated with more than one trait. Among these QTLs, eight QTLs associated with days to 50% flowering and days to 75% maturity were located on CcLG07. One SNP "S7_14185076" marker in BP-2 population has been found associated with four traits, namely days to 50% flowering, days to 75% maturity, primary branches per plant and secondary branches per plant with positive additive effect. Hence, the present study has not only identified QTLs for yield-related traits, but also discovered novel alleles from wild species, which can be used for improvement of traits through genomics-assisted breeding.

Construction of three half-sib SSR linkage maps derived from overwintering cultivated rice and segregation distortion loci mapping

Segregation distortion is a common phenomenon that has been observed in genetics and plant breeding; however, the mechanism of segregation distortion is unknown. In the present study, three half-sib F₂ populations derived from three japonica overwinter (perennial) rice varieties (W1, W2, and W3) crossed to the indica rice variety Minghui725 (MH725) were developed to construct three half-sib linkage maps. We established linkage map lengths of 2032.8, 2317.4, and 2108.7 cM with average intervals of 20.1, 20.5, and 19.7 cM using 101, 113, and 107 SSR markers in W1/MH725, W2/MH725, and W3/MH725, respectively. Discrepancies in marker order and genetic linkage distance occurred in the three half-sib linkage maps due to segregation distortion. A total of 88 markers exhibited segregation distortion across the three linkage maps at P < 0.01 level, 42 segregation distortion loci (SDLs) were detected across the three half-sib populations and exhibited variable LOD value that ranged from 3.2 (SDL2f) to 30.1 (SDL5d), and 13 of the 42 SDLs were repeatedly located at the same chromosomal regions of the previously published hybrid sterility quantitative trait loci. Data from this study provide an extensive archive for investigating the genetic characteristic of overwintering cultivated rice and the future exploration and innovation of overwintering rice breeding.

Identification and Molecular Mapping of a Gummy Stem Blight Resistance Gene in Wild Watermelon (Citrullus amarus) Germplasm PI 189225

Gummy stem blight (GSB), caused by Stagonosporopsis cucurbitacearum (syn. Didymella bryoniae), is a destructive foliar disease of watermelon in areas with hot and humid climates. The wild watermelon germplasm PI 189225 is a known source of resistance to GSB. The identification and use of molecular markers linked to resistance genes in the wild-type germplasm will speed up the introgression of GSB resistance into new watermelon varieties. An F₂ segregating population was obtained from a cross between the resistant wild watermelon genotype PI 189225 and the susceptible genotype K3. The F₂-derived F₃ families were inoculated with a single isolate of S. cucurbitacearum (JS002) from Jiangsu Academy of Agricultural Sciences. The results of the genetic analysis demonstrated that GSB resistance in PI 189225 was controlled by a major quantitative trait locus (QTL), temporarily designated Qgsb8.1. Based on the results of bulk sergeant analysis and sequencing, one associated region spanning 5.7 Mb (10,358,659 to 16,101,517) on chromosome 8 was identified as responsible for the resistance to GSB using the Δ(single-nucleotide polymorphism [SNP]-index) method. The result of a QTL linkage analysis with Kompetitive allele-specific PCR (KASP) SNP markers further mapped the GSB resistance locus between the SNP markers KASP_JS9383 and KASP_JS9168 in a region of 571.27 kb on chromosome 8. According to the watermelon gene annotation database, the region contains approximately 19 annotated genes and, of these 19 genes, 2 are disease resistance gene analogs: Cla001017 (coiled-coil nucleotide-binding site leucine-rich repeat resistance protein) and Cla001019 (pathogenesis related). Reverse-transcription quantitative PCR demonstrated that the expression of the two genes changed following S. cucurbitacearum infection, suggesting that they play important roles in GSB resistance in watermelon. This result will facilitate fine mapping and cloning of the Qgsb8.1 locus, and the linked markers will further provide a useful tool for marker-assisted selection of this locus in watermelon breeding programs.

Analysis of Crossover Events and Allele Segregation Distortion in Interspecific Citrus Hybrids by Single Pollen Genotyping

In citrus, a classical method of studying crossovers and segregation distortion (SD) is the genetic analysis of progenies. A new strategy combining fluorescence-activated cell sorting and whole genome amplification of haploid pollen nuclei with a large set of molecular markers, offers the opportunity to efficiently determine the frequency of crossovers and the identification of SD without the need to generate segregating populations. Here we have analyzed meiotic crossover events in a pollen nuclei population from "Eureka" lemon and the allelic SD was evaluated in a pollen nuclei population from a clementine × sweet orange hybrid ("CSO"). Data obtained from the "CSO" pollen nuclei population were compared to those obtained from genotyping of a segregating population ("RTSO") arising from a hand-made sexual hybridization between diploid non apomictic selected tangor (mandarin × sweet orange; "RTO" tangor) as female parent pollinated with "CSO" tangor as male parent. The analysis of crossovers rates on chromosome 1 revealed the presence of up to five crossovers events on one arm and four on the corresponding other arm, with an average of 1.97 crossovers per chromosome while no crossover events were observed in five "Eureka" lemon pollen nuclei. The rate of SD observed in "CSO" pollen nuclei (13.8%) was slightly lower than that recovered in the "RTSO" population (20.7%). In the pollen nuclei population, SD was found on linkage group (LG) 2, while the "RTSO" population showed SD on LGs 2 and 7. Potential male gametic selection mechanisms were distinguished in pollen grains, while in the population, mechanisms of gametophytic selection and/or zygotic selection were observed. This methodology is a very useful tool to facilitate research focused on the reproductive biology of citrus and study the mechanisms that affect crossovers and SD.

A high density SLAF-SNP genetic map and QTL detection for fibre quality traits in Gossypium hirsutum

Background

Upland Cotton (Gossypium hirsutum) is a very important cash crop known for its high quality natural fiber. Recent advances in sequencing technologies provide powerful tools with which to explore the cotton genome for single nucleotide polymorphism marker identification and high density genetic map construction toward more reliable quantitative trait locus mapping.

Results

In the present study, a RIL population was developed by crossing a Chinese high fiber quality cultivar (Yumian 1) and an American high fiber quality line (CA3084), with distinct genetic backgrounds. Specific locus amplified fragment sequencing (SLAF-seq) technology was used to discover SNPs, and a genetic map containing 6254 SNPs was constructed, covering 3141.72 cM with an average distance of 0.5 cM between markers. A total of 95 QTL were detected for fiber quality traits in three environments, explaining 5.5-24.6% of the phenotypic variance. Fifty-five QTL found in multiple environments were considered stable QTL. Nine of the stable QTL were found in all three environments. We identified 14 QTL clusters on 13 chromosomes, each containing one or more stable QTL.

Conclusion

A high-density genetic map of Gossypium hirsutum developed by using specific locus amplified fragment sequencing technology provides detailed mapping of fiber quality QTL, and identification of ‘stable QTL’ found in multiple environments. A marker-rich genetic map provides a foundation for fine mapping, candidate gene identification and marker-assisted selection of favorable alleles at stable QTL in breeding programs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5294-5) contains supplementary material, which is available to authorized users.

Anther culture in rice proportionally rescues microspores according to gametophytic gene effect and enhances genetic study of hybrid sterility

BACKGROUND

To investigate plant hybrid sterility, we studied interspecific hybrids of two cultivated rice species, Asian rice (Oryza sativa) and African rice (O. glaberrima). Male gametes of these hybrids display complete sterility owing to a dozen of hybrid sterility loci, termed HS loci, but this complicated genetic system remains poorly understood.

RESULTS

Microspores from these interspecific hybrids form sterile pollen but are viable at the immature stage. Application of the anther culture (AC) method caused these immature microspores to induce callus. The segregation distortion of 11 among 13 known HS loci was assessed in the callus population. Using many individual calli, fine mapping of the HS loci was attempted based on heterozygotes produced from chromosome segment substitution lines (CSSLs). Transmission ratio distortion (TRD) from microspores was detected at 6 of 11 HS loci in the callus population. The fine mapping of S ₁ and S ₁₉ loci using CSSLs revealed precise distances of markers from the positions of HS loci exhibiting excessive TRD.

CONCLUSIONS

We demonstrated that AC to generate callus populations derived from immature microspores is a useful methodology for genetic study. The callus population facilitated detection of TRD at multiple HS loci and dramatically shortened the process for mapping hybrid sterility genes.

SSR-Linkage map of interspecific populations derived from Gossypium trilobum and Gossypium thurberi and determination of genes harbored within the segregating distortion regions

Wild cotton species have significant agronomic traits that can be introgressed into elite cultivated varieties. The use of a genetic map is important in exploring, identification and mining genes which carry significant traits. In this study, 188 F₂mapping individuals were developed from Gossypium thurberi (female) and Gossypium trilobum (male), and were genotyped by using simple sequence repeat (SSR) markers. A total of 12,560 simple sequence repeat (SSR) markers, developed by Southwest University, thus coded SWU were screened out of which only 994 were found to be polymorphic, and 849 markers were linked in all the 13 chromosomes. The map had a length of 1,012.458 cM with an average marker distance of 1.193 cM. Segregation distortion regions (SDRs) were observed on Chr01, Chr02, Chr06, Chr07 Chr09, Chr10 and Chr11 with a large proportion of the SDR regions segregating towards the heterozygous allele. There was good syntenic block formation that revealed good collinearity between the genetic and physical map of G. raimondii, compared to the Dt_sub genome of the G. hirsutum and G. barbadense. A total of 2,496 genes were mined within the SSR related regions. The proteins encoding the mined genes within the SDR had varied physiochemical properties; their molecular weights ranged from 6.586 to 252.737 kDa, charge range of -39.5 to 52, grand hydropathy value (GRAVY) of -1.177 to 0.936 and isoelectric (pI) value of 4.087 to 12.206. The low GRAVY values detected showed that the proteins encoding these genes were hydrophilic in nature, a property common among the stress responsive genes. The RNA sequence analysis revealed more of the genes were highly upregulated in various stages of fiber development for instance; Gorai.002G241300 was highly up regulated at 5, 10, 20 and 25 day post anthesis (DPA). Validation through RT-qPCR further revealed that these genes mined within the SDR regions might be playing a significant role under fiber development stages, therefore we infer that Gorai.007G347600 (TFCA), Gorai.012G141600 (FOLB1), Gorai.006G024500 (NMD3), Gorai.002G229900 (LST8) and Gorai.002G235200 (NSA2) are significantly important in fiber development and in turn the quality, and further researches needed to be done to elucidate their exact roles in the fiber development process. The construction of the genetic map between the two wild species paves away for the mapping of quantitative trait loci (QTLs) since the average distance between the markers is small, and mining of genes on the SSR regions will provide an insight in identifying key genes that can be introgressed into the cultivated cotton cultivars.

Adaptive selection of founder segments and epistatic control of plant height in the MAGIC winter wheat population WM-800

BACKGROUND

Multi-parent advanced generation intercross (MAGIC) populations are a newly established tool to dissect quantitative traits. We developed the high resolution MAGIC wheat population WM-800, consisting of 910 F4:6 lines derived from intercrossing eight recently released European winter wheat cultivars.

RESULTS

Genotyping WM-800 with 7849 SNPs revealed a low mean genetic similarity of 59.7% between MAGIC lines. WM-800 harbours distinct genomic regions exposed to segregation distortion. These are mainly located on chromosomes 2 to 6 of the wheat B genome where founder specific DNA segments were positively or negatively selected. This suggests adaptive selection of individual founder alleles during population development. The application of a genome-wide association study identified 14 quantitative trait loci (QTL) controlling plant height in WM-800, including the known semi-dwarf genes Rht-B1 and Rht-D1 and a potentially novel QTL on chromosome 5A. Additionally, epistatic effects controlled plant height. For example, two loci on chromosomes 2B and 7B gave rise to an additive epistatic effect of 13.7 cm.

CONCLUSION

The present study demonstrates that plant height in the MAGIC-WHEAT population WM-800 is mainly determined by large-effect QTL and di-genic epistatic interactions. As a proof of concept, our study confirms that WM-800 is a valuable tool to dissect the genetic architecture of important agronomic traits.

Mapping Grain Iron and Zinc Content Quantitative Trait Loci in an Iniadi-Derived Immortal Population of Pearl Millet

Pearl millet is a climate-resilient nutritious crop requiring low inputs and is capable of giving economic returns in marginal agro-ecologies. In this study, we report large-effect iron (Fe) and zinc (Zn) content quantitative trait loci (QTLs) using diversity array technology (DArT) and simple sequence repeats (SSRs) markers to generate a genetic linkage map using 317 recombinant inbred line (RIL) population derived from the (ICMS 8511-S1-17-2-1-1-B-P03 × AIMP 92901-S1-183-2-2-B-08) cross. The base map [seven linkage groups (LGs)] of 196 loci was 964.2 cM in length (Haldane). AIMP 92901-S1-183-2-2-B-08 is an Iniadi line with high grain Fe and Zn, tracing its origin to the Togolese Republic, West Africa. The content of grain Fe in the RIL population ranged between 20 and 131 ppm (parts per million), and that of Zn from 18 to 110 ppm. QTL analysis revealed a large number of QTLs for high grain iron (Fe) and zinc (Zn) content. A total of 19 QTLs for Fe and Zn were detected, of which 11 were for Fe and eight were for Zn. The portion of the observed phenotypic variance explained by different QTLs for grain Fe and Zn content varied from 9.0 to 31.9% (cumulative 74%) and from 9.4 to 30.4% (cumulative 65%), respectively. Three large-effect QTLs for both minerals were co-mapped in this population, one on LG1 and two on LG7. The favorable QTL alleles of both mineral micronutrients were contributed by the male parent (AIMP 92901-deriv-08). Three putative epistasis interactions were observed for Fe content, while a single digenic interaction was found for Zn content. The reported QTLs may be useful in marker-assisted selection (MAS) programs, in genomic selection (GS) breeding pipelines for seed and restorer parents, and in population improvement programs for pearl millet.

First genetic linkage map of Lathyrus cicera based on RNA sequencing-derived markers: Key tool for genetic mapping of disease resistance

The Lathyrus cicera transcriptome was analysed in response to rust (Uromyces pisi) infection to develop novel molecular breeding tools with potential for genetic mapping of resistance in this robust orphan legume species. One RNA-seq library each was generated from control and rust-inoculated leaves from two L. cicera genotypes with contrasting quantitative resistance, de novo assembled into contigs and sequence polymorphisms were identified. In toto, 19,224 SNPs differentiate the susceptible from the partially resistant genotype's transcriptome. In addition, we developed and tested 341 expressed E-SSR markers from the contigs, of which 60.7% varied between the two L. cicera genotypes. A first L. cicera linkage map was created using part of the developed markers in a RIL population from the cross of the two genotypes. This map contains 307 markers, covered 724.2 cM and is organised in 7 major and 2 minor linkage groups, with an average mapping interval of 2.4 cM. The genic markers also enabled us to compare their position in L. cicera map with the physical position of the same markers mapped on Medicago truncatula genome, highlighting a high macrosyntenic conservation between both species. This study provides a large new set of genic polymorphic molecular markers with potential for mapping rust resistances. It represents the first step towards genomics-assisted precision breeding in L. cicera.

Genetic Determinants of Crop Timing and Quality Traits in Two Interspecific Petunia Recombinant Inbred Line Populations

The rate at which plants develop new nodes (development rate) is a major determinant of crop production time, yet the genetic control of this process, including genetic interactions with crop quality parameters, is poorly understood. We employed a modified genotyping-by-sequencing approach and generated genetic linkage maps with 6,291 and 3,297 single nucleotide polymorphisms (SNPs) for the interspecific Petunia recombinant inbred line (RIL) population - P. axillaris × P. exserta (AE) and P. integrifolia × P. axillaris (IA), respectively. Comparative mapping between the populations revealed perfect collinearity of marker order but different recombination frequency at the corresponding linkage groups (LGs). Quantitative trait loci (QTL) mapping conducted for development traits and other important quality traits indicated QTL clustered on chromosome 1, 2, 4 and 6 for the AE population and chromosome 1, 2, 5 and 6 for the IA population. Additionally, 209 differentially expressed unique transcripts were identified in shoot apex tissue between fast- and slow-developing RILs, 13 of which mapped to within 1 cM of a development rate QTL. These results will facilitate the identification of novel genes controlling crop timing and quality traits in Petunia and highlight the power of using multiple interspecific populations to elucidate genetic determinants of natural variation.

Genetic Architecture of a Rice Nested Association Mapping Population

Describing the genetic diversity in the gene pool of crops will provide breeders with novel resources for varietal improvement. Nested Association Mapping (NAM) populations are uniquely suited for characterizing parental diversity through the shuffling and fixation of parental haplotypes. Here, we describe a set of 1879 rice NAM lines created through the selfing and single-seed descent of F1 hybrids derived from elite IR64 indica crossed with 10 diverse tropical japonica lines. Genotyping data indicated tropical japonica alleles were captured at every queried locus despite the presence of segregation distortion factors. Several distortion loci were mapped, both shared and unique, among the 10 populations. Using two-point and multi-point genetic map calculations, our datasets achieved the ∼1500 cM expected map size in rice. Finally, we highlighted the utility of the NAM lines for QTL mapping, including joint analysis across the 10 populations, by confirming known QTL locations for the trait days to heading.

Enhancing genetic gain in the era of molecular breeding

As one of the important concepts in conventional quantitative genetics and breeding, genetic gain can be defined as the amount of increase in performance that is achieved annually through artificial selection. To develop pro ducts that meet the increasing demand of mankind, especially for food and feed, in addition to various industrial uses, breeders are challenged to enhance the potential of genetic gain continuously, at ever higher rates, while they close the gaps that remain between the yield potential in breeders' demonstration trials and the actual yield in farmers' fields. Factors affecting genetic gain include genetic variation available in breeding materials, heritability for traits of interest, selection intensity, and the time required to complete a breeding cycle. Genetic gain can be improved through enhancing the potential and closing the gaps, which has been evolving and complemented with modern breeding techniques and platforms, mainly driven by molecular and genomic tools, combined with improved agronomic practice. Several key strategies are reviewed in this article. Favorable genetic variation can be unlocked and created through molecular and genomic approaches including mutation, gene mapping and discovery, and transgene and genome editing. Estimation of heritability can be improved by refining field experiments through well-controlled and precisely assayed environmental factors or envirotyping, particularly for understanding and controlling spatial heterogeneity at the field level. Selection intensity can be significantly heightened through improvements in the scale and precision of genotyping and phenotyping. The breeding cycle time can be shortened by accelerating breeding procedures through integrated breeding approaches such as marker-assisted selection and doubled haploid development. All the strategies can be integrated with other widely used conventional approaches in breeding programs to enhance genetic gain. More transdisciplinary approaches, team breeding, will be required to address the challenge of maintaining a plentiful and safe food supply for future generations. New opportunities for enhancing genetic gain, a high efficiency breeding pipeline, and broad-sense genetic gain are also discussed prospectively.

SSR Marker Development, Linkage Mapping, and QTL Analysis for Establishment Rate in Common Bermudagrass

Common bermudagrass has been widely used as a major warm-season turf, forage, and soil stabilization grass in the southern United States. However, codominant marker development, linkage, and quantitative trait loci (QTL) mapping resources are limited in the important taxon. Accordingly, the objectives of this study were to develop simple sequence repeat (SSR) markers, construct a genetic map, and identify genomic regions associated with establishment rate. Five genomic SSR libraries were constructed, sequenced, and used in the development of 1003 validated SSR primer pairs (PPs). A linkage map was constructed using a first-generation selfed population derived from a genotype A12359 (2 = 4 = 36). A total of 249 polymorphic SSR PPs were mapped to 18 linkage groups (LGs). The total length of the map is 1094.7 cM, with an average marker interval of 4.3 cM. Ninety-eight out of 252 mapped loci (39%) were found to be distorted from the Mendelian 1:2:1 segregation ratio. Among the other 154 nondistorted loci, 88 coupling vs. 66 repulsion linkage phases were observed to confirm the allopolyploid origin of the parent. Ground coverage (GCR) phenotypic data in the establishment stage were collected in two replicated field trials. Quantitative trait loci mapping identified five genomic regions significantly related to the trait. The findings of this study provide valuable genetic tools and resources for genomic research, genetic improvement, and breeding new cultivars in the species.

Sequencing of Single Pollen Nuclei Reveals Meiotic Recombination Events at Megabase Resolution and Circumvents Segregation Distortion Caused by Postmeiotic Processes

Meiotic recombination is a fundamental mechanism to generate novel allelic combinations which can be harnessed by breeders to achieve crop improvement. The recombination landscape of many crop species, including the major crop barley, is characterized by a dearth of recombination in 65% of the genome. In addition, segregation distortion caused by selection on genetically linked loci is a frequent and undesirable phenomenon in double haploid populations which hampers genetic mapping and breeding. Here, we present an approach to directly investigate recombination at the DNA sequence level by combining flow-sorting of haploid pollen nuclei of barley with single-cell genome sequencing. We confirm the skewed distribution of recombination events toward distal chromosomal regions at megabase resolution and show that segregation distortion is almost absent if directly measured in pollen. Furthermore, we show a bimodal distribution of inter-crossover distances, which supports the existence of two classes of crossovers which are sensitive or less sensitive to physical interference. We conclude that single pollen nuclei sequencing is an approach capable of revealing recombination patterns in the absence of segregation distortion.

Development of genome-wide PCR-based markers from insertion, deletion and single nucleotide polymorphisms for closely related Japanese rice cultivars and identification of QTLs for the appearance of cooked rice and polished rice

Appearance of rice grain is an important property, affecting its acceptance by consumers. Moreover, appearance is a complex characteristic involving many components, including glossiness and whiteness. The genetic bases for the glossiness of cooked rice and the whiteness of polished rice (WPR) were determined using 133 recombinant inbred lines (RILs) derived from a cross between two closely related cultivars from Hokkaido, Joiku462, with high glossiness and whiteness, and Yukihikari, an ancestor of Joiku462 with low glossiness and whiteness. Analyses identified 167 genome-wide InDel markers, five cleaved amplified polymorphic sequences (CAPS) and eight derived CAPS markers differentiating the parental lines. The glossiness area (GLA) and glossiness strength (GLS) of cooked rice and WPR were determined for RILs in two locations, Pippu and Sapporo, Hokkaido. Four QTLs were detected. qGLA10 and qGLS9 were detected on chromosomes 10 and 9, respectively, with both being significant at both geographic locations. qWPR1 on chromosome 1 was significant at Pippu, and qWPR4 on chromosome 4 was significant at Sapporo. The Joiku462 alleles at all QTLs increased each trait. The PCR-based markers flanking these four QTLs may be useful for improvement of GLA, GLS and WPR.

Mapping quantitative trait loci responsible for resistance to Bakanae disease in rice

BACKGROUND

Bakanae or foot rot disease caused by Fusarium fujikuroi [teleomorph: Gibberella fujikuroi (Sawada) Ito] is emerging as a serious disease in rice. The disease causes both quantitative and qualitative losses to the grains under the field conditions. Breeding for resistance to Bakanae disease is a promising strategy to manage this emerging disease. In this study, we used a population of 168 F14 recombinant inbred lines (RILs) derived from two indica rice parents Pusa 1342, a highly resistant variety and Pusa Basmati 1121, a highly susceptible variety to map quantitative trait loci (QTLs) governing resistance against Bakanae disease.

RESULTS

The disease reaction of 168 F14 RILs were measured on the seedlings inoculated using Fusarium fujikuroi culture using high-throughput screening protocol under glasshouse conditions. Utilizing inclusive composite interval mapping, three QTLs governing resistance to Bakanae were identified, namely qBK1.1, qBK1.2 and qBK1.3 which accounted 4.76, 24.74 and 6.49 % of phenotypic variation, respectively. The major effect QTL designated qBK1.2 was mapped in 0.26 Mb region between RM5336 and RM10153. A total of 55 annotated genes were identified within the identified QTL region qBK1.2.

CONCLUSIONS

The novel QTLs identified in this study are useful resource for efficiently breeding rice cultivars resistant to Bakanae disease. This is the first report on identification of QTLs governing resistance against Bakanae in rice using inclusive composite interval mapping strategy in a RIL population.

Genome-Wide SNP Linkage Mapping and QTL Analysis for Fiber Quality and Yield Traits in the Upland Cotton Recombinant Inbred Lines Population

It is of significance to discover genes related to fiber quality and yield traits and tightly linked markers for marker-assisted selection (MAS) in cotton breeding. In this study, 188 F8 recombinant inbred lines (RILs), derived from a intraspecific cross between HS46 and MARCABUCAG8US-1-88 were genotyped by the cotton 63K single nucleotide polymorphism (SNP) assay. Field trials were conducted in Sanya, Hainan Province, during the 2014-2015 cropping seasons under standard conditions. Results revealed significant differences (P < 0.05) among RILs, environments and replications for fiber quality and yield traits. Broad-sense heritabilities of all traits including fiber length, fiber uniformity, micronaire, fiber elongation, fiber strength, boll weight, and lint percentage ranged from 0.26 to 0.66. A 1784.28 cM (centimorgans) linkage map, harboring 2618 polymorphic SNP markers, was constructed, which had 0.68 cM per marker density. Seventy-one quantitative trait locus (QTLs) for fiber quality and yield traits were detected on 21 chromosomes, explaining 4.70∼32.28% phenotypic variance, in which 16 were identified as stable QTLs across two environments. Meanwhile, 12 certain regions were investigated to be involved in the control of one (hotspot) or more (cluster) traits, mainly focused on Chr05, Chr09, Chr10, Chr14, Chr19, and Chr20. Nineteen pairs of epistatic QTLs (e-QTLs) were identified, of which two pairs involved in two additive QTLs. These additive QTLs, e-QTLs, and QTL clusters were tightly linked to SNP markers, which may serve as target regions for map-based cloning, gene discovery, and MAS in cotton breeding.

Development and molecular characterization of a doubled haploid population derived from a hybrid between japonica rice and wide compatible indica rice

Doubled haploid (DH) populations, particularly those from subspecies crosses possessing the wide compatible gene S5ⁿ , are important germplasm resources for rice genetic studies and breeding, but their feature and potential have not been fully assessed and explored. In the present study, we produced a DH population from the hybrid of japonica 668B and wide compatible indica T23. Genotyping of the S5 locus with allele-specific markers for ORF3, ORF4 and ORF5 revealed a potential recombination hot spot in the ORF3-ORF4 region. Haplotyping analysis revealed that 21/34 subspecies specific Indel markers segregated in distortion in the DH population, with a few lines having indica alleles either extremely low (1.7%) or high (98.3%), with little effect of the S5 allele. While DH lines with the S5ⁿ allele had higher frequency of indica alleles, no effect of the S5ⁿ allele was observed on all agronomic traits but flowering time. Taken together, the present study advanced understanding of the genetics of wide crosses in general, and DH production in particular between the two rice subspecies, and the new DH population generated will become a useful resource for rice genetic study and breeding in the future.

Comparative Analysis of Regions with Distorted Segregation in Three Diploid Populations of Potato

Genes associated with gametic and zygotic selection could underlie segregation distortion, observed as alterations of expected Mendelian genotypic frequencies in mapping populations. We studied highly dense genetic maps based on single nucleotide polymorphisms to elucidate the genetic nature of distorted segregation in potato. Three intra- and interspecific diploid segregating populations were used. DRH and D84 are crosses between the sequenced doubled monoploid DM 1-3 516 R44 Solanum tuberosum Group Phureja and either RH89-039-16 S. tuberosum or 84SD22, a S. tuberosum × S. chacoense hybrid. MSX902 is an interspecific cross between 84SD22 and Ber83 S. berthaultii × 2 × species mosaic. At the 0.05 significance level, 21%, 57%, and 51% of the total markers mapped in DRH, D84, and MSX902 exhibited distorted segregation, respectively. Segregation distortion regions for DRH were located on chromosomes 9 and 12; for D84 on chromosomes 2, 3, 4, 6, 7, and 8; and on chromosomes 1, 2, 7, 9, and 12 for MSX902. In general, each population had unique segregation distortion regions and directions of distortion. Interspecific crosses showed greater levels of distorted segregation and lower recombination rates as determined from the male parents. The different genomic regions where the segregation distortion regions occurred in the three populations likely reflect unique genetic combinations producing distorted segregation.

Development and Genetic Characterization of an Advanced Backcross-Nested Association Mapping (AB-NAM) Population of Wild × Cultivated Barley

The ability to access alleles from unadapted germplasm collections is a long-standing problem for geneticists and breeders. Here we developed, characterized, and demonstrated the utility of a wild barley advanced backcross-nested association mapping (AB-NAM) population. We developed this population by backcrossing 25 wild barley accessions to the six-rowed malting barley cultivar Rasmusson. The 25 wild barley parents were selected from the 318 accession Wild Barley Diversity Collection (WBDC) to maximize allelic diversity. The resulting 796 BC2F4:6 lines were genotyped with 384 SNP markers, and an additional 4022 SNPs and 263,531 sequence variants were imputed onto the population using 9K iSelect SNP genotypes and exome capture sequence of the parents, respectively. On average, 96% of each wild parent was introgressed into the Rasmusson background, and the population exhibited low population structure. While linkage disequilibrium (LD) decay (r(2) = 0.2) was lowest in the WBDC (0.36 cM), the AB-NAM (9.2 cM) exhibited more rapid LD decay than comparable advanced backcross (28.6 cM) and recombinant inbred line (32.3 cM) populations. Three qualitative traits: glossy spike, glossy sheath, and black hull color were mapped with high resolution to loci corresponding to known barley mutants for these traits. Additionally, a total of 10 QTL were identified for grain protein content. The combination of low LD, negligible population structure, and high diversity in an adapted background make the AB-NAM an important tool for high-resolution gene mapping and discovery of novel allelic variation using wild barley germplasm.

Construction of a high-density, high-quality genetic map of cultivated lotus (Nelumbo nucifera) using next-generation sequencing

BACKGROUND

The sacred lotus (Nelumbo nucifera) is widely cultivated in China for its edible rhizomes and seeds. Traditional plant breeding methods have been used to breed cultivars with increased yields and quality of rhizomes and seeds with limited success. Currently, the available genetic maps and molecular markers in lotus are too limited to be useful for molecular genetics based breeding programs. However, the development of next-generation sequencing (NGS) technologies has enabled large-scale identification of single-nucleotide polymorphisms (SNPs) for genetic map construction. In this study, we constructed an SNP-based high-density genetic map for cultivated lotus using double digest restriction site-associated DNA sequencing (ddRADseq).

RESULTS

An F2 population of 96 individuals was derived from a cross between the rhizome lotus cultivar 'Juwuba' (male parent) and the seed lotus cultivar 'Mantianxing' (female parent). Genomic DNAs from this population were digested with the restriction enzymes EcoRI and MspI and then sequenced. In total, 133.65 Gb of raw data containing 1,088,935,610 pair-end reads were obtained. The coverage of reads on a reference genome was 7.2 % for the female parent, 6.56 % for the male parent, and 1.46 % for F2 individuals. From these reads, 10,753 valid SNP markers were used for genetic map construction. Finally, 791 bin markers (so-segregated adjacent SNPs treated as a bin marker), consisting of 8,971 SNP markers, were sorted into 8 linkage groups (LGs) that spanned 581.3 cM, with an average marker interval of 0.74 cM. A total of 809 genome sequence scaffolds, covering about 565.9 cM of the wild sacred lotus genome, were anchored on the genetic map, accounting for 70.6 % of the genome assembly.

CONCLUSIONS

This study reports the large-scale discovery of SNPs between cultivars of rhizome and seed lotus using a ddRADseq library combined with NGS. These SNPs have been used to construct the first high-density genetic map for cultivated lotus that can serve as a genomic reference and will facilitate genetic mapping of important traits in the parental cultivars.

High-density linkage mapping and distribution of segregation distortion regions in the oak genome

We developed the densest single-nucleotide polymorphism (SNP)-based linkage genetic map to date for the genus Quercus An 8k gene-based SNP array was used to genotype more than 1,000 full-sibs from two intraspecific and two interspecific full-sib families of Quercus petraea and Quercus robur A high degree of collinearity was observed between the eight parental maps of the two species. A composite map was then established with 4,261 SNP markers spanning 742 cM over the 12 linkage groups (LGs) of the oak genome. Nine genomic regions from six LGs displayed highly significant distortions of segregation. Two main hypotheses concerning the mechanisms underlying segregation distortion are discussed: genetic load vs. reproductive barriers. Our findings suggest a predominance of pre-zygotic to post-zygotic barriers.

Genotyping-by-Sequencing on Pooled Samples and its Use in Measuring Segregation Bias during the Course of Androgenesis in Barley

Estimation of allelic frequencies is often required in breeding but genotyping many individuals at many loci can be expensive. We have developed a genotyping-by-sequencing (GBS) approach for estimating allelic frequencies on pooled samples (Pool-GBS) and used it to examine segregation distortion in doubled haploid (DH) populations of barley ( L.). In the first phase, we genotyped each line individually and exploited these data to explore a strategy to call single nucleotide polymorphisms (SNPs) on pooled reads. We measured both the number of SNPs called and the variance of the estimated allelic frequencies at various depths of coverage on a subset of reads containing 5 to 25 million reads. We show that allelic frequencies could be cost-effectively and accurately estimated at a depth of 50 reads per SNP using 15 million reads. This Pool-GBS approach yielded 1984 SNPs whose allelic frequency estimates were highly reproducible (CV = 10.4%) and correlated ( = 0.9167) with the "true" frequency derived from analysis of individual lines. In a second phase, we used Pool-GBS to investigate segregation bias throughout androgenesis from microspores to a population of regenerated plants. No strong bias was detected among the microspores resulting from the meiotic divisions, whereas significant biases could be shown to arise during embryo formation and plant regeneration. In summary, this methodology provides an approach to estimate allelic frequencies more efficiently and on materials that are unsuitable for individual analysis. In addition, it allowed us to shed light on the process of androgenesis in barley.

SNP-Based QTL Mapping of 15 Complex Traits in Barley under Rain-Fed and Well-Watered Conditions by a Mixed Modeling Approach

This study identified single nucleotide polymorphism (SNP) markers associated with 15 complex traits in a breeding population of barley (Hordeum vulgare L.) consisting of 137 recombinant chromosome substitution lines (RCSL), evaluated under contrasting water availability conditions in the Mediterranean climatic region of central Chile. Given that markers showed a very strong segregation distortion, a quantitative trait locus/loci (QTL) mapping mixed model was used to account for the heterogeneity in genetic relatedness between genotypes. Fifty-seven QTL were detected under rain-fed conditions, which accounted for 5-22% of the phenotypic variation. In full irrigation conditions, 84 SNPs were significantly associated with the traits studied, explaining 5-35% of phenotypic variation. Most of the QTL were co-localized on chromosomes 2H and 3H. Environment-specific genomic regions were detected for 12 of the 15 traits scored. Although most QTL-trait associations were environment and trait specific, some important and stable associations were also detected. In full irrigation conditions, a relatively major genomic region was found underlying hectoliter weight (HW), on chromosome 1H, which explained between 27% (SNP 2711-234) and 35% (SNP 1923-265) of the phenotypic variation. Interestingly, the locus 1923-265 was also detected for grain yield at both environmental conditions, accounting for 9 and 18%, in the rain-fed and irrigation conditions, respectively. Analysis of QTL in this breeding population identified significant genomic regions that can be used for marker-assisted selection (MAS) of barley in areas where drought is a significant constraint.

A High-Density SNP and SSR Consensus Map Reveals Segregation Distortion Regions in Wheat

Segregation distortion is a widespread phenomenon in plant and animal genomes and significantly affects linkage map construction and identification of quantitative trait loci (QTLs). To study segregation distortion in wheat, a high-density consensus map was constructed using single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) markers by merging two genetic maps developed from two recombinant-inbred line (RIL) populations, Ning7840 × Clark and Heyne × Lakin. Chromosome regions with obvious segregation distortion were identified in the map. A total of 3541 SNPs and 145 SSRs were mapped, and the map covered 3258.7 cM in genetic distance with an average interval of 0.88 cM. The number of markers that showed distorted segregation was 490 (18.5%) in the Ning7840 × Clark population and 225 (10.4%) in the Heyne × Lakin population. Most of the distorted markers (630) were mapped in the consensus map, which accounted for 17.1% of mapped markers. The majority of the distorted markers clustered in the segregation distortion regions (SDRs) on chromosomes 1B, 2A, 2B, 3A, 3B, 4B, 5A, 5B, 5D, 6B, 7A, and 7D. All of the markers in a given SDR skewed toward one of the parents, suggesting that gametophytic competition during zygote formation was most likely one of the causes for segregation distortion in the populations.

Multiplex SSR-PCR approaches for semi-automated genotyping and characterization of loci linked to blast disease resistance genes in rice

In the present study, 63 polymorphic microsatellite markers related to rice blast resistance genes were fluorescently labelled at the 5'-end with either 6-FAM or HEX using the G5 dye set and incorporated into a multiplex SSR-PCR for the detection of fragments using an automated system. For rice F3 families obtained from crosses between Pongsu Seribu 2 (Malaysian blast resistant cultivar) and Mahsuri (a susceptible rice cultivar), the genotypes for 13 designated multiplex SSR panels were determined. The genotyping assays were performed using a capillary-based ABIPRISM 3100 genetic analyser. The sizes of the SSRs alleles observed in the range from 79 to 324 bp. The observed marker segregation data were analysed using the Chi(2) test. A genetic linkage map covering ten chromosomes and comprising 63 polymorphic SSR markers was constructed, and the distorted loci were localised to linkage groups. The results indicated that distorted loci are presented on eight chromosomes.