Construction of a high-density composite map and comparative mapping of segregation distortion regions in barley

Identification of novel genes responsible for a pollen killer present in local natural populations of Arabidopsis thaliana

Gamete killers are genetic loci that distort segregation in the progeny of hybrids because the killer allele promotes the elimination of the gametes that carry the sensitive allele. They are widely distributed in eukaryotes and are important for understanding genome evolution and speciation. We had previously identified a pollen killer in hybrids between two distant natural accessions of Arabidopsis thaliana. This pollen killer involves three genetically linked genes, and we previously reported the identification of the gene encoding the antidote that protects pollen grains from the killer activity. In this study, we identified the two other genes of the pollen killer by using CRISPR-Cas9 induced mutants. These two genes are necessary for the killer activity that we demonstrated to be specific to pollen. The cellular localization of the pollen killer encoded proteins suggests that the pollen killer activity involves the mitochondria. Sequence analyses reveal predicted domains from the same families in the killer proteins. In addition, the C-terminal half of one of the killer proteins is identical to the antidote, and one amino acid, crucial for the antidote activity, is also essential for the killer function. Investigating more than 700 worldwide accessions of A. thaliana, we confirmed that the locus is subject to important structural rearrangements and copy number variation. By exploiting available de novo genomic sequences, we propose a scenario for the emergence of this pollen killer in A. thaliana. Furthermore, we report the co-occurrence and behavior of killer and sensitive genotypes in several local populations, a prerequisite for studying gamete killer evolution in the wild. This highlights the potential of the Arabidopsis model not only for functional studies of gamete killers but also for investigating their evolutionary trajectories at complementary geographical scales.

QTL mapping and transcriptome analysis of seed germination under PEG-induced water stress in Lactuca spp

The impact of limited water availability on lettuce growth has been well documented. However, the mechanisms by which lettuce controls seed germination under water stress remain unknown. Germination percentage was evaluated in the cv. Salinas (Lactuca sativa) (L. sativa) × US96UC23 (Lactuca serriola) (L. serriola) recombinant inbred line (RIL) population and USDA germplasm collection using 10% polyethylene glycol (PEG). About 50% of both populations displayed less than 90% germination. The average broad-sense heritability (H2) for germination percentage was 0.81 across both populations. Two quantitative trait loci (QTL) for germination percentage were identified on chromosomes 4 and 8 in the RIL population. The RNA-Seq and network analyses of wild lettuce, US96UC23, were performed using the control (distilled water, dH2O) and treatment (10% PEG) datasets. The number of differentially expressed genes (DEGs) was 4,095. The top 20 gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were assessed by enrichment analysis. The consensus network analysis captured 44 modules. Gene networks were constructed for the top 20 hub genes in 10 significant modules from each dataset. This study comprehensively explains QTL, GO terms, KEGG pathways, and gene networks associated with lettuce seed germination under osmotic stress.

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.

Adult resistance genes to barley powdery mildew confer basal penetration resistance associated with broad-spectrum resistance

Powdery mildew isa major disease of barley (Hordeum vulgare L.) for which breeders have traditionally relied on dominant, pathogen race-specific resistance genes for genetic control. Directional selection pressures in extensive monocultures invariably result in such genes being overcome as the pathogen mutates to evade recognition. This has led to a widespread reliance on fungicides and a single broad-spectrum recessive resistance provided by the mlo gene. The range of resistance genes and alleles found in wild crop relatives and landraces has been reduced in agricultural cultivars through an erosion of genetic diversity during domestication and selective breeding. Three novel major-effect adult plant resistance (APR) genes from landraces, designated Resistance to Blumeria graminis f. sp. hordei (Rbgh1 to Rbgh3), were identified in the terminal regions of barley chromosomes 5HL, 7HS, and 1HS, respectively. The phenotype of the new APR genes showed neither pronounced penetration resistance, nor the spontaneous necrosis and mesophyll cell death typical of mlo resistance, nor a whole epidermal cell hypersensitive response, typical of race-specific resistance. Instead, resistance was localized to the site of attempted penetration in an epidermal cell and was associated with cell wall appositions and cytosolic vesicle-like bodies, and lacked strong induction of reactive oxygen species. The APR genes exhibited differences in vesicle-like body sizes, their distribution, and the extent of localized 3,3-diaminobenzidine staining in individual doubled haploid lines. The results revealed a set of unique basal penetration resistance genes that offer opportunities for combining different resistance mechanisms in breeding programs for robust mildew resistance.

Efficient construction of a linkage map and haplotypes for Mentha suaveolens using sequence capture

The sustainability of many crops is hindered by the lack of genomic resources and a poor understanding of natural genetic diversity. Particularly, application of modern breeding requires high-density linkage maps that are integrated into a highly contiguous reference genome. Here, we present a rapid method for deriving haplotypes and developing linkage maps, and its application to Mentha suaveolens, one of the diploid progenitors of cultivated mints. Using sequence-capture via DNA hybridization to target single nucleotide polymorphisms (SNPs), we successfully genotyped ∼5000 SNPs within the genome of >400 individuals derived from a self cross. After stringent quality control, and identification of nonredundant SNPs, 1919 informative SNPs were retained for linkage map construction. The resulting linkage map defined a total genetic space of 942.17 cM divided among 12 linkage groups, ranging from 56.32 to 122.61 cM in length. The linkage map is in good agreement with pseudomolecules from our preliminary genome assembly, proving this resource effective for the correction and validation of the reference genome. We discuss the advantages of this method for the rapid creation of linkage maps.

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.

QTLs and candidate genes analyses for fruit size under domestication and differentiation in melon (Cucumis melo L.) based on high resolution maps

BACKGROUND

Melon is a very important horticultural crop produced worldwide with high phenotypic diversity. Fruit size is among the most important domestication and differentiation traits in melon. The molecular mechanisms of fruit size in melon are largely unknown.

RESULTS

Two high-density genetic maps were constructed by whole-genome resequencing with two F₂ segregating populations (WAP and MAP) derived from two crosses (cultivated agrestis × wild agrestis and cultivated melo × cultivated agrestis). We obtained 1,871,671 and 1,976,589 high quality SNPs that show differences between parents in WAP and MAP. A total of 5138 and 5839 recombination events generated 954 bins in WAP and 1027 bins in MAP with the average size of 321.3 Kb and 301.4 Kb respectively. All bins were mapped onto 12 linkage groups in WAP and MAP. The total lengths of two linkage maps were 904.4 cM (WAP) and 874.5 cM (MAP), covering 86.6% and 87.4% of the melon genome. Two loci for fruit size were identified on chromosome 11 in WAP and chromosome 5 in MAP, respectively. An auxin response factor and a YABBY transcription factor were inferred to be the candidate genes for both loci.

CONCLUSION

