Estimating the effects of population size and type on the accuracy of genetic maps

High-density linkage to physical mapping in a unique Tall × Dwarf coconut (Cocos nucifera L.) outbred F2 uncovers a major QTL for flowering time colocalized with the FLOWERING LOCUS T (FT)

Introduction

The coconut tree crop (Cocos nucifera L.) provides vital resources for millions of people worldwide. Coconut germplasm is largely classified into 'Tall' (Typica) and 'Dwarf' (Nana) types. While Tall coconuts are outcrossing, stress tolerant, and late flowering, Dwarf coconuts are inbred and flower early with a high rate of bunch emission. Precocity determines the earlier production of a plantation and facilitates management and harvest.

Methods

A unique outbred F2 population was used, generated by intercrossing F1 hybrids between Brazilian Green Dwarf from Jiqui (BGDJ) and West African Tall (WAT) cultivars. Single-nucleotide polymorphism (SNP) markers fixed for alternative alleles in the two varieties, segregating in an F2 configuration, were used to build a high-density linkage map with ~3,000 SNPs, anchored to the existing chromosome-level genome assemblies, and a quantitative trait locus (QTL) mapping analysis was carried out.

Results

The linkage map established the chromosome numbering correspondence between the two reference genome versions and the relationship between recombination rate, physical distance, and gene density in the coconut genomes. Leveraging the strong segregation for precocity inherited from the Dwarf cultivar in the F2, a major effect QTL with incomplete dominance was mapped for flowering time. FLOWERING LOCUS T (FT) gene homologs of coconut previously described as putatively involved in flowering time by alternative splice variant analysis were colocalized within a ~200-kb window of the major effect QTL [logarithm of the odds (LOD) = 11.86].

Discussion

Our work provides strong phenotype-based evidence for the role of the FT locus as the putative underlying functional variant for the flowering time difference between Dwarf and Tall coconuts. Major effect QTLs were also detected for developmental traits of the palm, plausibly suggesting pleiotropism of the FT locus for other precocity traits. Haplotypes of the two SNPs flanking the flowering time QTL inherited from the Dwarf parent BGDJ caused a reduction in the time to flower of approximately 400 days. These SNPs could be used for high-throughput marker-assisted selection of early-flowering and higher-productivity recombinant lines, providing innovative genetic material to the coconut industry.

Omics technologies towards sesame improvement: a review

Genetic improvement of sesame (Sesamum indicum L.), one of the most important oilseed crops providing edible oil, proteins, minerals, and vitamins, is important to ensure a balanced diet for the growing world population. Increasing yield, seed protein, oil, minerals, and vitamins is urgently needed to meet the global demand. The production and productivity of sesame is very low due to various biotic and abiotic stresses. Therefore, various efforts have been made to combat these constraints and increase the production and productivity of sesame through conventional breeding. However, less attention has been paid to the genetic improvement of the crop through modern biotechnological methods, leaving it lagging behind other oilseed crops. Recently, however, the scenario has changed as sesame research has entered the era of "omics" and has made significant progress. Therefore, the purpose of this paper is to provide an overview of the progress made by omics research in improving sesame. This review presents a number of efforts that have been made over past decade using omics technologies to improve various traits of sesame, including seed composition, yield, and biotic and abiotic resistant varieties. It summarizes the advances in genetic improvement of sesame using omics technologies, such as germplasm development (web-based functional databases and germplasm resources), gene discovery (molecular markers and genetic linkage map construction), proteomics, transcriptomics, and metabolomics that have been carried out in the last decade. In conclusion, this review highlights future directions that may be important for omics-assisted breeding in sesame genetic improvement.

Population size in QTL detection using quantile regression in genome-wide association studies

The aim of this study was to evaluate the performance of Quantile Regression (QR) in Genome-Wide Association Studies (GWAS) regarding the ability to detect QTLs (Quantitative Trait Locus) associated with phenotypic traits of interest, considering different population sizes. For this, simulated data was used, with traits of different levels of heritability (0.30 and 0.50), and controlled by 3 and 100 QTLs. Populations of 1,000 to 200 individuals were defined, with a random reduction of 100 individuals for each population. The power of detection of QTLs and the false positive rate were obtained by means of QR considering three different quantiles (0.10, 0.50 and 0.90) and also by means of the General Linear Model (GLM). In general, it was observed that the QR models showed greater power of detection of QTLs in all scenarios evaluated and a relatively low false positive rate in scenarios with a greater number of individuals. The models with the highest detection power of true QTLs at the extreme quantils (0.10 and 0.90) were the ones with the highest detection power of true QTLs. In contrast, the analysis based on the GLM detected few (scenarios with larger population size) or no QTLs in the evaluated scenarios. In the scenarios with low heritability, QR obtained a high detection power. Thus, it was verified that the use of QR in GWAS is effective, allowing the detection of QTLs associated with traits of interest even in scenarios with few genotyped and phenotyped individuals.

Genetic mapping of the powdery mildew resistance gene Pm7 on oat chromosome 5D

Three independent experiments with different genetic backgrounds mapped the resistance gene Pm7 in the oat genome to the distal part of the long arm of chromosome 5D. Resistance of oat to Blumeria graminis DC. f. sp. avenae is an important breeding goal in Central and Western Europe. In this study, the position of the effective and widely used resistance gene Pm7 in the oat genome was determined based on three independent experiments with different genetic backgrounds: genome-wide association mapping in a diverse set of inbred oat lines and binary phenotype mapping in two bi-parental populations. Powdery mildew resistance was assessed in the field as well as by detached leaf tests in the laboratory. Genotyping-by-sequencing was conducted to establish comprehensive genetic fingerprints for subsequent genetic mapping experiments. All three mapping approaches located the gene to the distal part of the long arm of chromosome 5D in the hexaploid oat genome sequences of OT3098 and 'Sang.' Markers from this region were homologous to a region of chromosome 2Ce of the C-genome species, Avena eriantha, the donor of Pm7, which appears to be the ancestral source of a translocated region on the hexaploid chromosome 5D.

Construction of a high density genetic linkage map to define the locus conferring seedlessness from Mukaku Kishu mandarin

Mukaku Kishu ('MK'), a small sized mandarin, is an important source of seedlessness in citrus breeding. Identification and mapping the gene(s) governing 'MK' seedlessness will expedite seedless cultivar development. In this study, two 'MK'-derived mapping populations- LB8-9 Sugar Belle^® ('SB') × 'MK' (N=97) and Daisy ('D') × 'MK' (N=68) were genotyped using an Axiom_Citrus56 Array encompassing 58,433 SNP probe sets, and population specific male and female parent linkage maps were constructed. The parental maps of each population were integrated to produce sub-composite maps, which were further merged to develop a consensus linkage map. All the parental maps (except 'MK_D') had nine major linkage groups, and contained 930 ('SB'), 810 ('MK_SB'), 776 ('D') and 707 ('MK_D') SNPs. The linkage maps displayed 96.9 ('MK_D') to 98.5% ('SB') chromosomal synteny with the reference Clementine genome. The consensus map was comprised of 2588 markers including a phenotypic seedless (Fs)-locus and spanned a genetic distance of 1406.84 cM, with an average marker distance of 0.54 cM, which is substantially lower than the reference Clementine map. For the phenotypic Fs-locus, the distribution of seedy and seedless progenies in both 'SB' × 'MK' (55:42, χ² = 1.74) and 'D' × 'MK' populations (33:35, χ² = 0.06) followed a test cross pattern. The Fs-locus mapped on chromosome 5 with SNP marker 'AX-160417325' at 7.4 cM in 'MK_SB' map and between two SNP markers 'AX-160536283' and 'AX-160906995' at a distance of 2.4 and 4.9 cM, respectively in 'MK_D' map. The SNPs 'AX-160417325' and 'AX-160536283' correctly predicted seedlessness of 25-91.9% progenies in this study. Based on the alignment of flanking SNP markers to the Clementine reference genome, the candidate gene for seedlessness hovered in a ~ 6.0 Mb region between 3.97 Mb (AX-160906995) to 10.00 Mb (AX-160536283). This region has 131 genes of which 13 genes (belonging to seven gene families) reportedly express in seed coat or developing embryo. The findings of the study will prove helpful in directing future research for fine mapping this region and eventually underpinning the exact causative gene governing seedlessness in 'MK'.

