acephala) using cross-species microsatellite markers

Kale is known for its exceptional nourishing and functional benefits to human body. However, it is an understudied species from genomic as well as agronomic aspects. It is important to characterize niche kale germplasms around the world to systematically conserve and utilize its genetic variability, especially for commercial traits in the interest of growers, consumers and industry. With this view, genomic and phenotypic characterizations of 62 Kashmiri kale accessions including popular landraces were done to estimate and partition genetic diversity, understand trait relationships, develop population structure and divulge marker-trait associations of economic significance. Sixty-six cross species microsatellite (SSR) markers within Brassica genus amplified 269 alleles in the germplasm. Their polymorphic information content (PIC) ranged from 0.00078 to 0.953 with an average of 0.407. The population structure analysis and neighbour joining tree clustering categorized the germplasm into three sub-populations. AMOVA revealed more within-population variance (67.73 %) than among-populations (32.27 %) variance. The principal component analysis (PCA) involving 24 agronomical traits revealed seven PCs (PC1 to PC7) having Eigen values more than 1, which explained a cumulative variation of 69.21 %. Association mapping with respect to these 24 agronomical traits using mixed linear model and general linear model revealed six overlapping significant marker-trait relationships with five being significant at probability value of 0.001/0.0001. The highly significant associations of two SSRs with economically important traits (siliqua length and seed weight) significantly correlated/related with leaf yield and seed yield were revealed for their possible utilization in marker assisted breeding for higher leaf and seed yields.

Association study of morpho-phenological traits in quinoa (Chenopodium quinoa Willd.) using SSR markers

In this study, the genetic and molecular diversity of 60 quinoa accessions was assessed using agronomically important traits related to grain yield as well as microsatellite (SSR) markers, and informative markers linked to the studied traits were identified using association study. The results showed that most of the studied traits had a relatively high diversity, but grain saponin and protein content showed the highest diversity. High diversity was also observed in all SSR markers, but KAAT023, KAAT027, KAAT036, and KCAA014 showed the highest values for most of the diversity indices and can be introduced as the informative markers to assess genetic diversity in quinoa. Population structure analysis showed that the studied population probably includes two subclusters, so that out of 60 quinoa accessions, 29 (48%) and 23 (38%) accessions were assigned to the first and second subclusters, respectively, and eight (13%) accessions were considered as the mixed genotypes. The study of the population structure using Structure software showed two possible subgroups (K = 2) in the studied population and the results of the bar plot confirmed it. Association study using the general linear model (GLM) and mixed linear model (MLM) identified the number of 35 and 32 significant marker-trait associations (MTAs) for the first year (2019) and 37 and 35 significant MTAs for the second year (2020), respectively. Among the significant MTAs identified for different traits, the highest number of significant MTAs were obtained for grain yield and 1000-grain weight with six and five MTAs, respectively.

Developing Chenopodium ficifolium as a potential B genome diploid model system for genetic characterization and improvement of allotetraploid quinoa (Chenopodium quinoa)

BACKGROUND

Quinoa (Chenopodium quinoa) is a high-value grain known for its excellent nutritional balance. It is an allotetraploid species (AABB, 2n = 4x = 36) formed by the hybridization between AA and BB genome diploid (2n = 2x = 18) species. This study reports genetic studies in Chenopodium ficifolium as a potential B genome diploid model system to simplify the genetic studies of quinoa including gene identification and marker-assisted breeding.

RESULTS

Portsmouth, New Hampshire and Quebec City, Quebec accessions of C. ficifolium were used to develop an F2 population segregating for agronomically relevant traits including flowering time, plant height, the number of branches, branch angle, and internode length. Marker-trait associations were identified for the FLOWERING LOCUS T-LIKE 1 (FTL1) marker gene, where the alternate alleles (A1/A2) were segregating among the F2 generation plants in association with flowering time, plant height, and the number of branches. There was a strong correlation of the flowering time trait with both plant height and the number of branches. Thus, a possible multifaceted functional role for FTL1 may be considered. The parental Portsmouth and Quebec City accessions were homozygous for the alternate FTL1 alleles, which were found to be substantially diverged. SNPs were identified in the FTL1 coding sequence that could have some functional significance in relation to the observed trait variation.