The high-resolution genetic maps and QTLs analyses for fruit size described here will provide a better understanding the genetic basis of domestication and differentiation, and provide a valuable tool for map-based cloning and molecular marker assisted breeding.

Distal Bias of Meiotic Crossovers in Hexaploid Bread Wheat Reflects Spatio-Temporal Asymmetry of the Meiotic Program

Meiotic recombination generates genetic variation and provides physical links between homologous chromosomes (crossovers) essential for accurate segregation. In cereals the distribution of crossovers, cytologically evident as chiasmata, is biased toward the distal regions of chromosomes. This creates a bottleneck for plant breeders in the development of varieties with improved agronomic traits, as genes situated in the interstitial and centromere proximal regions of chromosomes rarely recombine. Recent advances in wheat genomics and genome engineering combined with well-developed wheat cytogenetics offer new opportunities to manipulate recombination and unlock genetic variation. As a basis for these investigations we have carried out a detailed analysis of meiotic progression in hexaploid wheat (Triticum aestivum) using immunolocalization of chromosome axis, synaptonemal complex and recombination proteins. 5-Bromo-2'-deoxyuridine (BrdU) labeling was used to determine the chronology of key events in relation to DNA replication. Axis morphogenesis, synapsis and recombination initiation were found to be spatio-temporally coordinated, beginning in the gene-dense distal chromosomal regions and later occurring in the interstitial/proximal regions. Moreover, meiotic progression in the distal regions was coordinated with the conserved chromatin cycles that are a feature of meiosis. This mirroring of the chiasma bias was also evident in the distribution of the gene-associated histone marks, H3K4me3 and H3K27me3; the repeat-associated mark, H3K27me1; and H3K9me3. We believe that this study provides a cytogenetic framework for functional studies and ongoing initiatives to manipulate recombination in the wheat genome.

An Induced Mutation in HvRECQL4 Increases the Overall Recombination and Restores Fertility in a Barley HvMLH3 Mutant Background

Plant breeding relies on the meiotic recombination or crossing over to generate the new combinations of the alleles along and among the chromosomes. However, crossing over is constrained in the crops such as barley by a combination of the low frequency and biased distribution. In this study, we attempted to identify the genes that limit the recombination by performing a suppressor screen for the restoration of fertility to the semi-fertile barley mutant desynaptic10 (des10), carrying a mutation in the barley ortholog of MutL-Homolog 3 (HvMLH3), a member of the MutL-homolog (MLH) family of DNA mismatch repair genes. des10 mutants exhibit reduced recombination and fewer chiasmata, resulting in the loss of obligate crossovers (COs) leading to chromosome mis-segregation. We identified several candidate suppressor lines and confirmed their restored fertility in an Hvmlh3 background in the subsequent generations. We focus on one of the candidate suppressor lines, SuppLine2099, which showed the most complete restoration of fertility. We characterized this line by using a target-sequence enrichment and sequencing (TENSEQ) capture array representing barley orthologs of 46 meiotic genes. We found that SuppLine2099 contained a C/T change in the anti-CO gene RecQ-like helicase 4 (RECQL4) resulting in the substitution of a non-polar glycine to a polar aspartic acid (G700D) amino acid in the conserved helicase domain. Single nucleotide polymorphism (SNP) genotyping of F₃ populations revealed a significant increase in the recombination frequency in lines with Hvrecql4 in the Hvmlh3 background that was associated with the restoration of fertility. The genotyping also indicated that there was nearly double the recombination levels in homozygous Hvrecql4 lines compared to the wild type (WT). However, we did not observe any significant change in the distribution of CO events. Our results confirm the anti-CO role of RECQL4 in a large genome cereal and establish the possibility of testing the utility of increasing recombination in the context of traditional crop improvement.

Distal Bias of Meiotic Crossovers in Hexaploid Bread Wheat Reflects Spatio-Temporal Asymmetry of the Meiotic Program

Meiotic recombination generates genetic variation and provides physical links between homologous chromosomes (crossovers) essential for accurate segregation. In cereals the distribution of crossovers, cytologically evident as chiasmata, is biased toward the distal regions of chromosomes. This creates a bottleneck for plant breeders in the development of varieties with improved agronomic traits, as genes situated in the interstitial and centromere proximal regions of chromosomes rarely recombine. Recent advances in wheat genomics and genome engineering combined with well-developed wheat cytogenetics offer new opportunities to manipulate recombination and unlock genetic variation. As a basis for these investigations we have carried out a detailed analysis of meiotic progression in hexaploid wheat (Triticum aestivum) using immunolocalization of chromosome axis, synaptonemal complex and recombination proteins. 5-Bromo-2'-deoxyuridine (BrdU) labeling was used to determine the chronology of key events in relation to DNA replication. Axis morphogenesis, synapsis and recombination initiation were found to be spatio-temporally coordinated, beginning in the gene-dense distal chromosomal regions and later occurring in the interstitial/proximal regions. Moreover, meiotic progression in the distal regions was coordinated with the conserved chromatin cycles that are a feature of meiosis. This mirroring of the chiasma bias was also evident in the distribution of the gene-associated histone marks, H3K4me3 and H3K27me3; the repeat-associated mark, H3K27me1; and H3K9me3. We believe that this study provides a cytogenetic framework for functional studies and ongoing initiatives to manipulate recombination in the wheat genome.

Segregation distortion: Utilizing simulated genotyping data to evaluate statistical methods

Segregation distortion is the phenomenon in which genotypes deviate from expected Mendelian ratios in the progeny of a cross between two varieties or species. There is not currently a widely used consensus for the appropriate statistical test, or more specifically the multiple testing correction procedure, used to detect segregation distortion for high-density single-nucleotide polymorphism (SNP) data. Here we examine the efficacy of various multiple testing procedures, including chi-square test with no correction for multiple testing, false-discovery rate correction and Bonferroni correction using an in-silico simulation of a biparental mapping population. We find that the false discovery rate correction best approximates the traditional p-value threshold of 0.05 for high-density marker data. We also utilize this simulation to test the effect of segregation distortion on the genetic mapping process, specifically on the formation of linkage groups during marker clustering. Only extreme segregation distortion was found to effect genetic mapping. In addition, we utilize replicate empirical mapping populations of wheat varieties Avalon and Cadenza to assess how often segregation distortion conforms to the same pattern between closely related wheat varieties.

The first genetic linkage map for Fraxinus pennsylvanica and syntenic relationships with four related species

