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

A SLAF-based high-density genetic map construction and genetic architecture of thermotolerant traits in maize (Zea mays L.)

The leaf scorching trait at flowering is a crucial thermosensitive phenotype in maize under high temperature stress (HS), yet the genetic basis of this trait remains poorly understood. In this study, we genotyped a 254 RIL-F2:8 population, derived from the leaf scorch-free parental inbred line Abe2 and the leaf scorching maternal inbred line B73, using the specific-locus amplified fragment sequencing (SLAF-seq) method. A total of 10,112 polymorphic SLAF markers were developed, and a high-density genetic map with a total length of 1,475.88 cM was constructed. The average sequencing depth of the parents was 55.23X, and that of the progeny was 12.53X. Then, we identified a total of 16 QTLs associated with thermotolerant traits at flowering, of which four QTLs of leaf scorching damage (LS) were distributed on chromosomes 1 (qLS1), 2 (qLS2.1, qLS2.2) and 3 (qLS3), which could explain 19.73% of phenotypic variation. Combining one qLS1 locus with QTL-seq results led to the identification of 6 candidate genes. Expression experiments and sequence variation indicated that Zm00001d033328, encoding N-acetyl-gamma-glutamyl-phosphate reductase, was the most likely candidate gene controlling thermotolerant traits at flowering. In summary, the high-density genetic map and genetic basis of thermotolerant traits lay a critical foundation for mapping other complex traits and identifying the genes associated with thermotolerant traits in maize.

Genetic dissection of thousand-seed weight in linseed (Linum usitatissimum L.) using multi-locus genome-wide association study

Flaxseed/linseed is an important oilseed crop having applications in the food, nutraceutical, and paint industry. Seed weight is one of the most crucial determinants of seed yield in linseed. Here, quantitative trait nucleotides (QTNs) associated with thousand-seed weight (TSW) have been identified using multi-locus genome-wide association study (ML-GWAS). Field evaluation was carried out in five environments in multi-year-location trials. SNP genotyping information of the AM panel of 131 accessions comprising 68,925 SNPs was employed for ML-GWAS. From the six ML-GWAS methods employed, five methods helped identify a total of 84 unique significant QTNs for TSW. QTNs identified in ≥ 2 methods/environments were designated as stable QTNs. Accordingly, 30 stable QTNs have been identified for TSW accounting up to 38.65% trait variation. Alleles with positive effect on trait were analyzed for 12 strong QTNs with r ² ≥ 10.00%, which showed significant association of specific alleles with higher trait value in three or more environments. A total of 23 candidate genes have been identified for TSW, which included B3 domain-containing transcription factor, SUMO-activating enzyme, protein SCARECROW, shaggy-related protein kinase/BIN2, ANTIAUXIN-RESISTANT 3, RING-type E3 ubiquitin transferase E4, auxin response factors, WRKY transcription factor, and CBS domain-containing protein. In silico expression analysis of candidate genes was performed to validate their possible role in different stages of seed development process. The results from this study provide significant insight and elevate our understanding on genetic architecture of TSW trait in linseed.

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.

Development and Complex Application of Methods for the Identification of Mutations in the FAD3A and FAD3B Genes Resulting in the Reduced Content of Linolenic Acid in Flax Oil

Flax is grown worldwide for seed and fiber production. Linseed varieties differ in their oil composition and are used in pharmaceutical, food, feed, and industrial production. The field of application primarily depends on the content of linolenic (LIN) and linoleic (LIO) fatty acids. Inactivating mutations in the FAD3A and FAD3B genes lead to a decrease in the LIN content and an increase in the LIO content. For the identification of the three most common low-LIN mutations in flax varieties (G-to-A in exon 1 of FAD3A substituting tryptophan with a stop codon, C-to-T in exon 5 of FAD3A leading to arginine to a stop codon substitution, and C-to-T in exon 2 of FAD3B resulting in histidine to tyrosine substitution), three approaches were proposed: (1) targeted deep sequencing, (2) high resolution melting (HRM) analysis, (3) cleaved amplified polymorphic sequences (CAPS) markers. They were tested on more than a thousand flax samples of various types and showed promising results. The proposed approaches can be used in marker-assisted selection to choose parent pairs for crosses, separate heterogeneous varieties into biotypes, and select genotypes with desired homozygous alleles of the FAD3A and FAD3B genes at the early stages of breeding for the effective development of varieties with a particular LIN and LIO content, as well as in basic studies of the molecular mechanisms of fatty acid synthesis in flax seeds to select genotypes adequate to the tasks.