High-Density Genetic Linkage Map of the Southern Blue-ringed Octopus (Octopodidae: Hapalochlaena maculosa)

Genetic linkage maps provide a useful resource for non-model genomes and can aid in genome reassembly to form more contiguous pseudo-chromosomes. We present the first linkage map of any cephalopod, H. maculosa, composed of 47 linkage groups (LG). A total of 2166 single nucleotide polymorphisms and 2455 presence–absence variant loci were utilised by Lep-Map3 in linkage map construction. The map length spans 2016.62 cM with an average marker distance of 0.85 cM. Integration of the recent H. maculosa genome allowed 1151 scaffolds comprising 34% of the total genomic sequence to be orientated and/or placed using 1278 markers across all 47 LG. The linkage map generated provides a new perspective on HOX gene distribution in octopods. In the H. maculosa linkage map three (SCR, LOX4 and POST1) of six identified HOX genes (HOX1/LAB, SCR, LOX2, LOX4, LOX5, POST1) were located within the same LG (LG 9). The generation of a linkage map for H. maculosa has provided a valuable resource for understanding the evolution of cephalopod genomes and will provide a base for future work.

Genetic mapping of the Andean anthracnose resistance gene present in the common bean cultivar BRSMG Realce

The rajado seeded Andean bean (Phaseolus vulgaris L.) cultivar BRSMG Realce (striped seed coat) developed by Embrapa expressed a high level of anthracnose resistance, caused by Colletotrichum lindemuthianum, in field and greenhouse screenings. The main goal of this study was to evaluate the inheritance of anthracnose resistance in BRSMG Realce, map the resistance locus or major gene cluster previously named as Co-Realce, identify resistance-related positional genes, and analyze potential markers linked to the resistance allele. F₂ plants derived from the cross BRSMG Realce × BRS FC104 (Mesoamerican) and from the cross BRSMG Realce × BRS Notável (Mesoamerican) were inoculated with the C. lindemuthianum races 475 and 81, respectively. The BRSMG Realce × BRS FC104 F₂ population was also genotyped using the DArTseq technology. Crosses between BRSMG Realce and BAT 93 (Mesoamerican) were also conducted and resulting F₂ plants were inoculated with the C. lindemuthianum races 65 and 1609, individually. The results shown that anthracnose resistance in BRSMG Realce is controlled by a single locus with complete dominance. A genetic map including 1,118 SNP markers was built and shown 78% of the markers mapped at a distances less than 5.0 cM, with a total genetic length of 4,473.4 cM. A major locus (Co-Realce) explaining 54.6% of the phenotypic variation of symptoms caused by the race 475 was identified in Pv04, flanked by the markers snp1327 and snp12782 and 4.48 cM apart each other. These SNPs are useful for marker-assisted selection, due to an estimated selection efficiency of 99.2%. The identified resistance allele segregates independently of the resistance allele Co-3³ (Pv04) present in BAT 93. The mapped genomic region with 704,867 bp comprising 63 putative genes, 44 of which were related to the pathogen-host interaction. Based on all these results and evidence, anthracnose resistance in BRSMG Realce should be considered as monogenic, useful for breeding purpose. It is proposed that locus Co-Realce is unique and be provisionally designated as CoPv04^R until be officially nominated in accordance with the rules established by the Bean Improvement Cooperative Genetics Committee.

Identification of QTL for Stem Traits in Wheat (Triticum aestivum L.)

Lodging in wheat (Triticum aestivum L.) is a complicated phenomenon that is influenced by physiological, genetics, and external factors. It causes a great yield loss and reduces grain quality and mechanical harvesting efficiency. Lodging resistance is contributed by various traits, including increased stem strength. The aim of this study was to map quantitative trait loci (QTL) controlling stem strength-related features (the number of big vascular bundles, stem diameter, stem wall thickness) using a doubled haploid (DH) population derived from a cross between Baiqimai and Neixiang 5. Field experiments were conducted during 2020-2022, and glasshouse experiments were conducted during 2021-2022. Significant genetic variations were observed for all measured traits, and they were all highly heritable. Fifteen QTL for stem strength-related traits were identified on chromosomes 2D, 3A, 3B, 3D, 4B, 5A, 6B, 7A, and 7D, respectively, and 7 QTL for grain yield-related traits were identified on chromosomes 2B, 2D, 3D, 4B, 7A, and 7B, respectively. The superior allele of the major QTL for the number of big vascular bundle (VB) was independent of plant height (PH), making it possible to improve stem strength without a trade-off of PH, thus improving lodging resistance. VB also showed positive correlations with some of the yield components. The result will be useful for molecular marker-assisted selection (MAS) for high stem strength and high yield potential.

Genome properties of key oil palm (Elaeis guineensis Jacq.) breeding populations

A good knowledge of the genome properties of the populations makes it possible to optimize breeding methods, in particular genomic selection (GS). In oil palm (Elaeis guineensis Jacq), the world's main source of vegetable oil, this would provide insight into the promising GS results obtained so far. The present study considered two complex breeding populations, Deli and La Mé, with 943 individuals and 7324 single-nucleotide polymorphisms (SNPs) from genotyping-by-sequencing. Linkage disequilibrium (LD), haplotype sharing, effective size (N_e), and fixation index (F_st) were investigated. A genetic linkage map spanning 1778.52 cM and with a recombination rate of 2.85 cM/Mbp was constructed. The LD at r²=0.3, considered the minimum to get reliable GS results, spanned over 1.05 cM/0.22 Mbp in Deli and 0.9 cM/0.21 Mbp in La Mé. The significant degree of differentiation existing between Deli and La Mé was confirmed by the high F_st value (0.53), the pattern of correlation of SNP heterozygosity and allele frequency among populations, and the decrease of persistence of LD and of haplotype sharing among populations with increasing SNP distance. However, the level of resemblance between the two populations over short genomic distances (correlation of r values between populations >0.6 for SNPs separated by <0.5 cM/1 kbp and percentage of common haplotypes >40% for haplotypes <3600 bp/0.20 cM) likely explains the superiority of GS models ignoring the parental origin of marker alleles over models taking this information into account. The two populations had low N_e (<5). Population-specific genetic maps and reference genomes are recommended for future studies.

Characterization and QTL Mapping of a Major Field Resistance Locus for Bacterial Blight in Rice

Bacterial blight (BB) disease, caused by Xanthomonas oryzae pv. oryzae (Xoo), is among the major factors that can cause rice yields to decrease. To address BB disease, researchers have been looking for ways to change pesticides and cultivation methods, but developing resistant cultivars is the most effective method. However, the resistance and genetic factors of cultivars may be destroyed due to the emergence of new Xoo species caused by recent and rapid climate changes. Therefore, breeders need to identify resistance genes that can be sustained during unpredictable climate changes and utilized for breeding. Here, qBBR11, a quantitative trait locus (QTL) for resistance to BB disease, was detected in KJ (Korea Japonica varieties) 11_067 to KJ11_068 on chromosome 11 in a population derived by crossing JJ (Jeonju) 623 and HR(High resistant)27,195, which possess similar genetic backgrounds but different degrees of resistance to BB disease. qBBR11 was reduced from 18.49–18.69 Mbp of chromosome 11 to 200 kbp segment franked. In this region, 16 candidate genes were detected, and we identified 24 moderate-impact variations and four high-impact variations. In particular, high-impact variations were detected in Os11g0517800 which encode the domain region of GCN2 which is the eIF-2-alpha kinase associated with the resistance of abiotic/biotic stress in rice. In JJ623, which is moderately resistant to BB disease, a stop codon was created due to single nucleotide polymorphism (SNP). Therefore, compared with HR27195, JJ623 has weaker resistance to BB disease, though the two have similar genetic backgrounds. The results suggest that variation in the qBBR11 region regulates an important role in improving resistance to BB diseases, and qBBR11 is useful in providing an important resource for marker-assisted selection to improve mechanisms of resistance to BB disease.