The genetic linkage map for green ash (Fraxinus pennsylvanica) contains 1201 DNA markers in 23 linkage groups spanning 2008.87cM. The green ash map shows stronger synteny with coffee than tomato. Green ash (Fraxinus pennsylvanica) is an outcrossing, diploid (2n = 46) hardwood tree species, native to North America. Native ash species in North America are being threatened by the rapid spread of the emerald ash borer (EAB, Agrilus planipennis), an invasive pest from Asia. Green ash, the most widely distributed ash species, is severely affected by EAB infestation, yet few genomic resources for genetic studies and improvement of green ash are available. In this study, a total of 5712 high quality single nucleotide polymorphisms (SNPs) were discovered using a minimum allele frequency of 1% across the entire genome through genotyping-by-sequencing. We also screened hundreds of genomic- and EST-based microsatellite markers (SSRs) from previous de novo assemblies (Staton et al., PLoS ONE 10:e0145031, 2015; Lane et al., BMC Genom 17:702, 2016). A first genetic linkage map of green ash was constructed from 90 individuals in a full-sib family, combining 2719 SNP and 84 SSR segregating markers among the parental maps. The consensus SNP and SSR map contains a total of 1201 markers in 23 linkage groups spanning 2008.87 cM, at an average inter-marker distance of 1.67 cM with a minimum logarithm of odds of 6 and maximum recombination fraction of 0.40. Comparisons of the organization the green ash map with the genomes of asterid species coffee and tomato, and genomes of the rosid species poplar and peach, showed areas of conserved gene order, with overall synteny strongest with coffee.

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.

A regulator of early flowering in barley (Hordeum vulgare L.)

Heading date (HD) of cereals is an important trait for adaptation to diverse environments and is critical for determining yield and quality and the number of genes and gene combinations that confer earliness in barley under short days is limited. In our study, a QTL for early flowering was identified from the cross between an Australian malting barley cultivar and a Chinese landrace. Four sets of near isogenic lines (NILs) were developed with a QTL located on chromosome 5H at the interval of 122.0-129.0 cM. Further experiments were conducted to investigate how this gene was regulated by photoperiod using the NILs with three sowing dates from autumn to summer. The NILs carrying the earliness allele were significantly earlier than the late genotype at all sowing dates. This gene was different from previously reported vernalisation genes that are located at a similar position as no vernalisation was required for all the NILs. The difference between this gene and Eam5 (HvPHYC) locus which also located between two co-segregated markers (3398516S5, 122.5 cM, and 4014046D5, 126.1 cM), is that with the existence of Ppd-H1 (Eam1), Eam5 has no effect on ear emergence under long days while the gene from TX9425 still reduced the time to ear emergency. The locus showed no pleiotropic effects on grain pasting properties and agronomic traits except for spike length and number of spikelets per spike, and thus can be effectively used in breeding programs. The array of early heading dates caused by interactions of Eam5 gene with other maturity genes provides an opportunity to better fine tune heading dates with production environments, which can be critical factor in barley breeding.

QTL mapping uncovers a semi-dwarf 1 (sdw1) allele in the barley (Hordeum vulgare) ND23049 line

In barley, semi-dwarf varieties are attractive for their superior harvest index and lodging resistance, but many semi-dwarf barley genotypes suffer from poor spike emergence. We performed a genetic characterization of a semi-dwarf line (ND23049) that combines short stature, strong stiff culms, and adequate spike emergence. We developed a doubled haploid (DH) population derived by crossing ND23049 and the cultivar CLE253. A subset of 88 DH lines and parents were characterized for plant height in 2013 and 2014 and genotyped. In total, 1984 SNPs (345 unique loci) were used to produce a linkage map of 1127.1 cM. Three QTLs for plant height were detected in this population and coincided with the HvGA20ox2/Sdw1, HvBRI1/Uzu1, and HvPRR95 gene loci. The phenotypic variation explained by each QTL was 75.8%, 7.7%, and 4.1%, respectively, and jointly explained 83.3% (2013) and 87.7% (2014) of plant height. Our results suggest that ND23049 contributed the “short” allele at the HvGA20ox2/sdw1 locus while CLE253 provided “short” alleles at the HvBRI1/uzu1 and HvPRR95 loci. We identified a large deletion (at least 92.7 Kb), including HvGA20ox2 (Sdw1), as the causal mutation in ND23049. A set of tightly flanked SNP markers will help breeders to develop improved semi-dwarf varieties.

High-Density Genetic Map Construction and Identification of QTLs Controlling Oleic and Linoleic Acid in Peanut using SLAF-seq and SSRs

The cultivated peanut, A. hypogaea L., is an important oil and food crop globally.High-density genetic linkage mapping is a valuable and effective method for exploring complex quantitative traits. In this context, a recombinant inbred line (RIL) of 146 lines was developed by crossing Huayu28 and P76. We developed 433,679 high-quality SLAFs, of which 29,075 were polymorphic. 4,817 SLAFs were encoded and grouped into different segregation patterns. A high-resolution genetic map containing 2,334 markers (68 SSRs and 2,266 SNPs) on 20 linkage groups (LGs) spanning 2586.37 cM was constructed for peanut. The average distance between adjacent markers was 2.25 cM. Based on phenotyping in seven environments, QTLs for oleic acid (C18:1), linoleic acid (C18:2) and the ratio of oleic acid to linoleic acid (O/L) were identified and positioned on linkage groups A03, A04, A09, B09 and B10. Marker2575339 and Marker2379598 in B09 were associated with C18:1, C18:2 and O/L in seven environments, Marker4391589 and Marker4463600 in A09 were associated with C18:1, C18:2 and O/L in six environments. This map exhibits high resolution and accuracy, which will facilitate QTL discovery for essential agronomic traits in peanut.

Genotyping-by-Sequencing Derived High-Density Linkage Map and its Application to QTL Mapping of Flag Leaf Traits in Bread Wheat

Winter wheat parents ‘Harry’ (drought tolerant) and ‘Wesley’ (drought susceptible) were used to develop a recombinant inbred population with future goals of identifying genomic regions associated with drought tolerance. To precisely map genomic regions, high-density linkage maps are a prerequisite. In this study genotyping-by- sequencing (GBS) was used to construct the high-density linkage map. The map contained 3,641 markers distributed on 21 chromosomes and spanned 1,959 cM with an average distance of 1.8 cM between markers. The constructed linkage map revealed strong collinearity in marker order across 21 chromosomes with POPSEQ-v2.0, which was based on a high-density linkage map. The reliability of the linkage map for QTL mapping was demonstrated by co-localizing the genes to previously mapped genomic regions for two highly heritable traits, chaff color, and leaf cuticular wax. Applicability of linkage map for QTL mapping of three quantitative traits, flag leaf length, width, and area, identified 21 QTLs in four environments, and QTL expression varied across the environments. Two major stable QTLs, one each for flag leaf length (Qfll.hww-7A) and flag leaf width (Qflw.hww-5A) were identified. The map constructed will facilitate QTL and fine mapping of quantitative traits, map-based cloning, comparative mapping, and in marker-assisted wheat breeding endeavors.

The First Genetic Map in Sweet Osmanthus (Osmanthus fragrans Lour.) Using Specific Locus Amplified Fragment Sequencing

