Rapid identification of lettuce seed germination mutants by bulked segregant analysis and whole genome sequencing

Advance in the Thermoinhibition of Lettuce (Lactuca sativa L.) Seed Germination

Thermoinhibition refers to the inability of seeds to germinate when inhibited by high temperatures, but when environmental conditions return to normal, the seeds are able to germinate rapidly again, which is different from thermodormancy. Meanwhile, with global warming, the effect of the thermoinhibition phenomenon on the yield and quality of crops in agricultural production is becoming common. Lettuce, as a horticultural crop sensitive to high temperature, is particularly susceptible to the effects of thermoinhibition, resulting in yield reduction. Therefore, it is crucial to elucidate the intrinsic mechanism of action of thermoinhibition in lettuce seeds. This review mainly outlines several factors affecting thermoinhibition of lettuce seed germination, including endosperm hardening, alteration of endogenous or exogenous phytohormone concentrations, action of photosensitizing pigments, production and inhibition of metabolites, maternal effects, genetic expression, and other physical and chemical factors. Finally, we also discuss the challenges and potential of lettuce seed germination thermoinhibition research. The purpose of this study is to provide theoretical support for future research on lettuce seed germination thermoinhibition, and with the aim of revealing the mechanisms and effects behind lettuce seed thermoinhibition. This will enable the identification of more methods to alleviate seed thermoinhibition or the development of superior heat-tolerant lettuce seeds.

Advancements in genetic techniques and functional genomics for enhancing crop traits and agricultural sustainability

Genetic variability is essential for the development of new crop varieties with economically beneficial traits. The traits can be inherited from wild relatives or induced through mutagenesis. Novel genetic elements can then be identified and new gene functions can be predicted. In this study, forward and reverse genetics approaches were described, in addition to their applications in modern crop improvement programs and functional genomics. By using heritable phenotypes and linked genetic markers, forward genetics searches for genes by using traditional genetic mapping and allele frequency estimation. Despite recent advances in sequencing technology, omics and computation, genetic redundancy remains a major challenge in forward genetics. By analyzing close-related genes, we will be able to dissect their functional redundancy and predict possible traits and gene activity patterns. In addition to these predictions, sophisticated reverse gene editing tools can be used to verify them, including TILLING, targeted insertional mutagenesis, gene silencing, gene targeting and genome editing. By using gene knock-down, knock-up and knock-out strategies, these tools are able to detect genetic changes in cells. In addition, epigenome analysis and editing enable the development of novel traits in existing crop cultivars without affecting their genetic makeup by increasing epiallelic variants. Our understanding of gene functions and molecular dynamics of various biological phenomena has been revised by all of these findings. The study also identifies novel genetic targets in crop species to improve yields and stress tolerances through conventional and non-conventional methods. In this article, genetic techniques and functional genomics are specifically discussed and assessed for their potential in crop improvement.

Mapping and Screening of Candidate Gene Regulating the Biomass Yield of Sorghum (Sorghum bicolor L.)

Biomass yield is one of the important traits of sorghum, which is greatly affected by leaf morphology. In this study, a lobed-leaf mutant (sblob) was screened and identified, and its F2 inbred segregating line was constructed. Subsequently, MutMap and whole-genome sequencing were employed to identify the candidate gene (sblob1), the locus of which is Sobic.003G010300. Pfam and homologous analysis indicated that sblob1 encodes a Cytochrome P450 protein and plays a crucial role in the plant serotonin/melatonin biosynthesis pathway. Structural and functional changes in the sblob1 protein were elucidated. Hormone measurements revealed that sblob1 regulates both leaf morphology and sorghum biomass through regulation of the melatonin metabolic pathway. These findings provide valuable insights for further research and the enhancement of breeding programs, emphasizing the potential to optimize biomass yield in sorghum cultivation.

Seed color in lettuce is determined by the LsTT2, LsCHS, and Ls2OGD genes from the flavonoid biosynthesis pathway

KEY MESSAGE

The mutated LsTT2 and Ls2OGD genes are responsible for white seeds and yellow seeds in lettuce, respectively. Three LsCHS genes are involved in the biosynthesis of flavonoid in seed coats. Lettuce seeds have several different colors, including black, yellow, and white. The genetic mechanisms underlying color variations of lettuce seeds remain unknown. We used genome-wide association studies (GWAS) and map-based cloning approaches to clone genes controlling the color of lettuce seeds. LsTT2, which encodes an R2R3-MYB transcription factor and is homologous to the TT2 gene in Arabidopsis, was shown to be the causal gene for the variation of black and white seeds in lettuce. A point mutation leads to the lack of stop codon in the LsTT2 transcript, resulting in white seeds. Knockout of the LsTT2 gene converted black seeds to white seeds. The locus controlling yellow seeds was mapped to Chromosome 2. Knockout of two 2-oxoglutarate-dependent dioxygenases (2OGD) genes from the candidate region converted black seeds to yellow seeds, suggesting that these two 2OGD proteins catalyze the conversion of yellow metabolites to black metabolites. We also showed that three LsCHS genes from the candidate region are associated with flavonoid biosynthesis in seeds. Knockout mutants of the three LsCHS genes decreased color intensity. This study provides new insights into the regulation of flavonoid biosynthesis in plants.