Five plants per RIL for phenotyping traits of high or moderate heritability ensure the power of QTL mapping in a rice MAGIC population

Currently, the power of QTL mapping is mainly dependent on the quality of phenotypic data in a given population, regardless of the statistical method, as the quality of genotypic data is easily guaranteed in the laboratory. Increasing the sample size per line used for phenotyping is a good way to improve the quality of phenotypic data. However, accommodating a large-scale mapping population takes a large area of rice field, which frequently results in high costs and extra environmental noises. To acquire a reasonable small sample size without a penalty in mapping power, we conducted three experiments with a 4-way MAGIC population and measured phenotypes of 5, 10, and 20 plants per RIL. Three traits including heading date, plant height, and tillers per plant were focused. With SNP- and bin-based QTL mapping, 3 major and 3 minor QTLs for heading date with high heritability and 2 major QTLs for plant height with moderate heritability were commonly detected across the three experiments, but no QTL for tillers per plant with low heritability were commonly identified. In addition, bin-based QTL mapping was more powerful than SNP-based mapping and able to rank the genetic effects of parental alleles. Thus, 5 plants per RIL for phenotyping ensure the power of QTL mapping for traits of high or moderate heritability, and bin-based QTL mapping is recommended for multiparent populations.

Construction of Consensus Genetic Map With Applications in Gene Mapping of Wheat (Triticum aestivum L.) Using 90K SNP Array

Wheat is one of the most important cereal crops worldwide. A consensus map combines genetic information from multiple populations, providing an effective alternative to improve the genome coverage and marker density. In this study, we constructed a consensus map from three populations of recombinant inbred lines (RILs) of wheat using a 90K single nucleotide polymorphism (SNP) array. Phenotypic data on plant height (PH), spike length (SL), and thousand-kernel weight (TKW) was collected in six, four, and four environments in the three populations, and then used for quantitative trait locus (QTL) mapping. The mapping results obtained using the constructed consensus map were compared with previous results obtained using individual maps and previous studies on other populations. A simulation experiment was also conducted to assess the performance of QTL mapping with the consensus map. The constructed consensus map from the three populations spanned 4558.55 cM in length, with 25,667 SNPs, having high collinearity with physical map and individual maps. Based on the consensus map, 21, 27, and 19 stable QTLs were identified for PH, SL, and TKW, much more than those detected with individual maps. Four PH QTLs and six SL QTLs were likely to be novel. A putative gene called TraesCS4D02G076400 encoding gibberellin-regulated protein was identified to be the candidate gene for one major PH QTL located on 4DS, which may enrich genetic resources in wheat semi-dwarfing breeding. The simulation results indicated that the length of the confidence interval and standard errors of the QTLs detected using the consensus map were much smaller than those detected using individual maps. The consensus map constructed in this study provides the underlying genetic information for systematic mapping, comparison, and clustering of QTL, and gene discovery in wheat genetic study. The QTLs detected in this study had stable effects across environments and can be used to improve the wide adaptation of wheat cultivars through marker-assisted breeding.

GBS-SNP and SSR based genetic mapping and QTL analysis for drought tolerance in upland cotton

A recombinant inbred line mapping population of intra-species upland cotton was generated from a cross between the drought-tolerant female parent (AS2) and the susceptible male parent (MCU13). A linkage map was constructed deploying 1,116 GBS-based SNPs and public domain-based 782 SSRs spanning a total genetic distance of 28,083.03 cM with an average chromosomal span length of 1,080.12 cM with inter-marker distance of 10.19 cM.A total of 19 quantitative trait loci (QTLs) were identified in nine chromosomes for field drought tolerance traits. Chromosomes 3 and 8 harbored important drought tolerant QTLs for chlorophyll stability index trait while for relative water content trait, three QTLs on chromosome 8 and one QTL each on chromosome 4, 12 were identified. One QTL on each chromosome 8, 5, and 7, and two QTLs on chromosome 15 linking to proline content were identified. For the nitrate reductase activity trait, two QTLs were identified on chromosome 3 and one on each chromosome 8, 13, and 26. To complement our QTL study, a meta-analysis was conducted along with the public domain database and resulted in a consensus map for chromosome 8. Under field drought stress, chromosome 8 harbored a drought tolerance QTL hotspot with two in-house QTLs for chlorophyll stability index (qCSI01, qCSI02) and three public domain QTLs (qLP.FDT_1, qLP.FDT_2, qCC.ST_3). Identified QTL hotspot on chromosome 8 could play a crucial role in exploring abiotic stress-associated genes/alleles for drought trait improvement.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01041-y.

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.

QTL Analysis of Five Morpho-Physiological Traits in Bread Wheat Using Two Mapping Populations Derived from Common Parents

Traits such as plant height (PH), juvenile growth habit (GH), heading date (HD), and tiller number are important for both increasing yield potential and improving crop adaptation to climate change. In the present study, these traits were investigated by using the same bi-parental population at early (F₂ and F₂-derived F₃ families) and late (F₆ and F₇, recombinant inbred lines, RILs) generations to detect quantitative trait loci (QTLs) and search for candidate genes. A total of 176 and 178 lines were genotyped by the wheat Illumina 25K Infinium SNP array. The two genetic maps spanned 2486.97 cM and 3732.84 cM in length, for the F₂ and RILs, respectively. QTLs explaining the highest phenotypic variation were found on chromosomes 2B, 2D, 5A, and 7D for HD and GH, whereas those for PH were found on chromosomes 4B and 4D. Several QTL detected in the early generations (i.e., PH and tiller number) were not detected in the late generations as they were due to dominance effects. Some of the identified QTLs co-mapped to well-known adaptive genes (i.e., Ppd-1, Vrn-1, and Rht-1). Other putative candidate genes were identified for each trait, of which PINE1 and PIF4 may be considered new for GH and TTN in wheat. The use of a large F₂ mapping population combined with NGS-based genotyping techniques could improve map resolution and allow closer QTL tagging.

Quantile regression in genomic selection for oligogenic traits in autogamous plants: A simulation study

This study assessed the efficiency of Genomic selection (GS) or genome-wide selection (GWS), based on Regularized Quantile Regression (RQR), in the selection of genotypes to breed autogamous plant populations with oligogenic traits. To this end, simulated data of an F2 population were used, with traits with different heritability levels (0.10, 0.20 and 0.40), controlled by four genes. The generations were advanced (up to F6) at two selection intensities (10% and 20%). The genomic genetic value was computed by RQR for different quantiles (0.10, 0.50 and 0.90), and by the traditional GWS methods, specifically RR-BLUP and BLASSO. A second objective was to find the statistical methodology that allows the fastest fixation of favorable alleles. In general, the results of the RQR model were better than or equal to those of traditional GWS methodologies, achieving the fixation of favorable alleles in most of the evaluated scenarios. At a heritability level of 0.40 and a selection intensity of 10%, RQR (0.50) was the only methodology that fixed the alleles quickly, i.e., in the fourth generation. Thus, it was concluded that the application of RQR in plant breeding, to simulated autogamous plant populations with oligogenic traits, could reduce time and consequently costs, due to the reduction of selfing generations to fix alleles in the evaluated scenarios.

Identification of a quantitative trait loci (QTL) associated with ammonia tolerance in the Pacific white shrimp (Litopenaeus vannamei)