Osmanthus fragrans is an ornamental plant of substantial commercial value, and no genetic linkage maps of this species have previously been reported. Specific-locus amplified fragment sequencing (SLAF-seq) is a recently developed technology that allows massive single nucleotide polymorphisms (SNPs) to be identified and high-resolution genotyping. In our current research, we generated the first genetic map of O. fragrans using SLAF-seq, which is composed with 206.92 M paired-end reads and 173,537 SLAF markers. Among total 90,715 polymorphic SLAF markers, 15,317 polymorphic SLAFs could be used for genetic map construction. The integrated map contained 14,189 high quality SLAFs that were grouped in 23 genetic linkage groups, with a total length of 2962.46 cM and an average distance of 0.21 cM between two adjacent markers. In addition, 23,664 SNPs were identified from the mapped markers. As far as we know, this is the first of the genetic map of O. fragrans. Our results are further demonstrate that SLAF-seq is a very effective method for developing markers and constructing high-density linkage maps. The SNP markers and the genetic map reported in this study should be valuable resource in future research.

Genotyping-by-Sequencing Facilitates a High-Density Consensus Linkage Map for Aegilops umbellulata, a Wild Relative of Cultivated Wheat

High-density genetic maps are useful to precisely localize QTL or genes that might be used to improve traits of nutritional and/or economical importance in crops. However, high-density genetic maps are lacking for most wild relatives of crop species, including wheat. Aegilops umbellulata is a wild relative of wheat known for its potential as a source of biotic and abiotic stress resistance genes. In this work, we have developed a framework consensus genetic map using two biparental populations derived from accessions PI 298905, PI 542369, PI 5422375, and PI 554395. The framework map comprised 3009 genotype-by-sequence SNPs with a total map size of 948.72 cM. On average, there were three SNPs per centimorgan for each chromosome. Chromosome 1U was the shortest (66.5 cM), with only 81 SNPs, whereas the remaining chromosomes had between 391 and 591 SNP markers. A total of 2395 unmapped SNPs were added to the linkage maps through a recombination frequency approach, and increased the number of SNPs placed on the consensus map to a total of 5404 markers. Segregation distortion was disproportionally high for chromosome 1U for both populations used to construct component linkage maps, and thus segregation distortion could be one of the probable reasons for the exceptionally reduced linkage size for chromosome 1U. From comparative analysis, Ae. umbellulata chromosomes except 4U showed moderate to strong collinearity with corresponding homeologous chromosomes of hexaploid wheat and barley. The present consensus map may serve as a reference map in QTL mapping and validation projects, and also in genome assembly to develop a reference genome sequence for Ae. umbellulata.

Construction of a high-density genetic map using specific-locus amplified fragments in sorghum

BACKGROUND

Sorghum is mainly used as a human food and beverage source, playing an important role in the production of ethanol and other bio-industrial products. Thus it is regarded as a model crop for energy plants. Genetic map construction is the foundation for marker-assisted selection and gene cloning. So far several sorghum linkage maps have been reported using different kinds of molecular markers. However marker numbers and chromosome coverage are limited. As a result, it is difficult to get consistent results and the maps are hard to unify. In the present study, the genomes of 130 individuals consisting an F₂ population together with their parents were surveyed using a high-throughput sequencing technique. A high-density linkage map was constructed using specific-locus amplified fragments (SLAF) markers. This map can provide information and serve as a reference for effective gene exploration, and for marker assisted-breeding program.

RESULTS

A high-throughput sequencing method was adopted to screen SLAF markers with 130 F₂ individuals from a cross between a grain sorghum variety, J204, and a sweet sorghum variety, Keter. In the present study, 52,928 suitable SLAF markers out of 43,528,021 pair-end reads were chosen to conduct genetic map construction, 12.0% of which were polymorphic. Among the 6353 polymorphic SLAF markers, 5829 (91.8%) were successfully genotyped in the F₂ mapping population. Finally 2246 SLAF markers were obtained to construct a high-density genetic linkage map. The total distance of linkage map covering all 10 chromosomes was 2158.1 cM. The largest gap on each chromosome was 10.2 cM on average. The proportion of gaps less than and/or equal to 5.0 cM was averagely 98.1%. The markers on each chromosome ranged from 123 (chromosome 9) to 315 (chromosome 4) with a mean value of 224.6, the distance between adjacent markers ranged from 0.6 (chromosome 10) to 1.3 cM (chromosome 9) with an average distance of only 0.98 cM.

CONCLUSION

A high density sorghum genetic map was constructed in this study. The total length was 2158.1 cM covering all 10 chromosomes with a total number of 2246 SLAF markers. The construction of this map can provide detailed information for accurate gene localization and cloning and application of marker-assisted breeding.

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.

Extent and overlap of segregation distortion regions in 12 barley crosses determined via a Pool-GBS approach

Extent and overlap of segregation distortion regions in 12 barley crosses determined via a Pool-GBS approach. Segregation distortion is undesirable as it alters the frequency of alleles and can reduce the chances of obtaining a particular combination of alleles. In this work, we have used a pooled genotyping-by-sequencing (Pool-GBS) approach to estimate allelic frequencies and used it to examine segregation distortion in 12 segregating populations of barley derived from androgenesis. Thanks to the extensive genome-wide SNP coverage achieved (between 674 and 1744 markers), we determined that the proportion of distorted markers averaged 28.9 % while 25.3 % of the genetic map fell within segregation distortion regions (SDRs). These SDRs were characterized and identified based on the position of the marker showing the largest distortion and the span of each SDR. Summed across all 12 crosses, 36 different SDR peaks could be distinguished from a total of 50 SDRs and a majority of these SDRs (27 of 36) were observed in only one population. While most shared SDRs were common to only two crosses, two SDRs (SDR3.1 and SDR4.2) were exceptionally recurrent (seen in five and four crosses, respectively). Because of the broad span of most SDRs, an average of 30 % of crosses showed segregation distortion in any given chromosomal segment. In reciprocal crosses, although some SDRs were clearly shared, others were unique to a single direction. In summary, segregation distortion is highly variable in its extent and the number of loci underpinning these distortions seems to be quite large even in a narrow germplasm such as six-row spring barley.

Quantitative Trait Loci for Yield and Yield-Related Traits in Spring Barley Populations Derived from Crosses between European and Syrian Cultivars

In response to climatic changes, breeding programmes should be aimed at creating new cultivars with improved resistance to water scarcity. The objective of this study was to examine the yield potential of barley recombinant inbred lines (RILs) derived from three cross-combinations of European and Syrian spring cultivars, and to identify quantitative trait loci (QTLs) for yield-related traits in these populations. RILs were evaluated in field experiments over a period of three years (2011 to 2013) and genotyped with simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers; a genetic map for each population was constructed and then one consensus map was developed. Biological interpretation of identified QTLs was achieved by reference to Ensembl Plants barley gene space. Twelve regions in the genomes of studied RILs were distinguished after QTL analysis. Most of the QTLs were identified on the 2H chromosome, which was the hotspot region in all three populations. Syrian parental cultivars contributed alleles decreasing traits' values at majority of QTLs for grain weight, grain number, spike length and time to heading, and numerous alleles increasing stem length. The phenomic and molecular approaches distinguished the lines with an acceptable grain yield potential combining desirable features or alleles from their parents, that is, early heading from the Syrian breeding line (Cam/B1/CI08887//CI05761) and short plant stature from the European semidwarf cultivar (Maresi).

