Dissecting Key Adaptation Traits in the Polyploid Perennial Medicago sativa Using GBS-SNP Mapping

Morphological diversity variation of seed traits among 587 germplasm resources of Medicago Genus and 32 germplasm resources of Trigonella Genus

Germplasm resources within the Medicago genus are highly regarded for their value as forage crops and their critical roles in nitrogen cycling, ecosystem restoration, and soil structure improvement. Therefore, understanding the diversity of seed morphology in this genus is essential for advancing its development and utilization. This study analyzed seed samples from 587 germplasm accessions representing 77 species within Medicago genus, as well as 32 accessions from 21 species within the closely related genus Trigonella. A statistical analysis was conducted on twelve traits, including seven quantitative traits-straight length (SL), straight width (SW), width-to-length ratio (WL), perimeter (PE), radicle length (RL), hilum length (HL), and 100-seed weight (SY)-and five qualitative traits, including seed coat condition, radicle characteristics, seed size, shape, and color. The results revealed that: (1) there was significant diversity (P < 0.05) in SL, SW, WL, PE, RL, HL, and SY across Medicago species; (2) principal component analysis of the 587 Medicago accessions identified SL, SW, PE, HL, RL, and SY as the primary contributors to morphological diversity; and (3) high-resolution images of seeds from various accessions were captured for future research. This study provides a solid foundation for the establishment of seed banks and the enhancement of germplasm resources through the systematic analysis of these morphological traits.

Genomic characterization and gene bank curation of Aegilops: the wild relatives of wheat

Genetic diversity found in crop wild relatives is critical to preserve and utilize for crop improvement to achieve sustainable food production amid climate change and increased demand. We genetically characterized a large collection of 1,041 Aegilops accessions distributed among 23 different species using more than 45K single nucleotide polymorphisms identified by genotyping-by-sequencing. The Wheat Genetics Resource Center (WGRC) Aegilops germplasm collection was curated through the identification of misclassified and redundant accessions. There were 49 misclassified and 28 sets of redundant accessions within the four diploid species. The curated germplasm sets now have improved utility for genetic studies and wheat improvement. We constructed a phylogenetic tree and principal component analysis cluster for all Aegilops species together, giving one of the most comprehensive views of Aegilops. The Sitopsis section and the U genome Aegilops clade were further scrutinized with in-depth population analysis. The genetic relatedness among the pair of Aegilops species provided strong evidence for the species evolution, speciation, and diversification. We inferred genome symbols for two species Ae. neglecta and Ae. columnaris based on the sequence read mapping and the presence of segregating loci on the pertinent genomes as well as genetic clustering. The high genetic diversity observed among Aegilops species indicated that the genus could play an even greater role in providing the critical need for untapped genetic diversity for future wheat breeding and improvement. To fully characterize these Aegilops species, there is an urgent need to generate reference assemblies for these wild wheats, especially for the polyploid Aegilops.

Multispectral Phenotyping and Genetic Analyses of Spring Appearance in Greening Plant, Phedimus spp

The change in appearance during the seasonal transitions in ornamental greening plants is an important characteristic. In particular, the early onset of green leaf color is a desirable trait for a cultivar. In this study, we established a method for phenotyping leaf color change by multispectral imaging and performed genetic analysis based on the phenotypes to clarify the potential of the approach in breeding greening plants. We performed multispectral phenotyping and quantitative trait locus (QTL) analysis of an F₁ population derived from 2 parental lines of Phedimus takesimensis, known to be a drought and heat-tolerant rooftop plant species. The imaging was conducted in April of 2019 and 2020 when dormancy breakage occurs and growth extension begins. Principal component analysis of 9 different wavelength values showed a high contribution from the first principal component (PC1), which captured variation in the visible light range. The high interannual correlation in PC1 and in the intensity of visible light indicated that the multispectral phenotyping captured genetic variation in the color of leaves. We also performed restriction site-associated DNA sequencing and obtained the first genetic linkage map of Phedimus spp. QTL analysis revealed 2 QTLs related to early dormancy breakage. Based on the genotypes of the markers underlying these 2 QTLs, the F₁ phenotypes with early (late) dormancy break, green (red or brown) leaves, and a high (low) degree of vegetative growth were classified. The results suggest the potential of multispectral phenotyping in the genetic dissection of seasonal leaf color changes in greening plants.

