High-Density Genetic Map Construction and Gene Mapping of Basal Branching Habit and Flowers per Leaf Axil in Sesame

A good genetic map can provide the framework for quantitative trait loci (QTL) analysis, map-based gene cloning, and genome sequence assembling. The main objectives of this study were to develop a high-density genetic linkage map using specific length amplified fragment sequencing (SLAF-seq) in sesame. In the result, a high-resolution genetic map with 9,378 SLAF markers and 13 linkage groups (LGs) was constructed. The map spanned a total genetic distance of 1,974.23 cM, and the mean LG length was 151.86 cM, with an average genetic distance of 0.22 cM between adjacent markers. Based on the newly constructed genetic map, genes for basal branching habit (SiBH) and flowers per leaf axil (SiFA) were mapped to LG5 and LG11, respectively.

Molecular Mapping of Oil Content and Fatty Acids Using Dense Genetic Maps in Groundnut (Arachis hypogaea L.)

Enhancing seed oil content with desirable fatty acid composition is one of the most important objectives of groundnut breeding programs globally. Genomics-assisted breeding facilitates combining multiple traits faster, however, requires linked markers. In this context, we have developed two different F₂ mapping populations, one for oil content (OC-population, ICGV 07368 × ICGV 06420) and another for fatty acid composition (FA-population, ICGV 06420 × SunOleic 95R). These two populations were phenotyped for respective traits and genotyped using Diversity Array Technology (DArT) and DArTseq genotyping platforms. Two genetic maps were developed with 854 (OC-population) and 1,435 (FA-population) marker loci with total map distance of 3,526 and 1,869 cM, respectively. Quantitative trait locus (QTL) analysis using genotyping and phenotyping data identified eight QTLs for oil content including two major QTLs, qOc-A10 and qOc-A02, with 22.11 and 10.37% phenotypic variance explained (PVE), respectively. For seven different fatty acids, a total of 21 QTLs with 7.6-78.6% PVE were identified and 20 of these QTLs were of major effect. Two mutant alleles, ahFAD2B and ahFAD2A, also had 18.44 and 10.78% PVE for palmitic acid, in addition to oleic (33.8 and 17.4% PVE) and linoleic (41.0 and 19.5% PVE) acids. Furthermore, four QTL clusters harboring more than three QTLs for fatty acids were identified on the three LGs. The QTLs identified in this study could be further dissected for candidate gene discovery and development of diagnostic markers for breeding improved groundnut varieties with high oil content and desirable oil quality.

A gapless genome sequence of the fungus Botrytis cinerea

Following earlier incomplete and fragmented versions of a genome sequence for the grey mould Botrytis cinerea, a gapless, near-finished genome sequence for B. cinerea strain B05.10 is reported. The assembly comprised 18 chromosomes and was confirmed by an optical map and a genetic map based on approximately 75 000 single nucleotide polymorphism (SNP) markers. All chromosomes contained fully assembled centromeric regions, and 10 chromosomes had telomeres on both ends. The genetic map consisted of 4153 cM and a comparison of the genetic distances with the physical distances identified 40 recombination hotspots. The linkage map also identified two mutations, located in the previously described genes Bos1 and BcsdhB, that conferred resistance to the fungicides boscalid and iprodione. The genome was predicted to encode 11 701 proteins. RNAseq data from >20 different samples were used to validate and improve gene models. Manual curation of chromosome 1 revealed interesting features, such as the occurrence of a dicistronic transcript and fully overlapping genes in opposite orientations, as well as many spliced antisense transcripts. Manual curation also revealed that the untranslated regions (UTRs) of genes can be complex and long, with many UTRs exceeding lengths of 1 kb and possessing multiple introns. Community annotation is in progress.

Evidence of Genomic Exchanges between Homeologous Chromosomes in a Cross of Peanut with Newly Synthetized Allotetraploid Hybrids

