Integrated consensus genetic and physical maps of flax (Linum usitatissimum L.)

High-Quality Genome Assembly and Genome-Wide Association Study of Male Sterility Provide Resources for Flax Improvement

Flax is an economic crop with a long history. It is grown worldwide and is mainly used for edible oil, industry, and textiles. Here, we reported a high-quality genome assembly for "Neiya No. 9", a popular variety widely grown in China. Combining PacBio long reads, Hi-C sequencing, and a genetic map reported previously, a genome assembly of 473.55 Mb was constructed, which covers ~94.7% of the flax genome. These sequences were anchored onto 15 chromosomes. The N50 lengths of the contig and scaffold were 0.91 Mb and 31.72 Mb, respectively. A total of 32,786 protein-coding genes were annotated, and 95.9% of complete BUSCOs were found. Through morphological and cytological observation, the male sterility of flax was considered dominant nuclear sterility. Through GWAS analysis, the gene LUSG00017705 (cysteine synthase gene) was found to be closest to the most significant SNP, and the expression level of this gene was significantly lower in male sterile plants than in fertile plants. Among the significant SNPs identified in the GWAS analysis, only two were located in the coding region, and these two SNPs caused changes in the protein encoded by LUSG00017565 (cysteine protease gene). It was speculated that these two genes may be related to male sterility in flax. This is the first time the molecular mechanism of male sterility in flax has been reported. The high-quality genome assembly and the male sterility genes revealed, provided a solid foundation for flax breeding.

Molecular Advances to Combat Different Biotic and Abiotic Stresses in Linseed (Linum usitatissimum L.): A Comprehensive Review

Flax, or linseed, is considered a "superfood", which means that it is a food with diverse health benefits and potentially useful bioactive ingredients. It is a multi-purpose crop that is prized for its seed oil, fibre, nutraceutical, and probiotic qualities. It is suited to various habitats and agro-ecological conditions. Numerous abiotic and biotic stressors that can either have a direct or indirect impact on plant health are experienced by flax plants as a result of changing environmental circumstances. Research on the impact of various stresses and their possible ameliorators is prompted by such expectations. By inducing the loss of specific alleles and using a limited number of selected varieties, modern breeding techniques have decreased the overall genetic variability required for climate-smart agriculture. However, gene banks have well-managed collectionns of landraces, wild linseed accessions, and auxiliary Linum species that serve as an important source of novel alleles. In the past, flax-breeding techniques were prioritised, preserving high yield with other essential traits. Applications of molecular markers in modern breeding have made it easy to identify quantitative trait loci (QTLs) for various agronomic characteristics. The genetic diversity of linseed species and the evaluation of their tolerance to abiotic stresses, including drought, salinity, heavy metal tolerance, and temperature, as well as resistance to biotic stress factors, viz., rust, wilt, powdery mildew, and alternaria blight, despite addressing various morphotypes and the value of linseed as a supplement, are the primary topics of this review.

Identification of new potential molecular actors related to fiber quality in flax through Omics

One of the biggest challenges for a more widespread utilization of plant fibers is to better understand the different molecular factors underlying the variability in fineness and mechanical properties of both elementary and scutched fibers. Accordingly, we analyzed genome-wide transcription profiling from bast fiber bearing tissues of seven different flax varieties (4 spring, 2 winter fiber varieties and 1 winter linseed) and identified 1041 differentially expressed genes between varieties, of which 97 were related to cell wall metabolism. KEGG analysis highlighted a number of different enriched pathways. Subsequent statistical analysis using Partial Least-Squares Discriminant Analysis showed that 73% of the total variance was explained by the first 3 X-variates corresponding to 56 differentially expressed genes. Calculation of Pearson correlations identified 5 genes showing a strong correlation between expression and morphometric data. Two-dimensional gel proteomic analysis on the two varieties showing the most discriminant and significant differences in morphometrics revealed 1490 protein spots of which 108 showed significant differential abundance. Mass spectrometry analysis successfully identified 46 proteins representing 32 non-redundant proteins. Statistical clusterization based on the expression level of genes corresponding to the 32 proteins showed clear discrimination into three separate clusters, reflecting the variety type (spring-/winter-fiber/oil). Four of the 32 proteins were also highly correlated with morphometric features. Examination of predicted functions for the 9 (5 + 4) identified genes highlighted lipid metabolism and senescence process. Calculation of Pearson correlation coefficients between expression data and retted fiber mechanical measurements (strength and maximum force) identified 3 significantly correlated genes. The genes were predicted to be connected to cell wall dynamics, either directly (Expansin-like protein), or indirectly (NAD(P)-binding Rossmann-fold superfamily protein). Taken together, our results have allowed the identification of molecular actors potentially associated with the determination of both in-planta fiber morphometrics, as well as ex-planta fiber mechanical properties, both of which are key parameters for elementary fiber and scutched fiber quality in flax.

Diversity of Linum genetic resources in global genebanks: from agro-morphological characterisation to novel genomic technologies - a review

Linseed or flaxseed is a well-recognized nutritional food with nutraceutical properties owing to high omega-3 fatty acid (α-Linolenic acid), dietary fiber, quality protein, and lignan content. Currently, linseed enjoys the status of a 'superfood' and its integration in the food chain as a functional food is evolving continuously as seed constituents are associated with lowering the risk of chronic ailments, such as heart diseases, cancer, diabetes, and rheumatoid arthritis. This crop also receives much attention in the handloom and textile sectors as the world's coolest fabric linen is made up of its stem fibers which are endowed with unique qualities such as luster, tensile strength, density, bio-degradability, and non-hazardous nature. Worldwide, major linseed growing areas are facing erratic rainfall and temperature patterns affecting flax yield, quality, and response to biotic stresses. Amid such changing climatic regimes and associated future threats, diverse linseed genetic resources would be crucial for developing cultivars with a broad genetic base for sustainable production. Furthermore, linseed is grown across the world in varied agro-climatic conditions; therefore it is vital to develop niche-specific cultivars to cater to diverse needs and keep pace with rising demands globally. Linseed genetic diversity conserved in global genebanks in the form of germplasm collection from natural diversity rich areas is expected to harbor genetic variants and thus form crucial resources for breeding tailored crops to specific culinary and industrial uses. Global genebank collections thus potentially play an important role in supporting sustainable agriculture and food security. Currently, approximately 61,000 germplasm accessions of linseed including 1,127 wild accessions are conserved in genebanks/institutes worldwide. This review analyzes the current status of Linum genetic resources in global genebanks, evaluation for agro-morphological traits, stress tolerance, and nutritional profiling to promote their effective use for sustainable production and nutrition enhancement in our modern diets.

RNA-Seq combined with population-level analysis reveals important candidate genes related to seed size in flax (Linum usitatissimum L.)

Seed size is a key determinant of crop yields. Understanding the regulatory mechanisms of seed size is beneficial for improving flax seed yield. In this study, the development of large flax seeds lagged behind that of small seeds, and 1,751 protein-coding genes were differentially expressed in early seeds, torpedo-stage embryos, and endosperms of CIli2719 and Z11637 using RNA sequencing. Homologous alignment revealed that 129 differentially expressed genes (DEGs) in flax were homologous with 71 known seed size-related genes in Arabidopsis thaliana and rice (Oryza sativa L.). These DEGs controlled seed size through multiple processes and factors, among which phytohormone pathways and transcription factors were the most important. Moreover, 54 DEGs were found to be associated with seed size and weight in a DEG-based association study. Nucleotide diversity (π) analysis of seed size-related candidate DEGs by homologous alignment and association analysis showed that the π values decreased significantly during flax acclimation from oil to fiber flax, suggesting that some seed size-related candidate genes were selected in this acclimation process. These results provide important resources and genetic foundation for further research on seed size regulation and seed improvement in flax.