Background

Ammonia is one of the most common toxicological environment factors affecting shrimp health. Although ammonia tolerance in shrimp is closely related to successful industrial production, few genetic studies of this trait are available.

Results

In this study, we constructed a high-density genetic map of the Pacific white shrimp (Litopenaeus vannamei) using specific length amplified fragment sequencing (SLAF-seq). The constructed genetic map contained 17,338 polymorphic markers spanning 44 linkage groups, with a total distance of 6360.12 centimorgans (cM) and an average distance of 0.37 cM. Using this genetic map, we identified a quantitative trait locus (QTL) that explained 7.41–8.46% of the phenotypic variance in L. vannamei survival time under acute ammonia stress. We then sequenced the transcriptomes of the most ammonia-tolerant and the most ammonia-sensitive individuals from each of four genetically distinct L. vannamei families. We found that 7546 genes were differentially expressed between the ammonia-tolerant and ammonia-sensitive individuals. Using QTL analysis and the transcriptomes, we identified one candidate gene (annotated as an ATP synthase g subunit) associated with ammonia tolerance.

Conclusions

In this study, we constructed a high-density genetic map of L. vannamei and identified a QTL for ammonia tolerance. By combining QTL and transcriptome analyses, we identified a candidate gene associated with ammonia tolerance. Our work provides the basis for future genetic studies focused on molecular marker-assisted selective breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07254-x.

How Machine Learning Methods Helped Find Putative Rye Wax Genes Among GBS Data

The standard approach to genetic mapping was supplemented by machine learning (ML) to establish the location of the rye gene associated with epicuticular wax formation (glaucous phenotype). Over 180 plants of the biparental F2 population were genotyped with the DArTseq (sequencing-based diversity array technology). A maximum likelihood (MLH) algorithm (JoinMap 5.0) and three ML algorithms: logistic regression (LR), random forest and extreme gradient boosted trees (XGBoost), were used to select markers closely linked to the gene encoding wax layer. The allele conditioning the nonglaucous appearance of plants, derived from the cultivar Karlikovaja Zelenostebelnaja, was mapped at the chromosome 2R, which is the first report on this localization. The DNA sequence of DArT-Silico 3585843, closely linked to wax segregation detected by using ML methods, was indicated as one of the candidates controlling the studied trait. The putative gene encodes the ABCG11 transporter.

Identifying a novel powdery mildew resistance gene in a barley landrace from Morocco

Powdery mildew is a barley foliar disease that causes great loss in yield. Because of the limited number of effective resistance genes, efforts to identify new sources of resistance are frequently focused on genetically diversified landraces. The goal of this study was to characterise the powdery mildew resistance gene in barley line 2553-3 selected from the Moroccan landrace. Phytopathological testing against a set of differential pathogen isolates revealed different pattern responses of this gene from those of other known resistance genes. F₂ and F_2:3 (2553-3 × Manchuria) mapping populations were employed to investigate resistance inheritance. Two approaches were applied for the linkage analysis: in the first approach, 22 resistant and 21 susceptible homozygous F₂ plants genotyped by the DArTseq platform (Diversity Arrays Technology, Pty. Ltd.) were used; in the second, 94 F₂ plants were genotyped by converted DArTseq markers and SSRs. Both analyses delineated a new resistance gene on the short arm of chromosome 2H. The authors propose MlMor as a gene symbol for newly characterized powdery mildew resistance genes in barley line 255-3-3. The results presented herein provide a good foundation for the development of closer linkage markers and MAS breeding.

Construction of a genetic linkage map in Pyropia yezoensis (Bangiales, Rhodophyta) and QTL analysis of several economic traits of blades

Pyropia yezoensis is an economically important seaweed but its molecular genetics is poorly understood. In the present study, we used a doubled haploid (DH) population that was established in our previous work to construct a genetic linkage map of P. yezoensis and analyze the quantitative trait loci (QTLs) of blades. The DH population was genotyped with fluorescent sequence-related amplified polymorphism (SRAP) markers. A chi-square test identified 301 loci with normal segregation (P ≥ 0.01) and 96 loci (24.18%) with low-level skewed segregation (0.001 ≤ P < 0.01). The genetic map was constructed after a total of 92 loci were assembled into three linkage groups (LGs). The map spanned 557.36 cM covering 93.71% of the estimated genome, with a mean interlocus space of 6.23 cM. Kolmogorov-Smirnov test (α = 5%) showed a uniform distribution of the markers along each LG. On the genetic map, 10 QTLs associated with five economic traits of blades were detected. One QTL was for length, one for width, two for fresh weight, two for specific growth rate of length and four for specific growth rate of fresh weight. These QTLs could explain 2.29–7.87% of the trait variations, indicating that their effects were all minor. The results may serve as a framework for future marker-assisted breeding in P. yezoensis.

High-density SNP linkage map construction and QTL mapping for flavonoid-related traits in a tea plant (Camellia sinensis) using 2b-RAD sequencing

BACKGROUND

Flavonoids are important components that confer upon tea plants a unique flavour and health functions. However, the traditional breeding method for selecting a cultivar with a high or unique flavonoid content is time consuming and labour intensive. High-density genetic map construction associated with quantitative trait locus (QTL) mapping provides an effective way to facilitate trait improvement in plant breeding. In this study, an F1 population (LJ43×BHZ) was genotyped using 2b-restriction site-associated DNA (2b-RAD) sequencing to obtain massive single nucleotide polymorphism (SNP) markers to construct a high-density genetic map for a tea plant. Furthermore, QTLs related to flavonoids were identified using our new genetic map.

RESULTS

A total of 13,446 polymorphic SNP markers were developed using 2b-RAD sequencing, and 4,463 of these markers were available for constructing the genetic linkage map. A 1,678.52-cM high-density map at an average interval of 0.40 cM with 4,217 markers, including 427 frameset simple sequence repeats (SSRs) and 3,800 novel SNPs, mapped into 15 linkage groups was successfully constructed. After QTL analysis, a total of 27 QTLs related to flavonoids or caffeine content (CAF) were mapped to 8 different linkage groups, LG01, LG03, LG06, LG08, LG10, LG11, LG12, and LG13, with an LOD from 3.14 to 39.54, constituting 7.5% to 42.8% of the phenotypic variation.

CONCLUSIONS

To our knowledge, the highest density genetic map ever reported was constructed since the largest mapping population of tea plants was adopted in present study. Moreover, novel QTLs related to flavonoids and CAF were identified based on the new high-density genetic map. In addition, two markers were located in candidate genes that may be involved in flavonoid metabolism. The present study provides valuable information for gene discovery, marker-assisted selection breeding and map-based cloning for functional genes that are related to flavonoid content in tea plants.

Quantitative Trait Loci Mapping for Theobromine and Caffeine Contents in Tea Plant ( Camellia sinensis)

Understanding the genetic basis of theobromine and caffeine accumulation in the tea plant is important due to their contribution to tea flavor. Quantitative trait loci (QTL) analyses were carried out to identify genetic variants associated with theobromine and caffeine contents and ratio using a pseudo-testcross population derived from an intervarietal cross between two varieties of Camellia sinensis. A total of 10 QTL controlling caffeine content (CAF), theobromine content (TBR), sum of caffeine and theobromine (SCT), and caffeine-to-theobromine ratio (CTR) were identified over four measurement years. The major QTL controlling CAF, qCAF1, was mapped onto LG01 and validated across years, explaining an average of 20.1% of the phenotypic variance. The other QTL were detected in 1 or 2 years, and of them there were four, two, and three for TBR, SCT, and CTR, respectively. The present results provide valuable information for further fine mapping and cloning functional genes and for genetic improvement in tea plant.

QTL mapping of melon fruit quality traits using a high-density GBS-based genetic map

BACKGROUND