Role of Brassica orphan gene BrLFM on leafy head formation in Chinese cabbage (Brassica rapa)

Brassica orphan gene BrFLM, identified by two allelic mutants, was involved in leafy head formation in Chinese cabbage. Leafy head formation is a unique agronomic trait of Chinese cabbage that determines its yield and quality. In our previous study, an EMS mutagenesis Chinese cabbage mutant library was constructed using the heading Chinese cabbage double haploid (DH) line FT as the wild-type. Here, we screened two extremely similar leafy head deficiency mutants lfm-1 and lfm-2 with geotropic growth leaves from the library to investigate the gene(s) related to leafy head formation. Reciprocal crossing results showed that these two mutants were allelic. We utilized lfm-1 to identify the mutant gene(s). Genetic analysis showed that the mutated trait was controlled by a single nuclear gene Brlfm. Mutmap analysis showed that Brlfm was located on chromosome A05, and BraA05g012440.3C or BraA05g021450.3C were the candidate gene. Kompetitive allele-specific PCR analysis eliminated BraA05g012440.3C from the candidates. Sanger sequencing identified an SNP from G to A at the 271st nucleotide on BraA05g021450.3C. The sequencing of lfm-2 detected another non-synonymous SNP (G to A) located at the 266st nucleotide on BraA05g021450.3C, which verified its function on leafy head formation. We blasted BraA05g021450.3C on database and found that it belongs to a Brassica orphan gene encoding an unknown 13.74 kDa protein, named BrLFM. Subcellular localization showed that BrLFM was located in the nucleus. These findings reveal that BrLFM is involved in leafy head formation in Chinese cabbage.

Identification of an epicuticular wax crystal deficiency gene Brwdm1 in Chinese cabbage (Brassica campestris L. ssp. pekinensis)

Introduction

The cuticle wax covering the plant surface is a whitish hydrophobic protective barrier in Chinese cabbage, and the epicuticular wax crystal deficiency normally has higher commodity value for a tender texture and glossy appearance. Herein, two allelic epicuticular wax crystal deficiency mutants, wdm1 and wdm7, were obtained from the EMS mutagenesis population of a Chinese cabbage DH line 'FT'.

Methods

The cuticle wax morphology was observed by Cryo-scanning electron microscopy (Cryo-SEM) and the composition of wax was determined by GC-MS. The candidate mutant gene was found by MutMap and validated by KASP. The function of candidate gene was verified by allelic variation.

Results

The mutants had fewer wax crystals and lower leaf primary alcohol and ester content. Genetic analysis revealed that the epicuticular wax crystal deficiency phenotype was controlled by a recessive nuclear gene, named Brwdm1. MutMap and KASP analyses indicated that BraA01g004350.3C, encoding an alcohol-forming fatty acyl-CoA reductase, was the candidate gene for Brwdm1. A SNP 2,113,772 (C to T) variation in the 6^th exon of Brwdm1 in wdm1 led to the 262^nd amino acid substitution from threonine (T) to isoleucine (I), which existed in a rather conserved site among the amino acid sequences from Brwdm1 and its homologs. Meanwhile, the substitution changed the three-dimensional structure of Brwdm1. The SNP 2,114,994 (G to A) in the 10^th exon of Brwdm1 in wdm7 resulted in the change of the 434^th amino acid from valine (V) to isoleucine (I), which occurred in the STERILE domain. KASP genotyping showed that SNP 2,114,994 was co-segregated with glossy phenotype. Compared with the wild type, the relative expression of Brwdm1 was significantly decreased in the leaves, flowers, buds and siliques of wdm1.

Discussion

These results indicated that Brwdm1 was indispensable for the wax crystals formation and its mutation resulted in glossy appearance in Chinese cabbage.

A Locus Controlling Leaf Rolling Degree in Wheat under Drought Stress Identified by Bulked Segregant Analysis