Integrated omics approaches for flax improvement under abiotic and biotic stress: Current status and future prospects

Flax (Linum usitatissimum L.) or linseed is one of the important industrial crops grown all over the world for seed oil and fiber. Besides oil and fiber, flax offers a wide range of nutritional and therapeutic applications as a feed and food source owing to high amount of α-linolenic acid (omega-3 fatty acid), lignans, protein, minerals, and vitamins. Periodic losses caused by unpredictable environmental stresses such as drought, heat, salinity-alkalinity, and diseases pose a threat to meet the rising market demand. Furthermore, these abiotic and biotic stressors have a negative impact on biological diversity and quality of oil/fiber. Therefore, understanding the interaction of genetic and environmental factors in stress tolerance mechanism and identification of underlying genes for economically important traits is critical for flax improvement and sustainability. In recent technological era, numerous omics techniques such as genomics, transcriptomics, metabolomics, proteomics, phenomics, and ionomics have evolved. The advancements in sequencing technologies accelerated development of genomic resources which facilitated finer genetic mapping, quantitative trait loci (QTL) mapping, genome-wide association studies (GWAS), and genomic selection in major cereal and oilseed crops including flax. Extensive studies in the area of genomics and transcriptomics have been conducted post flax genome sequencing. Interestingly, research has been focused more for abiotic stresses tolerance compared to disease resistance in flax through transcriptomics, while the other areas of omics such as metabolomics, proteomics, ionomics, and phenomics are in the initial stages in flax and several key questions remain unanswered. Little has been explored in the integration of omic-scale data to explain complex genetic, physiological and biochemical basis of stress tolerance in flax. In this review, the current status of various omics approaches for elucidation of molecular pathways underlying abiotic and biotic stress tolerance in flax have been presented and the importance of integrated omics technologies in future research and breeding have been emphasized to ensure sustainable yield in challenging environments.

Fine-mapping of a putative glutathione S-transferase (GST) gene responsible for yellow seed colour in flax (Linum usitatissimum)

Objective

The brown seed coat colour of flax (Linum ustiatissimum) results from proanthocyanidin synthesis and accumulation. Glutathione S-transferases (GSTs), such as the TT19 protein in Arabidopsis, have been implicated in the transport of anthocyanidins during the synthesis of the brown proanthocyanidins. This study fine mapped the g allele responsible for yellow seed colour in S95407 and identified it as a putative mutated GST.

Results

We developed a Recombinant Inbred Line population with 320 lines descended from a cross between CDC Bethune (brown seed coat) and S95407 (yellow seed) and used molecular markers to fine map the G gene on Chromosome 6 (Chr 6). We used Next Generation Sequencing (NGS) to identify a putative GST was identified in this region and Sanger sequenced the gene from CDC Bethune, S95407 and other yellow seeded genotypes. The putative GST from S95407 had 13 SNPs encoding, including four non-synonymous amino acid changes, compared to the CDC Bethune reference sequence and the other genotypes. The GST encoded by Lus10019895 is a lambda-GST in contrast to the Arabidopsis TT19 which is a phi-GST.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-022-05964-x.

Genome-wide analysis of genomic imprinting in the endosperm and allelic variation in flax

Genomic imprinting is an epigenetic phenomenon that causes biased expression of maternally and paternally inherited alleles. In flowering plants, genomic imprinting predominantly occurs in the triploid endosperm and plays a vital role in seed development. In this study, we identified 248 candidate imprinted genes including 114 maternally expressed imprinted genes (MEGs) and 134 paternally expressed imprinted genes (PEGs) in flax (Linum usitatissimum L.) endosperm using deep RNA sequencing. These imprinted genes were neither clustered in specific chromosomal regions nor well conserved among flax and other plant species. MEGs tended to be expressed specifically in the endosperm, whereas the expression of PEGs was not tissue-specific. Imprinted single nucleotide polymorphisms differentiated 200 flax cultivars into the oil flax, oil-fiber dual purpose flax and fiber flax subgroups, suggesting that genomic imprinting contributed to intraspecific variation in flax. The nucleotide diversity of imprinted genes in the oil flax subgroup was significantly higher than that in the fiber flax subgroup, indicating that some imprinted genes underwent positive selection during flax domestication from oil flax to fiber flax. Moreover, imprinted genes that underwent positive selection were related to flax functions. Thirteen imprinted genes related to flax seed size and weight were identified using a candidate gene-based association study. Therefore, our study provides information for further exploration of the function and genomic variation of imprinted genes in the flax population.

Mapping QTLs for Alternaria blight in Linseed (Linum usitatissimum L.)

A SSR-based linkage map of linseed constructed based on 154 individual lines of F 2 mapping population derived from JRF-4 (disease-resistant) and Chambal (disease susceptible) genotypes. QTLs for Alternaria blight and other yield related traits identified. Out of 1720 SSRs, 216 SSRs were found polymorphic among the parents but due to segregation distortion 18 SSRs could not be used for linkage map construction. Total 191 SSRs were used to construct the linkage map and distributed in 15 linkage groups covering genome length of 1802.4 cM. A total of 10 QTLs were identified for 4 phenotypic traits including 4 QTLs for capsules/plant, 2 for capsule weight/plant, 2 for seed weight/plant and 2 for Alternaria blight resistance. This study laid a foundation for further validation and fine mapping with more advance and large set of marker for different QTL identification and marker-assisted selection in linseed.

The Toolbox for Fiber Flax Breeding: A Pipeline From Gene Expression to Fiber Quality

The goal of any plant breeding program is to improve quality of a target crop. Crop quality is a comprehensive feature largely determined by biological background. To improve the quality parameters of crops grown for the production of fiber, a functional approach was used to search for genes suitable for the effective manipulation of technical fiber quality. A key step was to identify genes with tissue and stage-specific pattern of expression in the developing fibers. In the current study, we investigated the relationship between gene expression evaluated in bast fibers of developing flax plants and the quality parameters of technical fibers measured after plant harvesting. Based on previously published transcriptomic data, two sets of genes that are upregulated in fibers during intrusive growth and tertiary cell wall deposition were selected. The expression level of the selected genes and fiber quality parameters were measured in fiber flax, linseed (oil flax) cultivars, and wild species that differ in type of yield and fiber quality parameters. Based on gene expression data, linear regression models for technical stem length, fiber tensile strength, and fiber flexibility were constructed, resulting in the identification of genes that have high potential for manipulating fiber quality. Chromosomal localization and single nucleotide polymorphism distribution in the selected genes were characterized for the efficacy of their use in conventional breeding and genome editing programs. Transcriptome-based selection is a highly targeted functional approach that could be used during the development of new cultivars of various crops.

Microsatellite marker development in the crop wild relative Linum bienne using genome skimming

PREMISE

Nuclear microsatellite markers were developed for Linum bienne, the sister species of the crop L. usitatissimum, to provide molecular genetic tools for the investigation of L. bienne genetic diversity and structure.

METHODS AND RESULTS

Fifty microsatellite loci were identified in L. bienne by means of genome skimming, and 44 loci successfully amplified. Of these, 16 loci evenly spread across the L. usitatissimum reference nuclear genome were used for genotyping six L. bienne populations. Excluding one monomorphic locus, the number of alleles per locus ranged from two to 12. Four out of six populations harbored private alleles. The levels of expected and observed heterozygosity were 0.076 to 0.667 and 0.000 to 1.000, respectively. All 16 loci successfully cross-amplified in L. usitatissimum.

