Genome-Wide Association Studies of Seven Root Traits in Soybean (Glycine max L.) Landraces

Soybean [Glycine max (L.) Merr.], an important oilseed crop, is a low-cost source of protein and oil. In Southeast Asia and Africa, soybeans are widely cultivated for use as traditional food and feed and industrial purposes. Given the ongoing changes in global climate, developing crops that are resistant to climatic extremes and produce viable yields under predicted climatic conditions will be essential in the coming decades. To develop such crops, it will be necessary to gain a thorough understanding of the genetic basis of agronomic and plant root traits. As plant roots generally lie beneath the soil surface, detailed observations and phenotyping throughout plant development present several challenges, and thus the associated traits have tended to be ignored in genomics studies. In this study, we phenotyped 357 soybean landraces at the early vegetative (V2) growth stages and used a 180 K single-nucleotide polymorphism (SNP) soybean array in a genome-wide association study (GWAS) conducted to determine the phenotypic relationships among root traits, elucidate the genetic bases, and identify significant SNPs associated with root trait-controlling genomic regions/loci. A total of 112 significant SNP loci/regions were detected for seven root traits, and we identified 55 putative candidate genes considered to be the most promising. Our findings in this study indicate that a combined approach based on SNP array and GWAS analyses can be applied to unravel the genetic basis of complex root traits in soybean, and may provide an alternative high-resolution marker strategy to traditional bi-parental mapping. In addition, the identified SNPs, candidate genes, and diverse variations in the root traits of soybean landraces will serve as a valuable basis for further application in genetic studies and the breeding of climate-resilient soybeans characterized by improved root traits.     

Genetic Diversity and Population Structure of Broomcorn Millet (Panicum miliaceum L.) Cultivars and Landraces in China Based on Microsatellite Markers

Broomcorn millet (Panicum miliaceum L.), one of the first domesticated crops, has been grown in Northern China for at least 10,000 years. The species is presently a minor crop, and evaluation of its genetic diversity has been very limited. In this study, we analyzed the genetic diversity of 88 accessions of broomcorn millet collected from various provinces of China. Amplification with 67 simple sequence repeat (SSR) primers revealed moderate levels of diversity in the investigated accessions. A total of 179 alleles were detected, with an average of 2.7 alleles per locus. Polymorphism information content and expected heterozygosity ranged from 0.043 to 0.729 (mean = 0.376) and 0.045 to 0.771 (mean = 0.445), respectively. Cluster analysis based on the unweighted pair group method of mathematical averages separated the 88 accessions into four groups at a genetic similarity level of 0.633. A genetic structure assay indicated a close correlation between geographical regions and genetic diversity. The uncovered information will be valuable for defining gene pools and developing breeding programs for broomcorn millet. Furthermore, the millet-specific SSR markers developed in this study should serve as useful tools for assessment of genetic diversity and elucidation of population structure in broomcorn millet.     

Molecular Diversity Assessment Using Sequence Related Amplified Polymorphism (SRAP) Markers in Vicia faba L.

Sequence-related amplified polymorphism (SRAP) markers were used to assess the genetic diversity and relationship among 58 faba bean (Vicia faba L.) genotypes. Fourteen SRAP primer combinations amplified a total of 1036 differently sized well-resolved peaks (fragments), of which all were polymorphic with a 0.96 PIC value and discriminated all of the 58 faba bean genotypes. An average pairwise similarity of 21% was revealed among the genotypes ranging from 2% to 65%. At a similarity of 28%, UPGMA clustered the genotypes into three main groups comprising 78% of the genotypes. The local landraces and most of the Egyptian genotypes in addition to the Sudan genotypes were grouped in the first main cluster. The advanced breeding lines were scattered in the second and third main clusters with breeding lines from the ICARDA and genotypes introduced from Egypt. At a similarity of 47%, all the genotypes formed separated clusters with the exceptions of Hassawi 1 and Hassawi 2. Group analysis of the genotypes according to their geographic origin and type showed that the landraces were grouped according to their origin, while others were grouped according to their seed type. To our knowledge, this is the first application of SRAP markers for the assessment of genetic diversity in faba bean. Such information will be useful to determine optimal breeding strategies to allow continued progress in faba bean breeding.     

