Elymus nutans genes for seed shattering and candidate gene-derived EST-SSR markers for germplasm evaluation

RNA-Seq and Genome-Wide Association Studies Reveal Potential Genes for Rice Seed Shattering

The loss of the shattering ability is one of the key events in rice domestication. The strength of the seed shattering ability is closely related to the harvest yield and the adaptability of modern mechanical harvesting methods. In this study, using a population of 587 natural rice cultivars, quantitative trait loci associated with seed shattering were detected by genome-wide association studies (GWASs). We consider the quantitative trait loci (QTLs) qBTS1 and qBTS3 to be the key loci for seed shattering in rice. Additionally, the abscission zone (AZ) and nonabscission zone (NAZ) of materials with a loss of shattering (DZ129) and easy shattering (W517) were subjected to RNA-Seq, and high-quality differential expression profiles were obtained. The AZ-specific differentially expressed genes (DEGs) of W517 were significantly enriched in plant hormone signal transduction, while the AZ-specific DEGs of DZ129 were enriched in phenylpropanoid biosynthesis. We identified candidate genes for the lignin-associated laccase precursor protein (LOC_Os01g63180) and the glycoside hydrolase family (LOC_Os03g14210) in the QTLs qBTS1 (chromosome 1) and qBTS3 (chromosome 3), respectively. In summary, our findings lay the foundation for the further cloning of qBTS1 and qBTS3, which would provide new insights into seed shattering in rice.

Histological characteristics, cell wall hydrolytic enzyme activity, and transcriptome analysis with seed shattering of Stylosanthes accessions

Stylosanthes spp. (stylo) are annual or perennial legume forages that are widely grown as forage and cover crops in tropical and subtropical regions. However, the seed yield of stylo is very low due to serious seed shattering. In the present study, we found that, although seed shattering was common among the stylo accessions, the shattering rates were genetically different. Therefore, we first synthesized the morphological, histological characteristic, physiochemical, and transcriptome analyses to determine the seed shattering mechanism in stylo. In general, the stylo germplasm with shorter lobules and thicker stems had a lower seed shattering rate and a higher seed weight. The seed and seed stalk joint is the abscission zone in stylo. Multiplex histology and hydrolytic enzyme activity analysis showed that the tearing of the abscission zone occurs due to the intense enzymatic degradation of polygalacturonase and cellulase in the seed shattering-susceptible accession TF0275. cDNA libraries from the abscission zone tissues of TF0041 and TF0275 at 14, 21, and 28 days after flowering were constructed and sequenced. A total of 47,606 unigenes were annotated and 18,606 differentially expressed genes (DEGs) were detected, including 9,140 upregulated and 9,446 downregulated unigenes. Furthermore, the 26 candidate DEGs involved in lignin biosynthesis, cellulase synthesis, and plant hormone signal transduction were found at all three developmental stages. This study provides valuable insights for future mechanistic studies of seed shattering in stylo.

Altitudinal Patterns in Adaptive Evolution of Genome Size and Inter-Genome Hybridization Between Three Elymus Species From the Qinghai–Tibetan Plateau

Genome size variation and hybridization occur frequently within or between plant species under diverse environmental conditions, which enrich species diversification and drive the evolutionary process. Elymus L. is the largest genus in Triticeae with five recognized basic genomes (St, H, P, W, and Y). However, the data on population cytogenetics of Elymus species are sparse, especially whether genome hybridization and chromosomal structure can be affected by altitude are still unknown. In order to explore the relationship between genome sizes, we studied interspecific hybridization and altitude of Elymus species at population genetic and cytological levels. Twenty-seven populations at nine different altitudes (2,800–4,300 m) of three Elymus species, namely, hexaploid E. nutans (StHY, 2n = 6x = 42), tetraploid E. burchan-buddae (StY, 2n = 4x = 28), and E. sibiricus (StH, 2n = 4x = 28), were sampled from the Qinghai–Tibetan Plateau (QTP) to estimate whether intraspecific variation could affect the genomic relationships by genomic in situ hybridization (GISH), and quantify the genome size of Elymus among different altitude ecological groups by flow cytometry. The genome size of E. nutans, E. burchan-buddae, and E. sibiricus varied from 12.38 to 22.33, 8.81 to 18.93, and 11.46 to 20.96 pg/2C with the averages of 19.59, 12.39, and 16.85 pg/2C, respectively. The curve regression analysis revealed a strong correlation between altitude and nuclear DNA content in three Elymus species. In addition, the chromosomes of the St and Y genomes demonstrated higher polymorphism than that of the H genome. Larger genome size variations occurred in the mid-altitude populations (3,900–4,300 m) compared with other-altitude populations, suggesting a notable altitudinal pattern in genome size variation, which shaped genome evolution by altitude. This result supports our former hypothesis that genetic richness center at medium altitude is useful and valuable for species adaptation to highland environmental conditions, germplasm utilization, and conservation.

Comparative Transcriptomics and Metabolites Analysis of Two Closely Related Euphorbia Species Reveal Environmental Adaptation Mechanism and Active Ingredients Difference