CONCLUSIONS

The 16 microsatellite loci developed here can be used for population genetic studies in L. bienne, and 28 additional loci that successfully amplified are available for further testing.

High-throughput sequencing techniques to flax genetics and breeding

Flax (Linum usitatissimum L.) is an important oil and fiber crop. Using modern methods for flax breeding allows accelerating the introduction of some desired genes into the genotypes of future varieties. Today, an important condition for their creation is the development of research, that is based on next-generation sequencing (NGS). This review summarizes the results obtained using NGS in flax research. To date, a linkage map with a high marker density has been obtained for L. usitatissimum, which is already being used for a more efficient search for quantitative traits loci. Comparative studies of transcriptomes and miRNomes of flax under stress and in control conditions elucidated molecular-genetic mechanisms of abiotic and biotic stress responses. The very accurate model for genomic selection of flax resistant to pasmo was constructed. Based on NGS-sequencing also some details of the genus Linum evolution were clarified. The knowledge systematized in the review can be useful for researchers working in flax breeding and whereas fundamental interest for understanding the phylogenetic relationships within the genus Linum, the ontogenesis, and the mechanisms of the response of flax plants to various stress factors.

Mapping Quantitative Trait Loci onto Chromosome-Scale Pseudomolecules in Flax

Quantitative trait loci (QTL) are genomic regions associated with phenotype variation of quantitative traits. To date, a total of 313 QTL for 31 quantitative traits have been reported in 14 studies on flax. Of these, 200 QTL from 12 studies were identified based on genetic maps, the scaffold sequences, or the pre-released chromosome-scale pseudomolecules. Molecular markers for QTL identification differed across studies but the most used ones were simple sequence repeats (SSRs) or single nucleotide polymorphisms (SNPs). To uniquely map the SSR and SNP markers from different references onto the recently released chromosome-scale pseudomolecules, methods with several scripts and database files were developed to locate PCR- and SNP-based markers onto the same reference, co-locate QTL, and scan genome-wide candidate genes. Using these methods, 195 out of 200 QTL were successfully sorted onto the 15 flax chromosomes and grouped into 133 co-located QTL clusters; the candidate genes that co-located with these QTL clusters were also predicted. The methods and tools presented in this article facilitate marker re-mapping to a new reference, genome-wide QTL analysis, candidate gene scanning, and breeding applications in flax and other crops.

Genetic Mapping of Prince Rupprecht's Larch (Larix principis-rupprechtii Mayr) by Specific-Locus Amplified Fragment Sequencing

A high-density genetic linkage map is essential for plant genetics and genomics research. However, due to the deficiency of genomic data and high-quality molecular markers, no genetic map has been published for Prince Rupprecht’s larch (Larix principis-rupprechtii Mayr), a conifer species with high ecological and commercial value in northern China. In this study, 145 F1 progeny individuals from an intraspecific cross between two elite clones of L. principis-rupprechtii and their parents were employed to construct the first genetic map in this important tree species using specific-locus amplified fragment sequencing (SLAF-seq). After preprocessing, the procedure yielded 300.20 Gb of raw data containing 1501.22 M pair-end reads. A total of 324,352 SNP markers were detected and 122,785 of them were polymorphic, with a polymorphism rate of 37.86%. Ultimately, 6099 SNPs were organized into a genetic map containing 12 linkage groups, consistent with the haploid chromosome number of larch and most other species in the Pinaceae family. The linkage map spanned 2415.58 cM and covered 99.6% of the L. principis-rupprechtii genome with an average of 0.4 cM between adjacent markers. To the best of our knowledge, this map is the first reference map for L. principis-rupprechtii, as well as the densest one obtained in larch species thus far. The genome-wide SNPs and the high-resolution genetic map will provide a foundation for future quantitative trait loci mapping, map-based cloning, marker-assisted selection, comparative genomics, and genome sequence assembly for larch trees.

Flax tubulin and CesA superfamilies represent attractive and challenging targets for a variety of genome- and base-editing applications

Flax is both a valuable resource and an interesting model crop. Despite a long history of flax genetic transformation only one transgenic linseed cultivar has been so far registered in Canada. Implementation and use of the genome-editing technologies that allow site-directed modification of endogenous genes without the introduction of foreign genes might improve this situation. Besides its potential for boosting crop yields, genome editing is now one of the best tools for carrying out reverse genetics and it is emerging as an especially versatile tool for studying basic biology. A complex interplay between the flax tubulin family (6 α-, 14 β-, and 2 γ-tubulin genes), the building block of microtubules, and the CesA (15-16 genes), the subunit of the multimeric cellulose-synthesizing complex devoted to the oriented deposition of the cellulose microfibrils is fundamental for the biosynthesis of the cell wall. The role of the different members of each family in providing specificities to the assembled complexes in terms of structure, dynamics, activity, and interaction remains substantially obscure. Genome-editing strategies, recently shown to be successful in flax, can therefore be useful to unravel the issue of functional redundancy and provide evidence for specific interactions between different members of the tubulin and CesA gene families, in relation to different phase and mode of cell wall biosynthesis.

Resequencing 200 Flax Cultivated Accessions Identifies Candidate Genes Related to Seed Size and Weight and Reveals Signatures of Artificial Selection

Seed size and weight are key traits determining crop yield, which often undergo strongly artificial selection during crop domestication. Although seed sizes differ significantly between oil flax and fiber flax, the genetic basis of morphological differences and artificial selection characteristics in seed size remains largely unclear. Here we re-sequenced 200 flax cultivated accessions to generate a genome variation map based on chromosome assembly reference genomes. We provide evidence that oil flax group is the ancestor of cultivated flax, and the oil-fiber dual purpose group (OF) is the evolutionary intermediate transition state between oil and fiber flax. Genome-wide association studies (GWAS) were combined with LD Heatmap to identify candidate regions related to seed size and weight, then candidate genes were screened based on detailed functional annotations and estimation of nucleotide polymorphism effects. Using this strategy, we obtained 13 candidate genes related to seed size and weight. Selective sweeps analysis indicates human-involved selection of small seeds during the oil to fiber flax transition. Our study shows the existence of elite alleles for seed size and weight in flax germplasm and provides molecular insights into approaches for further improvement.

Genome-wide regulatory gene-derived SSRs reveal genetic differentiation and population structure in fiber flax genotypes

We designed a set of 580 simple sequence repeat markers; 506 from transcription factor-coding genes, and 74 from long non-coding RNAs and designated them as regulatory gene-derived simple sequence repeat (ReG-SSR) markers. From this set, we could anchor 559 ReG-SSR markers on 15 flax chromosomes with an average marker distance of 0.56 Mb. Thirty-one polymorphic ReG-SSR primers, amplifying SSR loci length of at least 20 bp were chosen from 134 screened primers. This primer set was used to characterize a diversity panel of 93 flax accessions. The panel included 33 accessions from India, including released varieties, dual-purpose lines and landraces, and 60 fiber flax accessions from the global core collection. Thirty-one ReG-SSR markers generated 76 alleles, with an average of 2.5 alleles per primer and a mean allele frequency of 0.77. These markers recorded 0.32 average gene diversity, 0.26 polymorphism information content and 1.35% null alleles. All the 31 ReG-SSR loci were found selectively neutral and showed no evidence of population reduction. A model-based clustering analysis separated the flax accessions into two sub-populations-Indian and global, with some accessions showing admixtures. The distinct clustering pattern of the Indian accessions compared to the global accessions, conforms to the principal coordinate analysis, genetic dissimilarity-based unweighted neighbor-joining tree and analysis of molecular variance. Fourteen flax accessions with 99.3% allelic richness were found optimum to adopt in breeding programs. In summary, the genome-wide ReG-SSR markers will serve as a functional marker resource for genetic and phenotypic relationship studies, marker-assisted selections, and provide a basis for selection of accessions from the Indian and global gene pool in fiber flax breeding programs.