On the usefulness of mock genomes to define heterotic pools, testers, and hybrid predictions in orphan crops

The advances in genomics in recent years have increased the accuracy and efficiency of breeding programs for many crops. Nevertheless, the adoption of genomic enhancement for several other crops essential in developing countries is still limited, especially for those that do not have a reference genome. These crops are more often called orphans. This is the first report to show how the results provided by different platforms, including the use of a simulated genome, called the mock genome, can generate in population structure and genetic diversity studies, especially when the intention is to use this information to support the formation of heterotic groups, choice of testers, and genomic prediction of single crosses. For that, we used a method to assemble a reference genome to perform the single-nucleotide polymorphism (SNP) calling without needing an external genome. Thus, we compared the analysis results using the mock genome with the standard approaches (array and genotyping-by-sequencing (GBS)). The results showed that the GBS-Mock presented similar results to the standard methods of genetic diversity studies, division of heterotic groups, the definition of testers, and genomic prediction. These results showed that a mock genome constructed from the population's intrinsic polymorphisms to perform the SNP calling is an effective alternative for conducting genomic studies of this nature in orphan crops, especially those that do not have a reference genome.

Transition of bioeconomy as a key concept for the agriculture and agribusiness development: An extensive review on ASEAN countries

Transition of bioeconomy plays an important role in the development of any economy. Therefore, the purpose of the current paper is to review the key concept of transition of the bioeconomy in the ASEAN economies that is considered to be an innovative collective method for observing the operation of national economies' numerous sectors that produce and process biological wealth. Based on current literature, it has been reviewed that the transition of bioeconomy is being presented with respect to the theoretical concept, and their essence and bioeconomy size, and the challenges and risks which are related to the bioeconomy transition. The various ways and the area of action along with present strategies which are supporting the development of bio economies and their relationship with the circular economies model had been presented. It was seen in the literature that ASEAN agricultures sectors had been contributing their major part in the development of value-added products and employment, except for Malaysia, Thailand, and Laos. Moreover, the analysis results had also shown that the transition of bioeconomy is considered to be an important concept in the growth of agriculture, forestry, agribusiness, and various sectors that produce and use bio-based raw materials. It was further found that for real bio-economy opportunities, it is essential that there should be a proper national and regional strategy for the bioeconomy development in place and for the development of a proper design and proper management coordination at every level of sectors, enterprises, and provincial systems.

OMICS in Fodder Crops: Applications, Challenges, and Prospects

Biomass yield and quality are the primary targets in forage crop improvement programs worldwide. Low-quality fodder reduces the quality of dairy products and affects cattle's health. In multipurpose crops, such as maize, sorghum, cowpea, alfalfa, and oat, a plethora of morphological and biochemical/nutritional quality studies have been conducted. However, the overall growth in fodder quality improvement is not on par with cereals or major food crops. The use of advanced technologies, such as multi-omics, has increased crop improvement programs manyfold. Traits such as stay-green, the number of tillers per plant, total biomass, and tolerance to biotic and/or abiotic stresses can be targeted in fodder crop improvement programs. Omic technologies, namely genomics, transcriptomics, proteomics, metabolomics, and phenomics, provide an efficient way to develop better cultivars. There is an abundance of scope for fodder quality improvement by improving the forage nutrition quality, edible quality, and digestibility. The present review includes a brief description of the established omics technologies for five major fodder crops, i.e., sorghum, cowpea, maize, oats, and alfalfa. Additionally, current improvements and future perspectives have been highlighted.

A high-throughput skim-sequencing approach for genotyping, dosage estimation and identifying translocations