Roots of Euphorbia fischeriana and Euphorbia ebracteolata are recorded as the source plant of traditional Chinese medicine "Langdu," containing active ingredients with anticancer and anti-AIDS activity. However, the two species have specific patterns in the graphic distribution. Compared with E. ehracteolata, E. fischeriana distributes in higher latitude and lower temperature areas and might have experienced cold stress adaptation. To reveal the molecular mechanism of environmental adaptation, RNA-seq was performed toward the roots, stems, and leaves of E. fischeriana and E. ehracteolata. A total of 6,830 pairs of putative orthologs between the two species were identified. Estimations of non-synonymous or synonymous substitution rate ratios for these orthologs indicated that 533 of the pairs may be under positive selection (Ka/Ks > 0.5). Functional enrichment analysis revealed that significant proportions of the orthologs were in the TCA cycle, fructose and mannose metabolism, starch and sucrose metabolism, fatty acid biosynthesis, and terpenoid biosynthesis providing insights into how the two closely related Euphorbia species adapted differentially to extreme environments. Consistent with the transcriptome, a higher content of soluble sugars and proline was obtained in E. fischeriana, reflecting the adaptation of plants to different environments. Additionally, 5 primary or secondary metabolites were screened as the biomarkers to distinguish the two species. Determination of 4 diterpenoids was established and performed, showing jolkinolide B as a representative component in E. fischeriana, whereas ingenol endemic to E. ebracteolate. To better study population genetics, EST-SSR markers were generated and tested in 9 species of Euphorbia. A total of 33 of the 68 pairs were screened out for producing clear fragments in at least four species, which will furthermore facilitate the studies on the genetic improvement and phylogenetics of this rapidly adapting taxon. In this study, transcriptome and metabolome analyses revealed the evolution of genes related to cold stress tolerance, biosynthesis of TCA cycle, soluble sugars, fatty acids, and amino acids, consistent with the molecular strategy that genotypes adapting to environment. The key active ingredients of the two species were quantitatively analyzed to reveal the difference in pharmacodynamic substance basis and molecular mechanism, providing insights into rational crude drug use.

Development and Validation of SNP and InDel Markers for Pod-Shattering Tolerance in Soybean

Pod-shattering causes a significant yield loss in many soybean cultivars. Shattering-tolerant cultivars provide the most effective approach to minimizing this loss. We developed molecular markers for pod-shattering and validated them in soybeans with diverse genetic backgrounds. The genes Glyma.16g141200, Glyma.16g141500, and Glyma.16g076600, identified in our previous study by quantitative trait locus (QTL) mapping and whole-genome resequencing, were selected for marker development. The whole-genome resequencing of three parental lines (one shattering-tolerant and two shattering-susceptible) identified single nucleotide polymorphism (SNP) and/or insertion/deletion (InDel) regions within or near the selected genes. Two SNPs and one InDel were converted to Kompetitive Allele-Specific PCR (KASP) and InDel markers, respectively. The accuracy of the markers was examined in the two recombinant inbred line populations used for the QTL mapping, as well as the 120 varieties and elite lines, through allelic discrimination and phenotyping by the oven-drying method. Both types of markers successfully discriminated the pod shattering-tolerant and shattering-susceptible genotypes. The prediction accuracy, which was as high as 90.9% for the RILs and was 100% for the varieties and elite lines, also supported the accuracy and usefulness of these markers. Thus, the markers can be used effectively for genetic and genomic studies and the marker-assisted selection for pod-shattering tolerance in soybean.

Seed Shattering: A Trait of Evolutionary Importance in Plants

Seed shattering refers to the natural shedding of seeds when they ripe, a phenomenon typically observed in wild and weedy plant species. The timing and extent of this phenomenon varies considerably among plant species. Seed shattering is primarily a genetically controlled trait; however, it is significantly influenced by environmental conditions, management practices and their interactions, especially in agro-ecosystems. This trait is undesirable in domesticated crops where consistent efforts have been made to minimize it through conventional and molecular breeding approaches. However, this evolutionary trait serves as an important fitness and survival mechanism for most weeds that utilize it to ensure efficient dispersal of their seeds, paving the way for persistent soil seedbank development and sustained future populations. Weeds have continuously evolved variations in seed shattering as an adaptation under changing management regimes. High seed retention is common in many cropping weeds where weed maturity coincides with crop harvest, facilitating seed dispersal through harvesting operations, though some weeds have notoriously high seed shattering before crop harvest. However, high seed retention in some of the most problematic agricultural weed species such as annual ryegrass (Lolium rigidum), wild radish (Raphanus raphanistrum), and weedy amaranths (Amaranthus spp.) provides an opportunity to implement innovative weed management approaches such as harvest weed seed control, which aims at capturing and destroying weed seeds retained at crop harvest. The integration of such management options with other practices is important to avoid the rapid evolution of high seed shattering in target weed species. Advances in genetics and molecular biology have shown promise for reducing seed shattering in important crops, which could be exploited for manipulating seed shattering in weed species. Future research should focus on developing a better understanding of various seed shattering mechanisms in plants in relation to changing climatic and management regimes.