Genome-Wide Association Analysis of Mucilage and Hull Content in Flax (Linum usitatissimum L.) Seeds

New flaxseed cultivars differing in seed mucilage content (MC) with low hull content (HC) represent an attractive option to simultaneously target the food and feed markets. Here, a genome-wide association study (GWAS) was conducted for MC and HC in 200 diverse flaxseed accessions genotyped with 1.7 million single nucleotide polymorphism (SNP) markers. The data obtained for MC and HC indicated a broad phenotypic variation and high (~70%) and a moderate (~49%) narrow sense heritability, respectively. MC and HC did not differ statistically between fiber and oil morphotypes, but yellow-seeded accessions had 2.7% less HC than brown-seeded ones. The genome-wide linkage disequilibrium (LD) decayed to r² = 0.1 at a physical distance of ~100 kb. Seven and four quantitative trait loci (QTL) were identified for MC and HC, respectively. Promising candidate genes identified include Linum usitatissimum orthologs of the Arabidopsis thaliana genes TRANSPARENT TESTA 8, SUBTILISIN-LIKE SERINE PROTEASE, GALACTUROSYL TRANSFERASE-LIKE 5, MUCILAGE-MODIFIED 4, AGAMOUS-LIKE MADS-BOX PROTEIN AGL62, GLYCOSYL HYDROLASE FAMILY 17, and UDP-GLUCOSE FLAVONOL 3-O-GLUCOSYLTRANSFERASE. These genes have been shown to play a role in mucilage synthesis and release, seed coat development and anthocyanin biosynthesis in A. thaliana. The favorable alleles will be useful in flaxseed breeding towards the goal of achieving the ideal MC and HC composition for food and feed by genomic-based breeding.

Consensus genetic linkage map construction and QTL mapping for plant height-related traits in linseed flax (Linum usitatissimum L.)

BACKGROUND

Flax is an important field crop that can be used for either oilseed or fiber production. Plant height and technical length are important characters for flax. For linseed flax, plants usually have a short technical length and plant height than those for fiber flax. As an important agronomical character for fiber and linseed flax, plant height is usually a selection target for breeding. However, because of limited technologies and methods available, there has been little research focused on discovering the molecular mechanism controlling plant height.

RESULTS

In this study, two related recombinant inbred line (RIL) populations developed from crosses of linseed and fiber parents were developed and phenotyped for plant height and technical length in four environments. A consensus linkage map based on two RIL populations was constructed using SNP markers generated by genotyping by sequencing (GBS) technology. A total of 4497 single nucleotide polymorphism (SNP) markers were included on 15 linkage groups with an average marker density of one marker every 2.71 cM. Quantitative trait locus (QTL) mapping analysis was performed for plant height and technical length using the two populations. A total of 19 QTLs were identified for plant height and technical length. For the MH population, eight plant height QTLs and seven technical length QTLs were identified, five of which were common QTLs for both traits. For the PH population, six plant height and three technical length QTLs were identified. By comparing the QTLs and candidate gene information in the two population, two common QTLs and three candidate genes were discovered.

CONCLUSIONS

This study provides a foundation for map-based cloning of QTLs and marker-assisted selection for plant height-related traits in linseed and fiber flax.

Consensus genetic linkage map construction and QTL mapping for plant height-related traits in linseed flax (Linum usitatissimum L.)

BACKGROUND

Flax is an important field crop that can be used for either oilseed or fiber production. Plant height and technical length are important characters for flax. For linseed flax, plants usually have a short technical length and plant height than those for fiber flax. As an important agronomical character for fiber and linseed flax, plant height is usually a selection target for breeding. However, because of limited technologies and methods available, there has been little research focused on discovering the molecular mechanism controlling plant height.

RESULTS

In this study, two related recombinant inbred line (RIL) populations developed from crosses of linseed and fiber parents were developed and phenotyped for plant height and technical length in four environments. A consensus linkage map based on two RIL populations was constructed using SNP markers generated by genotyping by sequencing (GBS) technology. A total of 4497 single nucleotide polymorphism (SNP) markers were included on 15 linkage groups with an average marker density of one marker every 2.71 cM. Quantitative trait locus (QTL) mapping analysis was performed for plant height and technical length using the two populations. A total of 19 QTLs were identified for plant height and technical length. For the MH population, eight plant height QTLs and seven technical length QTLs were identified, five of which were common QTLs for both traits. For the PH population, six plant height and three technical length QTLs were identified. By comparing the QTLs and candidate gene information in the two population, two common QTLs and three candidate genes were discovered.

CONCLUSIONS

This study provides a foundation for map-based cloning of QTLs and marker-assisted selection for plant height-related traits in linseed and fiber flax.

Genome-Wide Association Study and Selection Signatures Detect Genomic Regions Associated with Seed Yield and Oil Quality in Flax

A genome-wide association study (GWAS) was performed on a set of 260 lines which belong to three different bi-parental flax mapping populations. These lines were sequenced to an averaged genome coverage of 19× using the Illumina Hi-Seq platform. Phenotypic data for 11 seed yield and oil quality traits were collected in eight year/location environments. A total of 17,288 single nucleotide polymorphisms were identified, which explained more than 80% of the phenotypic variation for days to maturity (DTM), iodine value (IOD), palmitic (PAL), stearic, linoleic (LIO) and linolenic (LIN) acid contents. Twenty-three unique genomic regions associated with 33 quantitative trait loci (QTL) for the studied traits were detected, thereby validating four genomic regions previously identified. The 33 QTL explained 48–73% of the phenotypic variation for oil content, IOD, PAL, LIO and LIN but only 8–14% for plant height, DTM and seed yield. A genome-wide selective sweep scan for selection signatures detected 114 genomic regions that accounted for 7.82% of the flax pseudomolecule and overlapped with the 11 GWAS-detected genomic regions associated with 18 QTL for 11 traits. The results demonstrate the utility of GWAS combined with selection signatures for dissection of the genetic structure of traits and for pinpointing genomic regions for breeding improvement.

Chromosome-scale pseudomolecules refined by optical, physical and genetic maps in flax

Genomes of varying sizes have been sequenced with next-generation sequencing platforms. However, most reference sequences include draft unordered scaffolds containing chimeras caused by mis-scaffolding. A BioNano genome (BNG) optical map was constructed to improve the previously sequenced flax genome (Linum usitatissimum L., 2n = 30, about 373 Mb), which consisted of 3852 scaffolds larger than 1 kb and totalling 300.6 Mb. The high-resolution BNG map of cv. CDC Bethune totalled 317 Mb and consisted of 251 BNG contigs with an N50 of 2.15 Mb. A total of 622 scaffolds (286.6 Mb, 94.9%) aligned to 211 BNG contigs (298.6 Mb, 94.2%). Of those, 99 scaffolds, diagnosed to contain assembly errors, were refined into 225 new scaffolds. Using the newly refined scaffold sequences and the validated bacterial artificial chromosome-based physical map of CDC Bethune, the 211 BNG contigs were scaffolded into 94 super-BNG contigs (N50 of 6.64 Mb) that were further assigned to the 15 flax chromosomes using the genetic map. The pseudomolecules total about 316 Mb, with individual chromosomes of 15.6 to 29.4 Mb, and cover 97% of the annotated genes. Evidence from the chromosome-scale pseudomolecules suggests that flax has undergone palaeopolyploidization and mesopolyploidization events, followed by rearrangements and deletions or fusion of chromosome arms from an ancient progenitor with a haploid chromosome number of eight.