CONCLUSION

These results draw further attention to the possible functional roles of the FTL1 locus in Chenopodium and justify continued exploration of C. ficifolium as a potential diploid model system for the genetic study of quinoa. We expect our findings to aid in quinoa breeding as well as to any studies related to the Chenopodium genus.

Exploring natural diversity reveals alleles to enhance antioxidant system in barley under salt stress

Soil salinity stress causes osmotic/ionic imbalances and induces oxidative stress that causes cellular structure damage, perturbs metabolism, antioxidant system (comprising enzymatic and non-enzymatic components) and hence inhibits plant growth performance. In this study, we used genome-wide association scan (GWAS) in 174 diverse spring barley accessions which were exposed to salt stress under field conditions at the vegetative stage to uncover the genetic basis of antioxidant components and agronomic traits. High activities of enzymatic and content of non-enzymatic antioxidants were observed under salt stress compared to control conditions. Under salt stress, all the agronomic and yield-related traits were significantly reduced. Six genomic regions were associated with antioxidants and agronomic traits under salt stress conditions which were found to be linked with candidate genes. Several significant associations were physically located inside or near genes which are potentially involved in antioxidants. Two candidate genes at 2H (40,659,364 bp) and 7H (416,743,127 bp) were found to be involved in Dihydroflavonol 4-reductase/flavanone protein and Glyceraldehyde-3-phosphate dehydrogenase, respectively. The allelic variation at SNP of BK_07 at 7H inside the GAPDH gene demonstrates a negative selection of accessions carrying A allele. This allele appears in cultivars with lower activity of enzymatic antioxidants e.g. superoxide dismutase and catalases under salt stress conditions. These accessions are predominantly two-rowed, cultivars, originated from Europe, and carrying photoperiod sensitive alleles. The detected associated molecular markers in this work are considered as an important source for selection of increased amount of antioxidant compounds in barley under stress conditions.

Genome wide association mapping for heat tolerance in sub-tropical maize

Background

Heat tolerance is becoming increasingly important where maize is grown under spring season in India which coincide with grain filling stage of crop resulting in tassel blast, reduced pollen viability, pollination failure and barren ears that causes devastating yield losses. So, there is need to identify the genomic regions associated with heat tolerance component traits which could be further employed in maize breeding program.

Results

An association mapping panel, consisting of 662 doubled haploid (DH) lines, was evaluated for yield contributing traits under normal and natural heat stress conditions. Genome wide association studies (GWAS) carried out using 187,000 SNPs and 130 SNPs significantly associated for grain yield (GY), days to 50% anthesis (AD), days to 50% silking (SD), anthesis-silking interval (ASI), plant height (PH), ear height (EH) and ear position (EPO) were identified under normal conditions. A total of 46 SNPs strongly associated with GY, ASI, EH and EPO were detected under heat stress conditions. Fifteen of the SNPs was found to have common association with more than one trait such as two SNPs viz. S10_1,905,273 and S10_1,905,274 showed colocalization with GY, PH and EH whereas S10_7,132,845 SNP associated with GY, AD and SD under normal conditions. No such colocalization of SNP markers with multiple traits was observed under heat stress conditions. Haplotypes trend regression analysis revealed 122 and 85 haplotype blocks, out of which, 20 and 6 haplotype blocks were associated with more than one trait under normal and heat stress conditions, respectively. Based on SNP association and haplotype mapping, nine and seven candidate genes were identified respectively, which belongs to different gene models having different biological functions in stress biology.

Conclusions

The present study identified significant SNPs and haplotype blocks associated with yield contributing traits that help in selection of donor lines with favorable alleles for multiple traits. These results provided insights of genetics of heat stress tolerance. The genomic regions detected in the present study need further validation before being applied in the breeding pipelines.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07463-y.

Genome-wide association studies for agronomical traits in winter rice accessions of Assam