The development of next-generation sequencing (NGS) enabled a shift from array-based genotyping to directly sequencing genomic libraries for high-throughput genotyping. Even though whole-genome sequencing was initially too costly for routine analysis in large populations such as breeding or genetic studies, continued advancements in genome sequencing and bioinformatics have provided the opportunity to capitalize on whole-genome information. As new sequencing platforms can routinely provide high-quality sequencing data for sufficient genome coverage to genotype various breeding populations, a limitation comes in the time and cost of library construction when multiplexing a large number of samples. Here we describe a high-throughput whole-genome skim-sequencing (skim-seq) approach that can be utilized for a broad range of genotyping and genomic characterization. Using optimized low-volume Illumina Nextera chemistry, we developed a skim-seq method and combined up to 960 samples in one multiplex library using dual index barcoding. With the dual-index barcoding, the number of samples for multiplexing can be adjusted depending on the amount of data required, and could be extended to 3,072 samples or more. Panels of doubled haploid wheat lines (Triticum aestivum, CDC Stanley x CDC Landmark), wheat-barley (T. aestivum x Hordeum vulgare) and wheat-wheatgrass (Triticum durum x Thinopyrum intermedium) introgression lines as well as known monosomic wheat stocks were genotyped using the skim-seq approach. Bioinformatics pipelines were developed for various applications where sequencing coverage ranged from 1 × down to 0.01 × per sample. Using reference genomes, we detected chromosome dosage, identified aneuploidy, and karyotyped introgression lines from the skim-seq data. Leveraging the recent advancements in genome sequencing, skim-seq provides an effective and low-cost tool for routine genotyping and genetic analysis, which can track and identify introgressions and genomic regions of interest in genetics research and applied breeding programs.

A QTL approach in faba bean highlights the conservation of genetic control of frost tolerance among legume species

Frost is a major abiotic stress of winter type faba beans (Vica faba L.) and has adverse effects on crop yield. Climate change, far from reducing the incidence of frost events, is making these phenomena more and more common, severe, and prolonged. Despite the important interaction that the environment has in the tolerance of faba bean to frost, this trait seems to have good levels of heritability. Several QTLs for frost tolerance have already been reported, however, a more robust identification is needed to more precisely identify the genomic regions involved in faba bean tolerance to sub-zero temperatures. Several pea (Pisum sativum L.) and barrel medic (Medicago truncatula L.) frost tolerance QTLs appear to be conserved between these two species, furthering the hypothesis that the genetic control of frost tolerance in legume species might be more generally conserved. In this work, the QTL mapping in two faba bean recombinant inbred line (RIL) populations connected by a common winter-type parent has led to the identification of five genomic regions involved in the control of frost tolerance on linkage groups I, III, IV, and V. Among them, a major and robust QTL of great interest for marker-assisted selection was identified on the lower part of the long-arm of LGI. The synteny between the faba bean frost tolerance QTLs and those previously identified in other legume species such as barrel medic, pea or soybean highlighted at least partial conservation of the genetic control of frost tolerance among different faba bean genetic pools and legume species. Four novel RILs showing high and stable levels of tolerance and the ability to recover from freezing temperatures by accumulating frost tolerance QTLs are now available for breeding programs.

Combining QTL mapping and RNA-Seq Unravels candidate genes for Alfalfa (Medicago sativa L.) leaf development

BACKGROUND

Leaf size affects crop canopy morphology and photosynthetic efficiency, which can influence forage yield and quality. It is of great significance to mine the key genes controlling leaf development for breeding new alfalfa varieties. In this study, we mapped leaf length (LL), leaf width (LW), and leaf area (LA) in an F1 mapping population derived from a cultivar named ZhongmuNo.1 with larger leaf area and a landrace named Cangzhou with smaller leaf area.

RESULTS

This study showed that the larger LW was more conducive to increasing LA. A total of 24 significant quantitative trait loci (QTL) associated with leaf size were identified on both the paternal and maternal linkage maps. Among them, nine QTL explained about 11.50-22.45% phenotypic variation. RNA-seq analysis identified 2,443 leaf-specific genes and 3,770 differentially expressed genes. Combining QTL mapping, RNA-seq alalysis, and qRT-PCR, we identified seven candidate genes associated with leaf development in five major QTL regions.