A genome-wide analysis of the flax (Linum usitatissimum L.) dirigent protein family: from gene identification and evolution to differential regulation

Identification of DIR encoding genes in flax genome. Analysis of phylogeny, gene/protein structures and evolution. Identification of new conserved motifs linked to biochemical functions. Investigation of spatio-temporal gene expression and response to stress. Dirigent proteins (DIRs) were discovered during 8-8' lignan biosynthesis studies, through identification of stereoselective coupling to afford either (+)- or (-)-pinoresinols from E-coniferyl alcohol. DIRs are also involved or potentially involved in terpenoid, allyl/propenyl phenol lignan, pterocarpan and lignin biosynthesis. DIRs have very large multigene families in different vascular plants including flax, with most still of unknown function. DIR studies typically focus on a small subset of genes and identification of biochemical/physiological functions. Herein, a genome-wide analysis and characterization of the predicted flax DIR 44-membered multigene family was performed, this species being a rich natural grain source of 8-8' linked secoisolariciresinol-derived lignan oligomers. All predicted DIR sequences, including their promoters, were analyzed together with their public gene expression datasets. Expression patterns of selected DIRs were examined using qPCR, as well as through clustering analysis of DIR gene expression. These analyses further implicated roles for specific DIRs in (-)-pinoresinol formation in seed-coats, as well as (+)-pinoresinol in vegetative organs and/or specific responses to stress. Phylogeny and gene expression analysis segregated flax DIRs into six distinct clusters with new cluster-specific motifs identified. We propose that these findings can serve as a foundation to further systematically determine functions of DIRs, i.e. other than those already known in lignan biosynthesis in flax and other species. Given the differential expression profiles and inducibility of the flax DIR family, we provisionally propose that some DIR genes of unknown function could be involved in different aspects of secondary cell wall biosynthesis and plant defense.

QTL Mapping of Fiber-Related Traits Based on a High-Density Genetic Map in Flax (Linum usitatissimum L.)

A genetic map is an important and valuable tool for quantitative trait locus (QTL) mapping, marker-assisted selection (MAS)-based breeding, and reference-assisted chromosome assembly. In this study, 112 F₂ plants from a cross between Linum usitatissimum L. "DIANE" and "NY17" and parent plants were subjected to high-throughput sequencing and specific-locus amplified fragment (SLAF) library construction. After preprocessing, 61.64 Gb of raw data containing 253.71 Mb paired-end reads, each 101 bp in length, were obtained. A total of 192,797 SLAFs were identified, of which 23,115 were polymorphic, with a polymorphism rate of 11.99%. Finally, 2,339 SLAFs were organized into a linkage map consisting of 15 linkage groups (LGs). The total length of the genetic map was 1483.25 centimorgans (cM) and the average distance between adjacent markers was 0.63 cM. Combined with flax chromosome-scale pseudomolecules, 12 QTLs associating with 6 flax fiber-related traits were mapped on the chromosomal scaffolds. This high-density genetic map of flax should serve as a foundation for flax fine QTL mapping, draft genome assembly, and MAS-guided breeding. Ultimately, the genomic regions identified in this research could potentially be valuable for improving flax fiber cultivars, as well as for identification of candidate genes involved in flax fiber formation processes.

SIGNIFICANCE STATEMENT

A high-density genetic map of flax was constructed, and QTLs were identified on the sequence scaffolds to be interrelated with fiber-related traits. The results of this study will not only provide a platform for gene/QTL fine mapping, map-based gene isolation, and molecular breeding for flax, but also provide a reference to help position sequence scaffolds on the physical map and assist in the process of assembling the flax genome sequence.

Construction of an SNP-based high-density linkage map for flax (Linum usitatissimum L.) using specific length amplified fragment sequencing (SLAF-seq) technology

Flax is an important crop for oil and fiber, however, no high-density genetic maps have been reported for this species. Specific length amplified fragment sequencing (SLAF-seq) is a high-resolution strategy for large scale de novo discovery and genotyping of single nucleotide polymorphisms. In this study, SLAF-seq was employed to develop SNP markers in an F2 population to construct a high-density genetic map for flax. In total, 196.29 million paired-end reads were obtained. The average sequencing depth was 25.08 in male parent, 32.17 in the female parent, and 9.64 in each F2 progeny. In total, 389,288 polymorphic SLAFs were detected, from which 260,380 polymorphic SNPs were developed. After filtering, 4,638 SNPs were found suitable for genetic map construction. The final genetic map included 4,145 SNP markers on 15 linkage groups and was 2,632.94 cM in length, with an average distance of 0.64 cM between adjacent markers. To our knowledge, this map is the densest SNP-based genetic map for flax. The SNP markers and genetic map reported in here will serve as a foundation for the fine mapping of quantitative trait loci (QTLs), map-based gene cloning and marker assisted selection (MAS) for flax.

QTL mapping and molecular characterization of the classical D locus controlling seed and flower color in Linum usitatissimum (flax)

The flowers of flax (linseed) are blue-hued, ephemeral and self-pollinating, and the seeds are typically brown. A century-old interest in natural yellow seed variants and a historical model point to recessive alleles in B1, D and G loci being responsible, but the functional aspects had remained unknown. Here, we characterized the "D" locus by quantitative trait loci (QTL) mapping and identified a FLAVONOID 3'5' HYDROXYLASE (F3'5'H) gene therein. It does not belong to the F3'5'H clade, but resembles biochemically characterized F3'Hs (flavonoid 3' hydroxylase) but without F3'H activity. The genome lacks other F3'H or F3'H-like genes. The apparent neo-functionalization from F3'H is associated with a Thr498 → Ser498 substitution in a substrate recognition site (SRS). The yellow seed and white flower phenotypes of the classical d mutation was found to be due to one nucleotide deletion that would truncate the deduced product and remove three of the six potential SRS, negatively impacting delphinidin synthesis. Delphinidin is sporadic in angiosperms, and flax has no known pollination syndrome(s) with functional pollinator group(s) that are attracted to blue flowers, raising questions on the acquisition of F3'5'H. The appearance of d allele is suggestive of the beginning of the loss of F3'5'H in this species.

Integration of Physical, Genetic, and Cytogenetic Mapping Data for Cellulose Synthase (CesA) Genes in Flax (Linum usitatissimum L.)