The 297 winter rice accessions of Assam, North East India were genotyped by sequencing (GBS). The 50,985 high-quality SNPs were filtered and assigned to 12 rice chromosomes. The population structure analysis revealed three major subgroups SG1, SG2, and SG3 consisting of 30, 8, and 143 accessions respectively. The remaining 116 accessions were grouped as admixture population. Phenotypic data were recorded on13 agronomical traits for genome-wide association studies (GWAS). The 60 significant marker-trait associations (MTAs) were identified for 11 agronomical traits, which explained 0 to 15% of phenotypic variance (PV). A QTL 'hot spot' was detected near the centromeric region on chromosome 6. The identified QTLs may be validated and utilized in 'genomics assisted breeding' for improvement of existing rice cultivars of Assam and North East India.

Genome-wide association analysis of stripe rust resistance in modern Chinese wheat

BACKGROUND

Stripe rust (yellow rust) is a significant disease for bread wheat (Triticum aestivum L.) worldwide. A genome-wide association study was conducted on 240 Chinese wheat cultivars and elite lines genotyped with the wheat 90 K single nucleotide polymorphism (SNP) arrays to decipher the genetic architecture of stripe rust resistance in Chinese germplasm.

RESULTS

Stripe rust resistance was evaluated at the adult plant stage in Pixian and Xindu in Sichuan province in the 2015-2016 cropping season, and in Wuhan in Hubei province in the 2013-2014, 2016-2017 and 2018-2019 cropping seasons. Twelve stable loci for stripe rust resistance were identified by GWAS using TASSEL and GAPIT software. These loci were distributed on chromosomes 1B, 1D, 2A, 2B, 3A, 3B, 4B (3), 4D, 6D, and 7B and explained 3.6 to 10.3% of the phenotypic variation. Six of the loci corresponded with previously reported genes/QTLs, including Sr2/Yr30/Lr27, while the other six (QYr.hbaas-1BS, QYr.hbaas-2BL, QYr.hbaas-3AL, QYr.hbaas-4BL.3, QYr.hbaas-4DL, and QYr.hbaas-6DS) are probably novel. The results suggest high genetic diversity for stripe rust resistance in this population. The resistance alleles of QYr.hbaas-2AS, QYr.hbaas-3BS, QYr.hbaas-4DL, and QYr.hbaas-7BL were rare in the present panel, indicating their potential use in breeding for stripe rust resistance in China. Eleven penta-primer amplification refractory mutation system (PARMS) markers were developed from SNPs significantly associated with seven mapped QTLs. Twenty-seven genes were predicted for mapped QTLs. Six of them were considered as candidates for their high relative expression levels post-inoculation.

CONCLUSION

The resistant germplasm, mapped QTLs, and PARMS markers developed in this study are resources for enhancing stripe rust resistance in wheat breeding.

Association of molecular markers with physio-biochemical traits related to seed vigour in rice

Eighteen physio-biochemical traits influencing seed vigour were studied for their association with molecular markers using a mini core set constituted from 120 germplasm lines. High genetic variation was detected in the parameters namely chlrophyll a, Chlrophyll b, total chlorophyll, carotenoids, total anthocyanin content, gamma-oryzanols, total phenolics content, superoxide dismutase, catalase, guaicol peroxidase, total soluble sugar, total protein, seed vigour index -I and seed vigour index -II. Strong positive correlation of seed vigour index II was observed with amylose content, total anthocyanin content, catalase, total phenolic content and total flavonoid content while a negative association was observed for gamma-oryzanol content. High gene diversity (0.7169) and informative markers value (0.6789) were estimated from the investigation. Three genetic structure groups were observed in the panel population and genotypes were grouped in the subpopulations based on the seed vigour trait. Differences in the fixation indices of the three sub populations indicated existence of linkage disequilibrium in the studied panel population. Association of the traits namely total flavonoids, superoxide dismutase, catalase, chlorophyll a, Chlorophyll b, total chlorophyll, carotenoids, starch, amylose, total anthocyanin, gamma-oryzanol, total phenolics with the molecular markers were detected by Generalized Linear Model and Mixed Linear Model showing > 0.10 R² value. Association of the trait, total flavonoids with marker RM7364 located on chromosome 8 reported in earlier study was validated in this investigation. The validated markers and the novel markers detected showing higher R² value will be useful for improvement of seed vigour in rice.

Association mapping for mungbean yellow mosaic India virus resistance in mungbean (Vigna radiata L. Wilczek)