Genetic Determinants of Fiber-Associated Traits in Flax Identified by Omics Data Integration

In this paper, we explore potential genetic factors in control of flax phenotypes associated with fiber by mining a collection of 306 flax accessions from the Federal Research Centre of the Bast Fiber Crops, Torzhok, Russia. In total, 11 traits were assessed in the course of 3 successive years. A genome-wide association study was performed for each phenotype independently using six different single-locus models implemented in the GAPIT3 R package. Moreover, we applied a multivariate linear mixed model implemented in the GEMMA package to account for trait correlations and potential pleiotropic effects of polymorphisms. The analyses revealed a number of genomic variants associated with different fiber traits, implying the complex and polygenic control. All stable variants demonstrate a statistically significant allelic effect across all 3 years of the experiment. We tested the validity of the predicted variants using gene expression data available for the flax fiber studies. The results shed new light on the processes and pathways associated with the complex fiber traits, while the pinpointed candidate genes may be further used for marker-assisted selection.

Identification of QTNs Associated With Flowering Time, Maturity, and Plant Height Traits in Linum usitatissimum L. Using Genome-Wide Association Study

Early flowering, maturity, and plant height are important traits for linseed to fit in rice fallows, for rainfed agriculture, and for economically viable cultivation. Here, Multi-Locus Genome-Wide Association Study (ML-GWAS) was undertaken in an association mapping panel of 131 accessions, genotyped using 68,925 SNPs identified by genotyping by sequencing approach. Phenotypic evaluation data of five environments comprising 3 years and two locations were used. GWAS was performed for three flowering time traits including days to 5%, 50%, and 95% flowering, days to maturity, and plant height by employing five ML-GWAS methods: FASTmrEMMA, FASTmrMLM, ISIS EM-BLASSO, mrMLM, and pLARmEB. A total of 335 unique QTNs have been identified for five traits across five environments. 109 QTNs were stable as observed in ≥2 methods and/or environments, explaining up to 36.6% phenotypic variance. For three flowering time traits, days to maturity, and plant height, 53, 30, and 27 stable QTNs, respectively, were identified. Candidate genes having roles in flower, pollen, embryo, seed and fruit development, and xylem/phloem histogenesis have been identified. Gene expression of candidate genes for flowering and plant height were studied using transcriptome of an early maturing variety Sharda (IC0523807). The present study unravels QTNs/candidate genes underlying complex flowering, days to maturity, and plant height traits in linseed.

Comparative Phenotyping of Two Commonly Used Chlamydomonas reinhardtii Background Strains: CC-1690 (21gr) and CC-5325 (The CLiP Mutant Library Background)

The unicellular green alga Chlamydomonas reinhardtii is an excellent model organism to investigate many essential cellular processes in photosynthetic eukaryotes. Two commonly used background strains of Chlamydomonas are CC-1690 and CC-5325. CC-1690, also called 21gr, has been used for the Chlamydomonas genome project and several transcriptome analyses. CC-5325 is the background strain for the Chlamydomonas Library Project (CLiP). Photosynthetic performance in CC-5325 has not been evaluated in comparison with CC-1690. Additionally, CC-5325 is often considered to be cell-wall deficient, although detailed analysis is missing. The circadian rhythms in CC-5325 are also unclear. To fill these knowledge gaps and facilitate the use of the CLiP mutant library for various screens, we performed phenotypic comparisons between CC-1690 and CC-5325. Our results showed that CC-5325 grew faster heterotrophically in dark and equally well in mixotrophic liquid medium as compared to CC-1690. CC-5325 had lower photosynthetic efficiency and was more heat-sensitive than CC-1690. Furthermore, CC-5325 had an intact cell wall which had comparable integrity to that in CC-1690 but appeared to have reduced thickness. Additionally, CC-5325 could perform phototaxis, but could not maintain a sustained circadian rhythm of phototaxis as CC1690 did. Finally, in comparison to CC-1690, CC-5325 had longer cilia in the medium with acetate but slower swimming speed in the medium without nitrogen and acetate. Our results will be useful for researchers in the Chlamydomonas community to choose suitable background strains for mutant analysis and employ the CLiP mutant library for genome-wide mutant screens under appropriate conditions, especially in the areas of photosynthesis, thermotolerance, cell wall, and circadian rhythms.