Genome Wide SSR High Density Genetic Map Construction from an Interspecific Cross of Gossypium hirsutum × Gossypium tomentosum

A high density genetic map was constructed using F2 population derived from an interspecific cross of G. hirsutum × G. tomentosum. The map consisted of 3093 marker loci distributed across all the 26 chromosomes and covered 4365.3 cM of cotton genome with an average inter-marker distance of 1.48 cM. The maximum length of chromosome was 218.38 cM and the minimum was 122.09 cM with an average length of 167.90 cM. A sub-genome covers more genetic distance (2189.01 cM) with an average inter loci distance of 1.53 cM than D sub-genome which covers a length of 2176.29 cM with an average distance of 1.43 cM. There were 716 distorted loci in the map accounting for 23.14% and most distorted loci were distributed on D sub-genome (25.06%), which were more than on A sub-genome (21.23%). In our map 49 segregation hotspots (SDR) were distributed across the genome with more on D sub-genome as compared to A genome. Two post-polyploidization reciprocal translocations of "A2/A3 and A4/A5" were suggested by seven pairs of duplicate loci. The map constructed through these studies is one of the three densest genetic maps in cotton however; this is the first dense genome wide SSR interspecific genetic map between G. hirsutum and G. tomentosum.

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.

Identification and Characterization of Segregation Distortion Loci on Cotton Chromosome 18

Segregation distortion is commonly detected via genetic mapping and this phenomenon has been reported in many species. However, the genetic causes of the segregation distortion regions in a majority of species are still unclear. To genetically dissect the SD on chromosome 18 in cotton, eight reciprocal backcross populations and two F₂ populations were developed. Eleven segregation distortion loci (SDL) were detected in these ten populations. Comparative analyses among populations revealed that SDL18.1 and SDL18.9 were consistent with male gametic competition; whereas SDL18.4 and SDL18.11 reflected female gametic selection. Similarly, other SDL could reflect zygotic selection. The surprising finding was that SDL18.8 was detected in all populations, and the direction was skewed towards heterozygotes. Consequently, zygotic selection or heterosis could represent the underlying genetic mechanism for SDL18.8. Among developed introgression lines, SDL18.8 was introgressed as a heterozygote, further substantiating that a heterozygote state was preferred under competition. Six out of 11 SDL on chromosome 18 were dependent on the cytoplasmic environment. These results indicated that different SDL showed varying responses to the cytoplasmic environment. Overall, the results provided a novel strategy to analyze the molecular mechanisms, which could be further exploited in cotton interspecific breeding programs.

An ultra-high-density map as a community resource for discerning the genetic basis of quantitative traits in maize

BACKGROUND

To safeguard the food supply for the growing human population, it is important to understand and exploit the genetic basis of quantitative traits. Next-generation sequencing technology performs advantageously and effectively in genetic mapping and genome analysis of diverse genetic resources. Hence, we combined re-sequencing technology and a bin map strategy to construct an ultra-high-density bin map with thousands of bin markers to precisely map a quantitative trait locus.

RESULTS

In this study, we generated a linkage map containing 1,151,856 high quality SNPs between Mo17 and B73, which were verified in the maize intermated B73 × Mo17 (IBM) Syn10 population. This resource is an excellent complement to existing maize genetic maps available in an online database  (iPlant, http://data.maizecode.org/maize/qtl/syn10/ ). Moreover, in this population combined with the IBM Syn4 RIL population, we detected 135 QTLs for flowering time and plant height traits across the two populations. Eighteen known functional genes and twenty-five candidate genes for flowering time and plant height trait were fine-mapped into a 2.21-4.96 Mb interval. Map expansion and segregation distortion were also analyzed, and evidence for inadvertent selection of early flowering time in the process of mapping population development was observed. Furthermore, an updated integrated map with 1,151,856 high-quality SNPs, 2,916 traditional markers and 6,618 bin markers was constructed. The data were deposited into the iPlant Discovery Environment (DE), which provides a fundamental resource of genetic data for the maize genetic research community.

CONCLUSIONS

Our findings provide basic essential genetic data for the maize genetic research community. An updated IBM Syn10 population and a reliable, verified high-quality SNP set between Mo17 and B73 will aid in future molecular breeding efforts.

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.

A high-density SSR genetic map constructed from a F2 population of Gossypium hirsutum and Gossypium darwinii

The cultivated allotetraploid species Gossypium hirsutum, accounts for 90% of the world cotton production, has narrow genetic basis that's why its yield, quality or stress resistance breeding is stagnant. It is therefore, essential to explore desirable genes from Gossypium darwinii which has enviable traits such as high fiber fineness, drought tolerance, fusarium and verticillium resistance. We used G. darwinii as primary plant materials in this study not only to enrich the genetic diversity of exiting germplasm but also to better understand its genome structure. An interspecific high density linkage map of allotetraploid cotton was constructed using F2 population (G. hirsutum×G. darwinii). The map was based entirely on genome-wide simple sequence repeat (SSR) markers. A total of 2763 markers were mapped in 26 linkage groups (chromosomes) covering a genome length of 4176.7cM with an average inter-locus distance of 1.5cM. The length of the chromosomes ranged from 84.7 to 238.5cM with an average length of 160.6cM. At subgenome length was 2160.7cM with an average distance of 1.6cM, where as Dt genome length was 2016cM with an average distance of 1.4cM. There were 601 distorted SSR loci. Less number of segregation distortion loci were located in At subgenome than in Dt subgenome. Two post-polyploidization reciprocal translocations of "A2/A3 and A4/A5" were suggested by 44 pairs of duplicate loci.

A new allele of acid soil tolerance gene from a malting barley variety

Background

Acid soil is a serious limitation to crop production all over the world. Toxic aluminium (Al) cations in acid soil inhibit root growth and reduce yield. Although a gene tolerant to acid soil has been identified, it has not been used in malting barley breeding, which is partly due to the acid soil tolerance gene being linked to unfavorable malting quality traits.

Results

A Brazilian malting barley variety Br2 was identified as tolerant to acid soil. A doubled haploid (DH) population was developed from a cross between Br2 and the Australian acid-sensitive cultivar Hamelin. The DH population was tested for acid soil tolerance in native acid soil and a hydroponic system with pH 4.2, pH 4.2 + Al or pH 6.5, and genotyped using SSR, DArT and gene-specific markers. A single QTL was detected for all parameters related to acid soil tolerance. The QTL was mapped to the known HvMATE location on chromosome 4H. Sequence alignment of the HvMATE gene identified 13 INDELs and 87 SNPs, where one SNP coded for a single amino acid difference between the two varieties. A gene-specific marker was developed to detect the single nucleotide polymorphism between Hamelin and Br2. This marker co-segregated with aluminium tolerance and accounted for 79 % of phenotypic variation for acid soil tolerance.

Conclusion

The present study identified a novel source of acid soil/Al tolerance from a Brazilian malting barley cultivar Br2. This variety tolerated Al toxicity but was sensitive to low pH which is similar to most other Al-tolerant varieties. A gene-specific marker Cit7 was developed based on the HvMATE gene sequence. Cit7 will improve the efficiency of molecular-assisted selection of new barley varieties with tolerance to acid soil. Multiple alleles exist for the acid soil tolerance gene on chromosome 4H, so a malting barley variety that tolerates acid soil could be developed by selecting suitable tolerant alleles. Tolerance to low pH may play an important role for barley to adapt to acid soils.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0254-4) contains supplementary material, which is available to authorized users.