CONCLUSION

Our study will provide a theoretical basis for marker-assisted breeding and lay a foundation for further revealing molecular mechanism of leaf development in alfalfa.

Application of machine learning to explore the genomic prediction accuracy of fall dormancy in autotetraploid alfalfa

Fall dormancy (FD) is an essential trait to overcome winter damage and for alfalfa (Medicago sativa) cultivar selection. The plant regrowth height after autumn clipping is an indirect way to evaluate FD. Transcriptomics, proteomics, and quantitative trait locus mapping have revealed crucial genes correlated with FD; however, these genes cannot predict alfalfa FD very well. Here, we conducted genomic prediction of FD using whole-genome SNP markers based on machine learning-related methods, including support vector machine (SVM) regression, and regularization-related methods, such as Lasso and ridge regression. The results showed that using SVM regression with linear kernel and the top 3000 genome-wide association study (GWAS)-associated markers achieved the highest prediction accuracy for FD of 64.1%. For plant regrowth height, the prediction accuracy was 59.0% using the 3000 GWAS-associated markers and the SVM linear model. This was better than the results using whole-genome markers (25.0%). Therefore, the method we explored for alfalfa FD prediction outperformed the other models, such as Lasso and ElasticNet. The study suggests the feasibility of using machine learning to predict FD with GWAS-associated markers, and the GWAS-associated markers combined with machine learning would benefit FD-related traits as well. Application of the methodology may provide potential targets for FD selection, which would accelerate genetic research and molecular breeding of alfalfa with optimized FD.

Genomic prediction for canopy height and dry matter yield in alfalfa using family bulks

Genomic selection (GS) has proven to be an effective method to increase genetic gain rates and accelerate breeding cycles in many crop species. However, its implementation requires large investments to phenotype of the training population and for routine genotyping. Alfalfa (Medicago sativa L.) is one of the major cultivated forage legumes, showing high-quality nutritional value. Alfalfa breeding is usually carried out by phenotypic recurrent selection and is commonly done at the family level. The application of GS in alfalfa could be simplified and less costly by genotyping and phenotyping families in bulks. For this study, an alfalfa reference population composed of 142 full-sib and 35 half-sib families was bulk-genotyped using target enrichment sequencing and phenotyped for dry matter yield (DMY) and canopy height (CH) in Florida, USA. Genotyping of the family bulks with 17,707 targeted probes resulted in 114,945 single-nucleotide polymorphisms. The markers revealed a population structure that matched the mating design, and the linkage disequilibrium slowly decayed in this breeding population. After exploring multiple prediction scenarios, a strategy was proposed including data from multiple harvests and accounting for the G×E in the training population, which led to a higher predictive ability of up to 38 and 24% for DMY and CH, respectively. Although this study focused on the implementation of GS in alfalfa families, the bulk methodology and the prediction schemes used herein could guide future studies in alfalfa and other crops bred in bulks.

Genome assembly of alfalfa cultivar zhongmu-4 and identification of SNPs associated with agronomic traits

Alfalfa (Medicago sativa L.) is the most important legume forage crop worldwide with high nutritional value and yield. For a long time, the breeding of alfalfa was hampered by lacking reliable information on the autotetraploid genome and molecular markers linked to important agronomic traits. We herein reported the de novo assembly of the allele-aware chromosome-level genome of Zhongmu-4, a cultivar widely cultivated in China, and a comprehensive database of genomic variations based on resequencing of 220 germplasms. Approximate 2.74 Gb contigs (N50 of 2.06 Mb), accounting for 88.39% of the estimated genome, were assembled, and 2.56 Gb contigs were anchored to 32 pseudo-chromosomes. A total of 34,922 allelic genes were identified from the allele-aware genome. We observed the expansion of gene families, especially those related to the nitrogen metabolism, and the increase of repetitive elements including transposable elements, which probably resulted in the increase of Zhongmu-4 genome compared with Medicago truncatula. Population structure analysis revealed that the accessions from Asia and South America had relatively lower genetic diversity than those from Europe, suggesting that geography may influence alfalfa genetic divergence during local adaption. Genome-wide association studies identified 101 single nucleotide polymorphisms (SNPs) associated with 27 agronomic traits. Two candidate genes were predicted to be correlated with fall dormancy and salt response. We believe that the allele-aware chromosome-level genome sequence of Zhongmu-4 combined with the resequencing data of the diverse alfalfa germplasms will facilitate genetic research and genomics-assisted breeding in variety improvement of alfalfa.