QTL mapping for growth-related traits by constructing the first genetic linkage map in Simao pine

BACKGROUND

Simao pine is one of the primary economic tree species for resin and timber production in southwest China. The exploitation and utilization of Simao pine are constrained by the relatively lacking of genetic information. Construction a fine genetic linkage map and detecting quantitative trait locis (QTLs) for growth-related traits is a prerequisite section of Simao Pine's molecular breeding program.

RESULTS

In our study, a high-resolution Simao pine genetic map employed specific locus amplified fragment sequencing (SLAF-seq) technology and based on an F₁ pseudo-testcross population has been constructed. There were 11,544 SNPs assigned to 12 linkage groups (LGs), and the total length of the map was 2,062.85 cM with a mean distance of 0.37 cM between markers. According to the phenotypic variation analysis for three consecutive years, a total of seventeen QTLs for four traits were detected. Among 17 QTLs, there were six for plant height (Dh.16.1, Dh16.2, Dh17.1, Dh18.1-3), five for basal diameter (Dbd.17.1-5), four for needle length (Dnl17.1-3, Dnl18.1) and two for needle diameter (Dnd17.1 and Dnd18.1) respectively. These QTLs individually explained phenotypic variance from 11.0-16.3%, and the logarithm of odds (LOD) value ranged from 2.52 to 3.87.

CONCLUSIONS

In our study, a fine genetic map of Simao pine applied the technology of SLAF-seq has been constructed for the first time. Based on the map, a total of 17 QTLs for four growth-related traits were identified. It provides helpful information for genomic studies and marker-assisted selection (MAS) in Simao pine.

Genes Associated with the Flax Plant Type (Oil or Fiber) Identified Based on Genome and Transcriptome Sequencing Data

As a result of the breeding process, there are two main types of flax (Linum usitatissimum L.) plants. Linseed is used for obtaining seeds, while fiber flax is used for fiber production. We aimed to identify the genes associated with the flax plant type, which could be important for the formation of agronomically valuable traits. A search for polymorphisms was performed in genes involved in the biosynthesis of cell wall components, lignans, fatty acids, and ion transport based on genome sequencing data for 191 flax varieties. For 143 of the 424 studied genes (4CL, C3'H, C4H, CAD, CCR, CCoAOMT, COMT, F5H, HCT, PAL, CTL, BGAL, ABC, HMA, DIR, PLR, UGT, TUB, CESA, RGL, FAD, SAD, and ACT families), one or more polymorphisms had a strong correlation with the flax type. Based on the transcriptome sequencing data, we evaluated the expression levels for each flax type-associated gene in a wide range of tissues and suggested genes that are important for the formation of linseed or fiber flax traits. Such genes were probably subjected to the selection press and can determine not only the traits of seeds and stems but also the characteristics of the root system or resistance to stresses at a particular stage of development, which indirectly affects the ability of flax plants to produce seeds or fiber.

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.

Regression based fast multi-trait genome-wide QTL analysis

Multivariate simple interval mapping (SIM) is one of the most popular approaches for multiple quantitative trait locus (QTL) analysis. Both maximum likelihood (ML) and least squares (LS) multivariate regression (MVR) are widely used methods for multi-trait SIM. ML-based MVR (MVR-ML) is an expectation maximization (EM) algorithm based iterative and complex time-consuming approach. Although the LS-based MVR (MVR-LS) approach is not an iterative process, the calculation of likelihood ratio (LR) statistic in MVR-LS is also a time-consuming complex process. We have introduced a new approach (called FastMtQTL) for multi-trait QTL analysis based on the assumption of multivariate normal distribution of phenotypic observations. Our proposed method can identify almost the same QTL positions as those identified by the existing methods. Moreover, the proposed method takes comparatively less computation time because of the simplicity in the calculation of LR statistic by this method. In the proposed method, LR statistic is calculated only using the sample variance-covariance matrix of phenotypes and the conditional probability of QTL genotype given the marker genotypes. This improvement in computation time is advantageous when the numbers of phenotypes and individuals are larger, and the markers are very dense resulting in a QTL mapping with a bigger dataset.

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.