Using QTL mapping to investigate the relationships between abiotic stress tolerance (drought and salinity) and agronomic and physiological traits

Background

Drought and salinity are two major abiotic stresses that severely limit barley production worldwide. Physiological and genetic complexity of these tolerance traits has significantly slowed the progress of developing stress-tolerant cultivars. Marker-assisted selection (MAS) may potentially overcome this problem. In the current research, seventy two double haploid (DH) lines from a cross between TX9425 (a Chinese landrace variety with superior drought and salinity tolerance) and a sensitive variety, Franklin were used to identify quantitative trait loci (QTL) for drought and salinity tolerance, based on a range of developmental and physiological traits.

Results

Two QTL for drought tolerance (leaf wilting under drought stress) and one QTL for salinity tolerance (plant survival under salt stress) were identified from this population. The QTL on 2H for drought tolerance determined 42% of phenotypic variation, based on three independent experiments. This QTL was closely linked with a gene controlling ear emergency. The QTL on 5H for drought tolerance was less affected by agronomic traits and can be effectively used in breeding programs. A candidate gene for this QTL on 5H was identified based on the draft barley genome sequence. The QTL for proline accumulation, under both drought and salinity stresses, were located on different positions to those for drought and salinity tolerance, indicating no relationship with plant tolerance to either of these stresses.

Conclusions

Using QTL mapping, the relationships between QTL for agronomic and physiological traits and plant drought and salinity tolerance were studied. A new QTL for drought tolerance which was not linked to any of the studied traits was identified. This QTL can be effectively used in breeding programs. It was also shown that proline accumulation under stresses was not necessarily linked with drought or salinity tolerance based on methods of phenotyping used in this experiment. The use of proline content in breeding programs can also be limited by the accuracy of phenotyping.

Electronic supplementary material

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

A high-density SNP Map of sunflower derived from RAD-sequencing facilitating fine-mapping of the rust resistance gene R12

A high-resolution genetic map of sunflower was constructed by integrating SNP data from three F₂ mapping populations (HA 89/RHA 464, B-line/RHA 464, and CR 29/RHA 468). The consensus map spanned a total length of 1443.84 cM, and consisted of 5,019 SNP markers derived from RAD tag sequencing and 118 publicly available SSR markers distributed in 17 linkage groups, corresponding to the haploid chromosome number of sunflower. The maximum interval between markers in the consensus map is 12.37 cM and the average distance is 0.28 cM between adjacent markers. Despite a few short-distance inversions in marker order, the consensus map showed high levels of collinearity among individual maps with an average Spearman's rank correlation coefficient of 0.972 across the genome. The order of the SSR markers on the consensus map was also in agreement with the order of the individual map and with previously published sunflower maps. Three individual and one consensus maps revealed the uneven distribution of markers across the genome. Additionally, we performed fine mapping and marker validation of the rust resistance gene R₁₂, providing closely linked SNP markers for marker-assisted selection of this gene in sunflower breeding programs. This high resolution consensus map will serve as a valuable tool to the sunflower community for studying marker-trait association of important agronomic traits, marker assisted breeding, map-based gene cloning, and comparative mapping.

Construction of a SNP-based genetic linkage map in cultivated peanut based on large scale marker development using next-generation double-digest restriction-site-associated DNA sequencing (ddRADseq)

BACKGROUND

Cultivated peanut, or groundnut (Arachis hypogaea L.), is an important oilseed crop with an allotetraploid genome (AABB, 2n=4x=40). In recent years, many efforts have been made to construct linkage maps in cultivated peanut, but almost all of these maps were constructed using low-throughput molecular markers, and most show a low density, directly influencing the value of their applications. With advances in next-generation sequencing (NGS) technology, the construction of high-density genetic maps has become more achievable in a cost-effective and rapid manner. The objective of this study was to establish a high-density single nucleotide polymorphism (SNP)-based genetic map for cultivated peanut by analyzing next-generation double-digest restriction-site-associated DNA sequencing (ddRADseq) reads.

RESULTS

We constructed reduced representation libraries (RRLs) for two A. hypogaea lines and 166 of their recombinant inbred line (RIL) progenies using the ddRADseq technique. Approximately 175 gigabases of data containing 952,679,665 paired-end reads were obtained following Solexa sequencing. Mining this dataset, 53,257 SNPs were detected between the parents, of which 14,663 SNPs were also detected in the population, and 1,765 of the obtained polymorphic markers met the requirements for use in the construction of a genetic map. Among 50 randomly selected in silico SNPs, 47 were able to be successfully validated. One linkage map was constructed, which was comprised of 1,685 marker loci, including 1,621 SNPs and 64 simple sequence repeat (SSR) markers. The map displayed a distribution of the markers into 20 linkage groups (LGs A01-A10 and B01-B10), spanning a distance of 1,446.7 cM. The alignment of the LGs from this map was shown in comparison with a previously integrated consensus map from peanut.

CONCLUSIONS

This study showed that the ddRAD library combined with NGS allowed the rapid discovery of a large number of SNPs in the cultivated peanut. The first high density SNP-based linkage map for A. hypogaea was generated that can serve as a reference map for cultivated Arachis species and will be useful in genetic mapping. Our results contribute to the available molecular marker resources and to the assembly of a reference genome sequence for the peanut.

High density SNP and SSR-based genetic maps of two independent oil palm hybrids

BACKGROUND

Oil palm is an important perennial oil crop with an extremely long selection cycle of 10 to 12 years. As such, any tool that speeds up its genetic improvement process, such as marker-assisted breeding is invaluable. Previously, genetic linkage maps based on AFLP, RFLP and SSR markers were developed and QTLs for fatty acid composition and yield components identified. High density genetic maps of crosses of different genetic backgrounds are indispensable tools for investigating oil palm genetics. They are also useful for comparative mapping analyses to identify markers closely linked to traits of interest.