Flax, Linum usitatissimum L., is a valuable multi-purpose plant, and currently, its genome is being extensively investigated. Nevertheless, mapping of genes in flax genome is still remaining a challenging task. The cellulose synthase (CesA) multigene family involving in the process of cellulose synthesis is especially important for metabolism of this fiber crop. For the first time, fluorescent in situ hybridization (FISH)-based chromosomal localization of the CesA conserved fragment (KF011584.1), 5S, and 26S rRNA genes was performed in landrace, oilseed, and fiber varieties of L. usitatissimum. Intraspecific polymorphism in chromosomal distribution of KF011584.1 and 5S DNA loci was revealed, and the generalized chromosome ideogram was constructed. Using BLAST analysis, available data on physical/genetic mapping and also whole-genome sequencing of flax, localization of KF011584.1, 45S, and 5S rRNA sequences on genomic scaffolds, and their anchoring to the genetic map were conducted. The alignment of the results of FISH and BLAST analyses indicated that KF011584.1 fragment revealed on chromosome 3 could be anchored to linkage group (LG) 11. The common LG for 45S and 5S rDNA was not found probably due to the polymorphic localization of 5S rDNA on chromosome 1. Our findings indicate the complexity of integration of physical, genetic, and cytogenetic mapping data for multicopy gene families in plants. Nevertheless, the obtained results can be useful for future progress in constructing of integrated physical/genetic/cytological maps in L. usitatissimum which are essential for flax breeding.

Development of Novel SSR Markers for Flax (Linum usitatissimum L.) Using Reduced-Representation Genome Sequencing

Flax (Linum usitatissimum L.) is a major fiber and oil yielding crop grown in northeastern China. Identification of flax molecular markers is a key step toward improving flax yield and quality via marker-assisted breeding. Simple sequence repeat (SSR) markers, which are based on genomic structural variation, are considered the most valuable type of genetic marker for this purpose. In this study, we screened 1574 microsatellites from Linum usitatissimum L. obtained using reduced representation genome sequencing (RRGS) to systematically identify SSR markers. The resulting set of microsatellites consisted mainly of trinucleotide (56.10%) and dinucleotide (35.23%) repeats, with each motif consisting of 5-8 repeats. We then evaluated marker sensitivity and specificity based on samples of 48 flax isolates obtained from northeastern China. Using the new SSR panel, the results demonstrated that fiber flax and oilseed flax varieties clustered into two well separated groups. The novel SSR markers developed in this study show potential value for selection of varieties for use in flax breeding programs.

Genetic diversity, population structure and association analysis in linseed (Linum usitatissimum L.)

The present investigation aimed to explore the level of genetic diversity, determine the population structure in a larger set of germplasm of linseed using microsatellite marker and identify linked markers through association mapping. A total of 168 accessions of linseed were evaluated for major agro-economic traits and SSRs markers deployed for diversity assessment. A total of 337 alleles were amplified by 50 SSRs ranging from 2 to 13 with an average of 6.74 ± 2.8 alleles per loci. The neighbor joining based clustering grouped all the accessions into three major clusters that were also confirmed by scatter plot of PCoA. While model based clustering determined four sub-populations (K = 4). Further, analysis of molecular variance analysis considering three population showed that maximum variation (79%) was within the population. We identified one putative SSR marker (Lu_3043) linked with days to 50% flowering through both GLM and MLM analysis of association mapping. The results of this preliminary study revealed genetic diversity, population structure in linseed and linked marker which could be utilized in future breeding program.

QTL for fatty acid composition and yield in linseed (Linum usitatissimum L.)

The combined SSR-SNP map and 20 QTL for agronomic and quality traits will assist in marker assisted breeding as well as map-based cloning of key genes in linseed. Flax is an important nutraceutical crop mostly because it is a rich source of omega-3 fatty acids and antioxidant compounds. Canada is the largest producer and exporter of oilseed flax (or linseed), creating a growing need to improve crop productivity and quality. In this study, a genetic map was constructed based on selected 329 single nucleotide polymorphic markers and 362 simple sequence repeat markers using a recombinant inbred line population of 243 individuals from a cross between the Canadian varieties CDC Bethune and Macbeth. The genetic map consisted of 15 linkage groups comprising 691 markers with an average marker density of one marker every 1.9 cM. A total of 20 quantitative trait loci (QTL) were identified corresponding to 14 traits. Three QTL each for oleic acid and stearic acid, two QTL each for linoleic acid and iodine value and one each for palmitic acid, linolenic acid, oil content, seed protein, cell wall, straw weight, thousand seed weight, seeds per boll, yield and days to maturity were identified. The QTL for cell wall, straw weight, seeds per boll, yield and days to maturity all co-located on linkage group 4. Analysis of the candidate gene regions underlying the QTL identified proteins involved in cell wall and fibre synthesis, fatty acid biosynthesis as well as their metabolism and yield component traits. This study provides the foundation for assisting in map-based cloning of the QTL and marker assisted selection of a wide range of quality and agronomic traits in linseed and potentially fibre flax.

QTL for fatty acid composition and yield in linseed (Linum usitatissimum L.)

The combined SSR-SNP map and 20 QTL for agronomic and quality traits will assist in marker assisted breeding as well as map-based cloning of key genes in linseed. Flax is an important nutraceutical crop mostly because it is a rich source of omega-3 fatty acids and antioxidant compounds. Canada is the largest producer and exporter of oilseed flax (or linseed), creating a growing need to improve crop productivity and quality. In this study, a genetic map was constructed based on selected 329 single nucleotide polymorphic markers and 362 simple sequence repeat markers using a recombinant inbred line population of 243 individuals from a cross between the Canadian varieties CDC Bethune and Macbeth. The genetic map consisted of 15 linkage groups comprising 691 markers with an average marker density of one marker every 1.9 cM. A total of 20 quantitative trait loci (QTL) were identified corresponding to 14 traits. Three QTL each for oleic acid and stearic acid, two QTL each for linoleic acid and iodine value and one each for palmitic acid, linolenic acid, oil content, seed protein, cell wall, straw weight, thousand seed weight, seeds per boll, yield and days to maturity were identified. The QTL for cell wall, straw weight, seeds per boll, yield and days to maturity all co-located on linkage group 4. Analysis of the candidate gene regions underlying the QTL identified proteins involved in cell wall and fibre synthesis, fatty acid biosynthesis as well as their metabolism and yield component traits. This study provides the foundation for assisting in map-based cloning of the QTL and marker assisted selection of a wide range of quality and agronomic traits in linseed and potentially fibre flax.

Resolution of genetic map expansion caused by excess heterozygosity in plant recombinant inbred populations

Recombinant inbred populations of many plant species exhibit more heterozygosity than expected under the Mendelian model of segregation. This segregation distortion causes the overestimation of recombination frequencies and consequent genetic map expansion. Here we build upon existing genetic models of differential zygotic viability to model a heterozygote fitness term and calculate expected genotypic proportions in recombinant inbred populations propagated by selfing. We implement this model using the existing open-source genetic map construction code base for R/qtl to estimate recombination fractions. Finally, we show that accounting for excess heterozygosity in a sorghum recombinant inbred mapping population shrinks the genetic map by 213 cM (a 13% decrease corresponding to 4.26 fewer recombinations per meiosis). More accurate estimates of linkage benefit linkage-based analyses used in the identification and utilization of causal genetic variation.

The potential of pale flax as a source of useful genetic variation for cultivated flax revealed through molecular diversity and association analyses