Drought stress frequently occurs, which seriously restricts the production of wheat (Triticum aestivum L.). Leaf rolling is a typical physiological phenomenon of plants during drought stress. To understand the genetic mechanism of wheat leaf rolling, we constructed an F₂ segregating population by crossing the slight-rolling wheat cultivar “Aikang 58” (AK58) with the serious-rolling wheat cultivar ″Zhongmai 36″ (ZM36). A combination of bulked segregant analysis (BSA) with Wheat 660K SNP Array was used to identify molecular markers linked to leaf rolling degree. A major locus for leaf rolling degree under drought stress was detected on chromosome 7A. We named this locus LEAF ROLLING DEGREE 1 (LERD1), which was ultimately mapped to a region between 717.82 and 720.18 Mb. Twenty-one genes were predicted in this region, among which the basic helix-loop-helix (bHLH) transcription factor TraesCS7A01G543300 was considered to be the most likely candidate gene for LERD1. The TraesCS7A01G543300 is highly homologous to the Arabidopsis ICE1 family proteins ICE/SCREAM, SCREAM2 and bHLH093, which control stomatal initiation and development. Two nucleotide variation sites were detected in the promoter region of TraesCS7A01G543300 between the two wheat cultivars. Gene expression assays indicated that TraesCS7A01G543300 was higher expressed in AK58 seedlings than that of ZM36. This research discovered a candidate gene related to wheat leaf rolling under drought stress, which may be helpful for understanding the leaf rolling mechanism and molecular breeding in wheat.

BrCPS1 Function in Leafy Head Formation Was Verified by Two Allelic Mutations in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

The formation of the leafy heads of Chinese cabbage is an important agricultural factor because it directly affects yield. In this study, we identified two allelic non-heading mutants, nhm4-1 and nhm4-2, from an ethyl methanesulfonate mutagenic population of a heading Chinese cabbage double haploid line "FT." Using MutMap, Kompetitive Allele-Specific PCR genotyping, and map-based cloning, we found that BraA09g001440.3C was the causal gene for the mutants. BraA09g001440.3C encodes an ent-copalyl diphosphate synthase 1 involved in gibberellin biosynthesis. A single non-synonymous SNP in the seventh and fourth exons of BraA09g001440.3C was responsible for the nhm4-1 and nhm4-2 mutant phenotypes, respectively. Compared with the wild-type "FT," the gibberellin content in the mutant leaves was significantly reduced. Both mutants showed a tendency to form leafy heads after exogenous GA₃ treatment. The two non-heading mutants and the work presented herein demonstrate that gibberellin is related to leafy head formation in Chinese cabbage.

Detection of a major QTL and development of KASP markers for seed weight by combining QTL-seq, QTL-mapping and RNA-seq in peanut

Combining QTL-seq, QTL-mapping and RNA-seq identified a major QTL and candidate genes, which contributed to the development of KASP markers and understanding of molecular mechanisms associated with seed weight in peanut. Seed weight, as an important component of seed yield, is a significant target of peanut breeding. However, relatively little is known about the quantitative trait loci (QTLs) and candidate genes associated with seed weight in peanut. In this study, three major QTLs on chromosomes A05, B02, and B06 were determined by applying the QTL-seq approach in a recombinant inbred line (RIL) population. Based on conventional QTL-mapping, these three QTL regions were successfully narrowed down through newly developed single nucleotide polymorphism (SNP) and simple sequence repeat markers. Among these three QTL regions, qSWB06.3 exhibited stable expression, contributing mainly to phenotypic variance across environments. Furthermore, differentially expressed genes (DEGs) were identified at the three seed developmental stages between the two parents of the RIL population. It was found that the DEGs were widely distributed in the ubiquitin-proteasome pathway, the serine/threonine-protein pathway, signal transduction of hormones and transcription factors. Notably, DEGs at the early stage were mostly involved in regulating cell division, whereas DEGs at the middle and late stages were primarily involved in cell expansion during seed development. The expression patterns of candidate genes related to seed weight in qSWB06.3 were investigated using quantitative real-time PCR. In addition, the allelic diversity of qSWB06.3 was investigated in peanut germplasm accessions. The marker Ah011475 has higher efficiency for discriminating accessions with different seed weights, and it would be useful as a diagnostic marker in marker-assisted breeding. This study provided insights into the genetic and molecular mechanisms of seed weight in peanut.

PpMYB36 Encodes a MYB-Type Transcription Factor That Is Involved in Russet Skin Coloration in Pear (Pyrus pyrifolia)