Mapping freezing tolerance QTL in alfalfa: based on indoor phenotyping

BACKGROUND

Winter freezing temperature impacts alfalfa (Medicago sativa L.) persistence and seasonal yield and can lead to the death of the plant. Understanding the genetic mechanisms of alfalfa freezing tolerance (FT) using high-throughput phenotyping and genotyping is crucial to select suitable germplasm and develop winter-hardy cultivars. Several clones of an alfalfa F1 mapping population (3010 x CW 1010) were tested for FT using a cold chamber. The population was genotyped with SNP markers identified using genotyping-by-sequencing (GBS) and the quantitative trait loci (QTL) associated with FT were mapped on the parent-specific linkage maps. The ultimate goal is to develop non-dormant and winter-hardy alfalfa cultivars that can produce extended growth in the areas where winters are often mild.

RESULTS

Alfalfa FT screening method optimized in this experiment comprises three major steps: clone preparation, acclimation, and freezing test. Twenty clones of each genotype were tested, where 10 samples were treated with freezing temperature, and 10 were used as controls. A moderate positive correlation (r ~ 0.36, P < 0.01) was observed between indoor FT and field-based winter hardiness (WH), suggesting that the indoor FT test is a useful indirect selection method for winter hardiness of alfalfa germplasm. We detected a total of 20 QTL associated with four traits; nine for visual rating-based FT, five for percentage survival (PS), four for treated to control regrowth ratio (RR), and two for treated to control biomass ratio (BR). Some QTL positions overlapped with WH QTL reported previously, suggesting a genetic relationship between FT and WH. Some favorable QTL from the winter-hardy parent (3010) were from the potential genic region for a cold tolerance gene CBF. The BLAST alignment of a CBF sequence of M. truncatula, a close relative of alfalfa, against the alfalfa reference showed that the gene's ortholog resides around 75 Mb on chromosome 6.

CONCLUSIONS

The indoor freezing tolerance selection method reported is useful for alfalfa breeders to accelerate breeding cycles through indirect selection. The QTL and associated markers add to the genomic resources for the research community and can be used in marker-assisted selection (MAS) for alfalfa cold tolerance improvement.

Effect of a fall cut on dry matter yield, nutritive value, and stand persistence of alfalfa

Information is lacking about the effect of date of a fall cut of alfalfa (Medicago sativa L.) on dry matter yield (DMY), forage nutritive value, and stands persistence. The objective of this study was to determine the effect of timing of a fall cut on DMY, forage nutritive value and stand persistence of three alfalfa varieties: low-lignin Hi–Gest 360, Roundup Ready Tonica, and conventional Gunner in Northeastern Kansas in the United States. The field study was carried out by splitting plot in randomized complete block design with four replications. The harvesting data of different maturity stages were collected in each year from 2015 to 2018. Three cuts were harvested based on the stage of maturity, and the last (fourth) cut was done on September 15, September 30, October 15, and October 30 of each year. The persistence of the alfalfa stands was determined each fall after the last cut, and each spring after the first cut, by counting the number of live plants in a randomly placed quadrat in each plot. Alfalfa cut on September 15 and September 30 had a higher stand persistence compared to alfalfa cut on October 15 and October 30. The DMY of the first cut in 2016 was significantly higher in roundup ready than the low-lignin alfalfa variety. In the second cut, DMY was significantly higher in conventional alfalfa than the roundup ready. There were no significant differences in DMY between alfalfa varieties in the rest of seasonal cuttings in 2016 and 2017 and annual total yield in both years. In general, low lignin alfalfa variety had higher crude protein and relative feed value and lower acid detergent fiber and neutral detergent fiber contents than those in roundup ready and conventional alfalfa varieties. On average, nutritive value of alfalfa was generally affected by last cutting dates in 2017. Based on 3-year data the last cutting of alfalfa in the fall could be done by September 30–October 15 without harmful effect on DMY.