Melon shows a broad diversity in fruit morphology and quality, which is still underexploited in breeding programs. The knowledge of the genetic basis of fruit quality traits is important for identifying new alleles that may be introduced in elite material by highly efficient molecular breeding tools.

RESULTS

In order to identify QTLs controlling fruit quality, a recombinant inbred line population was developed using two commercial cultivars as parental lines: "Védrantais", from the cantalupensis group, and "Piel de Sapo", from the inodorus group. Both have desirable quality traits for the market, but their fruits differ in traits such as rind and flesh color, sugar content, ripening behavior, size and shape. We used a genotyping-by-sequencing strategy to construct a dense genetic map, which included around five thousand variants distributed in 824 bins. The RIL population was phenotyped for quality and morphology traits, and we mapped 33 stable QTLs involved in sugar and carotenoid content, fruit and seed morphology and major loci controlling external color of immature fruit and mottled rind. The median confidence interval of the QTLs was 942 kb, suggesting that the high density of the genetic map helped in increasing the mapping resolution. Some of these intervals contained less than a hundred annotated genes, and an integrative strategy combining gene expression and resequencing data enabled identification of candidate genes for some of these traits.

CONCLUSION

Several QTLs controlling fruit quality traits in melon were identified and delimited to narrow genomic intervals, using a RIL population and a GBS-based genetic map.

Linkage mapping combined with association analysis reveals QTL and candidate genes for three husk traits in maize

Key message Combined linkage and association mapping analyses facilitate the emphasis on the candidate genes putatively involved in maize husk growth. The maize (Zea mays L.) husk consists of multiple leafy layers and plays important roles in protecting the ear from pathogen infection and in preventing grain dehydration. Although husk morphology varies widely among different maize inbred lines, the genetic basis of such variation is poorly understood. In this study, we used three maize recombinant inbred line (RIL) populations to dissect the genetic basis of three husk traits: i.e., husk length (HL), husk width (HW), and the number of husk layers (HN). Three husk traits in all three RIL populations showed wide phenotypic variation and high heritability. The HL showed stronger correlations with ear traits than did HW and HN. A total of 21 quantitative trait loci (QTL) were identified for the three traits in three RIL populations, and some of them were commonly observed for the same trait in different populations. The proportions of total phenotypic variation explained by QTL in three RIL populations were 31.8, 35.3, and 44.5% for HL, HW, and HN, respectively. The highest proportions of phenotypic variation explained by a single QTL were 14.7% for HL in the By815/K22 RIL population (BYK), 13.5% for HW in the By815/DE3 RIL population (BYD), and 19.4% for HN in the BYD population. A combined analysis of linkage mapping with a previous genome-wide association study revealed five candidate genes related to husk morphology situated within three QTL loci. These five genes were related to metabolism, gene expression regulation, and signal transduction.

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.

Improving Nelumbo nucifera genome assemblies using high-resolution genetic maps and BioNano genome mapping reveals ancient chromosome rearrangements

Genetic and physical maps are powerful tools to anchor fragmented draft genome assemblies generated from next-generation sequencing. Currently, two draft assemblies of Nelumbo nucifera, the genomes of 'China Antique' and 'Chinese Tai-zi', have been released. However, there is presently no information on how the sequences are assembled into chromosomes in N. nucifera. The lack of physical maps and inadequate resolution of available genetic maps hindered the assembly of N. nucifera chromosomes. Here, a linkage map of N. nucifera containing 2371 bin markers [217 577 single nucleotide polymorphisms (SNPs)] was constructed using restriction-site associated DNA sequencing data of 181 F₂ individuals and validated by adding 197 simple sequence repeat (SSR) markers. Additionally, a BioNano optical map covering 86.20% of the 'Chinese Tai-zi' genome was constructed. The draft assembly of 'Chinese Tai-zi' was improved based on the BioNano optical map, showing an increase of the scaffold N50 from 0.989 to 1.48 Mb. Using a combination of multiple maps, 97.9% of the scaffolds in the 'Chinese Tai-zi' draft assembly and 97.6% of the scaffolds in the 'China Antique' draft assembly were anchored into pseudo-chromosomes, and the centromere regions along the pseudo-chromosomes were identified. An evolutionary scenario was proposed to reach the modern N. nucifera karyotype from the seven ancestral eudicot chromosomes. The present study provides the highest-resolution linkage map, the optical map and chromosome level genome assemblies for N. nucifera, which are valuable for the breeding and cultivation of N. nucifera and future studies of comparative and evolutionary genomics in angiosperms.

A consensus linkage map of lentil based on DArT markers from three RIL mapping populations

Background

Lentil (Lens culinaris ssp. culinaris Medikus) is a diploid (2n = 2x = 14), self-pollinating grain legume with a haploid genome size of about 4 Gbp and is grown throughout the world with current annual production of 4.9 million tonnes.

Materials and methods

A consensus map of lentil (Lens culinaris ssp. culinaris Medikus) was constructed using three different lentils recombinant inbred line (RIL) populations, including “CDC Redberry” x “ILL7502” (LR8), “ILL8006” x “CDC Milestone” (LR11) and “PI320937” x “Eston” (LR39).

Results

The lentil consensus map was composed of 9,793 DArT markers, covered a total of 977.47 cM with an average distance of 0.10 cM between adjacent markers and constructed 7 linkage groups representing 7 chromosomes of the lentil genome. The consensus map had no gap larger than 12.67 cM and only 5 gaps were found to be between 12.67 cM and 6.0 cM (on LG3 and LG4). The localization of the SNP markers on the lentil consensus map were in general consistent with their localization on the three individual genetic linkage maps and the lentil consensus map has longer map length, higher marker density and shorter average distance between the adjacent markers compared to the component linkage maps.

Conclusion

This high-density consensus map could provide insight into the lentil genome. The consensus map could also help to construct a physical map using a Bacterial Artificial Chromosome library and map based cloning studies. Sequence information of DArT may help localization of orientation scaffolds from Next Generation Sequencing data.

High throughput SNP discovery and genotyping in hexaploid wheat

Because of their abundance and their amenability to high-throughput genotyping techniques, Single Nucleotide Polymorphisms (SNPs) are powerful tools for efficient genetics and genomics studies, including characterization of genetic resources, genome-wide association studies and genomic selection. In wheat, most of the previous SNP discovery initiatives targeted the coding fraction, leaving almost 98% of the wheat genome largely unexploited. Here we report on the use of whole-genome resequencing data from eight wheat lines to mine for SNPs in the genic, the repetitive and non-repetitive intergenic fractions of the wheat genome. Eventually, we identified 3.3 million SNPs, 49% being located on the B-genome, 41% on the A-genome and 10% on the D-genome. We also describe the development of the TaBW280K high-throughput genotyping array containing 280,226 SNPs. Performance of this chip was examined by genotyping a set of 96 wheat accessions representing the worldwide diversity. Sixty-nine percent of the SNPs can be efficiently scored, half of them showing a diploid-like clustering. The TaBW280K was proven to be a very efficient tool for diversity analyses, as well as for breeding as it can discriminate between closely related elite varieties. Finally, the TaBW280K array was used to genotype a population derived from a cross between Chinese Spring and Renan, leading to the construction a dense genetic map comprising 83,721 markers. The results described here will provide the wheat community with powerful tools for both basic and applied research.

Utilization of a Wheat660K SNP array-derived high-density genetic map for high-resolution mapping of a major QTL for kernel number