Fruit color is one of the most important external qualities of pear (Pyrus pyrifolia) fruits. However, the mechanisms that control russet skin coloration in pear have not been well characterized. Here, we explored the molecular mechanisms that determine the russet skin trait in pear using the F₁ population derived from a cross between russet skin ('Niitaka') and non-russet skin ('Dangshansu') cultivars. Pigment measurements indicated that the lignin content in the skin of the russet pear fruits was greater than that in the non-russet pear skin. Genetic analysis revealed that the phenotype of the russet skin pear is associated with an allele of the PpRus gene. Using bulked segregant analysis combined with the genome sequencing (BSA-seq), we identified two simple sequence repeat (SSR) marker loci linked with the russet-colored skin trait in pear. Linkage analysis showed that the PpRus locus maps to the scaffold NW_008988489.1: 53297-211921 on chromosome 8 in the pear genome. In the mapped region, the expression level of LOC103929640 was significantly increased in the russet skin pear and showed a correlation with the increase of lignin content during the ripening period. Genotyping results demonstrated that LOC103929640 encoding the transcription factor MYB36 is the causal gene for the russet skin trait in pear. Particularly, a W-box insertion at the PpMYB36 promoter of russet skin pears is essential for PpMYB36-mediated regulation of lignin accumulation and russet coloration in pear. Overall, these results show that PpMYB36 is involved in the regulation of russet skin trait in pear.

Mapping of a Pale Green Mutant Gene and Its Functional Verification by Allelic Mutations in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Leaves are the main organ for photosynthesis, and variations in leaf color affect photosynthesis and plant biomass formation. We created two similar whole-plant pale green mutants (pem1 and pem2) from the double haploid (DH) Chinese cabbage line "FT" through ethyl methanesulfonate (EMS) mutagenesis of seeds. Photosynthetic pigment contents and net photosynthetic rates were significantly lower in the mutants than in the wild-type "FT," and the chloroplast thylakoid endomembrane system was poor. Genetic analysis showed that the mutated phenotypes of pem1 and pem2 were caused by a single nuclear gene. Allelism tests showed that pem1 and pem2 were alleles. We mapped Brpem1 to a 64.25 kb region on chromosome A10, using BSR-Seq and map-based cloning of 979 F₂ recessive individuals. Whole-genome re-sequencing revealed a single nucleotide polymorphism (SNP) transition from guanine to adenosine on BraA10g021490.3C in pem1, causing an amino acid shift from glycine to glutamic acid (G to E); in addition, BraA10g021490.3C in pem2 was found to have a single nucleotide substitution from guanine to adenosine, causing an amino acid change from E to lysine (K). BraA10g021490.3C is a homolog of the Arabidopsisdivinyl chlorophyllide a 8-vinyl-reductase (DVR) gene that encodes 3,8-divinyl protochlorophyllide a 8-vinyl reductase, which is a key enzyme in chlorophyll biosynthesis. Enzyme activity assay and chlorophyll composition analysis demonstrated that impaired DVR had partial loss of function. These results provide a basis to understand chlorophyll metabolism and explore the mechanism of a pale green phenotype in Chinese cabbage.

Effects of Maternal Environment on Seed Germination and Seedling Vigor of Petunia × hybrida under Different Abiotic Stresses

Seed germination and seedling vigor can be affected by environmental cues experienced by the mother plant. However, information about how the maternal environment affects seed quality is scarce in ornamental plants. This study aimed to investigate the effects of two different maternal environments on the seed germination and seedling vigor of Petunia × hybrida under a variety of abiotic stresses. Petunia mother plants were grown in either a greenhouse during the summer months or an indoor controlled-temperature-and-light environment. Collected seeds were subjected to external stressors, including polyethylene glycol (PEG), sodium chloride (NaCl), high temperature, and abscisic acid (ABA), to determine seed germination percentage and seedling vigor. Results indicated that seeds harvested from the mother plants grown in a controlled environment germinated better than seeds harvested from the mother plants grown in the greenhouse when suboptimal germination conditions were applied. Additionally, the seedlings from the controlled maternal environment performed better in both ABA and salinity stress tests than the greenhouse seedlings. Interestingly, the greenhouse seedlings displayed less reactive oxygen species (ROS) damage and lower electrolyte leakage than the controlled environment seedlings under dehydration stress. The difference in germination and seedling vigor of seeds from the two different maternal environments might be due to the epigenetic memory inherited from the mother plants. This study highlighted the strong impact of the maternal environment on seed germination and seedling vigor in Petunia and may assist in high-quality seed production in ornamental plants.

Identification of two QTLs associated with high fruit acidity in apple using pooled genome sequencing analysis