Identification of Genetic Loci Associated With Crude Protein Content and Fiber Composition in Alfalfa (Medicago sativa L.) Using QTL Mapping

Forage quality determined mainly by protein content and fiber composition has a crucial influence on digestibility and nutrition intake for animal feeding. To explore the genetic basis of quality traits, we conducted QTL mapping based on the phenotypic data of crude protein (CP), neutral detergent fiber (NDF), acid detergent fiber (ADF), and lignin of an F1 alfalfa population generated by crossing of two alfalfa parents with significant difference in quality. In total, 83 QTLs were identified with contribution to the phenotypic variation (PVE) ranging from 1.45 to 14.35%. Among them, 47 QTLs interacted significantly with environment and 12 QTLs were associated with more than one trait. Epistatic effect was also detected for 73 pairs of QTLs with PVE of 1.08-14.06%. The results suggested that the inheritance of quality-related traits was jointly affected by additive, epistasis and environment. In addition, 83.33% of the co-localized QTLs were shared by ADF and NDF with the same genetic direction, while the additive effect of crude protein-associated QTLs was opposite to that fiber composition on the same locus, suggesting that the loci may antagonistically contribute to protein content and fiber composition. Further analysis of a QTL related to all the three traits of fiber composition (qNDF1C, qADF1C-2, and qlignin1C-2) showed that five candidate genes were homologs of cellulose synthase-like protein A1 in Medicago truncatula, indicating the potential role in fiber synthesis. For the protein-associated loci we identified, qCP4C-1 was located in the shortest region (chr 4.3 39.3-39.4 Mb), and two of the seven corresponding genes in this region were predicted to be E3 ubiquitin-protein ligase in protein metabolism. Therefore, our results provide some reliable regions significantly associated with alfalfa quality, and identification of the key genes would facilitate marker-assisted selection for favorable alleles in breeding program of alfalfa quality improvement.

Construction of high-density genetic linkage map and mapping quantitative trait loci (QTL) for flowering time in autotetraploid alfalfa (Medicago sativa L.) using genotyping by sequencing

Flowering time is an important agronomic trait of alfalfa (Medicago sativa L.). Managing flowering time can produce economic benefits for farmers. To understand the genetic basis of this trait, quantitative trait loci (QTL) mapping was conducted in a full-sib population that consisted of 392 individuals segregating based on flowering time. High density linkage maps were constructed using single nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing (GBS). The linkage maps contained 3,818 SNP markers on 64 linkage groups in two parents. The average marker density was 4.33 cM for Parent 1 (P1) and 1.47 cM for Parent 2 (P2). Phenotypic data for flowering time was collected for three years at one location. Twenty-eight QTLs were identified associated with flowering time. Eleven QTLs explained more than 10% of the phenotypic variation. Among them, five main effect QTLs located on linkage group (LG) 7D of P1 and five main effect QTLs located on LG 6D of P2 were identified. Three QTLs were co-located with other QTLs. The identified linked markers to QTLs could be used for marker-assisted selection in breeding programs to develop new alfalfa varieties to modulate flowering time.

Expanded functionality, increased accuracy, and enhanced speed in the de novo genotyping-by-sequencing pipeline GBS-SNP-CROP

Summary