RESULTS

A 4.5 K customized oil palm SNP array was developed using the Illumina Infinium platform. The SNPs and 252 SSRs were genotyped on two mapping populations, an intraspecific cross with 87 palms and an interspecific cross with 108 palms. Parental maps with 16 linkage groups (LGs), were constructed for the three fruit forms of E. guineensis (dura, pisifera and tenera). Map resolution was further increased by integrating the dura and pisifera maps into an intraspecific integrated map with 1,331 markers spanning 1,867 cM. We also report the first map of a Colombian E. oleifera, comprising 10 LGs with 65 markers spanning 471 cM. Although not very dense due to the high level of homozygosity in E. oleifera, the LGs were successfully integrated with the LGs of the tenera map. Direct comparison between the parental maps identified 603 transferable markers polymorphic in at least two of the parents. Further analysis revealed a high degree of marker transferability covering 1,075 cM, between the intra- and interspecific integrated maps. The interspecific cross displayed higher segregation distortion than the intraspecific cross. However, inclusion of distorted markers in the genetic maps did not disrupt the marker order and no map expansion was observed.

CONCLUSIONS

The high density SNP and SSR-based genetic maps reported in this paper have greatly improved marker density and genome coverage in comparison with the first reference map based on AFLP and SSR markers. Therefore, it is foreseen that they will be more useful for fine mapping of QTLs and whole genome association mapping studies in oil palm.

Genome-wide distribution of genetic diversity and linkage disequilibrium in a mass-selected population of maritime pine

Background

The accessibility of high-throughput genotyping technologies has contributed greatly to the development of genomic resources in non-model organisms. High-density genotyping arrays have only recently been developed for some economically important species such as conifers. The potential for using genomic technologies in association mapping and breeding depends largely on the genome wide patterns of diversity and linkage disequilibrium in current breeding populations. This study aims to deepen our knowledge regarding these issues in maritime pine, the first species used for reforestation in south western Europe.

Results

Using a new map merging algorithm, we first established a 1,712 cM composite linkage map (comprising 1,838 SNP markers in 12 linkage groups) by bringing together three already available genetic maps. Using rigorous statistical testing based on kernel density estimation and resampling we identified cold and hot spots of recombination. In parallel, 186 unrelated trees of a mass-selected population were genotyped using a 12k-SNP array. A total of 2,600 informative SNPs allowed to describe historical recombination, genetic diversity and genetic structure of this recently domesticated breeding pool that forms the basis of much of the current and future breeding of this species. We observe very low levels of population genetic structure and find no evidence that artificial selection has caused a reduction in genetic diversity. By combining these two pieces of information, we provided the map position of 1,671 SNPs corresponding to 1,192 different loci. This made it possible to analyze the spatial pattern of genetic diversity (H_e) and long distance linkage disequilibrium (LD) along the chromosomes. We found no particular pattern in the empirical variogram of H_e across the 12 linkage groups and, as expected for an outcrossing species with large effective population size, we observed an almost complete lack of long distance LD.

Conclusions

These results are a stepping stone for the development of strategies for studies in population genomics, association mapping and genomic prediction in this economical and ecologically important forest tree species.

A new QTL for plant height in barley (Hordeum vulgare L.) showing no negative effects on grain yield

INTRODUCTION

Reducing plant height has played an important role in improving crop yields. The success of a breeding program relies on the source of dwarfing genes. For a dwarfing or semi-dwarfing gene to be successfully used in a breeding program, the gene should have minimal negative effects on yield and perform consistently in different environments.

METHODS

In this study, 182 doubled haploid lines, generated from a cross between TX9425 and Naso Nijo, were grown in six different environments to identify quantitative trait loci (QTL) controlling plant height and investigate QTL × environments interaction.

RESULTS

A QTL for plant was identified on 7H. This QTL showed no significant effects on other agronomic traits and yield components and consistently expressed in the six environments. A sufficient allelic effect makes it possible for this QTL to be successfully used in breeding programs.

Constructing a new integrated genetic linkage map and mapping quantitative trait loci for vegetative mycelium growth rate in Lentinula edodes

The most saturated linkage map for Lentinula edodes to date was constructed based on a monokaryotic population of 146 single spore isolates (SSIs) using sequence-related amplified polymorphism (SRAP), target region amplification polymorphism (TRAP), insertion-deletion (InDel) markers, and the mating-type loci. Five hundred and twenty-four markers were located on 13 linkage groups (LGs). The map spanned a total length of 1006.1 cM, with an average marker spacing of 2.0 cM. Quantitative trait loci (QTLs) mapping was utilized to uncover the loci regulating and controlling the vegetative mycelium growth rate on various synthetic media, and complex medium for commercial cultivation of L. edodes. Two and 13 putative QTLs, identified respectively in the monokaryotic population and two testcross dikaryotic populations, were mapped on seven different LGs. Several vegetative mycelium growth rate-related QTLs uncovered here were clustered on LG4 (Qmgr1, Qdgr1, Qdgr2 and Qdgr9) and LG6 (Qdgr3, Qdgr4 and Qdgr5), implying the presence of main genomic areas responsible for growth rate regulation and control. The QTL hotspot region on LG4 was found to be in close proximity to the region containing the mating-type A (MAT-A) locus. Moreover, Qdgr2 on LG4 was detected on different media, contributing 8.07 %-23.71 % of the phenotypic variation. The present study provides essential information for QTL mapping and marker-assisted selection (MAS) in L. edodes.

Quantitative high resolution mapping of HvMLH3 foci in barley pachytene nuclei reveals a strong distal bias and weak interference

In barley (Hordeum vulgare L.), chiasmata (the physical sites of genetic crossovers) are skewed towards the distal ends of chromosomes, effectively consigning a large proportion of genes to recombination coldspots. This has the effect of limiting potential genetic variability, and of reducing the efficiency of map-based cloning and breeding approaches for this crop. Shifting the sites of recombination to more proximal chromosome regions by forward and reverse genetic means may be profitable in terms of realizing the genetic potential of the species, but is predicated upon a better understanding of the mechanisms governing the sites of these events, and upon the ability to recognize real changes in recombination patterns. The barley MutL Homologue (HvMLH3), a marker for class I interfering crossovers, has been isolated and a specific antibody has been raised. Immunolocalization of HvMLH3 along with the synaptonemal complex transverse filament protein ZYP1, used in conjunction with fluorescence in situ hybridization (FISH) tagging of specific barley chromosomes, has enabled access to the physical recombination landscape of the barley cultivars Morex and Bowman. Consistent distal localization of HvMLH3 foci throughout the genome, and similar patterns of HvMLH3 foci within bivalents 2H and 3H have been observed. A difference in total numbers of HvMLH3 foci between these two cultivars has been quantified, which is interpreted as representing genotypic variation in class I crossover frequency. Discrepancies between the frequencies of HvMLH3 foci and crossover frequencies derived from linkage analysis point to the existence of at least two crossover pathways in barley. It is also shown that interference of HvMLH3 foci is relatively weak compared with other plant species.