Pale flax (Linum bienne Mill.) is the wild progenitor of cultivated flax (Linum usitatissimum L.) and represents the primary gene pool to broaden its genetic base. Here, a collection of 125 pale flax accessions and the Canadian flax core collection of 407 accessions were genotyped using 112 genome-wide simple sequence repeat markers and phenotyped for nine traits with the aim of conducting population structure, molecular diversity and association mapping analyses. The combined population structure analysis identified two well-supported major groups corresponding to pale and cultivated flax. The L. usitatissimum convar. crepitans accessions most closely resembled its wild progenitor, both having dehiscent capsules. The unbiased Nei's genetic distance (0.65) confirmed the strong genetic differentiation between cultivated and pale flax. Similar levels of genetic diversity were observed in both species, albeit 430 (48 %) of pale flax alleles were unique, in agreement with their high genetic differentiation. Significant associations were identified for seven and four traits in pale and cultivated flax, respectively. Favorable alleles with potentially positive effect to improve yield through yield components were identified in pale flax. The allelic frequencies of markers associated with domestication-related traits such as capsular dehiscence indicated directional selection with the most common alleles in pale flax being absent or rare in cultivated flax and vice versa. Our results demonstrated that pale flax is a potential source of novel variation to improve multiple traits in cultivated flax and that association mapping is a suitable approach to screening pale flax germplasm to identify favorable quantitative trait locus alleles.

Diversity analysis in Indian genotypes of linseed (Linum usitatissimum L.) using AFLP markers

AFLP fingerprinting of 45 Indian genotypes of linseed was carried out to determine the genetic relationship among them. Sixteen primer combinations produced 1142 fragments with 1129 as polymorphic and 13 as monomorphic fragments. Polymorphic fragments varied from 44 (E-ACA/M-CTA) to 94 (E-AGC/M-CAC) with an average of 70.6 fragments per primer combination. The frequency of polymorphism varied from 93.7% to 100% with an average of 98.8% across all the genotypes. The PIC value ranged from 0.19 to 0.31 with an average of 0.23 per primer combination. The primer pair E-AGC/M-CAC showed the maximum PIC value (0.31) followed by E-AGC/M-CAG (0.29), E-AAC/M-CAG (0.26) and E-AGC/M-CTA (0.25). Resolving power (RP) and marker index (MI) varied from 13.73 to 43.50 and 8.81 to 28.91 respectively. The Jaccard's similarity coefficient varied from 0.16 to 0.57 with an average of 0.26 ± 0.05. The maximum genetic similarities (57%) were detected between genotypes Him Alsi-1 and Him Alsi-2, followed by Him Alsi-1 and GS41 and GS41 and LC-54. The genotypes R-552, Himani, RKY-14, Meera, Indira Alsi-32 and Suyog were found to be more divergent genotypes. The NJ clustering grouped all the 45 genotypes into three major clusters. In general the genotypes of cluster III had high oil content and those of cluster I had low oil content. At the population level, within population variance was much higher than between populations variance.

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.

Association mapping of seed quality traits using the Canadian flax (Linum usitatissimum L.) core collection

The identification of stable QTL for seed quality traits by association mapping of a diverse panel of linseed accessions establishes the foundation for assisted breeding and future fine mapping in linseed. Linseed oil is valued for its food and non-food applications. Modifying its oil content and fatty acid (FA) profiles to meet market needs in a timely manner requires clear understanding of their quantitative trait loci (QTL) architectures, which have received little attention to date. Association mapping is an efficient approach to identify QTL in germplasm collections. In this study, we explored the quantitative nature of seed quality traits including oil content (OIL), palmitic acid, stearic acid, oleic acid, linoleic acid (LIO) linolenic acid (LIN) and iodine value in a flax core collection of 390 accessions assayed with 460 microsatellite markers. The core collection was grown in a modified augmented design at two locations over 3 years and phenotypic data for all seven traits were obtained from all six environments. Significant phenotypic diversity and moderate to high heritability for each trait (0.73-0.99) were observed. Most of the candidate QTL were stable as revealed by multivariate analyses. Nine candidate QTL were identified, varying from one for OIL to three for LIO and LIN. Candidate QTL for LIO and LIN co-localized with QTL previously identified in bi-parental populations and some mapped nearby genes known to be involved in the FA biosynthesis pathway. Fifty-eight percent of the QTL alleles were absent (private) in the Canadian cultivars suggesting that the core collection possesses QTL alleles potentially useful to improve seed quality traits. The candidate QTL identified herein will establish the foundation for future marker-assisted breeding in linseed.

Genomic regions underlying agronomic traits in linseed (Linum usitatissimum L.) as revealed by association mapping

The extreme climate of the Canadian Prairies poses a major challenge to improve yield. Although it is possible to breed for yield per se, focusing on yield-related traits could be advantageous because of their simpler genetic architecture. The Canadian flax core collection of 390 accessions was genotyped with 464 simple sequence repeat markers, and phenotypic data for nine agronomic traits including yield, bolls per area, 1,000 seed weight, seeds per boll, start of flowering, end of flowering, plant height, plant branching, and lodging collected from up to eight environments was used for association mapping. Based on a mixed model (principal component analysis (PCA) + kinship matrix (K)), 12 significant marker-trait associations for six agronomic traits were identified. Most of the associations were stable across environments as revealed by multivariate analyses. Statistical simulation for five markers associated with 1000 seed weight indicated that the favorable alleles have additive effects. None of the modern cultivars carried the five favorable alleles and the maximum number of four observed in any accessions was mostly in breeding lines. Our results confirmed the complex genetic architecture of yield-related traits and the inherent difficulties associated with their identification while illustrating the potential for improvement through marker-assisted selection.

PT-Flax (phenotyping and TILLinG of flax): development of a flax (Linum usitatissimum L.) mutant population and TILLinG platform for forward and reverse genetics

BACKGROUND

Flax (Linum usitatissimum L.) is an economically important fiber and oil crop that has been grown for thousands of years. The genome has been recently sequenced and transcriptomics are providing information on candidate genes potentially related to agronomically-important traits. In order to accelerate functional characterization of these genes we have generated a flax EMS mutant population that can be used as a TILLinG (Targeting Induced Local Lesions in Genomes) platform for forward and reverse genetics.

RESULTS