Acidity is a critical component determining apple fruit quality. Previous studies reported two major acidity quantitative trait loci (QTLs) on linkage groups (LGs) 16 (Ma) and 8 (Ma3), respectively, and their homozygous genotypes mama and ma3ma3 usually confer low titratable acidity (TA) (<3.0 mg ml-1) to apple fruit. However, apples of genotypes Ma- (MaMa and Mama) or Ma3- (Ma3Ma3 and Ma3ma3) frequently show an acidity range spanning both regular (TA 3.0-10.0 mg ml-1) and high (TA > 10 mg ml-1) acidity levels. To date, the genetic control for high-acidity apples remains essentially unknown. In order to map QTLs associated with high acidity, two genomic DNA pools, one for high acidity and the other for regular acidity, were created in an interspecific F1 population Royal Gala (Malus domestica) × PI 613988 (M. sieversii) of 191 fruit-bearing progenies. By Illumina paired-end sequencing of the high and regular acidity pools, 1,261,640 single-nucleotide variants (SNVs) commonly present in both pools were detected. Using allele frequency directional difference and density (AFDDD) mapping approach, one region on chromosome 4 and another on chromosome 6 were identified to be putatively associated with high acidity, and were named Ma6 and Ma4, respectively. Trait association analysis of DNA markers independently developed from the Ma6 and Ma4 regions confirmed the mapping of Ma6 and Ma4. In the background of MaMa, 20.6% of acidity variation could be explained by Ma6, 28.5% by Ma4, and 50.7% by the combination of both. The effects of Ma6 and Ma4 in the background of Mama were also significant, but lower. These findings provide important genetic insight into high acidity in apple.

Role of BrSDG8 on bolting in Chinese cabbage (Brassica rapa)

Mapping and resequencing of two allelic early bolting mutants ebm5-1 and ebm5-2 revealed that the BrSDG8 gene is related to bolting in Chinese cabbage (Brassica rapa ssp. pekinensis). Bolting influences the leafy head formation and seed yield of Chinese cabbage therefore being an important agronomic trait. Herein, two allelic early bolting mutants, ebm5-1 and ebm5-2, stably inherited in Chinese cabbage were obtained from wild-type 'FT' seeds by ethyl methane sulfonate mutagenesis. Both mutants flowered significantly earlier than 'FT,' and genetic analysis revealed that the early bolting of the two mutants was controlled by one recessive nuclear gene. With BSR-seq, the mutations originating lines ebm5-1 and ebm5-2 were located to the same region in chromosome A07. Using the 1741 F₂ individuals with the ebm5-1 phenotype as the mapping population, this region was narrowed to 56.24 kb between markers InDel18 and InDel45. A single-nucleotide polymorphism (SNP) was aligned to the BraA07g040740.3C (BrSDG8) region by whole-genome resequencing of ebm5-1 mutant and 'FT.' BrSDG8 is a homolog of Arabidopsis thaliana SDG8 encoding a histone methyltransferase affecting H3K4 trimethylation in FLOWERING LOCUS C chromatin. Comparative sequencing established that the SNP occurred on BrSDG8 17th exon in ebm5-1. Genotype analysis showed full co-segregation of the early bolting phenotype with this SNP. Cloning of allelic mutant ebm5-2 indicated that it harbors a deletion mutation on the 12th exon of BrSDG8. Quantitative real-time PCR analysis indicated that BrSDG8 expression level was observably lower in mutant ebm5-1 than in 'FT.' Overall, the present results provide strong evidence that BrSDG8 mutation leads to early bolting in Chinese cabbage, thereby providing a basis to understand the molecular mechanisms underlying this phenotype.

The cloning and characterization of hypersensitive to salt stress mutant, affected in quinolinate synthase, highlights the involvement of NAD in stress-induced accumulation of ABA and proline

Nicotinamide adenine dinucleotide (NAD), a ubiquitous coenzyme, is required for many physiological reactions and processes. However, it remains largely unknown how NAD affects plant response to salt stress. We isolated a salt-sensitive mutant named hypersensitive to salt stress (hss) from an ethyl methanesulfonate-induced mutation population. A point mutation was identified by MutMap in the encoding region of Quinolinate Synthase (QS) gene required for the de novo synthesis of NAD. This point mutation caused a substitution of amino acid in the highly-conserved NadA domain of QS, resulting in an impairment of NAD biosynthesis in the mutant. Molecular and chemical complementation have restored the response of the hss mutant to salt stress, indicating that the decreased NAD contents in the mutant were responsible for its hypersensitivity to salt stress. Furthermore, the endogenous levels of abscisic acid (ABA) and proline were also reduced in stress-treated hss mutant. The application of ABA or proline could alleviate stress-induced oxidative damage of the mutant and partially rescue its hypersensitivity to salt stress, but not affect NAD concentration. Taken together, our results demonstrated that the NadA domain of QS is important for NAD biosynthesis, and NAD participates in plant response to salt stress by affecting stress-induced accumulation of ABA and proline.

Characterization of four polymorphic genes controlling red leaf colour in lettuce that have undergone disruptive selection since domestication