Cultivated peanut and synthetics are allotetraploids (2n = 4x = 40) with two homeologous sets of chromosomes. Meiosis in allotetraploid peanut is generally thought to show diploid-like behavior. However, a recent study pointed out the occurrence of recombination between homeologous chromosomes, especially when synthetic allotetraploids are used, challenging the view of disomic inheritance in peanut. In this study, we investigated the meiotic behavior of allotetraploid peanut using 380 SSR markers and 90 F2 progeny derived from the cross between Arachis hypogaea cv Fleur 11 (AABB) and ISATGR278-18 (AAKK), a synthetic allotetraploid that harbors a K-genome that was reported to pair with the cultivated B-genome during meiosis. Segregation analysis of SSR markers showed 42 codominant SSRs with unexpected null bands among some progeny. Chi-square tests for these loci deviate from the expected 1:2:1 Mendelian ratio under disomic inheritance. A linkage map of 357 codominant loci aligned on 20 linkage groups (LGs) with a total length of 1728 cM, averaging 5.1 cM between markers, was developed. Among the 10 homeologous sets of LGs, one set consisted of markers that all segregated in a polysomic-like pattern, six in a likely disomic pattern and the three remaining in a mixed pattern with disomic and polysomic loci clustered on the same LG. Moreover, we reported a substitution of homeologous chromosomes in some progeny. Our results suggest that the homeologous recombination events occurred between the A and K genomes in the newly synthesized allotetraploid and have been highlighted in the progeny. Homeologous exchanges are rarely observed in tetraploid peanut and have not yet been reported for AAKK and AABB genomes. The implications of these results on peanut breeding are discussed.

High-resolution genetic maps of Lotus japonicus and L. burttii based on re-sequencing of recombinant inbred lines

Recombinant inbred lines (RILs) derived from bi-parental populations are stable genetic resources, which are widely used for constructing genetic linkage maps. These genetic maps are essential for QTL mapping and can aid contig and scaffold anchoring in the final stages of genome assembly. In this study, two Lotus sp. RIL populations, Lotus japonicus MG-20 × Gifu and Gifu × L. burttii, were characterized by Illumina re-sequencing. Genotyping of 187 MG-20 × Gifu RILs at 87,140 marker positions and 96 Gifu × L. burttii RILs at 357,973 marker positions allowed us to accurately identify 1,929 recombination breakpoints in the MG-20 × Gifu RILs and 1,044 breakpoints in the Gifu × L. burttii population. The resulting high-density genetic maps now facilitate high-accuracy QTL mapping, identification of reference genome mis-assemblies, and characterization of structural variants.

Genetic Mapping of Millions of SNPs in Safflower (Carthamus tinctorius L.) via Whole-Genome Resequencing

Accurate assembly of complete genomes is facilitated by very high density genetic maps. We performed low-coverage, whole-genome shotgun sequencing on 96 F₆ recombinant inbred lines (RILs) of a cross between safflower (Carthamus tinctorius L.) and its wild progenitor (C. palaestinus Eig). We also produced a draft genome assembly of C. tinctorius covering 866 million bp (∼two-thirds) of the expected 1.35 Gbp genome after sequencing a single, short insert library to ∼21 × depth. Sequence reads from the RILs were mapped to this genome assembly to facilitate SNP identification, and the resulting polymorphisms were used to construct a genetic map. The resulting map included 2,008,196 genetically located SNPs in 1178 unique positions. A total of 57,270 scaffolds, each containing five or more mapped SNPs, were anchored to the map. This resulted in the assignment of sequence covering 14% of the expected genome length to a genetic position. Comparison of this safflower map to genetic maps of sunflower and lettuce revealed numerous chromosomal rearrangements, and the resulting patterns were consistent with a whole-genome duplication event in the lineage leading to sunflower. This sequence-based genetic map provides a powerful tool for the assembly of a low-cost draft genome of safflower, and the same general approach is expected to work for other species.

A gene-based map of the Nod factor-independent Aeschynomene evenia sheds new light on the evolution of nodulation and legume genomes

Aeschynomene evenia has emerged as a new model legume for the deciphering of the molecular mechanisms of an alternative symbiotic process that is independent of the Nod factors. Whereas most of the research on nitrogen-fixing symbiosis, legume genetics and genomics has so far focused on Galegoid and Phaseolid legumes, A. evenia falls in the more basal and understudied Dalbergioid clade along with peanut (Arachis hypogaea). To provide insights into the symbiotic genes content and the structure of the A. evenia genome, we established a gene-based genetic map for this species. Firstly, an RNAseq analysis was performed on the two parental lines selected to generate a F₂ mapping population. The transcriptomic data were used to develop molecular markers and they allowed the identification of most symbiotic genes. The resulting map comprised 364 markers arranged in 10 linkage groups (2n = 20). A comparative analysis with the sequenced genomes of Arachis duranensis and A. ipaensis, the diploid ancestors of peanut, indicated blocks of conserved macrosynteny. Altogether, these results provided important clues regarding the evolution of symbiotic genes in a Nod factor-independent context. They provide a basis for a genome sequencing project and pave the way for forward genetic analysis of symbiosis in A. evenia.

Large-scale development of SSR markers in tobacco and construction of a linkage map in flue-cured tobacco