Genetic Variation in Safflower (Carthamus tinctorious L.) for Seed Quality-Related Traits and Inter-Simple Sequence Repeat (ISSR) Markers

Safflower (Carthamus tinctorious L.) is an oilseed crop that is valued as a source of high quality vegetable oil. The genetic diversity of 16 safflower genotypes originated from different geographical regions of Iran and some with exotic origin were evaluated. Eight different seed quality-related traits including fatty acid composition of seed oil (stearic acid, palmitic acid, oleic acid and linoleic acid), the contents of, oil, protein, fiber and ash in its seeds, as well as 20 inter-simple sequence repeat (ISSR) polymorphic primers were used in this study. Analysis of variance showed significant variation in genotypes for the seed quality-related traits. Based on ISSR markers, a total of 204 bands were amplified and 149 bands (about 70%) of these were polymorphic. Cluster analysis based on either biochemical or molecular markers classified the genotypes into four groups, showing some similarities between molecular and biochemical markers for evaluated genotypes. A logical similarity between the genotype clusters based on molecular data with their geographical origins was observed.     

Isolation and Characterization of New Microsatellite Markers for the Invasive Softshell Clam, Mya arenaria (L.) (Bivalvia: Myidae)

The invasive softshell clam (Mya arenaria Linnaeus, 1758) is native to the northwestern region of the Atlantic Ocean. This species has been introduced in the northeast Pacific and along the European coasts, due to intense naval transports and aquaculture, and it is now present in all the European seas. In this paper we describe seven new microsatellite loci for Mya arenaria. The isolated loci are polymorphic with a number of alleles per locus between 6 and 14. The observed and expected heterozygosities ranged from 0.417 to 0.951, and from 0.643 to 0.895, with an average of 0.716 and 0.775, respectively. These microsatellite markers should be useful in analyzing this species' genetic diversity, which could explain various processes of its invasion history.     

Genome-Wide Identification and Characterization of the CC-NBS-LRR Gene Family in Cucumber (Cucumis sativus L.)

The NBS-LRR (NLR) gene family plays a pivotal role in regulating disease defense response in plants. Cucumber is one of the most important vegetable crops in the world, and various plant diseases, including powdery mildew (PM), cause severe losses in both cucumber productivity and quality annually. To characterize and understand the role of the CC-NBS-LRR(CNL) family of genes in disease defense response in cucumber plants, we performed bioinformatical analysis to characterize these genes systematically. We identified 33 members of the CNL gene family in cucumber plants, and they are distributed on each chromosome with chromosome 4 harboring the largest cluster of five different genes. The corresponding CNL family member varies in the number of amino acids and exons, molecular weight, theoretical isoelectric point (pI) and subcellular localization. Cis-acting element analysis of the CNL genes reveals the presence of multiple phytohormone, abiotic and biotic responsive elements in their promoters, suggesting that these genes might be responsive to plant hormones and stress. Phylogenetic and synteny analysis indicated that the CNL proteins are conserved evolutionarily in different plant species, and they can be divided into four subfamilies based on their conserved domains. MEME analysis and multiple sequence alignment showed that conserved motifs exist in the sequence of CNLs. Further DNA sequence analysis suggests that CsCNL genes might be subject to the regulation of different miRNAs upon PM infection. By mining available RNA-seq data followed by real-time quantitative PCR (qRT-PCR) analysis, we characterized expression patterns of the CNL genes, and found that those genes exhibit a temporospatial expression pattern, and their expression is also responsive to PM infection, ethylene, salicylic acid, and methyl jasmonate treatment in cucumber plants. Finally, the CNL genes targeted by miRNAs were predicted in cucumber plants. Our results in this study provided some basic information for further study of the functions of the CNL gene family in cucumber plants.     

Development of Microsatellite Markers in the Branched Broomrape Phelipanche ramosa L. (Pomel) and Evidence for Host-Associated Genetic Divergence