Anthocyanins protect plants from biotic and abiotic stressors and provide great health benefits to consumers. In this study, we cloned four genes (Red Lettuce Leaves 1 to 4: RLL1 to RLL4) that contribute to colour variations in lettuce. The RLL1 gene encodes a bHLH transcription factor, and a 5-bp deletion in some cultivars abolishes its function to activate the anthocyanin biosynthesis pathway. The RLL2 gene encodes an R2R3-MYB transcription factor, which was derived from a duplication followed by mutations in its promoter region. The RLL3 gene encodes an R2-MYB transcription factor, which down-regulates anthocyanin biosynthesis through competing with RLL2 for interaction with RLL1; a mis-sense mutation compromises the capacity of RLL3 to bind RLL1. The RLL4 gene encodes a WD-40 transcription factor, homologous to the RUP genes suppressing the UV-B signal transduction pathway in Arabidopsis; a mis-sense mutation in rll4 attenuates its suppressing function, leading to a high concentration of anthocyanins. Sequence analysis of the RLL1-RLL4 genes from wild and cultivated lettuce showed that their function-changing mutations occurred after domestication. The mutations in rll1 disrupt anthocyanin biosynthesis, while the mutations in RLL2, rll3 and rll4 activate anthocyanin biosynthesis, showing disruptive selection for leaf colour during domestication of lettuce. The characterization of multiple polymorphic genes in this study provides the necessary molecular resources for the rational breeding of lettuce cultivars with distinct levels of red pigments and green cultivars with high levels of health-promoting flavonoids.

Functional variants of DOG1 control seed chilling responses and variation in seasonal life-history strategies in Arabidopsis thaliana

The seasonal timing of seed germination is critical for plant fitness in different climates. To germinate at the right time of year, seeds respond to seasonal environmental cues, such as cold temperatures. We characterized genetic variation in seed dormancy responses to cold across the geographic range of a widespread annual plant. Induction of secondary seed dormancy during winter conditions (which restricts germination to autumn) was positively correlated with flowering time, constructing winter and spring seasonal life-history strategies. Variation in seed chilling responses was strongly associated with functional variants of a known dormancy gene. These variants showed evidence of ancient diversification associated with Pleistocene glacial cycles, and were associated with climate gradients across the species’ geographical range.

The seasonal timing of seed germination determines a plant’s realized environmental niche, and is important for adaptation to climate. The timing of seasonal germination depends on patterns of seed dormancy release or induction by cold and interacts with flowering-time variation to construct different seasonal life histories. To characterize the genetic basis and climatic associations of natural variation in seed chilling responses and associated life-history syndromes, we selected 559 fully sequenced accessions of the model annual species Arabidopsis thaliana from across a wide climate range and scored each for seed germination across a range of 13 cold stratification treatments, as well as the timing of flowering and senescence. Germination strategies varied continuously along 2 major axes: 1) Overall germination fraction and 2) induction vs. release of dormancy by cold. Natural variation in seed responses to chilling was correlated with flowering time and senescence to create a range of seasonal life-history syndromes. Genome-wide association identified several loci associated with natural variation in seed chilling responses, including a known functional polymorphism in the self-binding domain of the candidate gene DOG1. A phylogeny of DOG1 haplotypes revealed ancient divergence of these functional variants associated with periods of Pleistocene climate change, and Gradient Forest analysis showed that allele turnover of candidate SNPs was significantly associated with climate gradients. These results provide evidence that A. thaliana’s germination niche and correlated life-history syndromes are shaped by past climate cycles, as well as local adaptation to contemporary climate.

Whole-genome resequencing from bulked-segregant analysis reveals gene set based association analyses for the Vibrio anguillarum resistance of turbot (Scophthalmus maximus)

Many achievements have been made to develop quantitative trait loci (QTLs) and gene-associated single nucleotide polymorphisms (SNPs) to facilitate practical marker-assisted selection (MAS) in aquatic animals. However, the systematic studies of SNPs associated with extreme threshold traits were poor in populations lacking of parental genomic information. Coupling next generation sequencing with bulked segregant analysis (BSA) should allow identification of numerous associated SNPs with extreme phenotypes. In the present study, using combination of SNP frequency difference and Euclidean distance, we conducted linkage analysis of SNPs located in genes involved in immune responses, and identified markers associated with Vibrio anguillarum resistance in turbot (Scophthalmus maximus). A total of 221 SNPs was found as candidate SNPs between resistant and susceptible individuals. Among these SNPs, 35 loci located in immune related genes were genotyped in verification population and 7 of them showed significant association with V. anguillarum resistance in both alleles and genotypes (P < 0.05). Among these 7 genes, PIK3CA-like, CYLD, VCAM1, RhoB and RhoGEF are involved in PI3K/Akt/mTOR pathway and NF-κB pathway, which influence the efficiency of bacteria entering the host and inflammation. SNP-SNP interaction analysis was performed by generalized multifactor dimensionality reduction (GMDR). The combination of SNP loci in RhoB, PIK3CA-like and ADCY3 showed a significant effect on V. anguillarum resistance with the verification rate in the sequencing population up to 70.8%. Taken all, our findings demonstrated the feasibility of BSA-seq approach in identifying genes responsible for the extreme phenotypes and will aid in performing MAS in turbot.