The present study aimed to detect the marker-trait association of a selected diverse panel of 127 mungbean genotypes against mungbean yellow mosaic India virus (MYMIV). Virus-specific primers pairs viz., AC-abut/AV-abut and BC-abut/BV-abut confirmed the involvement of MYMIV in yellow mosaic disease development and the same was validated through restriction digestion analysis. 256 genome-wide microsatellite markers were screened on a test panel in which 93 polymorphic markers were used in association studies. Population structure analysis led to formation of six distinct subpopulations. 1097 alleles were detected among 127 test genotypes whereas number of alleles ranged 2-22 and PIC values ranged 0.27-0.92%, indicating ample amount of variation at genome level. 15 microsatellite markers were detected as associated with MYMIV resistance, among them three microsatellites explained 11-14% phenotypic variation. The specific regions close to CEDG293, DMB-SSR008 and DMB-SSR059 associated with MYMIV resistance were detected, located on linkage group 2, 4 and 9 and may prove useful in marker-assisted mungbean improvement programme for enhancing MYMIV resistance.

Genome-wide association mapping for adult resistance to powdery mildew in common wheat

Blumeria graminis f. sp. tritici, the causal agent of wheat powdery mildew disease, can occur at all stages of the crop and constantly threatens wheat production. To identify candidate resistance genes for powdery mildew, we performed GWAS (genome-wide association studies) on a total set of 329 wheat varieties obtained from different origins. These wheat materials were genotyped using wheat 90K SNP array and evaluated for their resistance in either field or glasshouse condition from 2016 to 2018. Using a mixed linear model, 33 SNP markers of which 14 QTL (quantitative trait loci) were later defined were observed to associate with powdery mildew resistance. Among these, QTL on chromosome 3A, 3B, 6D and 7D were concluded as potentially new QTL. Exploration of candidate genes for new QTL suggested roles of these genes involved in encoding disease resistance and defence-related proteins, and regulating early immune response to the pathogen. Overall, the results reveal that GWAS can be an effective means of identifying marker-trait associations, though further functional validation and fine-mapping of gene candidates are required before creating opportunities for developing new resistant genotypes.

Genome-wide association study reveals significant genomic regions for improving yield, adaptability of rice under dry direct seeded cultivation condition

BACKGROUND

Puddled transplanted system of rice cultivation despite having several benefits, is a highly labor, water and energy intensive system. In the face of changing climatic conditions, a successful transition from puddled to dry direct seeded rice (DDSR) cultivation system looks must in future. Genome-wide association study was performed for traits including, roots and nutrient uptake (14 traits), plant-morphological (5 traits), lodging-resistance (4 traits) and yield and yield attributing traits (7 traits) with the aim to identify significant marker-trait associations (MTAs) for traits enhancing rice adaptability to dry direct-seeded rice (DDSR) system.

RESULTS

Study identified a total of 37 highly significant MTAs for 20 traits. The false discovery rate (FDR) ranged from 0.264 to 3.69 × 10- 4, 0.0330 to 1.25 × 10- 4, and 0.0534 to 4.60 × 10- 6 in 2015WS, 2016DS and combined analysis, respectively. The percent phenotypic variance (PV) explained by SNPs ranged from 9 to 92%. Among the identified significant MTAs, 15 MTAs associated with the traits including nodal root, root hair length, root length density, stem and culm diameter, plant height and grain yield were reported to be located in the proximity of earlier identified candidate gene. The significant positive correlation of grain-yield with seedling establishment traits, root morphological and nutrient-uptake related traits and grain yield attributing traits pointing towards combining target traits to increase rice yield and adaptability under DDSR. Seven promising progenies with better root morphology, nutrient-uptake and higher grain yield were identified that can further be used in genomics assisted breeding for DDSR varietal development.

CONCLUSIONS

Once validated, the identified MTAs and the SNPs linked with trait of interest could be of direct use in genomic assisted breeding (GAB) to improve grain yield and adaptability of rice under DDSR.

Marker-trait association for low-light intensity tolerance in rice genotypes from Eastern India