Phelipanche ramosa is a parasitic plant that infects numerous crops worldwide. In Western Europe it recently expanded to a new host crop, oilseed rape, in which it can cause severe yield losses. We developed 13 microsatellite markers for P. ramosa using next-generation 454 sequencing data. The polymorphism at each locus was assessed in a sample of 96 individuals collected in France within 6 fields cultivated with tobacco, hemp or oilseed rape. Two loci were monomorphic. At the other 11 loci, the number of alleles and the expected heterozygosity ranged from 3 to 6 and from 0.31 to 0.60, respectively. Genetic diversity within each cultivated field was very low. The host crop from which individuals were collected was the key factor structuring genetic variation. Individuals collected on oilseed rape were strongly differentiated from individuals collected on hemp or tobacco, which suggests that P. ramosa infecting oilseed rape forms a genetically diverged race. The microsatellites we developed will be useful for population genetics studies and for elucidating host-associated genetic divergence in P. ramosa.     

Isolation and Characterization of 13 New Polymorphic Microsatellite Markers in the Phaseolus vulgaris L. (Common Bean) Genome

In this study, 13 polymorphic microsatellite markers were isolated from the Phaseolus vulgaris L. (common bean) by using the Fast Isolation by AFLP of Sequence COntaining Repeats (FIASCO) protocol. These markers revealed two to seven alleles, with an average of 3.64 alleles per locus. The polymorphic information content (PIC) values ranged from 0.055 to 0.721 over 13 loci, with a mean value of 0.492, and 7 loci having PIC greater than 0.5. The expected heterozygosity (H_E) and observed heterozygosity (H_O) levels ranged from 0.057 to 0.814 and from 0.026 to 0.531, respectively. Cross-species amplification of the 13 prime pairs was performed in its related specie of Vigna unguiculata L. Seven out of all these markers showed cross-species transferability. These markers will be useful for future genetic diversity and population genetics studies for this agricultural specie and its related species.     

Genome-Wide Association Mapping Identifies New Candidate Genes for Cold Stress and Chilling Acclimation at Seedling Stage in Rice (Oryza sativa L.)

Rice (Oryza sativa L.) is a chilling-sensitive staple food crop, and thus, low temperature significantly affects rice growth and yield. Many studies have focused on the cold shock of rice although chilling acclimation is more likely to happen in the field. In this paper, a genome-wide association study (GWAS) was used to identify the genes that participated in cold stress and chilling accumulation. A total of 235 significantly associated single-nucleotide polymorphisms (SNPs) were identified. Among them, we detected 120 and 88 SNPs for the relative shoot fresh weight under cold stress and chilling acclimation, respectively. Furthermore, 11 and 12 quantitative trait loci (QTLs) were identified for cold stress and chilling acclimation, respectively, by integrating the co-localized SNPs. Interestingly, we identified 10 and 15 candidate genes in 11 and 12 QTLs involved in cold stress and chilling acclimation, respectively, and two new candidate genes (LOC_Os01g62410, LOC_Os12g24490) were obviously up-regulated under chilling acclimation. Furthermore, OsMYB3R-2 (LOC_Os01g62410) that encodes a R1R2R3 MYB gene was associated with cold tolerance, while a new C3HC4-type zinc finger protein-encoding gene LOC_Os12g24490 was found to function as a putative E3 ubiquitin-protein ligase in rice. Moreover, haplotype, distribution, and Wright’s fixation index (FST) of both genes showed that haplotype 3 of LOC_Os12g24490 is more stable in chilling acclimation, and the SNP (A > T) showed a difference in latitudinal distribution. FST analysis of SNPs in OsMYB3R-2 (LOC_Os01g62410) and LOC_Os12g24490 indicated that several SNPs were under selection in rice indica and japonica subspecies. This study provided new candidate genes in genetic improvement of chilling acclimation response in rice.     

Isolation and Characterisation of 11 Polymorphic Microsatellite Markers in Papaver rhoeas L. (Corn Poppy), a Major Annual Plant Species from Cultivated Areas