GBS-SNP-CROP is a bioinformatics pipeline originally developed to support the cost-effective genome-wide characterization of plant genetic resources through paired-end genotyping-by-sequencing (GBS), particularly in the absence of a reference genome. Since its 2016 release, the pipeline’s functionality has greatly expanded, its computational efficiency has improved, and its applicability to a broad set of genomic studies for both plants and animals has been demonstrated. This note details the suite of improvements to date, as realized in GBS-SNP-CROP v.4.0, with specific attention paid to a new integrated metric that facilitates reliable variant identification despite the complications of homologs. Using the new de novo GBS read simulator GBS-Pacecar, also introduced in this note, results show an improvement in overall pipeline accuracy from 66% (v.1.0) to 84% (v.4.0), with a time saving of ∼70%. Both GBS-SNP-CROP versions significantly outperform TASSEL-UNEAK; and v.4.0 resolves the issue of non-overlapping variant calls observed between UNEAK and v.1.0.

Availability and implementation

GBS-SNP-CROP source code and user manual are available at https://github.com/halelab/GBS-SNP-CROP. The GBS read simulator GBS-Pacecar is available at https://github.com/halelab/GBS-Pacecar.

Supplementary information

Supplementary data are available at Bioinformatics online.

Identification of Exogenous Nitric Oxide-Responsive miRNAs from Alfalfa (Medicago sativa L.) under Drought Stress by High-Throughput Sequencing

Alfalfa (Medicago sativa L.) is a high quality leguminous forage. Drought stress is one of the main factors that restrict the development of the alfalfa industry. High-throughput sequencing was used to analyze the microRNA (miRNA) profiles of alfalfa plants treated with CK (normal water), PEG (polyethylene glycol-6000; drought stress), and PEG + SNP (sodium nitroprusside; nitric oxide (NO) sprayed externally under drought stress). We identified 90 known miRNAs belonging to 46 families and predicted 177 new miRNAs. Real-time quantitative fluorescent PCR (qRT-PCR) was used to validate high-throughput expression analysis data. A total of 32 (14 known miRNAs and 18 new miRNAs) and 55 (24 known miRNAs and 31 new miRNAs) differentially expressed miRNAs were identified in PEG and PEG + SNP samples. This suggested that exogenous NO can induce more new miRNAs. The differentially expressed miRNA maturation sequences in the two treatment groups were targeted by 86 and 157 potential target genes, separately. The function of target genes was annotated by gene ontology (GO) enrichment and kyoto encyclopedia of genes and genomes (KEGG) analysis. The expression profiles of nine selected miRNAs and their target genes verified that their expression patterns were opposite. This study has documented that analysis of miRNA under PEG and PEG + SNP conditions provides important insights into the improvement of drought resistance of alfalfa by exogenous NO at the molecular level. This has important scientific value and practical significance for the improvement of plant drought resistance by exogenous NO.

QTL mapping of flowering time and biomass yield in tetraploid alfalfa (Medicago sativa L.)

BACKGROUND

The genetic and genomic basis of flowering time and biomass yield in alfalfa (Medicago sativa L.) remains poorly understood mainly due to the autopolyploid nature of the species and the lack of adequate genomic resources. We constructed linkage maps using genotyping-by-sequencing (GBS) based single dose allele (SDA) SNP and mapped alfalfa timing of flowering (TOF), spring yield (SY), and cumulative summer biomass (CSB) in a pseudo-testcross F1 population derived from a fall dormant (3010) and a non-dormant (CW 1010) cultivars. We analyzed the quantitative trait loci (QTL) to identify conserved genomic regions and detected molecular markers and potential candidate genes associated with the traits to improve alfalfa and provide genomic resources for the future studies.

RESULTS