Mapping Cucumber Vein Yellowing Virus Resistance in Cucumber (Cucumis sativus L.) by Using BSA-seq Analysis

Cucumber vein yellowing virus (CVYV) causes severe yield losses in cucurbit crops across Mediterranean countries. The control of this virus is based on cultural practices to prevent the presence of its vector (Bemisia tabaci) and breeding for natural resistance, which requires the identification of the loci involved and the development of molecular markers for linkage analysis. In this work, we mapped a monogenic locus for resistance to CVYV in cucumber by using a Bulked Segregant Analysis (BSA) strategy coupled with whole-genome resequencing. We phenotyped 135 F₃ families from a segregating population between a susceptible pickling cucumber and a resistant Long Dutch type cucumber for CVYV resistance. Phenotypic analysis determined the monogenic and incomplete dominance inheritance of the resistance. We named the locus CsCvy-1. For mapping this locus, 15 resistant and 15 susceptible homozygous F₂ individuals were selected for whole genome resequencing. By using a customized bioinformatics pipeline, we identified a unique region in chromosome 5 associated to resistance to CVYV, explaining more than 80% of the variability. The resequencing data provided us with additional SNP markers to decrease the interval of CsCvy-1 to 625 kb, containing 24 annotated genes. Markers flanking CsCvy-1 in a 5.3 cM interval were developed for marker-assisted selection (MAS) in breeding programs and will be useful for the identification of the target gene in future studies.

Apple fruit acidity is genetically diversified by natural variations in three hierarchical epistatic genes: MdSAUR37, MdPP2CH and MdALMTII

Many efforts have been made to map quantitative trait loci (QTLs) to facilitate practical marker-assisted selection (MAS) in plants. In the present study, using MapQTL and BSA-seq (bulk segregant analysis using next generation sequencing) with two independent pedigree-based populations, we identified four major genome-wide QTLs responsible for apple fruit acidity. Candidate genes were screened in major QTL regions, and three functional gene markers, including a non-synonymous A/G single-nucleotide polymorphism (SNP) in the coding region of MdPP2CH, a 36-bp insertion in the promoter of MdSAUR37 and a previously reported SNP in MdALMTII, were validated to influence the malate content of apple fruits. In addition, MdPP2CH inactivated three vacuolar H⁺ -ATPases (MdVHA-A3, MdVHA-B2 and MdVHA-D2) and one aluminium-activated malate transporter (MdALMTII) via dephosphorylation and negatively influenced fruit malate accumulation. The dephosphotase activity of MdPP2CH was suppressed by MdSAUR37, which implied a higher hierarchy of genetic interaction. Therefore, the MdSAUR37/MdPP2CH/MdALMTII chain cascaded hierarchical epistatic genetic effects to precisely determine apple fruit malate content. An A/G SNP (-1010) on the MdMYB44 promoter region from a major QTL (qtl08.1) was closely associated with fruit malate content. The predicted phenotype values (PPVs) were estimated using the tentative genotype values of the gene markers, and the PPVs were significantly correlated with the observed phenotype values. Our findings provide an insight into plant genome-based selection in apples and will aid in conducting research to understand the fundamental physiological basis of quantitative genetics.

Exploring DNA variant segregation types in pooled genome sequencing enables effective mapping of weeping trait in Malus

To unlock the power of next generation sequencing-based bulked segregant analysis in allele discovery in out-crossing woody species, and to understand the genetic control of the weeping trait, an F1 population from the cross 'Cheal's Weeping' × 'Evereste' was used to create two genomic DNA pools 'weeping' (17 progeny) and 'standard' (16 progeny). Illumina pair-end (2 × 151 bp) sequencing of the pools to a 27.1× (weeping) and a 30.4× (standard) genome (742.3 Mb) coverage allowed detection of 84562 DNA variants specific to 'weeping', 92148 specific to 'standard', and 173169 common to both pools. A detailed analysis of the DNA variant genotypes in the pools predicted three informative segregation types of variants: (type I) in weeping pool-specific variants, and (type II) and (type III) in variants common to both pools, where the first allele is assumed to be weeping linked and the allele shown in bold is a variant in relation to the reference genome. Conducting variant allele frequency and density-based mappings revealed four genomic regions with a significant association with weeping: a major locus, Weeping (W), on chromosome 13 and others on chromosomes 10 (W2), 16 (W3), and 5 (W4). The results from type I variants were noisier and less certain than those from type II and type III variants, demonstrating that although type I variants are often the first choice, type II and type III variants represent an important source of DNA variants that can be exploited for genetic mapping in out-crossing woody species. Confirmation of the mapping of W and W2, investigation into their genetic interactions, and identification of expressed genes in the W and W2 regions provided insight into the genetic control of weeping and its expressivity in Malus.