In crop plants, a high-density genetic linkage map is essential for both genetic and genomic researches. The complexity and the large size of wheat genome have hampered the acquisition of a high-resolution genetic map. In this study, we report a high-density genetic map based on an individual mapping population using the Affymetrix Wheat660K single-nucleotide polymorphism (SNP) array as a probe in hexaploid wheat. The resultant genetic map consisted of 119 566 loci spanning 4424.4 cM, and 119 001 of those loci were SNP markers. This genetic map showed good collinearity with the 90 K and 820 K consensus genetic maps and was also in accordance with the recently released wheat whole genome assembly. The high-density wheat genetic map will provide a major resource for future genetic and genomic research in wheat. Moreover, a comparative genomics analysis among gramineous plant genomes was conducted based on the high-density wheat genetic map, providing an overview of the structural relationships among theses gramineous plant genomes. A major stable quantitative trait locus (QTL) for kernel number per spike was characterized, providing a solid foundation for the future high-resolution mapping and map-based cloning of the targeted QTL.

Utilization of a Wheat660K SNP array-derived high-density genetic map for high-resolution mapping of a major QTL for kernel number

In crop plants, a high-density genetic linkage map is essential for both genetic and genomic researches. The complexity and the large size of wheat genome have hampered the acquisition of a high-resolution genetic map. In this study, we report a high-density genetic map based on an individual mapping population using the Affymetrix Wheat660K single-nucleotide polymorphism (SNP) array as a probe in hexaploid wheat. The resultant genetic map consisted of 119 566 loci spanning 4424.4 cM, and 119 001 of those loci were SNP markers. This genetic map showed good collinearity with the 90 K and 820 K consensus genetic maps and was also in accordance with the recently released wheat whole genome assembly. The high-density wheat genetic map will provide a major resource for future genetic and genomic research in wheat. Moreover, a comparative genomics analysis among gramineous plant genomes was conducted based on the high-density wheat genetic map, providing an overview of the structural relationships among theses gramineous plant genomes. A major stable quantitative trait locus (QTL) for kernel number per spike was characterized, providing a solid foundation for the future high-resolution mapping and map-based cloning of the targeted QTL.

Construction of a high-density genetic map and QTL mapping for pearl quality-related traits in Hyriopsis cumingii

A high-density genetic map is essential for quantitative trait locus (QTL) fine mapping. In this study, 4,508 effective single nucleotide polymorphism markers (detected using specific-locus amplified fragment sequencing) and 475 microsatellites were mapped to 19 linkage groups (LGs) using a family with 157 individuals. The map spanned 2,713 cM, with an average of 259 markers and 79 loci per LG and an average inter-marker distance of 1.81 cM. To identify QTLs for pearl quality traits, 26 putatively significant QTLs were detected for 10 traits, including, three for shell width, seven for body weight, two for shell weight, two for margin mantle weight, five for inner mantle weight, and seven for shell nacre colour. Among them, five QTLs associated with shell nacre colour were mapped to LG17 and explained 19.7% to 22.8% of the trait variation; this suggests that some important genes or loci determine shell nacre colour in LG17. The linkage map and mapped QTLs for shell nacre colour would be useful for improving the quality of Hyriopsis cumingii via marker-assisted selection.

Exploiting genotyping by sequencing to characterize the genomic structure of the American cranberry through high-density linkage mapping

Background

The application of genotyping by sequencing (GBS) approaches, combined with data imputation methodologies, is narrowing the genetic knowledge gap between major and understudied, minor crops. GBS is an excellent tool to characterize the genomic structure of recently domesticated (~200 years) and understudied species, such as cranberry (Vaccinium macrocarpon Ait.), by generating large numbers of markers for genomic studies such as genetic mapping.

Results

We identified 10842 potentially mappable single nucleotide polymorphisms (SNPs) in a cranberry pseudo-testcross population wherein 5477 SNPs and 211 short sequence repeats (SSRs) were used to construct a high density linkage map in cranberry of which a total of 4849 markers were mapped. Recombination frequency, linkage disequilibrium (LD), and segregation distortion at the genomic level in the parental and integrated linkage maps were characterized for first time in cranberry. SSR markers, used as the backbone in the map, revealed high collinearity with previously published linkage maps. The 4849 point map consisted of twelve linkage groups spanning 1112 cM, which anchored 2381 nuclear scaffolds accounting for ~13 Mb of the estimated 470 Mb cranberry genome. Bin mapping identified 592 and 672 unique bins in the parentals and a total of 1676 unique marker positions in the integrated map. Synteny analyses comparing the order of anchored cranberry scaffolds to their homologous positions in kiwifruit, grape, and coffee genomes provided initial evidence of homology between cranberry and closely related species.

Conclusions

GBS data was used to rapidly saturate the cranberry genome with markers in a pseudo-testcross population. Collinearity between the present saturated genetic map and previous cranberry SSR maps suggests that the SNP locations represent accurate marker order and chromosome structure of the cranberry genome. SNPs greatly improved current marker genome coverage, which allowed for genome-wide structure investigations such as segregation distortion, recombination, linkage disequilibrium, and synteny analyses. In the future, GBS can be used to accelerate cranberry molecular breeding through QTL mapping and genome-wide association studies (GWAS).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2802-3) contains supplementary material, which is available to authorized users.

Construction of high resolution genetic linkage maps to improve the soybean genome sequence assembly Glyma1.01

BACKGROUND

A landmark in soybean research, Glyma1.01, the first whole genome sequence of variety Williams 82 (Glycine max L. Merr.) was completed in 2010 and is widely used. However, because the assembly was primarily built based on the linkage maps constructed with a limited number of markers and recombinant inbred lines (RILs), the assembled sequence, especially in some genomic regions with sparse numbers of anchoring markers, needs to be improved. Molecular markers are being used by researchers in the soybean community, however, with the updating of the Glyma1.01 build based on the high-resolution linkage maps resulting from this research, the genome positions of these markers need to be mapped.

RESULTS

Two high density genetic linkage maps were constructed based on 21,478 single nucleotide polymorphism loci mapped in the Williams 82 x G. soja (Sieb. & Zucc.) PI479752 population with 1083 RILs and 11,922 loci mapped in the Essex x Williams 82 population with 922 RILs. There were 37 regions or single markers where marker order in the two populations was in agreement but was not consistent with the physical position in the Glyma1.01 build. In addition, 28 previously unanchored scaffolds were positioned. Map data were used to identify false joins in the Glyma1.01 assembly and the corresponding scaffolds were broken and reassembled to the new assembly, Wm82.a2.v1. Based upon the plots of the genetic on physical distance of the loci, the euchromatic and heterochromatic regions along each chromosome in the new assembly were delimited. Genomic positions of the commonly used markers contained in BARCSOYSSR_1.0 database and the SoySNP50K BeadChip were updated based upon the Wm82.a2.v1 assembly.

CONCLUSIONS

The information will facilitate the study of recombination hot spots in the soybean genome, identification of genes or quantitative trait loci controlling yield, seed quality and resistance to biotic or abiotic stresses as well as other genetic or genomic research.

DArT Markers Effectively Target Gene Space in the Rye Genome

Large genome size and complexity hamper considerably the genomics research in relevant species. Rye (Secale cereale L.) has one of the largest genomes among cereal crops and repetitive sequences account for over 90% of its length. Diversity Arrays Technology is a high-throughput genotyping method, in which a preferential sampling of gene-rich regions is achieved through the use of methylation sensitive restriction enzymes. We obtained sequences of 6,177 rye DArT markers and following a redundancy analysis assembled them into 3,737 non-redundant sequences, which were then used in homology searches against five Pooideae sequence sets. In total 515 DArT sequences could be incorporated into publicly available rye genome zippers providing a starting point for the integration of DArT- and transcript-based genomics resources in rye. Using Blast2Go pipeline we attributed putative gene functions to 1101 (29.4%) of the non-redundant DArT marker sequences, including 132 sequences with putative disease resistance-related functions, which were found to be preferentially located in the 4RL and 6RL chromosomes. Comparative analysis based on the DArT sequences revealed obvious inconsistencies between two recently published high density consensus maps of rye. Furthermore we demonstrated that DArT marker sequences can be a source of SSR polymorphisms. Obtained data demonstrate that DArT markers effectively target gene space in the large, complex, and repetitive rye genome. Through the annotation of putative gene functions and the alignment of DArT sequences relative to reference genomes we obtained information, that will complement the results of the studies, where DArT genotyping was deployed, by simplifying the gene ontology and microcolinearity based identification of candidate genes.