Papaver rhoeas, an annual plant species in the Papaveraceae family, is part of the biodiversity of agricultural ecosystems and also a noxious agronomic weed. We developed microsatellite markers to study the genetic diversity of P. rhoeas, using an enriched microsatellite library coupled with 454 next-generation sequencing. A total of 13,825 sequences were obtained that yielded 1795 microsatellite loci. After discarding loci with less than six repeats of the microsatellite motif, automated primer design was successful for 598 loci. We tested 74 of these loci for amplification with a total of 97 primer pairs. Thirty loci passed our tests and were subsequently tested for polymorphism using 384 P. rhoeas plants originating from 12 populations from France. Of the 30 loci, 11 showed reliable polymorphism not affected by the presence of null alleles. The number of alleles and the expected heterozygosity ranged from 3 to 7.4 and from 0.27 to 0.73, respectively. A low but significant genetic differentiation among populations was observed (F_ST = 0.04; p < 0.001). The 11 validated polymorphic microsatellite markers developed in this work will be useful in studies of genetic diversity and population structure of P. rhoeas, assisting in designing management strategies for the control or the conservation of this species.     

木本油料植物山核桃新品种的选育

通过遗传多样性分析,确立山核桃选优策略.对山核桃进行大范围的选优,采取常规选育和分子标记相结合,初步筛选了山核桃优株143株,复选出高产稳产型、早熟型、大果型3种类型的50个优良单株,并对临安、安吉、淳安3地区山核桃复选优系的鲜果质量、干果质量、出籽率、出仁率、粗脂肪、粗蛋白6个经济性状的变异作了比较分析.结果表明: 1)所试山核桃优系有进一步选育品种的潜力,鲜果质量的变异系数最大,粗脂肪的变异系数最小; 2)经济性状在地区间的变化没有一致性,但地域分布对山核桃出籽率、粗蛋白含量有显著影响.     

Genome-Wide Identification,Evolution and Expressional Analysis of OSCA Gene Family in Barley (Hordeum vulgare L.)

The hyperosmolality-gated calcium-permeable channel gene family (OSCA) is one kind of conserved osmosensors, playing a crucial role in maintaining ion and water homeostasis and protecting cellular stability from the damage of hypertonic stress. Although it has been systematically characterized in diverse plants, it is necessary to explore the role of the OSCA family in barley, especially its importance in regulating abiotic stress response. In this study, a total of 13 OSCA genes (HvOSCAs) were identified in barley through an in silico genome search method, which were clustered into 4 clades based on phylogenetic relationships with members in the same clade showing similar protein structures and conserved motif compositions. These HvOSCAs had many cis-regulatory elements related to various abiotic stress, such as MBS and ARE, indicating their potential roles in abiotic stress regulation. Furthermore, their expression patterns were systematically detected under diverse stresses using RNA-seq data and qRT-PCR methods. All of these 13 HvOSCAs were significantly induced by drought, cold, salt and ABA treatment, demonstrating their functions in osmotic regulation. Finally, the genetic variations of the HvOSCAs were investigated using the re-sequencing data, and their nucleotide diversity in wild barley and landrace populations were 0.4966 × 10⁻³ and 0.391 × 10⁻³, respectively, indicating that a genetic bottleneck has occurred in the OSCA family during the barley evolution process. This study evaluated the genomic organization, evolutionary relationship and genetic expression of the OSCA family in barley, which not only provides potential candidates for further functional genomic study, but also contributes to genetically improving stress tolerance in barley and other crops.     

High Genetic Diversity and Low Differentiation of Michelia coriacea (Magnoliaceae), a Critically Endangered Endemic in Southeast Yunnan, China

Michelia coriacea, a critically endangered tree, has a restricted and fragmented distribution in Southeast Yunnan Province, China. The genetic diversity, genetic structure and gene flow in the three extant populations of this species were detected by 10 inter-simple sequence repeat (ISSR) markers and 11 simple sequence repeat (SSR) markers. Examination of genetic diversity revealed that the species maintained a relatively high level of genetic diversity at the species level (percentage of polymorphic bands) PPB = 96.36% from ISSRs; PPL (percentage of polymorphic loci) = 95.56% from SSRs, despite several fragmental populations. Low levels of genetic differentiation among the populations of M. coriacea were detected by Nei's G(st) = 0.187 for ISSR and Wright's F(st) = 0.090 for SSR markers, which is further confirmed by Bayesian model-based STRUCTURE and PCoA analysis that could not reveal a clear separation between populations, although YKP was differentiated to other two populations by ISSR markers. Meanwhile, AMOVA analysis also indicated that 22.84% and 13.90% of genetic variation existed among populations for ISSRs and SSRs, respectively. The high level of genetic diversity, low genetic differentiation, and the population, structure imply that the fragmented habitat and the isolated population of M. coriacea may be due to recent over-exploitation. Conservation and management of M. coriacea should concentrate on maintaining the high level of genetic variability through both in and ex-situ conservation actions.     