A population of 4,894 M2 mutant seed families was generated using 3 different EMS concentrations (0.3%, 0.6% and 0.75%) and used to produce M2 plants for subsequent phenotyping and DNA extraction. 10,839 viable M2 plants (4,033 families) were obtained and 1,552 families (38.5%) showed a visual developmental phenotype (stem size and diameter, plant architecture, flower-related). The majority of these families showed more than one phenotype. Mutant phenotype data are organised in a database and can be accessed and searched at UTILLdb (http://urgv.evry.inra.fr/UTILLdb). Preliminary screens were also performed for atypical fiber and seed phenotypes. Genomic DNA was extracted from 3,515 M2 families and eight-fold pooled for subsequent mutant detection by ENDO1 nuclease mis-match cleavage. In order to validate the collection for reverse genetics, DNA pools were screened for two genes coding enzymes of the lignin biosynthesis pathway: Coumarate-3-Hydroxylase (C3H) and Cinnamyl Alcohol Dehydrogenase (CAD). We identified 79 and 76 mutations in the C3H and CAD genes, respectively. The average mutation rate was calculated as 1/41 Kb giving rise to approximately 9,000 mutations per genome. Thirty-five out of the 52 flax cad mutant families containing missense or codon stop mutations showed the typical orange-brown xylem phenotype observed in CAD down-regulated/mutant plants in other species.

CONCLUSIONS

We have developed a flax mutant population that can be used as an efficient forward and reverse genetics tool. The collection has an extremely high mutation rate that enables the detection of large numbers of independant mutant families by screening a comparatively low number of M2 families. The population will prove to be a valuable resource for both fundamental research and the identification of agronomically-important genes for crop improvement in flax.

The LuWD40-1 gene encoding WD repeat protein regulates growth and pollen viability in flax (Linum Usitatissimum L.)

As a crop, flax holds significant commercial value for its omega-3 rich oilseeds and stem fibres. Canada is the largest producer of linseed but there exists scope for significant yield improvements. Implementation of mechanisms such as male sterility can permit the development of hybrids to assist in achieving this goal. Temperature sensitive male sterility has been reported in flax but the leakiness of this system in field conditions limits the production of quality hybrid seeds. Here, we characterized a 2,588 bp transcript differentially expressed in male sterile lines of flax. The twelve intron gene predicted to encode a 368 amino acid protein has five WD40 repeats which, in silico, form a propeller structure with putative nucleic acid and histone binding capabilities. The LuWD40-1 protein localized to the nucleus and its expression increased during the transition and continued through the vegetative stages (seed, etiolated seedling, stem) while the transcript levels declined during reproductive development (ovary, anthers) and embryonic morphogenesis of male fertile plants. Knockout lines for LuWD40-1 in flax failed to develop shoots while overexpression lines showed delayed growth phenotype and were male sterile. The non-viable flowers failed to open and the pollen grains from these flowers were empty. Three independent transgenic lines overexpressing the LuWD40-1 gene had ∼80% non-viable pollen, reduced branching, delayed flowering and maturity compared to male fertile genotypes. The present study provides new insights into a male sterility mechanism present in flax.

Genetic characterization of a core collection of flax (Linum usitatissimum L.) suitable for association mapping studies and evidence of divergent selection between fiber and linseed types

BACKGROUND

Flax is valued for its fiber, seed oil and nutraceuticals. Recently, the fiber industry has invested in the development of products made from linseed stems, making it a dual purpose crop. Simultaneous targeting of genomic regions controlling stem fiber and seed quality traits could enable the development of dual purpose cultivars. However, the genetic diversity, population structure and linkage disequilibrium (LD) patterns necessary for association mapping (AM) have not yet been assessed in flax because genomic resources have only recently been developed. We characterized 407 globally distributed flax accessions using 448 microsatellite markers. The data was analyzed to assess the suitability of this core collection for AM. Genomic scans to identify candidate genes selected during the divergent breeding process of fiber flax and linseed were conducted using the whole genome shotgun sequence of flax.

RESULTS

Combined genetic structure analysis assigned all accessions to two major groups with six sub-groups. Population differentiation was weak between the major groups (F(ST) = 0.094) and for most of the pairwise comparisons among sub-groups. The molecular coancestry analysis indicated weak relatedness (mean = 0.287) for most individual pairs. Abundant genetic diversity was observed in the total panel (5.32 alleles per locus), and some sub-groups showed a high proportion of private alleles. The average genome-wide LD (r²) was 0.036, with a relatively fast decay of 1.5 cM. Genomic scans between fiber flax and linseed identified candidate genes involved in cell-wall biogenesis/modification, xylem identity and fatty acid biosynthesis congruent with genes previously identified in flax and other plant species.

CONCLUSIONS

Based on the abundant genetic diversity, weak population structure and relatedness and relatively fast LD decay, we concluded that this core collection is suitable for AM studies targeting multiple agronomic and quality traits aiming at the improvement of flax as a true dual purpose crop. Our genomic scans provide the first insights into candidate regions affected by divergent selection in flax. In combination with AM, genomic scans have the ability to increase the power to detect loci influencing complex traits.

Outlier Loci and Selection Signatures of Simple Sequence Repeats (SSRs) in Flax (Linum usitatissimum L.)

Genomic microsatellites (gSSRs) and expressed sequence tag-derived SSRs (EST-SSRs) have gained wide application for elucidating genetic diversity and population structure in plants. Both marker systems are assumed to be selectively neutral when making demographic inferences, but this assumption is rarely tested. In this study, three neutrality tests were assessed for identifying outlier loci among 150 SSRs (85 gSSRs and 65 EST-SSRs) that likely influence estimates of population structure in three differentiated flax sub-populations (F_ST = 0.19). Moreover, the utility of gSSRs, EST-SSRs, and the combined sets of SSRs was also evaluated in assessing genetic diversity and population structure in flax. Six outlier loci were identified by at least two neutrality tests showing footprints of balancing selection. After removing the outlier loci, the STRUCTURE analysis and the dendrogram topology of EST-SSRs improved. Conversely, gSSRs and combined SSRs results did not change significantly, possibly as a consequence of the higher number of neutral loci assessed. Taken together, the genetic structure analyses established the superiority of gSSRs to determine the genetic relationships among flax accessions, although the combined SSRs produced the best results. Genetic diversity parameters did not differ statistically (P > 0.05) between gSSRs and EST-SSRs, an observation partially explained by the similar number of repeat motifs. Our study provides new insights into the ability of gSSRs and EST-SSRs to measure genetic diversity and structure in flax and confirms the importance of testing for the occurrence of outlier loci to properly assess natural and breeding populations, particularly in studies considering only few loci.

Tinkering with meiosis

Meiosis is at the heart of Mendelian heredity. Recently, much progress has been made in the understanding of this process, in various organisms. In the last 15 years, the functional characterization of numerous genes involved in meiosis has dramatically deepened our knowledge of key events, including recombination, the cell cycle, and chromosome distribution. Through a constantly advancing tool set and knowledge base, a number of advances have been made that will allow manipulation of meiosis from a plant breeding perspective. This review focuses on the aspects of meiosis that can be tinkered with to create and propagate new varieties. We would like to dedicate this review to the memory of Simon W. Chan (1974-2012) (http://www.plb.ucdavis.edu/labs/srchan/).

Genome wide SNP discovery in flax through next generation sequencing of reduced representation libraries

BACKGROUND

Flax (Linum usitatissimum L.) is a significant fibre and oilseed crop. Current flax molecular markers, including isozymes, RAPDs, AFLPs and SSRs are of limited use in the construction of high density linkage maps and for association mapping applications due to factors such as low reproducibility, intense labour requirements and/or limited numbers. We report here on the use of a reduced representation library strategy combined with next generation Illumina sequencing for rapid and large scale discovery of SNPs in eight flax genotypes. SNP discovery was performed through in silico analysis of the sequencing data against the whole genome shotgun sequence assembly of flax genotype CDC Bethune. Genotyping-by-sequencing of an F6-derived recombinant inbred line population provided validation of the SNPs.

RESULTS

Reduced representation libraries of eight flax genotypes were sequenced on the Illumina sequencing platform resulting in sequence coverage ranging from 4.33 to 15.64X (genome equivalents). Depending on the relatedness of the genotypes and the number and length of the reads, between 78% and 93% of the reads mapped onto the CDC Bethune whole genome shotgun sequence assembly. A total of 55,465 SNPs were discovered with the largest number of SNPs belonging to the genotypes with the highest mapping coverage percentage. Approximately 84% of the SNPs discovered were identified in a single genotype, 13% were shared between any two genotypes and the remaining 3% in three or more. Nearly a quarter of the SNPs were found in genic regions. A total of 4,706 out of 4,863 SNPs discovered in Macbeth were validated using genotyping-by-sequencing of 96 F6 individuals from a recombinant inbred line population derived from a cross between CDC Bethune and Macbeth, corresponding to a validation rate of 96.8%.

CONCLUSIONS

Next generation sequencing of reduced representation libraries was successfully implemented for genome-wide SNP discovery from flax. The genotyping-by-sequencing approach proved to be efficient for validation. The SNP resources generated in this work will assist in generating high density maps of flax and facilitate QTL discovery, marker-assisted selection, phylogenetic analyses, association mapping and anchoring of the whole genome shotgun sequence.