Systems biology and genome-wide approaches to unveil the molecular players involved in the pre-germinative metabolism: implications on seed technology traits

The pre-germinative metabolism is among the most fascinating aspects of seed biology. The early seed germination phase, or pre-germination, is characterized by rapid water uptake (imbibition), which directs a series of dynamic biochemical events. Among those are enzyme activation, DNA damage and repair, and use of reserve storage compounds, such as lipids, carbohydrates and proteins. Industrial seedling production and intensive agricultural production systems require seed stocks with high rate of synchronized germination and low dormancy. Consequently, seed dormancy, a quantitative trait related to the activation of the pre-germinative metabolism, is probably the most studied seed trait in model species and crops. Single omics, systems biology, QTLs and GWAS mapping approaches have unveiled a list of molecules and regulatory mechanisms acting at transcriptional, post-transcriptional and post-translational levels. Most of the identified candidate genes encode for regulatory proteins targeting ROS, phytohormone and primary metabolisms, corroborating the data obtained from simple molecular biology approaches. Emerging evidences show that epigenetic regulation plays a crucial role in the regulation of these mentioned processes, constituting a still unexploited strategy to modulate seed traits. The present review will provide an up-date of the current knowledge on seed pre-germinative metabolism, gathering the most relevant results from physiological, genetics, and omics studies conducted in model and crop plants. The effects exerted by the biotic and abiotic stresses and priming are also addressed. The possible implications derived from the modulation of pre-germinative metabolism will be discussed from the point of view of seed quality and technology.

Efficient Identification of Causal Mutations through Sequencing of Bulked F 2 from Two Allelic Bloomless Mutants of Sorghum bicolor

Sorghum (Sorghum bicolor Moench, L.) plant accumulates copious layers of epi-cuticular wax (EW) on its aerial surfaces, to a greater extent than most other crops. EW provides a vapor barrier that reduces water loss, and is therefore considered to be a major determinant of sorghum's drought tolerance. However, little is known about the genes responsible for wax accumulation in sorghum. We isolated two allelic mutants, bloomless40-1 (bm40-1) and bm40-2, from a mutant library constructed from ethyl methane sulfonate (EMS) treated seeds of an inbred, BTx623. Both mutants were nearly devoid of the EW layer. Each bm mutant was crossed to the un-mutated BTx623 to generated F₂ populations that segregated for the bm phenotype. Genomic DNA from 20 bm F₂ plants from each population was bulked for whole genome sequencing. A single gene, Sobic.001G228100, encoding a GDSL-like lipase/acylhydrolase, had unique homozygous mutations in each bulked F₂ population. Mutant bm40-1 harbored a missense mutation in the gene, whereas bm40-2 had a splice donor site mutation. Our findings thus provide strong evidence that mutation in this GDSL-like lipase gene causes the bm phenotype, and further demonstrate that this approach of sequencing two independent allelic mutant populations is an efficient method for identifying causal mutations. Combined with allelic mutants, MutMap provides powerful method to identify all causal genes for the large collection of bm mutants in sorghum, which will provide insight into how sorghum plants accumulate such abundant EW on their aerial surface. This knowledge may facilitate the development of tools for engineering drought-tolerant crops with reduced water loss.

Next-Generation Sequencing from Bulked-Segregant Analysis Accelerates the Simultaneous Identification of Two Qualitative Genes in Soybean

Next-generation sequencing (NGS)-based bulked-segregant analysis (BSA) approaches have been proven successful for rapidly mapping genes in plant species. However, most such methods are based on mutants and usually only one gene controlling the mutant phenotype is identified. In this study, NGS-based BSA was employed to map simultaneously two qualitative genes controlling cotyledon color of seed in soybean. Yellow-cotyledon (YC) and green-cotyledon (GC) bulks from progenies of a biparental population (Zhonghuang 30 × Jiyu 102) were sequenced. The SNP-index of each SNP locus in YC and GC bulks was calculated and two genomic regions on chromosomes 1 and 11 harboring, respectively, loci qCC1 and qCC2 were identified by Δ(SNP-index) analysis. These two BSA-seq-derived loci were further validated with SSR markers and fine-mapped. qCC1 was mapped to a 30.7-kb region containing four annotated genes and qCC2 was mapped to a 67.7-kb region with nine genes. These two regions contained, respectively, genes D1 and D2, which had previously been identified by homology-based cloning as being associated with cotyledon color. Sequence analysis of the NGS data also identified a frameshift deletion in the coding region of D1. These results suggested that BSA-seq could accelerate the mapping of loci controlling qualitative traits, even if a trait is controlled by more than one locus.