Molecular Genetic Variability of Commercial and Wild Accessions of Passion Fruit (Passiflora spp.) Targeting ex Situ Conservation and Breeding

Passiflora species are distributed throughout Latin America, and Brazil and Colombia serve as the centers of diversity for this genus. We performed cross-species amplification to evaluate 109 microsatellite loci in 14 Passiflora species and estimated the diversity and genetic structure of Passiflora cincinnata, Passiflora setaceae and Passiflora edulis. A total of 127 accessions, including 85 accessions of P. edulis, a commercial species, and 42 accessions of 13 wild species, were examined. The cross-species amplification was effective for obtaining microsatellite loci (average cross-amplification of 70%). The average number of alleles per locus (five) was relatively low, and the average diversity ranged from 0.52 in P. cincinnata to 0.32 in P. setacea. The Bayesian analyses indicated that the P. cincinnata and P. setacea accessions were distributed into two groups, and the P. edulis accessions were distributed into five groups. Private alleles were identified, and suggestions for core collections are presented. Further collections are necessary, and the information generated may be useful for breeding and conservation.     

Identification of Candidate Genes for Root Traits Using Genotype–Phenotype Association Analysis of Near-Isogenic Lines in Hexaploid Wheat (Triticum aestivum L.)

Global wheat (Triticum aestivum L.) production is constrained by different biotic and abiotic stresses, which are increasing with climate change. An improved root system is essential for adaptability and sustainable wheat production. In this study, 10 pairs of near-isogenic lines (NILs)—targeting four genomic regions (GRs) on chromosome arms 4BS, 4BL, 4AS, and 7AL of hexaploid wheat—were used to phenotype root traits in a semi-hydroponic system. Seven of the 10 NIL pairs significantly differed between their isolines for 11 root traits. The NIL pairs targeting qDSI.4B.1 GR varied the most, followed by the NIL pair targeting qDT.4A.1 and QHtscc.ksu-7A GRs. For pairs 5–7 targeting qDT.4A.1 GR, pair 6 significantly differed in the most root traits. Of the 4 NIL pairs targeting qDSI.4B.1 GR, pairs 2 and 4 significantly differed in 3 and 4 root traits, respectively. Pairs 9 and 10 targeting QHtscc.ksu-7A GR significantly differed in 1 and 4 root traits, respectively. Using the wheat 90K Illumina iSelect array, we identified 15 putative candidate genes associated with different root traits in the contrasting isolines, in which two UDP-glycosyltransferase (UGT)-encoding genes, TraesCS4A02G185300 and TraesCS4A02G442700, and a leucine-rich repeat receptor-like protein kinase (LRR-RLK)-encoding gene, TraesCS4A02G330900, also showed important functions for root trait control in other crops. This study characterized, for the first time, that these GRs control root traits in wheat, and identified candidate genes, although the candidate genes will need further confirmation and validation for marker-assisted wheat breeding.     

Genetic Diversity and Population Structure of Siberian apricot (Prunus sibirica L.) in China