Light intensity is a crucial environmental factor that affects photosynthesis and ultimately, grain yield in rice. However, no gene or marker directly associated with improved performance under low-light intensity under field conditions has been identified till date. With an aim of identifying genes and markers associated with improved performance (measured in terms of better yields) under low-light intensity, an integrated field screening, in silico and wet lab validation analysis was performed. Field-based screening of a diverse set of 110 genotypes led to the identification of two physiological and three morphological parameters critical for low-light tolerance in rice. In silico analysis using information available in public databases led to the identification of a set of 90 potential candidate genes which were narrowed to thirteen genic targets for possible marker-trait association. Marker-trait association on the panel of 48 diverse rice genotypes varying in their response to low-light intensity led to the identification of six markers [HvSSR02-44 (biological yield), HvSSR02-52 (spikelet fertility), HvSSR02-54 (grain yield), HvSSR06-56 (spikelet fertility), HvSSR06-69 (spikelet fertility; biological yield), HvSSR09-45 (spikelet fertility)] lying on chromosomes 2, 6 and 9 showing significant association (R² > 0.1) for traits like grain yield/plant, biological yield and spikelet fertility under low light. Eight rice genes [including member of BBX (B-box) family] lying within 10 kb distance of these identified markers already reported for their role in response to stress or change in plant architecture in rice were also identified. The eight rice genotypes, five traits, eight genes and six markers identified in the current study will help in devising strategies to increase yield under low light intensity and pave way for future application in marker-assisted breeding.

Identification of SNP markers associated with tolerance to ammonia toxicity by selective genotyping from de novo assembled transcriptome in Litopenaeus vannamei

The high concentration of ammonia from deteriorated aquaculture environments and the intensive culture system could increase the susceptibility to pathogens and even cause high mortality in Litopenaeus vannamei. In addition, we have revealed that the ammonia-tolerant shrimp also have high disease resistance in L. vannamei. In the present study, in order to identify SNP markers associated with tolerance to ammonia toxicity, we developed and characterized SNPs from our previous transcriptome sequencing data of ammonia-stressed and control groups, and a marker-trait association analysis was performed for marker-assisted selection (MAS) to increase production in L. vannamei. A total of 318,919 SNPs were identified from the transcriptome sequences, and 25,772 SNPs were found from the 1826 ammonia-responsive genes with functional annotation. We selected 49 SNPs from 26 ammonia-responsive genes that had strong homologies to known genes in the shrimp and probably involved in immune function as candidate markers for genotyping, among which 39 SNPs were polymorphic for further marker-trait association analysis with the ammonia-tolerant (AT) and ammonia-sensitive (AS) groups. Finally, 12 out of the 49 SNP markers were identified to be associated with ammonia tolerance, containing 10 loci with significantly different allele frequencies and 10 loci with significantly different genotyping frequencies between the AT and AS groups. Among the associated markers, the G allele of TSP-1 (the first locus from the thrombospondin gene), the A allele of TSP-3, and the C allele of XBP1-5 (the fifth locus from X-box binding protein 1) only presented in the AT groups, but they were absent from the AS groups, which would be the preference of the MAS for the ammonia-tolerant shrimp. In addition, when the 12 associated SNP markers were used for analysis, the genetic diversity of the AT groups were significantly higher than that of the AS groups, but when the 39 loci were used there was no difference. This is the first report for the markers associated with ammonia tolerance in this species, indirectly with disease resistance, which provided important potential for genetic selection to increase survival rate and production in shrimp farming.

New Method Application for Marker-Trait Association Studies in Plants: Partial Least Square Regression Aids Detection of Simultaneous Correlations