Tobacco (Nicotiana tabacum L.), particularly flue-cured tobacco, is one of the most economically important nonfood crops and is also an important model system in plant biotechnology. Despite its importance, only limited molecular marker resources are available for genome analysis, genetic mapping, and breeding. Simple sequence repeats (SSR) are one of the most widely-used molecular markers, having significant advantages including that they are generally co-dominant, easy to use, abundant in eukaryotic organisms, and produce highly reproducible results. In this study, based on the genome sequence data of flue-cured tobacco (K326), we developed a total of 13,645 mostly novel SSR markers, which were working in a set of eighteen tobacco varieties of four different types. A mapping population of 213 backcross (BC1) individuals, which were derived from an intra-type cross between two flue-cured tobacco varieties, Y3 and K326, was selected for mapping. Based on the newly developed SSR markers as well as published SSR markers, we constructed a genetic map consisting of 626 SSR loci distributed across 24 linkage groups and covering a total length of 1120.45 cM with an average distance of 1.79 cM between adjacent markers, which is the highest density map of flue-cured tobacco till date.

New Software for the Fast Estimation of Population Recombination Rates (FastEPRR) in the Genomic Era

Genetic recombination is a very important evolutionary mechanism that mixes parental haplotypes and produces new raw material for organismal evolution. As a result, information on recombination rates is critical for biological research. In this paper, we introduce a new extremely fast open-source software package (FastEPRR) that uses machine learning to estimate recombination rate ρ (=4Ner) from intraspecific DNA polymorphism data. When ρ>10 and the number of sampled diploid individuals is large enough (≥50), the variance of ρFastEPRR remains slightly smaller than that of ρLDhat. The new estimate ρcomb (calculated by averaging ρFastEPRR and ρLDhat) has the smallest variance of all cases. When estimating ρFastEPRR, the finite-site model was employed to analyze cases with a high rate of recurrent mutations, and an additional method is proposed to consider the effect of variable recombination rates within windows. Simulations encompassing a wide range of parameters demonstrate that different evolutionary factors, such as demography and selection, may not increase the false positive rate of recombination hotspots. Overall, accuracy of FastEPRR is similar to the well-known method, LDhat, but requires far less computation time. Genetic maps for each human population (YRI, CEU, and CHB) extracted from the 1000 Genomes OMNI data set were obtained in less than 3 d using just a single CPU core. The Pearson Pairwise correlation coefficient between the ρFastEPRR and ρLDhat maps is very high, ranging between 0.929 and 0.987 at a 5-Mb scale. Considering that sample sizes for these kinds of data are increasing dramatically with advances in next-generation sequencing technologies, FastEPRR (freely available at http://www.picb.ac.cn/evolgen/) is expected to become a widely used tool for establishing genetic maps and studying recombination hotspots in the population genomic era.

A High-Resolution SNP Array-Based Linkage Map Anchors a New Domestic Cat Draft Genome Assembly and Provides Detailed Patterns of Recombination

High-resolution genetic and physical maps are invaluable tools for building accurate genome assemblies, and interpreting results of genome-wide association studies (GWAS). Previous genetic and physical maps anchored good quality draft assemblies of the domestic cat genome, enabling the discovery of numerous genes underlying hereditary disease and phenotypes of interest to the biomedical science and breeding communities. However, these maps lacked sufficient marker density to order thousands of shorter scaffolds in earlier assemblies, which instead relied heavily on comparative mapping with related species. A high-resolution map would aid in validating and ordering chromosome scaffolds from existing and new genome assemblies. Here, we describe a high-resolution genetic linkage map of the domestic cat genome based on genotyping 453 domestic cats from several multi-generational pedigrees on the Illumina 63K SNP array. The final maps include 58,055 SNP markers placed relative to 6637 markers with unique positions, distributed across all autosomes and the X chromosome. Our final sex-averaged maps span a total autosomal length of 4464 cM, the longest described linkage map for any mammal, confirming length estimates from a previous microsatellite-based map. The linkage map was used to order and orient the scaffolds from a substantially more contiguous domestic cat genome assembly (Felis catus v8.0), which incorporated ∼20 × coverage of Illumina fragment reads. The new genome assembly shows substantial improvements in contiguity, with a nearly fourfold increase in N50 scaffold size to 18 Mb. We use this map to report probable structural errors in previous maps and assemblies, and to describe features of the recombination landscape, including a massive (∼50 Mb) recombination desert (of virtually zero recombination) on the X chromosome that parallels a similar desert on the porcine X chromosome in both size and physical location.

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.

Genome-Wide Single-Nucleotide Polymorphisms Discovery and High-Density Genetic Map Construction in Cauliflower Using Specific-Locus Amplified Fragment Sequencing

Molecular markers and genetic maps play an important role in plant genomics and breeding studies. Cauliflower is an important and distinctive vegetable; however, very few molecular resources have been reported for this species. In this study, a novel, specific-locus amplified fragment (SLAF) sequencing strategy was employed for large-scale single nucleotide polymorphism (SNP) discovery and high-density genetic map construction in a double-haploid, segregating population of cauliflower. A total of 12.47 Gb raw data containing 77.92 M pair-end reads were obtained after processing and 6815 polymorphic SLAFs between the two parents were detected. The average sequencing depths reached 52.66-fold for the female parent and 49.35-fold for the male parent. Subsequently, these polymorphic SLAFs were used to genotype the population and further filtered based on several criteria to construct a genetic linkage map of cauliflower. Finally, 1776 high-quality SLAF markers, including 2741 SNPs, constituted the linkage map with average data integrity of 95.68%. The final map spanned a total genetic length of 890.01 cM with an average marker interval of 0.50 cM, and covered 364.9 Mb of the reference genome. The markers and genetic map developed in this study could provide an important foundation not only for comparative genomics studies within Brassica oleracea species but also for quantitative trait loci identification and molecular breeding of cauliflower.

Whole genome sequencing of turbot (Scophthalmus maximus; Pleuronectiformes): a fish adapted to demersal life

The turbot is a flatfish (Pleuronectiformes) with increasing commercial value, which has prompted active genomic research aimed at more efficient selection. Here we present the sequence and annotation of the turbot genome, which represents a milestone for both boosting breeding programmes and ascertaining the origin and diversification of flatfish. We compare the turbot genome with model fish genomes to investigate teleost chromosome evolution. We observe a conserved macrosyntenic pattern within Percomorpha and identify large syntenic blocks within the turbot genome related to the teleost genome duplication. We identify gene family expansions and positive selection of genes associated with vision and metabolism of membrane lipids, which suggests adaptation to demersal lifestyle and to cold temperatures, respectively. Our data indicate a quick evolution and diversification of flatfish to adapt to benthic life and provide clues for understanding their controversial origin. Moreover, we investigate the genomic architecture of growth, sex determination and disease resistance, key traits for understanding local adaptation and boosting turbot production, by mapping candidate genes and previously reported quantitative trait loci. The genomic architecture of these productive traits has allowed the identification of candidate genes and enriched pathways that may represent useful information for future marker-assisted selection in turbot.

New Insights on Eggplant/Tomato/Pepper Synteny and Identification of Eggplant and Pepper Orthologous QTL

Eggplant, pepper, and tomato are the most exploited berry-producing vegetables within the Solanaceae family. Their genomes differ in size, but each has 12 chromosomes which have undergone rearrangements causing a redistribution of loci. The genome sequences of all three species are available but differ in coverage, assembly quality and percentage of anchorage. Determining their syntenic relationship and QTL orthology will contribute to exploit genomic resources and genetic data for key agronomic traits. The syntenic analysis between tomato and pepper based on the alignment of 34,727 tomato CDS to the pepper genome sequence, identified 19,734 unique hits. The resulting synteny map confirmed the 14 inversions and 10 translocations previously documented, but also highlighted 3 new translocations and 4 major new inversions. Furthermore, each of the 12 chromosomes exhibited a number of rearrangements involving small regions of 0.5-0.7 Mbp. Due to high fragmentation of the publicly available eggplant genome sequence, physical localization of most eggplant QTL was not possible, thus, we compared the organization of the eggplant genetic map with the genome sequence of both tomato and pepper. The eggplant/tomato syntenic map confirmed all the 10 translocations but only 9 of the 14 known inversions; on the other hand, a newly detected inversion was recognized while another one was not confirmed. The eggplant/pepper syntenic map confirmed 10 translocations and 8 inversions already detected and suggested a putative new translocation. In order to perform the assessment of eggplant and pepper QTL orthology, the eggplant and pepper sequence-based markers located in their respective genetic map were aligned onto the pepper genome. GBrowse in pepper was used as reference platform for QTL positioning. A set of 151 pepper QTL were located as well as 212 eggplant QTL, including 76 major QTL (PVE ≥ 10%) affecting key agronomic traits. Most were confirmed to cluster in orthologous chromosomal regions. Our results highlight that the availability of genome sequences for an increasing number of crop species and the development of "ultra-dense" physical maps provide new and key tools for detailed syntenic and orthology studies between related plant species.

Mapping QTLs for Fusarium Head Blight Resistance in an Interspecific Wheat Population

Fusarium head blight (scab) is one of the most widespread and damaging diseases of wheat, causing grain yield and quality losses and production of harmful mycotoxins. Development of resistant varieties is hampered by lack of effective resistance sources in the tetraploid wheat primary gene pool. Here we dissected the genetic basis of resistance in a new durum wheat (Triticum turgidum ssp. durum) Recombinant inbred lines (RILs) population obtained by crossing an hexaploid resistant line and a durum susceptible cultivar. A total of 135 RILs were used for constituting a genetic linkage map and mapping loci for head blight incidence, severity, and disease-related plant morphological traits (plant height, spike compactness, and awn length). The new genetic map accounted for 4,366 single nucleotide polymorphism markers assembled in 52 linkage groups covering a total length of 4,227.37 cM. Major quantitative trait loci (QTL) for scab incidence and severity were mapped on chromosomes 2AS, 3AL, and 2AS, 2BS, 4BL, respectively. Plant height loci were identified on 3A, 3B, and 4B, while major QTL for ear compactness were found on 4A, 5A, 5B, 6A, and 7A. In this work, resistance to Fusarium was transferred from hexaploid to durum wheat, and correlations between the disease and morphological traits were assessed.

Mapping QTLs for Fusarium Head Blight Resistance in an Interspecific Wheat Population

Fusarium head blight (scab) is one of the most widespread and damaging diseases of wheat, causing grain yield and quality losses and production of harmful mycotoxins. Development of resistant varieties is hampered by lack of effective resistance sources in the tetraploid wheat primary gene pool. Here we dissected the genetic basis of resistance in a new durum wheat (Triticum turgidum ssp. durum) Recombinant inbred lines (RILs) population obtained by crossing an hexaploid resistant line and a durum susceptible cultivar. A total of 135 RILs were used for constituting a genetic linkage map and mapping loci for head blight incidence, severity, and disease-related plant morphological traits (plant height, spike compactness, and awn length). The new genetic map accounted for 4,366 single nucleotide polymorphism markers assembled in 52 linkage groups covering a total length of 4,227.37 cM. Major quantitative trait loci (QTL) for scab incidence and severity were mapped on chromosomes 2AS, 3AL, and 2AS, 2BS, 4BL, respectively. Plant height loci were identified on 3A, 3B, and 4B, while major QTL for ear compactness were found on 4A, 5A, 5B, 6A, and 7A. In this work, resistance to Fusarium was transferred from hexaploid to durum wheat, and correlations between the disease and morphological traits were assessed.

A High-Resolution InDel (Insertion-Deletion) Markers-Anchored Consensus Genetic Map Identifies Major QTLs Governing Pod Number and Seed Yield in Chickpea

Development and large-scale genotyping of user-friendly informative genome/gene-derived InDel markers in natural and mapping populations is vital for accelerating genomics-assisted breeding applications of chickpea with minimal resource expenses. The present investigation employed a high-throughput whole genome next-generation resequencing strategy in low and high pod number parental accessions and homozygous individuals constituting the bulks from each of two inter-specific mapping populations [(Pusa 1103 × ILWC 46) and (Pusa 256 × ILWC 46)] to develop non-erroneous InDel markers at a genome-wide scale. Comparing these high-quality genomic sequences, 82,360 InDel markers with reference to kabuli genome and 13,891 InDel markers exhibiting differentiation between low and high pod number parental accessions and bulks of aforementioned mapping populations were developed. These informative markers were structurally and functionally annotated in diverse coding and non-coding sequence components of genome/genes of kabuli chickpea. The functional significance of regulatory and coding (frameshift and large-effect mutations) InDel markers for establishing marker-trait linkages through association/genetic mapping was apparent. The markers detected a greater amplification (97%) and intra-specific polymorphic potential (58-87%) among a diverse panel of cultivated desi, kabuli, and wild accessions even by using a simpler cost-efficient agarose gel-based assay implicating their utility in large-scale genetic analysis especially in domesticated chickpea with narrow genetic base. Two high-density inter-specific genetic linkage maps generated using aforesaid mapping populations were integrated to construct a consensus 1479 InDel markers-anchored high-resolution (inter-marker distance: 0.66 cM) genetic map for efficient molecular mapping of major QTLs governing pod number and seed yield per plant in chickpea. Utilizing these high-density genetic maps as anchors, three major genomic regions harboring each of pod number and seed yield robust QTLs (15-28% phenotypic variation explained) were identified on chromosomes 2, 4, and 6. The integration of genetic and physical maps at these QTLs mapped on chromosomes scaled-down the long major QTL intervals into high-resolution short pod number and seed yield robust QTL physical intervals (0.89-2.94 Mb) which were essentially got validated in multiple genetic backgrounds of two chickpea mapping populations. The genome-wide InDel markers including natural allelic variants and genomic loci/genes delineated at major six especially in one colocalized novel congruent robust pod number and seed yield robust QTLs mapped on a high-density consensus genetic map were found most promising in chickpea. These functionally relevant molecular tags can drive marker-assisted genetic enhancement to develop high-yielding cultivars with increased seed/pod number and yield in chickpea.

Molecular Mapping and QTL for Expression Profiles of Flavonoid Genes in Brassica napus

Flavonoids are secondary metabolites that are extensively distributed in the plant kingdom and contribute to seed coat color formation in rapeseed. To decipher the genetic networks underlying flavonoid biosynthesis in rapeseed, we constructed a high-density genetic linkage map with 1089 polymorphic loci (including 464 SSR loci, 97 RAPD loci, 451 SRAP loci, and 75 IBP loci) using recombinant inbred lines (RILs). The map consists of 19 linkage groups and covers 2775 cM of the B. napus genome with an average distance of 2.54 cM between adjacent markers. We then performed expression quantitative trait locus (eQTL) analysis to detect transcript-level variation of 18 flavonoid biosynthesis pathway genes in the seeds of the 94 RILs. In total, 72 eQTLs were detected and found to be distributed among 15 different linkage groups that account for 4.11% to 52.70% of the phenotypic variance atrributed to each eQTL. Using a genetical genomics approach, four eQTL hotspots together harboring 28 eQTLs associated with 18 genes were found on chromosomes A03, A09, and C08 and had high levels of synteny with genome sequences of A. thaliana and Brassica species. Associated with the trans-eQTL hotspots on chromosomes A03, A09, and C08 were 5, 17, and 1 genes encoding transcription factors, suggesting that these genes have essential roles in the flavonoid biosynthesis pathway. Importantly, bZIP25, which is expressed specifically in seeds, MYC1, which controls flavonoid biosynthesis, and the R2R3-type gene MYB51, which is involved in the synthesis of secondary metabolites, were associated with the eQTL hotspots, and these genes might thus be involved in different flavonoid biosynthesis pathways in rapeseed. Hence, further studies of the functions of these genes will provide insight into the regulatory mechanism underlying flavonoid biosynthesis, and lay the foundation for elaborating the molecular mechanism of seed coat color formation in B. napus.

Construction of a high-density genetic map and mapping of a sex-linked locus for the brown alga Undaria pinnatifida (Phaeophyceae) based on large scale marker development by specific length amplified fragment (SLAF) sequencing

BACKGROUND

Undaria pinnatifida is an important economic brown alga in East Asian countries. However, its genetic and genomic information is very scarce, which hinders further research in this species. A high-density genetic map is a basic tool for fundamental and applied research such as discovery of functional genes and mapping of quantitative trait loci (QTL). In this study the recently developed specific length amplified fragment sequencing (SLAF-seq) technology was employed to construct a high-density genetic linkage map and locate a sex determining locus for U. pinnatifida.

RESULTS

A total of 28.06 Gb data including 140.31 M pair-end reads was obtained. After linkage analysis 4626 SLAF markers were mapped onto the genetic map. After adding the sex linked simple sequence repeat (SSR) marker [GenBank:AY738602.1], the final genetic map was 1816.28 cM long, consisting of 30 linkage groups with an average distance of 0.39 cM between adjacent markers. The length of LGs ranged from 20.12 to 106.95 cM. A major sex associated QTL was mapped to LG22 within a window starting at 29.01 cM and ending at 68.81 cM with a total of 68 SLAF markers. The SSR marker and five SLAF markers (Marker6556, 19020, 43089, 60771 and 26359) were identified as tightly sex-linked markers, as indicated by the absence of recombination between them and the sex phenotype. These markers were located at the position of 59.50 cM, which was supposed to be the sex determining region.

CONCLUSIONS

A high-density genetic linkage map was constructed using SLAF-seq technique and F1 gametophyte population for the first time in the economically important brown alga U. pinnatifida. For the first time, a major sex associated QTL suggesting a sex determining region was mapped to a single LG. This map will facilitate the further fundamental and applied research such as QTL mapping and map-based gene clone in U. pinnatifida and provide a reference for studies in other kelp species.

Integration of genetic and physical maps of the Primula vulgaris S locus and localization by chromosome in situ hybridization

Heteromorphic flower development in Primula is controlled by the S locus. The S locus genes, which control anther position, pistil length and pollen size in pin and thrum flowers, have not yet been characterized. We have integrated S-linked genes, marker sequences and mutant phenotypes to create a map of the P. vulgaris S locus region that will facilitate the identification of key S locus genes. We have generated, sequenced and annotated BAC sequences spanning the S locus, and identified its chromosomal location. We have employed a combination of classical genetics and three-point crosses with molecular genetic analysis of recombinants to generate the map. We have characterized this region by Illumina sequencing and bioinformatic analysis, together with chromosome in situ hybridization. We present an integrated genetic and physical map across the P. vulgaris S locus flanked by phenotypic and DNA sequence markers. BAC contigs encompass a 1.5-Mb genomic region with 1 Mb of sequence containing 82 S-linked genes anchored to overlapping BACs. The S locus is located close to the centromere of the largest metacentric chromosome pair. These data will facilitate the identification of the genes that orchestrate heterostyly in Primula and enable evolutionary analyses of the S locus.

Genes with Restricted Introgression in a Field Cricket (Gryllus firmus/Gryllus pennsylvanicus) Hybrid Zone Are Concentrated on the X Chromosome and a Single Autosome

Characterizing the extent of genomic differentiation between recently diverged lineages provides an important context for understanding the early stages of speciation. When such lineages form discrete hybrid zones, patterns of differential introgression allow direct estimates of which genome regions are likely involved in speciation and local adaptation. Here we use a backcross experimental design to construct a genetic linkage map for the field crickets Gryllus firmus and Gryllus pennsylvanicus, which interact in a well-characterized hybrid zone in eastern North America. We demonstrate that loci with major allele frequency differences between allopatric populations are not randomly distributed across the genome. Instead, most are either X-linked or map to a few small autosomal regions. Furthermore, the subset of those highly differentiated markers that exhibit restricted introgression across the cricket hybrid zone are also concentrated on the X chromosome (39 of 50 loci) and in a single 7-cM region of one autosome. Although the accumulation on the sex chromosome of genes responsible for postzygotic barriers is a well-known phenomenon, less attention has been given to the genomic distribution of genes responsible for prezygotic barriers. We discuss the implications of our results for speciation, both in the context of the role of sex chromosomes and also with respect to the likely causes of heterogeneous genomic divergence. Although we do not yet have direct evidence for the accumulation of ecological, behavioral, or fertilization prezygotic barrier genes on the X chromosome, faster-X evolution could make these barriers more likely to be X-linked.

Construction of a high-density genetic map using specific length amplified fragment markers and identification of a quantitative trait locus for anthracnose resistance in walnut (Juglans regia L.)

BACKGROUND

Walnut (Juglans regia, 2n = 32, approximately 606 Mb per 1C genome) is an economically important tree crop. Resistance to anthracnose, caused by Colletotrichum gloeosporioides, is a major objective of walnut genetic improvement in China. The recently developed specific length amplified fragment sequencing (SLAF-seq) is an efficient strategy that can obtain large numbers of markers with sufficient sequence information to construct high-density genetic maps and permits detection of quantitative trait loci (QTLs) for molecular breeding.

RESULTS

SLAF-seq generated 161.64 M paired-end reads. 153,820 SLAF markers were obtained, of which 49,174 were polymorphic. 13,635 polymorphic markers were sorted into five segregation types and 2,577 markers of them were used to construct genetic linkage maps: 2,395 of these fell into 16 linkage groups (LGs) for the female map, 448 markers for the male map, and 2,577 markers for the integrated map. Taking into account the size of all LGs, the marker coverage was 2,664.36 cM for the female map, 1,305.58 cM for the male map, and 2,457.82 cM for the integrated map. The average intervals between two adjacent mapped markers were 1.11 cM, 2.91 cM and 0.95 cM for three maps, respectively. 'SNP_only' markers accounted for 89.25% of the markers on the integrated map. Mapping markers contained 5,043 single nucleotide polymorphisms (SNPs) loci, which corresponded to two SNP loci per SLAF marker. According to the integrated map, we used interval mapping (Logarithm of odds, LOD > 3.0) to detect our quantitative trait. One QTL was detected for anthracnose resistance. The interval of this QTL ranged from 165.51 cM to 176.33 cM on LG14, and ten markers in this interval that were above the threshold value were considered to be linked markers to the anthracnose resistance trait. The phenotypic variance explained by each marker ranged from 16.2 to 19.9%, and their LOD scores varied from 3.22 to 4.04.

CONCLUSIONS

High-density genetic maps for walnut containing 16 LGs were constructed using the SLAF-seq method with an F1 population. One QTL for walnut anthracnose resistance was identified based on the map. The results will aid molecular marker-assisted breeding and walnut resistance genes identification.

Degeneration of the nonrecombining regions in the mating-type chromosomes of the anther-smut fungi

Dimorphic mating-type chromosomes in fungi are excellent models for understanding the genomic consequences of recombination suppression. Their suppressed recombination and reduced effective population size are expected to limit the efficacy of natural selection, leading to genomic degeneration. Our aim was to identify the sequences of the mating-type chromosomes (a1 and a2) of the anther-smut fungi and to investigate degeneration in their nonrecombining regions. We used the haploid a1 Microbotryum lychnidis-dioicae reference genome sequence. The a1 and a2 mating-type chromosomes were both isolated electrophoretically and sequenced. Integration with restriction-digest optical maps identified regions of recombination and nonrecombination in the mating-type chromosomes. Genome sequence data were also obtained for 12 other Microbotryum species. We found strong evidence of degeneration across the genus in the nonrecombining regions of the mating-type chromosomes, with significantly higher rates of nonsynonymous substitution (dN/dS) than in nonmating-type chromosomes or in recombining regions of the mating-type chromosomes. The nonrecombining regions of the mating-type chromosomes also showed high transposable element content, weak gene expression, and gene losses. The levels of degeneration did not differ between the a1 and a2 mating-type chromosomes, consistent with the lack of homogametic/heterogametic asymmetry between them, and contrasting with X/Y or Z/W sex chromosomes.

Fine mapping and characterization of Sr21, a temperature-sensitive diploid wheat resistance gene effective against the Puccinia graminis f. sp. tritici Ug99 race group

KEY MESSAGE

The diploid wheat stem rust resistance gene Sr21 confers temperature-sensitive resistance to isolates of the Ug99 group and maps to the middle of the long arm of chromosome 2A (m). A race of Puccinia graminis f. sp. tritici, the causal pathogen of stem rust of wheat, known as Ug99, and its variants, are virulent to plants carrying stem rust resistance genes currently deployed in most wheat cultivars worldwide. Therefore, identification, mapping and deployment of effective resistance genes are critical to reduce this threat. Resistance gene Sr21 identified in diploid wheat T. monococcum can be effective against races from the Ug99 race group, but both susceptible and partial resistant reactions have been reported in previous studies. To clarify this conflicting information we screened four monogenic lines with Sr21 and four susceptible controls with 16 Pgt isolates including five isolates of the Ug99 race group under three different temperatures and three different photoperiods. We observed that, temperature influences the interaction between monogenic lines with Sr21 and Ug99 race group isolates, and may be one source of previous inconsistencies. This result indicates that, although Sr21 confers partial resistance against Ug99, its effectiveness can be modulated by environmental conditions and should not be deployed alone. Using two large diploid wheat-mapping populations (total 3,788 F2 plants) we mapped Sr21 approximately 50 cM from the centromere on the long arm of chromosome 2A(m) within a 0.20 cM interval flanked by sequence-based markers FD527726 and EX594406. The closely linked markers identified in this study will be useful to reduce the T. monococcum segments introgressed into common wheat, accelerate Sr21 deployment in wheat breeding programs, and facilitate the map-based cloning of this gene.

Enrichment of an intraspecific genetic map of upland cotton by developing markers using parental RAD sequencing

RAD sequencing was performed using DH962 and Jimian5 as upland cotton mapping parents. Sequencing data for DH962 and Jimian5 were assembled into the genome sequences of ≈55.27 and ≈57.06 Mb, respectively. Analysing genome sequences of the two parents, 1,323 SSR, 3,838 insertion/deletion (InDel), and 9,366 single-nucleotide polymorphism (SNP) primer pairs were developed. All of the SSRs, 121 InDels, 441 SNPs, and other 6,747 primer pairs were screened in the two parents, and a total of 535 new polymorphic loci were identified. A genetic map including 1,013 loci was constructed using these results and 506 loci previously published for this population. Twenty-seven new QTLs for yield and fibre quality were identified, indicating that the efficiency of QTL detection was greatly improved by the increase in map density. Comparative genomics showed there to be considerable homology and collinearity between the A_T and A₂ genomes and between the D_T and D₅ genomes, although there were a few exchanges and introgressions among the chromosomes of the A₂ genome. Here, the development of markers using parental RAD sequencing was effective, and a high-density intraspecific genetic map was constructed. This map can be used for molecular marker-assisted selection in cotton.