This study showed that both fall dormant and non-dormant alfalfa cultivars harbored QTL for early and late flowering, suggesting that flowering time in alfalfa is not an indicator of its fall dormancy (FD) levels. A weak phenotypic correlation between the flowering time and fall dormancy (FD) in F1 and checks also corroborated that alfalfa FD and TOF are not the predictors of one another. The relationship between flowering time and alfalfa biomass yield was not strong, but the non-dormant had relatively more SY than dormant. Therefore, selecting superior alfalfa cultivars that are non-dormant, winter-hardy, and early flowering would allow for an early spring harvest with enhanced biomass. In this study, we found 25 QTL for TOF, 17 for SY and six QTL for CSB. Three TOF related QTL were stable and four TOF QTL were detected in the corresponding genomic locations of the flowering QTL of M. truncatula, an indication of possible evolutionarily conserved regions. The potential candidate genes for the SNP sequences of QTL regions were identified for all three traits and these genes would be potential targets for further molecular studies.

CONCLUSIONS

This research showed that variation in alfalfa flowering time after spring green up has no association with dormancy levels. Here we reported QTL, markers, and potential candidate genes associated with spring flowering time and biomass yield of alfalfa, which constitute valuable genomic resources for improving these traits via marker-assisted selection (MAS).

High-density linkage map construction and mapping QTL for yield and yield components in autotetraploid alfalfa using RAD-seq

BACKGROUND

Alfalfa (Medicago sativa L.) is an important forage crop grown worldwide. Alfalfa is called the "queen of forage crops" due to its high forage yield and nutritional characteristics. The aim of this study was to undertake quantitative trait loci (QTL) mapping of yield and yield-related traits in an F₁ population of two alfalfa varieties that differ in their yield and yield-related traits.

RESULTS

We constructed a high-density linkage map using single nucleotide polymorphism (SNP) markers generated by restriction-site associated DNA sequencing (RAD-seq). The linkage map contains 4346 SNP and 119 simple sequence repeat (SSR) markers, with 32 linkage groups for each parent. The average marker distances were 3.00 and 1.32 cM, with coverages of 3455 cM and 4381 cM for paternal and maternal linkage maps, respectively. Using these maps and phenotypic data, we identified a total of 21 QTL for yield and yield components, including five for yield, five for plant height, five for branch number, and six for shoot diameter. Among them, six QTL were co-located for more than one trait. Five QTL explained more than 10% of the phenotypic variation.

CONCLUSIONS

We used RAD-seq to construct a linkage map for alfalfa that greatly enhanced marker density compared to previous studies. This high-density linkage map of alfalfa is a useful reference for mapping yield-related traits. Identified yield-related loci could be used to validate their usefulness in developing markers for maker-assisted selection in breeding populations to improve yield potential in alfalfa.

Chromosomal Locations and Interactions of Four Loci Associated With Seed Coat Color in Watermelon

Different species of edible seed watermelons (Citrullus spp.) are cultivated in Asia and Africa for their colorful nutritious seeds. Consumer preference varies for watermelon seed coat color. Therefore, it is an important consideration for watermelon breeders. In 1940s, a genetic model of four genes, R, T, W and D, was proposed to elucidate the inheritance of seed coat color in watermelon. In this study, we developed three segregating F₂ populations: Sugar Baby (dotted black seed, RRTTWW) × plant introduction (PI) 482379 (green seed, rrTTWW), Charleston Gray (dotted black seed, RRTTWW) × PI 189225 (red seed, rrttWW), and Charleston Gray (dotted black seed, RRTTWWdd) × UGA147 (clump seed, RRTTwwDD) to re-examine the four-gene model and to map the four genes. In the dotted black × green population, the dotted black seed coat color (R_) is dominant to green seed coat color (rr). In the dotted black × red population, the dominant dotted black seed coat color and the recessive red seed coat color segregate for the R and T genes, where the R gene is dominantly epistatic to the T gene. However, the inheritance of the T locus did not fit the four-gene model, thus we named it T¹ . In the dotted black × clump population, the clump seed coat color and the dotted black seed coat color segregate for W and D, where D is recessively epistatic to W. The R, T¹ , W, and D loci were mapped on chromosomes 3, 5, 6, and 8, respectively, using QTL-seq and genotyping-by-sequencing (GBS). Kompetitive Allele Specific PCR (KASP™) assays and SNP markers linked to the four loci were developed to facilitate maker-assisted selection (MAS) for watermelon seed coat color.