Large-Scale SNP Discovery and Genotyping for Constructing a High-Density Genetic Map of Tea Plant Using Specific-Locus Amplified Fragment Sequencing (SLAF-seq)

Genetic maps are important tools in plant genomics and breeding. The present study reports the large-scale discovery of single nucleotide polymorphisms (SNPs) for genetic map construction in tea plant. We developed a total of 6,042 valid SNP markers using specific-locus amplified fragment sequencing (SLAF-seq), and subsequently mapped them into the previous framework map. The final map contained 6,448 molecular markers, distributing on fifteen linkage groups corresponding to the number of tea plant chromosomes. The total map length was 3,965 cM, with an average inter-locus distance of 1.0 cM. This map is the first SNP-based reference map of tea plant, as well as the most saturated one developed to date. The SNP markers and map resources generated in this study provide a wealth of genetic information that can serve as a foundation for downstream genetic analyses, such as the fine mapping of quantitative trait loci (QTL), map-based cloning, marker-assisted selection, and anchoring of scaffolds to facilitate the process of whole genome sequencing projects for tea plant.

High-resolution genetic maps of Eucalyptus improve Eucalyptus grandis genome assembly

Genetic maps are key tools in genetic research as they constitute the framework for many applications, such as quantitative trait locus analysis, and support the assembly of genome sequences. The resequencing of the two parents of a cross between Eucalyptus urophylla and Eucalyptus grandis was used to design a single nucleotide polymorphism (SNP) array of 6000 markers evenly distributed along the E. grandis genome. The genotyping of 1025 offspring enabled the construction of two high-resolution genetic maps containing 1832 and 1773 markers with an average marker interval of 0.45 and 0.5 cM for E. grandis and E. urophylla, respectively. The comparison between genetic maps and the reference genome highlighted 85% of collinear regions. A total of 43 noncollinear regions and 13 nonsynthetic regions were detected and corrected in the new genome assembly. This improved version contains 4943 scaffolds totalling 691.3 Mb of which 88.6% were captured by the 11 chromosomes. The mapping data were also used to investigate the effect of population size and number of markers on linkage mapping accuracy. This study provides the most reliable linkage maps for Eucalyptus and version 2.0 of the E. grandis genome.

Chromosomal evolution and patterns of introgression in helianthus

Knowledge of the nature and extent of karyotypic differences between species provides insight into the evolutionary history of the genomes in question and, in the case of closely related species, the potential for genetic exchange between taxa. We constructed high-density genetic maps of the silverleaf sunflower (Helianthus argophyllus) and Algodones Dune sunflower (H. niveus ssp. tephrodes) genomes and compared them to a consensus map of cultivated sunflower (H. annuus) to identify chromosomal rearrangements between species. The genetic maps of H. argophyllus and H. niveus ssp. tephrodes included 17 linkage groups each and spanned 1337 and 1478 cM, respectively. Comparative analyses revealed greater divergence between H. annuus and H. niveus ssp. tephrodes (13 inverted segments, 18 translocated segments) than between H. annuus and H. argophyllus (10 inverted segments, 8 translocated segments), consistent with their known phylogenetic relationships. Marker order was conserved across much of the genome, with 83 and 64% of the H. argophyllus and H. niveus ssp. tephrodes genomes, respectively, being syntenic with H. annuus. Population genomic analyses between H. annuus and H. argophyllus, which are sympatric across a portion of the natural range of H. annuus, revealed significantly elevated genetic structure in rearranged portions of the genome, indicating that such rearrangements are associated with restricted gene flow between these two species.

Construction of a SSR-based genetic map and identification of QTLs for catechins content in tea plant (Camellia sinensis)

Catechins are the most important bioactive compounds in tea, and have been demonstrated to possess a wide variety of pharmacological activities. To characterize quantitative trait loci (QTLs) for catechins content in the tender shoots of tea plant, we constructed a moderately saturated genetic map using 406 simple sequence repeat (SSR) markers, based on a pseudo-testcross population of 183 individuals derived from an intraspecific cross of two Camellia sinensis varieties with diverse catechins composition. The map consisted of fifteen linkage groups (LGs), corresponding to the haploid chromosome number of tea plant (2n = 2x = 30). The total map length was 1,143.5 cM, with an average locus spacing of 2.9 cM. A total of 25 QTLs associated with catechins content were identified over two measurement years. Of these, nine stable QTLs were validated across years, and clustered into four main chromosome regions on LG03, LG11, LG12 and LG15. The population variability explained by each QTL was predominantly at moderate-to-high levels and ranged from 2.4% to 71.0%, with an average of 17.7%. The total number of QTL for each trait varied from four to eight, while the total population variability explained by all QTLs for a trait ranged between 38.4% and 79.7%. This is the first report on the identification of QTL for catechins content in tea plant. The results of this study provide a foundation for further cloning and functional characterization of catechin QTLs for utilization in improvement of tea plant.

AFLP-based genetic mapping of the "bud-flowering" trait in heather (Calluna vulgaris)

BACKGROUND

Calluna vulgaris is one of the most important landscaping plants produced in Germany. Its enormous economic success is due to the prolonged flower attractiveness of mutants in flower morphology, the so-called bud-bloomers. In this study, we present the first genetic linkage map of C. vulgaris in which we mapped a locus of the economically highly desired trait "flower type".

RESULTS

The map was constructed in JoinMap 4.1. using 535 AFLP markers from a single mapping population. A large fraction (40%) of markers showed distorted segregation. To test the effect of segregation distortion on linkage estimation, these markers were sorted regarding their segregation ratio and added in groups to the data set. The plausibility of group formation was evaluated by comparison of the "two-way pseudo-testcross" and the "integrated" mapping approach. Furthermore, regression mapping was compared to the multipoint-likelihood algorithm. The majority of maps constructed by different combinations of these methods consisted of eight linkage groups corresponding to the chromosome number of C. vulgaris.

CONCLUSIONS

All maps confirmed the independent inheritance of the most important horticultural traits "flower type", "flower colour", and "leaf colour". An AFLP marker for the most important breeding target "flower type" was identified. The presented genetic map of C. vulgaris can now serve as a basis for further molecular marker selection and map-based cloning of the candidate gene encoding the unique flower architecture of C. vulgaris bud-bloomers.

Genetic linkage mapping in fungi: current state, applications, and future trends

Genetic mapping is a basic tool for eukaryotic genomic research. Linkage maps provide insights into genome organization and can be used for genetic studies of traits of interest. A genetic linkage map is a suitable support for the anchoring of whole genome sequences. It allows the localization of genes of interest or quantitative trait loci (QTL) and map-based cloning. While genetic mapping has been extensively used in plant or animal models, this discipline is more recent in fungi. The present article reviews the current status of genetic linkage map research in fungal species. The process of linkage mapping is detailed, from the development of mapping populations to the construction of the final linkage map, and illustrated based on practical examples. The range of specific applications in fungi is browsed, such as the mapping of virulence genes in pathogenic species or the mapping of agronomically relevant QTL in cultivated edible mushrooms. Future prospects are finally discussed in the context of the most recent advances in molecular techniques and the release of numerous fungal genome sequences.

A molecular linkage map for Drosophila mediopunctata confirms synteny with Drosophila melanogaster and suggests a region that controls the variation in the number of abdominal spots