In this work, we investigated the suitability of performing partial least square regression (PLSR) on genotype-phenotype datasets to identify marker-trait associations. We utilized data collected on a cotton (Gossypium hirsutum L.) recombinant inbred line (RIL) mapping population that was evaluated under contrasting irrigation treatments, well-watered and water-limited conditions, in a hot, arid environment in 2012. Two phenotypic data sets were used in combination with the genetic data which consisted of 841 marker loci assigned to 117 linkage groups. The first dataset contained canopy traits that were gathered using a mobile, high-throughput phenotyping platform and included canopy temperature (CT), normalized difference vegetation index (NDVI), and canopy height (CHT) with leaf area index (LAI) being derived from NDVI and CHT measurements. The second phenotypic data set consisted of 14 elemental concentration measurements corresponding to the following elements: P, K, Ca, Mn, Fe, Zn, Ni, Cu, As, Co, Rb, Mo, S, and Mg. To conduct the PSLR analyses we used the "pls" and "pls depot" available in R statistical software version 3.2.4. The PLSR bi plot from the analysis of the first dataset showed that three (LAI, NDVI, and CHT) out of the four canopy traits were highly correlated, and by using multivariate analysis of variance (MANOVA), we detected 22 significant (p<0.01) marker-trait associations for the four traits. In contrast to the canopy trait analysis, our PLSR bi plot for the second dataset showed varying correlations for each of the 14 traits. Because of the lack of distinct trait similarities, MANOVA was not an ideal option to test for marker-trait associations so we implemented a jackknife re sampling technique. Jackknife re sampling failed to detect significant marker effects for several of the 14 elemental concentration traits. Thus, our future work aims to test other re sampling techniques such as boot straping for traits that do not exhibit high correlation. Overall, PLSR was a very informative way to comprehend data structure, displaying correlations within markers, within traits, and between marker and traits in one bi plot. Further studies are still needed to leverage detection of additional variance in correlated datasets and to prevent spurious results. To the best of our knowledge, this is the first time PLSR has been reported in such a context.

Genetic diversity, population structure and marker-trait associations for agronomic and grain traits in wild diploid wheat Triticum urartu

BACKGROUND

Wild diploid wheat, Triticum urartu (T. urartu) is the progenitor of bread wheat, and understanding its genetic diversity and genome function will provide considerable reference for dissecting genomic information of common wheat.

RESULTS

In this study, we investigated the morphological and genetic diversity and population structure of 238 T. urartu accessions collected from different geographic regions. This collection had 19.37 alleles per SSR locus and its polymorphic information content (PIC) value was 0.76, and the PIC and Nei's gene diversity (GD) of high-molecular-weight glutenin subunits (HMW-GSs) were 0.86 and 0.88, respectively. UPGMA clustering analysis indicated that the 238 T. urartu accessions could be classified into two subpopulations, of which Cluster I contained accessions from Eastern Mediterranean coast and those from Mesopotamia and Transcaucasia belonged to Cluster II. The wide range of genetic diversity along with the manageable number of accessions makes it one of the best collections for mining valuable genes based on marker-trait association. Significant associations were observed between simple sequence repeats (SSR) or HMW-GSs and six morphological traits: heading date (HD), plant height (PH), spike length (SPL), spikelet number per spike (SPLN), tiller angle (TA) and grain length (GL).

CONCLUSIONS

Our data demonstrated that SSRs and HMW-GSs were useful markers for identification of beneficial genes controlling important traits in T. urartu, and subsequently for their conservation and future utilization, which may be useful for genetic improvement of the cultivated hexaploid wheat.

Molecular polymorphism related to flowering trait variation in a Phaseolus vulgaris L. collection

The aim of this study was to investigate the flowering variation and the molecular polymorphism in key regulatory genes that control flowering in a Phaseolus vulgaris L. collection of 94 accessions from Europe and the Americas. The analysis of variance revealed that the difference in days-to-flowering between accessions was significant, with European accessions characterized by flowering precocity. Population structure analysis corroborated previous data on the genetic distinction between the Andean and Mesoamerican gene pools. A low level of admixture was detected. Genomic sequences of 15 gene fragments were obtained. About 7.0 kb per accession were sequenced and a total of 48 nucleotide substitutions identified. A Mixed Linear Model analysis, including population structure and kinship, was used to identify marker-trait associations. Haplotype tagging single nucleotide polymorphisms (htSNPs) associated with the studied traits were detected: in PvVRN1 and PvPHYB with days-to-flowering, in PvMYB29 with number of flower buds per inflorescence and in PvTFL1z and PvFCA with inflorescence length. The two genes associated with days-to-flowering control belong to the photoperiod and vernalization pathways. In particular, the PvVRN1 gene appears to play an important role in regulating the adaptation process of common bean.