An Ultra-High-Density, Transcript-Based, Genetic Map of Lettuce

We have generated an ultra-high-density genetic map for lettuce, an economically important member of the Compositae, consisting of 12,842 unigenes (13,943 markers) mapped in 3696 genetic bins distributed over nine chromosomal linkage groups. Genomic DNA was hybridized to a custom Affymetrix oligonucleotide array containing 6.4 million features representing 35,628 unigenes of Lactuca spp. Segregation of single-position polymorphisms was analyzed using 213 F7:8 recombinant inbred lines that had been generated by crossing cultivated Lactuca sativa cv. Salinas and L. serriola acc. US96UC23, the wild progenitor species of L. sativa The high level of replication of each allele in the recombinant inbred lines was exploited to identify single-position polymorphisms that were assigned to parental haplotypes. Marker information has been made available using GBrowse to facilitate access to the map. This map has been anchored to the previously published integrated map of lettuce providing candidate genes for multiple phenotypes. The high density of markers achieved in this ultradense map allowed syntenic studies between lettuce and Vitis vinifera as well as other plant species.

A consensus linkage map of common carp (Cyprinus carpio L.) to compare the distribution and variation of QTLs associated with growth traits

The ability to detect and identify quantitative trait loci (QTLs) in a single population is often limited. Analyzing multiple populations in QTL analysis improves the power of detecting QTLs and provides a better understanding of their functional allelic variation and distribution. In this study, a consensus map of the common carp was constructed, based on four populations, to compare the distribution and variation of QTLs. The consensus map spans 2371.6 cM across the 42 linkage groups and comprises 257 microsatellites and 421 SNPs, with a mean marker interval of 3.7 cM/marker. Sixty-seven QTLs affecting four growth traits from the four populations were mapped to the consensus map. Only one QTL was common to three populations, and nine QTLs were detected in two populations. However, no QTL was common to all four populations. The results of the QTL comparison suggest that the QTLs are responsible for the phenotypic variability observed for these traits in a broad array of common carp germplasms. The study also reveals the different genetic performances between major and minor genes in different populations.

Diversifying sunflower germplasm by integration and mapping of a novel male fertility restoration gene

The combination of a single cytoplasmic male-sterile (CMS) PET-1 and the corresponding fertility restoration (Rf) gene Rf1 is used for commercial hybrid sunflower (Helianthus annuus L., 2n = 34) seed production worldwide. A new CMS line 514A was recently developed with H. tuberosus cytoplasm. However, 33 maintainers and restorers for CMS PET-1 and 20 additional tester lines failed to restore the fertility of CMS 514A. Here, we report the discovery, characterization, and molecular mapping of a novel Rf gene for CMS 514A derived from an amphiploid (Amp H. angustifolius/P 21, 2n = 68). Progeny analysis of the male-fertile (MF) plants (2n = 35) suggested that this gene, designated Rf6, was located on a single alien chromosome. Genomic in situ hybridization (GISH) indicated that Rf6 was on a chromosome with a small segment translocation on the long arm in the MF progenies (2n = 34). Rf6 was mapped to linkage group (LG) 3 of the sunflower SSR map. Eight markers were identified to be linked to this gene, covering a distance of 10.8 cM. Two markers, ORS13 and ORS1114, were only 1.6 cM away from the gene. Severe segregation distortions were observed for both the fertility trait and the linked marker loci, suggesting the possibility of a low frequency of recombination or gamete selection in this region. This study discovered a new CMS/Rf gene system derived from wild species and provided significant insight into the genetic basis of this system. This will diversify the germplasm for sunflower breeding and facilitate understanding of the interaction between the cytoplasm and nuclear genes.

A reference linkage map for Eucalyptus

Background

Genetic linkage maps are invaluable resources in plant research. They provide a key tool for many genetic applications including: mapping quantitative trait loci (QTL); comparative mapping; identifying unlinked (i.e. independent) DNA markers for fingerprinting, population genetics and phylogenetics; assisting genome sequence assembly; relating physical and recombination distances along the genome and map-based cloning of genes. Eucalypts are the dominant tree species in most Australian ecosystems and of economic importance globally as plantation trees. The genome sequence of E. grandis has recently been released providing unprecedented opportunities for genetic and genomic research in the genus. A robust reference linkage map containing sequence-based molecular markers is needed to capitalise on this resource. Several high density linkage maps have recently been constructed for the main commercial forestry species in the genus (E. grandis, E. urophylla and E. globulus) using sequenced Diversity Arrays Technology (DArT) and microsatellite markers. To provide a single reference linkage map for eucalypts a composite map was produced through the integration of data from seven independent mapping experiments (1950 individuals) using a marker-merging method.

Results

The composite map totalled 1107 cM and contained 4101 markers; comprising 3880 DArT, 213 microsatellite and eight candidate genes. Eighty-one DArT markers were mapped to two or more linkage groups, resulting in the 4101 markers being mapped to 4191 map positions. Approximately 13% of DArT markers mapped to identical map positions, thus the composite map contained 3634 unique loci at an average interval of 0.31 cM.

Conclusion

The composite map represents the most saturated linkage map yet produced in Eucalyptus. As the majority of DArT markers contained on the map have been sequenced, the map provides a direct link to the E. grandis genome sequence and will serve as an important reference for progressing eucalypt research.

Development of wild barley (Hordeum chilense)-derived DArT markers and their use into genetic and physical mapping

Diversity arrays technology (DArT) genomic libraries were developed from H. chilense accessions to support robust genotyping of this species and a novel crop comprising H. chilense genome (e.g., tritordeums). Over 11,000 DArT clones were obtained using two complexity reduction methods. A subset of 2,209 DArT markers was identified on the arrays containing these clones as polymorphic between parents and segregating in a population of 92 recombinant inbred lines (RIL) developed from the cross between H. chilense accessions H1 and H7. Using the segregation data a high-density map of 1,503 cM was constructed with average inter-bin density of 2.33 cM. A subset of DArT markers was also mapped physically using a set of wheat-H. chilense chromosome addition lines. It allowed the unambiguous assignment of linkage groups to chromosomes. Four segregation distortion regions (SDRs) were found on the chromosomes 2H(ch), 3H(ch) and 5H(ch) in agreement with previous findings in barley. The new map improves the genome coverage of previous H. chilense maps. H. chilense-derived DArT markers will enable further genetic studies in ongoing projects on hybrid wheat, seed carotenoid content improvement or tritordeum breeding program. Besides, the genetic map reported here will be very useful as the basis to develop comparative genomics studies with barley and model species.