The classic approach to gene discovery relies on the construction of linkage maps. We report the first molecular-based linkage map for Drosophila mediopunctata, a neotropical species of the tripunctata group. Eight hundred F(2) individuals were genotyped at 49 microsatellite loci, resulting in a map that is ≈450 centimorgans long. Five linkage groups were detected, and the species' chromosomes were identified through cross-references to BLASTn searches and Müller elements. Strong synteny was observed when compared with the Drosophila melanogaster chromosome arms, but little conservation in the gene order was seen. The incorporation of morphological data corresponding to the number of central abdominal spots on the map was consistent with the expected location of a genomic region responsible for the phenotype on the second chromosome.

A high density consensus map of rye (Secale cereale L.) based on DArT markers

Background

Rye (Secale cereale L.) is an economically important crop, exhibiting unique features such as outstanding resistance to biotic and abiotic stresses and high nutrient use efficiency. This species presents a challenge to geneticists and breeders due to its large genome containing a high proportion of repetitive sequences, self incompatibility, severe inbreeding depression and tissue culture recalcitrance. The genomic resources currently available for rye are underdeveloped in comparison with other crops of similar economic importance. The aim of this study was to create a highly saturated, multilocus linkage map of rye via consensus mapping, based on Diversity Arrays Technology (DArT) markers.

Methodology/Principal Findings

Recombinant inbred lines (RILs) from 5 populations (564 in total) were genotyped using DArT markers and subjected to linkage analysis using Join Map 4.0 and Multipoint Consensus 2.2 software. A consensus map was constructed using a total of 9703 segregating markers. The average chromosome map length ranged from 199.9 cM (2R) to 251.4 cM (4R) and the average map density was 1.1 cM. The integrated map comprised 4048 loci with the number of markers per chromosome ranging from 454 for 7R to 805 for 4R. In comparison with previously published studies on rye, this represents an eight-fold increase in the number of loci placed on a consensus map and a more than two-fold increase in the number of genetically mapped DArT markers.

Conclusions/Significance

Through the careful choice of marker type, mapping populations and the use of software packages implementing powerful algorithms for map order optimization, we produced a valuable resource for rye and triticale genomics and breeding, which provides an excellent starting point for more in-depth studies on rye genome organization.

Genetic map of artichoke × wild cardoon: toward a consensus map for Cynara cardunculus

An integrated consensus linkage map is proposed for globe artichoke. Maternal and paternal genetic maps were constructed on the basis of an F(1) progeny derived from crossing an artichoke genotype (Mola) with its progenitor, the wild cardoon (Tolfa), using EST-derived SSRs, genomic SSRs, AFLPs, ten genes, and two morphological traits. For most genes, mainly belonging to the chlorogenic acid pathway, new markers were developed. Five of these were SNP markers analyzed through high-resolution melt technology. From the maternal (Mola) and paternal (Tolfa) maps, an integrated map was obtained, containing 337 molecular and one morphological markers ordered in 17 linkage groups (LGs), linked between Mola and Tolfa. The integrated map covers 1,488.8 cM, with an average distance of 4.4 cM between markers. The map was aligned with already existing maps for artichoke, and 12 LGs were linked via 31 bridge markers. LG numbering has been proposed. A total of 124 EST-SSRs and two genes were mapped here for the first time, providing a framework for the construction of a functional map in artichoke. The establishment of a consensus map represents a necessary condition to plan a complete sequencing of the globe artichoke genome.

Towards a unified genetic map for diploid roses

We have constructed the first integrated consensus map (ICM) for rose, based on the information of four diploid populations and more than 1,000 initial markers. The single population maps are linked via 59 bridge markers, on average 8.4 per linkage group (LG). The integrated map comprises 597 markers, 206 of which are sequence-based, distributed over a length of 530 cM on seven LGs. By using a larger effective population size and therefore higher marker density, the marker order in the ICM is more reliable than in the single population maps. This is supported by a more even marker distribution and a decrease in gap sizes in the consensus map as compared to the single population maps. This unified map establishes a standard nomenclature for rose LGs, and presents the location of important ornamental traits, such as self-incompatibility, black spot resistance (Rdr1), scent production and recurrent blooming. In total, the consensus map includes locations for 10 phenotypic single loci, QTLs for 7 different traits and 51 ESTs or gene-based molecular markers. This consensus map combines for the first time the information for traits with high relevance for rose variety development. It will serve as a tool for selective breeding and marker assisted selection. It will benefit future efforts of the rose community to sequence the whole rose genome and will be useful for synteny studies in the Rosaceae family and especially in the section Rosoideae.

Fine mapping QMi-C11 a major QTL controlling root-knot nematodes resistance in Upland cotton

The identification and utilization of a high-level of host plant resistance is the most effective and economical approach to control root-knot nematode (Meloidogyne incognita). In an earlier study, we identified a major quantitative trait locus (QTL) for resistance to root-knot nematode in the M-120 RNR Upland cotton line (Gossypium hirsutum L.) of the Auburn 623 RNR source. The QTL is located in a 12.9-cM interval flanked by the two SSR markers CIR069 and CIR316 on the distal segment of chromosome 11. To construct a fine map around the target region, a bulked segregation analysis was performed using two DNA pools consisting of five individuals, with each being homozygous for the two parental alleles. From a survey of 1,152 AFLP primer combinations, 9 AFLP markers closely linked to the target region were identified. By screening an additional 1,221 F(2) individuals developed from the initial mapping population, the Mi-C11 locus was delimited to a 3.6-cM interval flanked by the SSR marker CIR069 and the AFLP marker E14M27-375. These results further elucidate the genetic fine structure of the Mi-C11 locus and provide the basis for map-based isolation of the nematode resistance gene in M-120 RNR.

The first doubled haploid linkage map for cultivated oat

To date, all linkage maps of hexaploid oat (Avena sativa L.) have been constructed using recombinant inbred lines (RILs). Doubled haploids (DHs), however, have the advantage over RILs of their comprehensive homozygosity. DHs have been used for mapping in several cereal species, but in oats the production of large DH populations has only recently become an option. A linkage map of hexaploid oat was constructed using an anther culture-derived DH population (137 individuals) from the F1 individuals of a cross between the Finnish cultivar 'Aslak' and the Swedish cultivar 'Matilda'. The map is composed of 28 linkage groups containing 625 DNA markers: 375 AFLPs (amplified fragment length polymorphisms), 3 IRAPs (inter-retrotransposon amplified polymorphisms), 12 ISSRs (inter simple sequence repeats), 12 microsatellites, 57 RAPDs (random amplified polymorphic DNAs), 59 REMAPs (retrotransposon-microsatellite amplified polymorphisms), 105 SRAPs (sequence-related amplified polymorphisms), and 2 SNPs (single-nucleotide polymorphisms). The total map size is 1526 cM. Over half of the markers in the map showed distorted segregation, with alleles from 'Aslak' usually prevailing. This is explained by the better performance of 'Aslak' in anther culture. Quantitative trait loci affecting some important quality and agronomic traits are being localized on the map.

Predicting the size of the progeny mapping population required to positionally clone a gene

A key frustration during positional gene cloning (map-based cloning) is that the size of the progeny mapping population is difficult to predict, because the meiotic recombination frequency varies along chromosomes. We describe a detailed methodology to improve this prediction using rice (Oryza sativa L.) as a model system. We derived and/or validated, then fine-tuned, equations that estimate the mapping population size by comparing these theoretical estimates to 41 successful positional cloning attempts. We then used each validated equation to test whether neighborhood meiotic recombination frequencies extracted from a reference RFLP map can help researchers predict the mapping population size. We developed a meiotic recombination frequency map (MRFM) for approximately 1400 marker intervals in rice and anchored each published allele onto an interval on this map. We show that neighborhood recombination frequencies (R-map, >280-kb segments) extracted from the MRFM, in conjunction with the validated formulas, better predicted the mapping population size than the genome-wide average recombination frequency (R-avg), with improved results whether the recombination frequency was calculated as genes/cM or kb/cM. Our results offer a detailed road map for better predicting mapping population size in diverse eukaryotes, but useful predictions will require robust recombination frequency maps based on sampling more progeny.