Development of EST-SSR Markers Linked to Flowering Candidate Genes in Elymus sibiricus L. Based on RNA Sequencing

Elymus sibiricus L. is an important cold-season grass with excellent cold and drought tolerance, good palatability, and nutrition. Flowering time is a key trait that affects forage and seed yield. Development of EST-SSR (expressed sequence tag simple sequence repeat) markers based on flowering genes contributes to the improvement of flowering traits. In the study, we detected 155 candidate genes related to flowering traits from 10,591 unigenes via transcriptome sequencing in early- and late-flowering genotypes. These candidate genes were mainly involved in the photoperiodic pathway, vernalization pathway, central integrator, and gibberellin pathway. A total of 125 candidate gene-based EST-SSRs were developed. Further, 15 polymorphic EST-SSRs closely associated to 13 candidate genes were used for genetic diversity and population structure analysis among 20 E. sibiricus accessions, including two contrasting panels (early-flowering and late-flowering). Among them, primer 28366, designed from heading date 3a (HD3a), effectively distinguished early- and late-flowering genotypes using a specifically amplified band of 175 bp. The polymorphic information content (PIC) value ranged from 0.12 to 0.48, with an average of 0.25. The unweighted pair group method analysis (UPGMA) cluster and structure analysis showed that the 20 E. sibiricus genotypes with similar flowering times tended to group together. These newly developed EST-SSR markers have the potential to be used for molecular markers assisted selection and germplasm evaluation of flowering traits in E. sibiricus.

Development of SSR markers based on transcriptome data and association mapping analysis for fruit shell thickness associated traits in oil palm (Elaeis guineensis Jacq.)

Present study mainly aimed to ascertain the distribution characteristics of gene-based microsatellite loci and to develop polymorphic SSR markers from the already available transcriptome data of Elaeis guineensis Jacq, an important oil crop. From this study, we identified the sum of 5791 SSRs across 51,425 unigenes from the transcripts of oil palm. We were able to evaluate 331primer pairs and characterized 183 polymorphic gene-based SSR markers. We identified a total of 506 alleles from the 183 polymorphic SSR loci, with an average of 2.77 alleles per locus. The characterized gene-based SSR markers from the transcriptome data of oil palm exhibited moderate levels of polymorphism with a significant level of heterozygosity ranges from 0.096 to 0.594 (mean = 0.336 ± 0.11). Among the identified SSR markers, sixty polymorphic markers were used to analyze genotypes of 55 oil palm accessions selected from three different provinces of China. Association mapping analysis provided the information of four markers that are associated with fruit shell thickness trait of oil palm. Among the four markers identified from association analysis, one SSR marker obtained from Unigene17150 is strictly associated with the oil palm fruit shell thickness trait.

Construction of the first high-density genetic linkage map and identification of seed yield-related QTLs and candidate genes in Elymus sibiricus, an important forage grass in Qinghai-Tibet Plateau

Background

Elymus sibiricus is an ecologically and economically important perennial, self-pollinated, and allotetraploid (StStHH) grass, widely used for forage production and animal husbandry in Western and Northern China. However, it has low seed yield mainly caused by seed shattering, which makes seed production difficult for this species. The goals of this study were to construct the high-density genetic linkage map, and to identify QTLs and candidate genes for seed-yield related traits.

Results

An F₂ mapping population of 200 individuals was developed from a cross between single genotype from “Y1005” and “ZhN06”. Specific-locus amplified fragment sequencing (SLAF-seq) was applied to construct the first genetic linkage map. The final genetic map included 1971 markers on the 14 linkage groups (LGs) and was 1866.35 cM in total. The length of each linkage group varied from 87.67 cM (LG7) to 183.45 cM (LG1), with an average distance of 1.66 cM between adjacent markers. The marker sequences of E. sibiricus were compared to two grass genomes and showed 1556 (79%) markers mapped to wheat, 1380 (70%) to barley. Phenotypic data of eight seed-related traits (2016–2018) were used for QTL identification. A total of 29 QTLs were detected for eight seed-related traits on 14 linkage groups, of which 16 QTLs could be consistently detected for two or three years. A total of 6 QTLs were associated with seed shattering. Based on annotation with wheat and barley genome and transcriptome data of abscission zone in E. sibiricus, we identified 30 candidate genes for seed shattering, of which 15, 7, 6 and 2 genes were involved in plant hormone signal transcription, transcription factor, hydrolase activity and lignin biosynthetic pathway, respectively.

Conclusion

This study constructed the first high-density genetic linkage map and identified QTLs and candidate genes for seed-related traits in E. sibiricus. Results of this study will not only serve as genome-wide resources for gene/QTL fine mapping, but also provide a genetic framework for anchoring sequence scaffolds on chromosomes in future genome sequence assembly of E. sibiricus.