The genetic diversity and population genetic structure of 252 accessions from 21 Prunus sibirica L. populations were investigated using 10 ISSR, SSR, and SRAP markers. The results suggest that the entire population has a relatively high level of genetic diversity, with populations HR and MY showing very high diversity. A low level of inter-population genetic differentiation and a high level of intra-population genetic differentiation was found, which is supported by a moderate level of gene flow, and largely attributable to the cross-pollination and self-incompatibility reproductive system. A STRUCTURE (model-based program) analysis revealed that the 21 populations can be divided into two main groups, mainly based on geographic differences and genetic exchanges. The entire wild Siberia apricot population in China could be divided into two subgroups, including 107 accessions in subgroup (SG) 1 and 147 accessions in SG 2. A Mantel test revealed a significant positive correlation between genetic and geographic distance matrices, and there was a very significant positive correlation among three marker datasets. Overall, we recommend a combination of conservation measures, with ex situ and in situ conservation that includes the construction of a core germplasm repository and the implement of in situ conservation for populations HR, MY, and ZY.     

Genetic Analysis of Hexaploid Wheat (Triticum aestivum L.) Using the Complete Sequencing of Chloroplast DNA and Haplotype Analysis of the Wknox1 Gene

The aim of the presented study is a genetic characterization of the hexaploid wheat Triticum aestivum L. Two approaches were used for the genealogical study of hexaploid wheats—the complete sequencing of chloroplast DNA and PCR-based haplotype analysis of the fourth intron of Wknox1d and of the fifth-to-sixth-exon region of Wknox1b. The complete chloroplast DNA sequences of 13 hexaploid wheat samples were determined: Free-threshing—T. aestivum subsp. aestivum, one sample; T. aestivum subsp. compactum, two samples; T. aestivum subsp. sphaerococcum, one sample; T. aestivum subsp. carthlicoides, four samples. Hulled—T. aestivum subsp. spelta, three samples; T. aestivum subsp. vavilovii jakubz., two samples. The comparative analysis of complete cpDNA sequences of 20 hexaploid wheat samples (13 samples in this article plus 7 samples sequenced in this laboratory in 2018) was carried out. PCR-based haplotype analysis of the fourth intron of Wknox1d and of the fifth-to-sixth exon region of Wknox1b of all 20 hexaploid wheat samples was carried out. The 20 hexaploid wheat samples (13 samples in this article plus 7 samples in 2018) can be divided into two groups—T. aestivum subsp. spelta, three samples and T. aestivum subsp. vavilovii collected in Armenia, and the remaining 16 samples, including T. aestivum subsp. vavilovii collected in Europe (Sweden). If we take the cpDNA of Chinese Spring as a reference, 25 SNPs can be identified. Furthermore, 13–14 SNPs can be identified in T. aestivum subsp. spelta and subsp. vavilovii (Vav1). In the other samples up to 11 SNPs were detected. 22 SNPs are found in the intergenic regions, 2 found in introns, and 10 SNPs were found in the genes, of which seven are synonymous. PCR-based haplotype analysis of the fourth intron of Wknox1d and the fifth-to-sixth-exon region of Wknox1b provides an opportunity to make an assumption that hexaploid wheats T. aestivum subsp. macha var. palaeocolchicum and var. letshckumicum differ from other macha samples by the absence of a 42 bp insertion in the fourth intron of Wknox1d. One possible explanation for this observation would be that two Aegilops tauschii Coss. (A) and (B) participated in the formation of hexaploids through the D genome: Ae. tauschii (A)—macha (1–5, 7, 8, 10–12), and Ae. tauschii (B)—macha M6, M9, T. aestivum subsp. aestivum cv. ‘Chinese Spring’ and cv. ‘Red Doly’.     

Isolation and Characterization of the First Microsatellite Markers for the Endangered Relict Mussel Hypanis colorata (Mollusca: Bivalvia: Cardiidae)

Hypanis colorata (Eichwald, 1829) (Cardiidae: Lymnocardiinae) is a bivalve relict species with a Ponto-Caspian distribution and is under strict protection in Romania, according to national regulations. While the species is depressed in the western Black Sea lagoons from Romania and Ukraine, it is also a successful invader in the middle Dniepr and Volga regions. Establishing a conservation strategy for this species or studying its invasion process requires knowledge about the genetic structure of the species populations. We have isolated and characterized nine polymorphic microsatellite markers in H. colorata. The number of alleles per locus ranged from 4 to 28 and the observed heterozygosity ranged from 0.613 to 1.000. The microsatellites developed in the present study are highly polymorphic and they should be useful for the assessment of genetic variation within this species.