Molecular markers from the transcribed/expressed region of the genome in higher plants

Genetics, genomics and breeding of fennel

Fennel (Foeniculum vulgare Mill. or Anethum foeniculum) stands out as a versatile herb whose cultivation spans across various regions worldwide, thanks to its adaptability to diverse climatic conditions. Its economic importance is mainly due to its numerous pharmaceutical properties and its widespread use in culinary applications. In this review, we first reviewed the chemical composition of this species, stressing the importance of two volatile compounds: t-anethole and estragole. The few cytological and genetic information available in the scientific literature were summarized. Regarding this latter aspect, we pointed out the almost complete absence of classical genetic studies, the lack of a chromosome-level reference genome, and the shortage of adequate transcriptomic studies. We also reviewed the main agronomic practices, with particular emphasis on breeding schemes aimed at the production of F1 hybrids and synthetic varieties. The few available studies on biotic and abiotic stresses were discussed too. Subsequently, we summarized the main studies on genetic diversity conducted in fennel and the available germplasm collections. Finally, we outlined an overview of the main in vitro regeneration techniques successfully applied in this species.

Evaluation and characteristic analysis of SSRs from the transcriptomic sequences of Perilla crop (Perilla frutescens L.)

Perilla crop is a self-fertilizing annual plant, cultivated and used mainly in East Asia. Perilla frutescens var. frutescens seeds are rich in unsaturated fatty acids, which have health benefits, and Perilla frutescens var. crispa leaves are rich in anthocyanins. However, genomic analysis such as whole genome sequencing or genetic mapping has not been performed on Perilla crop. This current study confirms the abundance and diversity of 15,991 simple sequence repeats (SSRs) classified in previous studies in the Perilla genome, selects and designs 1,538 SSR primer sets, and confirms which SSR primer sets exhibit high polymorphism. Of the 15,991 SSRs classified, there were 9,910 (62%) dinucleotide repeats, 5,652 (35.3%) trinucleotide repeats, and 429 (2.7%) tetranucleotide repeats. Among these, the most identified was (CT)n with a total of 4,817. The 15,991 SSRs had 4 to 26 repeats. Four repeats were the most frequent with 11,084 (69.3%). A total of 1,538 SSR primers were selected and designed to confirm polymorphism, of which 157 showed persistent and clear polymorphism. Among these 157 SSR primer sets, 98 (62.4%) were dinucleotide repeats, 39 (24.8%) were trinucleotide repeats, and 20 (12.7%) were tetranucleotide repeats. Among 549 SSR primers that showed polymorphism, trinucleotide repeats showed persistent polymorphism at a high rate. Therefore, when developing SSR primer sets for Perilla crop in the future, it is recommended that trinucleotide repeats be selected first. These research results will be helpful in future genomic analysis and development of SSR primers in Perilla crop.

Population Genetics of Haliotis discus hannai in China Inferred Through EST-SSR Markers

BACKGROUND/OBJECTIVES

The Pacific abalone Haliotis discus hannai originated in cold waters and is an economically important aquaculture shellfish in China. Our goal was to clarify the current status of the genetic structure of Pacific abalone in China.

METHODS

In this study, eighteen polymorphic EST-SSR loci were successfully developed based on the hemolymph transcriptome data of Pacific abalone, and thirteen highly polymorphic EST-SSR loci were selected for the genetic variation analysis of the six populations collected.

RESULTS

The results showed that the average number of observed alleles was 8.0769 (RC)-11.3848 (DQ) in each population. The number of observed alleles in the DQ, NH, and TJ populations was significantly higher than that in the RC population. The cultivated population outside the Changshan Islands has experienced a 22.79% reduction in allele diversity compared to the wild population of DQ. The pairwise Fst values and analysis of molecular variance (AMOVA) revealed significant population differentiation among all populations except DQ and NH populations, with RC and ZZ cultured populations exhibiting the largest population differentiation (Fst = 0.1334). The phylogenetic tree and structural analysis divided the six populations into two groups (group 1: NH, DQ, and ZZ; group 2: DL, TJ, and RC), and there was no relationship between geographical distance and genetic distance.

CONCLUSIONS

These results may reflect the large-scale culture from different populations in China and the exchange of juveniles between hatcheries. Different breeding conditions have led to a higher degree of genetic differentiation between the RC and ZZ populations. This study enables a better understanding of the genetic diversity and structure of current Pacific abalone populations.

Metabolomic Profiling and DNA-Fingerprinting of Newly Recorded White-Flowered Populations of Salvia lanigera Poir. in Egypt

Salvia lanigera Poir. is a small herbaceous perennial species with violet flowers that grows in low-altitude deserts, and sandy loam. During the collection of S. lanigera, unusual populations with white flowers were found. Therefore, the two populations (violet- and white-flowered) were subjected to comparative investigations, including DNA fingerprinting, chemical composition, and biological evaluation. The two populations showed DNA variations, with 6.66 % polymorphism in ISSR and 25 % in SCoT markers. GC/MS and UHPLC/HRMS of aqueous methanol extracts, led to the tentative identification of 43 and 50 compounds in both populations. In addition, the structures of nine compounds, including four first-time reported compounds in the species, were confirmed by NMR. Furthermore, the total extracts exhibited weak radical scavenging activity against DPPH and a lower inhibitory effect towards acetylcholinesterase. In conclusion, the obtained data suggested that the white-colored flower could be an additional important character record for the Egyptian S. lanigera.

Transcriptome analysis and development of EST-SSR markers in the mushroom Auricularia heimuer

Auricularia heimuer, the third most frequently cultivated edible mushroom species worldwide, has high medicinal value. However, a shortage of molecular marker hinders the efficiency and accuracy of genetic breeding efforts for A. heimuer. High-throughput transcriptome sequencing data are essential for gene discovery and molecular markers development. This study aimed to clarify the distribution of SSR loci across the A. heimuer transcriptome and to develop highly informative EST-SSR markers. These tools can be used for phylogenetic analysis, functional gene mining, and molecular marker-assisted breeding of A. heimuer. This study used Illumina high-throughput sequencing technology to obtain A. heimuer transcriptome data. The results revealed 37,538 unigenes in the A. heimuer transcriptome. Of these unigenes, 24,777 (66.01%) were annotated via comparison with the COG, Pfam, and NR databases. Overall, 2510 SSRs were identified from the unigenes, including 6 types of SSRs. The most abundant type of repeats were trinucleotides (1425, 56.77%), followed by mononucleotides (391, 15.58%) and dinucleotides (456, 18.17%). Primer pairs for 102 SSR loci were randomly designed for validity confirmation and polymorphism identification; this process yielded 53 polymorphic EST-SSR markers. Finally, 13 pairs of highly polymorphic EST-SSR primers were used to analyze the genetic diversity and population structure of 52 wild A. heimuer germplasms, revealing that the 52 germplasms could be divided into three categories. These results indicated that SSR loci were abundant in types, numbers, and frequencies, providing a potential basis for germplasm resource identification, genetic diversity analysis, and molecular marker-assisted breeding of A. heimuer.

Transcriptomics reveal useful resources for examining fruit development and variation in fruit size in Coccinia grandis

Introduction

The Cucurbitaceae family comprises many agronomically important members, that bear nutritious fruits and vegetables of great economic importance. Coccinia grandis, commonly known as Ivy gourd, belongs to this family and is widely consumed as a vegetable. Members of this family are known to display an impressive range of variation in fruit morphology. Although there have been studies on flower development in Ivy gourd, fruit development remains unexplored in this crop.

Methods

In this study, comparative transcriptomics of two Ivy gourd cultivars namely "Arka Neelachal Kunkhi" (larger fruit size) and "Arka Neelachal Sabuja" (smaller fruit size) differing in their average fruit size was performed. A de novo transcriptome assembly for Ivy gourd was developed by collecting fruits at different stages of development (5, 10, 15, and 20 days after anthesis i.e. DAA) from these two varieties. The transcriptome was analyzed to identify differentially expressed genes, transcription factors, and molecular markers.

Results

The transcriptome of Ivy gourd consisted of 155205 unigenes having an average contig size of 1472bp. Unigenes were annotated on publicly available databases to categorize them into different biological functions. Out of these, 7635 unigenes were classified into 38 transcription factor (TF) families, of which Trihelix TFs were most abundant. A total of 11,165 unigenes were found to be differentially expressed in both the varieties and the in silico expression results were validated through real-time PCR. Also, 98768 simple sequence repeats (SSRs) were identified in the transcriptome of Ivy gourd.

Discussion

This study has identified a number of genes, including transcription factors, that could play a crucial role in the determination of fruit shape and size in Ivy gourd. The presence of polymorphic SSRs indicated a possibility for marker-assisted selection for crop breeding in Ivy gourd. The information obtained can help select candidate genes that may be implicated in regulating fruit development and size in other fruit crops.

Comparative analysis of VMT genes/proteins in selected plant species with emphasis on bread wheat (Triticum aestivum L.)

BACKGROUND

In recent years, the study of molecular basis of uptake, transport and utilization of grain Fe/Zn (GFe/GZn) in wheat has been an active area of research. As a result, it has been shown that a number of transporters are involved in uptake and transport of Fe. In a recent study, knockout of a transporter gene OsVMT (VACUOLAR MUGINEIC ACID TRANSPORTER) in rice was shown to be involved in Fe homoeostasis.

OBJECTIVE

In this study, we analysed VMT genes among six monocots and three dicots with major emphasis on wheat VMT genes (TaVMTs), taking OsVMT gene as a reference.

METHODS AND RESULTS

Using OsVMT gene as a reference, VMT genes were identified and sequence similarities were examined among six monocots and three dicots. Each VMT protein carried one functional domain and 7 to 10 distinct motifs (including 9 novel motifs). The qRT-PCR analysis showed differential expression by all the six TaVMT genes in pairs of contrasting wheat genotypes with high (FAR4 and WB02) and low (K8027 and HD3226) GFe/GZn at two different grain filling stages (14 DAA and 28 DAA). TaVMT1 genes showed up-regulation in high GFe/Zn genotypes relative to low GFe/Zn genotypes, whereas the TaVMT2 genes showed down-regulation or nonsignificant up-regulation in a few cases.

CONCLUSIONS

At 14 DAA, each of the six TaVMT genes exhibited higher expression in wheat genotypes with high GFe and GZn relative to those with low GFe and GZn, suggesting major role of VMT genes in improvement of grain Fe/Zn homoeostasis, thus making TaVMT genes useful for improvement in Fe/Zn in wheat grains.

Effect of foliar application of phosphorus nanoparticles on the performance and sustainable agriculture of sweet corn

Traditional phosphorus fertilizers are necessary for plant growth but about 80-90% are lost into the surrounding environment via irrigation, therefore nano-fertilizers have been developed as slow-release fertilizers to achieve sustainable agriculture. This trial investigated the impact of the foliar application of hydroxyapatite nanoparticles (HA-NPs) as a source of nano-phosphorus (P-NPs) on two cultivars of sweet corn (yellow and white) throughout two seasons. The morphology and structure of the prepared HA-NPs were characterized via transmission electron microscopy (TEM) and X-ray diffractometry (XRD). In addition, agro-morphological criteria, chemical contents (i.e., photosynthetic pigments, phenols, indoles, minerals, etc.), and genomic template stability percentage (GTS%) were evaluated in the produced sweet corn. The application of 50 mg/l HA-NPs improved the growth characteristics, yield per hectare, leaf pigments, and chemical content of yellow sweet corn, whereas the application of 100 mg/l of HA-NPs to white sweet corn enhanced the vegetative characteristics, production, photosynthetic pigments, phenols, and indoles. The difference in results may be due to the presence of a +ve unique band with SCoT-4 and SCot-2 primers at 1250 and 470 bp in yellow and white corn treated with 50 and 100 mg/l, respectively. The minimum GTS% was recorded at a concentration of 75 mg/l for both white and yellow corn. The HA-NPs can be applied as a foliar source of P-NPs to achieve agricultural sustainability.

Deciphering the genomic regions governing major agronomic traits and blast resistance using genome wide association mapping in finger millet

Finger millet (Eleusine coracana L.) is climate resilient minor millet of Asia and Africa with wide adaptation and unparallel nutritional profile. To date, genomic resources available in finger millet are scanty and genetic control of agronomic traits remains elusive. Here, a collection of eco-geographically diverse 186 genotypes was quantified for variation in 13 agronomic traits and reaction to blast to identify marker-trait associations (MTAs) using genotyping-by-sequencing (GBS) and genome-wide association study (GWAS). GBS generated 2977 high quality single nucleotide polymorphism (SNPs) markers and identified three subpopulations with varying admixture levels. General linear and mixed model approaches of GWAS to correct for population structure and genetic relatedness identified 132 common MTAs for agronomic traits across the years. The phenotypic variance explained by the makers varied from 4.8% (TP692389-flag leaf width) to 20% (TP714446-green fodder weight). Of these, 26 MTAs showed homology with candidate genes having role in plant growth, development and photosynthesis in the genomes of foxtail millet, rice, maize, wheat and barley. We also found 4 common MTAs for neck blast resistance, which explained 5.9-15.1% phenotypic variance. Three MTAs for neck blast resistance showed orthologues in related genera having putative functions in pathogen defense in plants. The results of this work lay a foundation for understanding the genetic architecture of agronomic traits and blast resistance in finger millet and provide a framework for genomics assisted breeding.

Phylogenomic analysis of 20S proteasome gene family reveals stress-responsive patterns in rapeseed (Brassica napus L.)

The core particle represents the catalytic portions of the 26S proteasomal complex. The genes encoding α- and β-subunits play a crucial role in protecting plants against various environmental stresses by controlling the quality of newly produced proteins. The 20S proteasome gene family has already been reported in model plants such as Arabidopsis and rice; however, they have not been studied in oilseed crops such as rapeseed (Brassica napus L.). In the present study, we identified 20S proteasome genes for α- (PA) and β-subunits (PB) in B. napus through systematically performed gene structure analysis, chromosomal location, conserved motif, phylogenetic relationship, and expression patterns. A total of 82 genes, comprising 35 BnPA and 47 BnPB of the 20S proteasome, were revealed in the B. napus genome. These genes were distributed on all 20 chromosomes of B. napus and most of these genes were duplicated on homoeologous chromosomes. The BnPA (α1-7) and BnPB (β1-7) genes were phylogenetically placed into seven clades. The pattern of expression of all the BnPA and BnPB genes was also studied using RNA-seq datasets under biotic and abiotic stress conditions. Out of 82 BnPA/PB genes, three exhibited high expression under abiotic stresses, whereas two genes were overexpressed in response to biotic stresses at both the seedling and flowering stages. Moreover, an additional eighteen genes were expressed under normal conditions. Overall, the current findings developed our understanding of the organization of the 20S proteasome genes in B. napus, and provided specific BnPA/PB genes for further functional research in response to abiotic and biotic stresses.

Seed coat colour of Indian mustard [Brassica juncea (L.) Czern. and Coss.] is associated with Bju.TT8 homologs identifiable by targeted functional markers

Seed coat colour is an important trait in Indian mustard. Breeding for seed coat colour needs precise knowledge of mode of inheritance and markers linked to it. The present study was focussed on genetics and development of functional markers for seed coat colour. F₁s (direct and reciprocal) and F₂ populations were developed by crossing two contrasting parents for seed coat colour (DRMRIJ-31, brown seeded and RLC-3, yellow seeded). Phenotypic results have shown that the seed coat colour trait was under the influence of maternal effect and controlled by digenic-duplicate gene action. Further, Bju.TT8 homologs of both parents (DRMRIJ-31 and RLC-3) were cloned and sequenced. Sequencing results of Bju.TT8 homologs revealed that in RLC-3, gene Bju.ATT8 had an insertion of 1279bp in the 7^th exon; whereas, gene Bju.BTT8 had an SNP (C→T) in the 7^th exon. These two mutations were found to be associated with yellow seed coat colour. Using sequence information, functional markers were developed for both Bju.TT8 homologs, validated on F₂ population and were found highly reliable with no recombination between the markers and the phenotype. Further, these markers were subjected to a germplasm assembly of Indian mustard, and their allelic combination for the seed coat colour genes has been elucidated. The comparative genomics of TT8 genes revealed high degree of similarity between and across the Brassica species, and the respective diploid progenitors in tetraploid Brassica species are the possible donors of TT8 homologs. This study will help in the marker-assisted breeding for seed coat colour, and aid in understanding seed coat colour genetics more precisely.

Measurement of Genetic Mobility Using a Transposon-Based Marker System in Sorghum

Transposable elements (TEs) are ubiquitous repetitive components of eukaryotic organisms that show mobility in the genome against diverse stresses. TEs contribute considerably to the size, structure, and plasticity of genomes and also play an active role in genome evolution by helping their hosts adapt to novel conditions by conferring useful characteristics. We developed a simple and rapid method for investigation of genetic mobility and diversity among TEs in combination with a target region amplification polymorphism (TE-TRAP) marker system in gamma-irradiated sorghum mutants. The TE-TRAP marker system reveals a high level of genetic diversity, which provides a useful marker resource for genetic mobility research.

GPR-Based Automatic Identification of Root Zones of Influence Using HDBSCAN

The belowground root zone of influence (ZOI) is fundamental to the study of the root–root and root–soil interaction mechanisms of plants and is vital for understanding changes in plant community compositions and ecosystem processes. However, traditional root research methods have a limited capacity to measure the actual ZOIs within plant communities without destroying them in the process. This study has developed a new approach to determining the ZOIs within natural plant communities. First, ground-penetrating radar (GPR), a non-invasive near-surface geophysical tool, was used to obtain a dataset of the actual spatial distribution of the coarse root system in a shrub quadrat. Second, the root dataset was automatically clustered and analyzed using the hierarchical density-based spatial clustering of applications with noise (HDBSCAN) algorithm to determine the ZOIs of different plants. Finally, the shape, size, and other characteristics of each ZOI were extracted based on the clustering results. The proposed method was validated using GPR-obtained root data collected in two field shrub plots and one simulation on a dataset from existing literature. The results show that the shrubs within the studied community exhibited either segregated and aggregated ZOIs, and the two types of ZOIs were distinctly in terms of shape and size, demonstrating the complexity of root growth in response to changes in the surrounding environment. The ZOIs extracted based on GPR survey data were highly consistent with the actual growth pattern of shrub roots and can thus be used to reveal the spatial competition strategies of plant roots responding to changes in the soil environment and the influence of neighboring plants.

Development and use of miRNA-derived SSR markers for the study of genetic diversity, population structure, and characterization of genotypes for breeding heat tolerant wheat varieties

Heat stress is an important abiotic factor that limits wheat production globally, including south-east Asia. The importance of micro (mi) RNAs in gene expression under various biotic and abiotic stresses is well documented. Molecular markers, specifically simple sequence repeats (SSRs), play an important role in the wheat improvement breeding programs. Given the role of miRNAs in heat stress-induced transcriptional regulation and acclimatization, the development of miRNA-derived SSRs would prove useful in studying the allelic diversity at the heat-responsive miRNA-genes in wheat. In the present study, efforts have been made to identify SSRs from 96 wheat heat-responsive miRNA-genes and their characterization using a panel of wheat genotypes with contrasting reactions (tolerance/susceptible) to heat stress. A set of 13 miRNA-derived SSR markers were successfully developed as an outcome. These miRNA-SSRs are located on 11 different common wheat chromosomes (2A, 3A, 3B, 3D, 4D, 5A, 5B, 5D, 6A, 6D, and 7A). Among 13 miRNA-SSRs, seven were polymorphic on a set of 37 selected wheat genotypes. Within these polymorphic SSRs, three makers, namely HT-169j, HT-160a, and HT-160b, were found promising as they could discriminate heat-tolerant and heat-susceptible genotypes. This is the first report of miRNA-SSR development in wheat and their deployment in genetic diversity and population structure studies and characterization of trait-specific germplasm. The study suggests that this new class of molecular makers has great potential in the marker-assisted breeding (MAB) programs targeted at improving heat tolerance and other adaptability or developmental traits in wheat and other crops.

Comprehensive functional analysis and mapping of SSR markers in the chickpea genome (Cicer arietinum L.)

Plant molecular breeding largely depends on the relationship between molecular markers and major traits. Herein, a total of 32,962 genomic simple sequence repeats (SSRs) were detected in the whole genome of chickpea with an average density of 94.93 SSRs/Mb. Chickpea chromosomes uniformity test indicated that the genomic SSRs (gSSRs) were steadily distributed across the genome. Moreover, 48,667 transcriptome sequences were analyzed and 1949 SSR-containing transcript assembly contigs (TACs) were identified. The analysis showed that di- and trinucleotide SSRs were the most frequent SSR motifs within the transcriptome sequences. Among them, AT and TTA and AG and TTC motifs within the transcriptome showed the highest frequencies among di- and trinucleotide repeat motifs, respectively. The SSRs-containing TACs were compared to the GenBank non-redundant database using BLASTX, and subsequently, gene ontology (GO) analysis was performed using QuickGO browser to reduce complexity and highlight biological processes associated with the SSRs-containing TACs. The identified SSRs-containing TACs were categorized into 35 enriched functional-related gene group. The mapping of characterized SSRs-containing TACs onto chickpea chromosomes was performed using BLASTN. The mapping result showed that, a total of 1798 SSRs-containing TACs were mapped onto the chickpea genome. Based on the functional analysis result, 249 and 242 of the mapped SSRs-containing TACs were found in the genes encoding for putative stress-related proteins and transcription factors, respectively. The results presented here can be applied to improve and speed up the chickpea breeding programs.

Discovery of powdery mildew resistance gene candidates from Aegilops biuncialis chromosome 2Mb based on transcriptome sequencing

Powdery mildew is one of the most widespread diseases of wheat. The development and deployment of resistant varieties are one of the most economical and effective methods to manage this disease. Our previous study showed that the gene(s) at 2Mb in Chinese Spring (CS)-Aegilops biuncialis 2Mb disomic addition line TA7733 conferred a high level of resistance to powdery mildew of wheat. In this study, resistance spectrum of TA7733 was assayed by using 15 Blumeria graminis f. sp. tritici (Bgt) isolates prevalent in different regions of China. The result indicated that TA7733 was highly resistant to all tested Bgt isolates and the gene(s) on chromosome 2Mb conferred broad-spectrum resistance to powdery mildew. In order to characterize mechanism of powdery mildew resistance by identifying candidates R-genes derived from Ae. biuncialis chromosome 2Mb and develop 2Mb-specific molecular markers, we performed RNA-seq analysis on TA7733 and CS. In total we identified 7,278 unigenes that showed specific expression in TA7733 pre and post Bgt-infection when compared to CS. Of these 7,278 unigenes, 295 were annotated as putative resistance (R) genes. Comparatively analysis of R-gene sequences from TA7733 and CS and integration CS Ref Seq v1.0 were used to develop R-gene specific primers. Of 295 R-genes we identified 53 R-genes were specific to 2Mb and could be involved in powdery mildew resistance. Functional annotation of majority of the 53 R-genes encoded nucleotide binding leucine rich repeat (NLR) protein. The broad-spectrum resistance to powdery mildew in TA7733 and availability of 2Mb-derived putative candidate R-gene specific molecular markers identified in this study will lay foundations for transferring powdery mildew resistance from 2Mb to common wheat by inducing CS-Ae. biuncialis homoeologous recombination. Our study also provides useful candidates for further isolation and cloning of powdery mildew resistance gene(s) from Ae. biuncialis chromosome 2Mb.

Whole Genome Diversity, Population Structure, and Linkage Disequilibrium Analysis of Chickpea (Cicer arietinum L.) Genotypes Using Genome-Wide DArTseq-Based SNP Markers

Characterization of genetic diversity, population structure, and linkage disequilibrium is a prerequisite for proper management of breeding programs and conservation of genetic resources. In this study, 186 chickpea genotypes, including advanced "Kabuli" breeding lines and Iranian landrace "Desi" chickpea genotypes, were genotyped using DArTseq-Based single nucleotide polymorphism (SNP) markers. Out of 3339 SNPs, 1152 markers with known chromosomal position were selected for genome diversity analysis. The number of mapped SNP markers varied from 52 (LG8) to 378 (LG4), with an average of 144 SNPs per linkage group. The chromosome size that was covered by SNPs varied from 16,236.36 kbp (LG8) to 67,923.99 kbp (LG5), while LG4 showed a higher number of SNPs, with an average of 6.56 SNPs per Mbp. Polymorphism information content (PIC) value of SNP markers ranged from 0.05 to 0.50, with an average of 0.32, while the markers on LG4, LG6, and LG8 showed higher mean PIC value than average. Unweighted neighbor joining cluster analysis and Bayesian-based model population structure grouped chickpea genotypes into four distinct clusters. Principal component analysis (PCoA) and discriminant analysis of principal component (DAPC) results were consistent with that of the cluster and population structure analysis. Linkage disequilibrium (LD) was extensive and LD decay in chickpea germplasm was relatively low. A few markers showed r2 ≥ 0.8, while 2961 pairs of markers showed complete LD (r2 = 1), and a huge LD block was observed on LG4. High genetic diversity and low kinship value between pairs of genotypes suggest the presence of a high genetic diversity among the studied chickpea genotypes. This study also demonstrates the efficiency of DArTseq-based SNP genotyping for large-scale genome analysis in chickpea. The genotypic markers provided in this study are useful for various association mapping studies when combined with phenotypic data of different traits, such as seed yield, abiotic, and biotic stresses, and therefore can be efficiently used in breeding programs to improve chickpea.

Conserved DNA-derived polymorphism, new markers for genetic diversity analysis of Tunisian Pistacia vera L

Pistachio trees (Pistacia vera L.) have been cultivated in Tunisia for decades and the plantation was extended mostly in the center of the country contributing to the economic growth of marginalized areas. Herein we used conserved DNA derived polymorphism (CDDP) technique, which target specifically conserved sequences of plant functional genes, to assess the genetic diversity and construct genetic relationships among 65 Tunisian pistachio trees. A set of nine primers were used and 157 CDDP markers were revealed with an average of 17.44 showing a high degree of polymorphism (99.37%). The average of polymorphism information content of CDDP markers was of 0.86, which indicates the efficiency of CDDP primers in the estimation of genetic diversity between pistachios. UPGMA dendrogram and the principal component analysis showed four clusters of analyzed pistachios trees. Our results showed that the genetic structure depends on: (1) the gene exchanges between groups, (2) the geographical origin and (3) the sex of the tree. The same result was revealed by the Bayesian analysis implemented in STRUCTURE at K = 4, in which the pistachio genotypes of El Guettar, Kasserine and Sfax were assigned with more than 80% of probability. Our results prove polymorphism and efficiency of CDDP markers in the characterization and genetic diversity analysis of P. vera L. genotypes to define conservation strategy.

Development and Validation of a Gene-Targeted dCAPS Marker for Marker-Assisted Selection of Low-Alkaloid Content in Seeds of Narrow-Leafed Lupin (Lupinus angustifolius L.)

Low-alkaloid content is an important breeding target to improve the quality of lupin seeds. An APETALA2/ethylene response transcription factor, RAP2-7, is likely a candidate gene for the major alkaloid locus iucundus, and plays a crucial role in regulation of seed alkaloid content in narrow-leafed lupin (NLL; Lupinus angustifolius L.). Here, we exploited a single-nucleotide polymorphism within RAP2-7 credibly associated with seed alkaloid content, to develop the co-dominant derived cleaved amplified polymorphic sequence (dCAPS) marker iuc_RAP2-7. Marker validation in 202 NLL accessions demonstrated that seed alkaloid content ≥0.9% of the seed dry weight was associated with the high-alkaloid marker band (Iucundus genotypes), whereas alkaloid content up to 0.5% of the seed dry weight was associated with the low-alkaloid marker band (iucundus genotypes). Within a given detection limit, iuc_RAP2-7 unambiguously identified all but three low-alkaloid accessions. The latter accessions apparently have a different regulatory mechanism for seed alkaloid content because the RAP2-7 gene/putative promoter sequence and expression of alkaloid-associated genes in the leaves of the three ambiguous accessions were similar to those of bitter Iucundus lines. We consider the iuc_RAP2-7 marker is a powerful tool that will facilitate NLL marker-assisted selection by rapid rejection of bitter Iucundus genotypes and thus accelerate development of new low-alkaloid cultivars.

Genome wide association mapping of agro-morphological traits among a diverse collection of finger millet (Eleusine coracana L.) genotypes using SNP markers

Finger millet (Eleusine coracana L.) is an important dry-land cereal in Asia and Africa because of its ability to provide assured harvest under extreme dry conditions and excellent nutritional properties. However, the genetic improvement of the crop is lacking in the absence of suitable genomic resources for reliable genotype-phenotype associations. Keeping this in view, a diverse global finger millet germplasm collection of 113 accessions was evaluated for 14 agro-morphological characters in two environments viz. ICAR-Vivekananda Institute of Hill Agriculture, Almora (E1) and Crop Research Centre (CRC), GBPUA&T, Pantnagar (E2), India. Principal component analysis and cluster analysis of phenotypic data separated the Indian and exotic accessions into two separate groups. Previously generated SNPs through genotyping by sequencing (GBS) were used for association mapping to identify reliable marker(s) linked to grain yield and its component traits. The marker trait associations were determined using single locus single trait (SLST), multi-locus mixed model (MLMM) and multi-trait mixed model (MTMM) approaches. SLST led to the identification of 20 marker-trait associations (MTAs) (p value<0.01 and <0.001) for 5 traits. While advanced models, MLMM and MTMM resulted in additional 36 and 53 MTAs, respectively. Nine MTAs were common out of total 109 associations in all the three mapping approaches (SLST, MLMM and MTMM). Among these nine SNPs, five SNP sequences showed homology to candidate genes of Oryza sativa (Rice) and Setaria italica (Foxtail millet), which play an important role in flowering, maturity and grain yield. In addition, 67 and 14 epistatic interactions were identified for 10 and 7 traits at E1 and E2 locations, respectively. Hence, the 109 novel SNPs associated with important agro-morphological traits, reported for the first time in this study could be precisely utilized in finger millet genetic improvement after validation.

The leaf transcriptome of fennel (Foeniculum vulgare Mill.) enables characterization of the t-anethole pathway and the discovery of microsatellites and single-nucleotide variants

Fennel is a plant species of both agronomic and pharmaceutical interest that is characterized by a shortage of genetic and molecular data. Taking advantage of NGS technology, we sequenced and annotated the first fennel leaf transcriptome using material from four different lines and two different bioinformatic approaches: de novo and genome-guided transcriptome assembly. A reference transcriptome for assembly was produced by combining these two approaches. Among the 79,263 transcripts obtained, 47,775 were annotated using BLASTX analysis performed against the NR protein database subset with 11,853 transcripts representing putative full-length CDS. Bioinformatic analyses revealed 1,011 transcripts encoding transcription factors, mainly from the BHLH, MYB-related, C2H2, MYB, and ERF families, and 6,411 EST-SSR regions. Single-nucleotide variants of SNPs and indels were identified among the 8 samples at a frequency of 0.5 and 0.04 variants per Kb, respectively. Finally, the assembled transcripts were screened to identify genes related to the biosynthesis of t-anethole, a compound well-known for its nutraceutical and medical properties. For each of the 11 genes encoding structural enzymes in the t-anethole biosynthetic pathway, we identified at least one transcript showing a significant match. Overall, our work represents a treasure trove of information exploitable both for marker-assisted breeding and for in-depth studies on thousands of genes, including those involved in t-anethole biosynthesis.

Comparative whole genome re-sequencing analysis in upland New Rice for Africa: insights into the breeding history and respective genome compositions

BACKGROUND

Increasing rice demand is one of the consequences of the steadily improving socio-economic status of the African countries. New Rice for Africa (NERICA), which are interspecific hybrids between Asian and African rice varieties, are one of successful breeding products utilizing biodiversity across the two different rice crop species. Upland NERICA varieties (NU) exhibit agronomic traits of value for the harsh eco-geography, including shorter duration, higher yield and stress tolerance, compared to local African varieties. However, the molecular basis of the traits in NU varieties is largely unknown.

RESULTS

Whole genome re-sequencing was performed of four NU lines (3, 4, 5, and 7) and for the parental Oryza sativa WAB56-104 and Oryza glaberrima CG14. The k-mer analysis predicted large genomes for the four NU lines, most likely inherited from WAB56-104. Approximately 3.1, 0.10, and 0.40 million single nucleotide polymorphisms, multi nucleotide polymorphisms, and short insertions/deletions were mined between the parental lines, respectively. Integrated analysis with another four NU lines (1, 2, 8, and 9) showed that the ratios of the donor CG14 allelic sites in the NU lines ranged from 1.3 to 9.8%. High resolution graphical genotype indicated genome-level similarities and common genetic events during the breeding process: five xyloglucan fucosyltransferase from O. glaberrima were introgressed in common. Segregation of genic segments revealed potential causal genes for some agronomic traits including grain shattering, awnness, susceptibility to bacterial leaf bright, and salt tolerance. Analysis of unmapped sequences against the reference cultivar Nipponbare indicated existence of unique genes for pathogen and abiotic stress resistance in the NU varieties.

CONCLUSIONS

The results provide understanding of NU genomes for rice improvement for Africa reinforcing local capacity for food security and insights into molecular events in breeding of interspecific hybrid crops.

De novo assembly, gene annotation, and marker development using Illumina paired-end transcriptome sequencing in the Crassadoma gigantea

Crassadoma gigantea is an important commercial marine bivalve species in Baja California and Mexico. In this study, we have applied RNA-Seq technology to profile the transcriptome of the C. gigantea for the first time. A total of 80,832,518 raw reads were produced from a Illumina HiSeq4000 platform, and 77,306,198 (95.64%) clean reads were generated after trimming the adaptor sequences. The transcriptome assembled into 158,855 transcripts with an N50 size of 1995 bp and an average size of 1008 bp. A number of DNA repair related genes, such as MSH3, EGF, TGF, IGF, FGF, encoding different groups of growth factors were found in the transcriptome data of C. gigantean. In addition, immune related genes Toll-like receptor (TLR) including TLR1, TLR2, TLR3, TLR4, TLR5, TLR6, TLR7, TLR8, and TLR9 was also observed in C. gigantean. A set of 12 polymorphic microsatellite loci were firstly developed and characterized in C. gigantea. The results show that the number of alleles and expected heterozygosity ranged from 3 to 9 and from 0.254 to 0.820, respectively. The average polymorphic information content was 0.790. These microsatellite loci will facilitate future studies of population structure and conservation genetics in this species.

Simple Sequence Repeat (SSR) Genetic Linkage Map of D Genome Diploid Cotton Derived from an Interspecific Cross between Gossypium davidsonii and Gossypium klotzschianum

The challenge in tetraploid cotton cultivars is the narrow genetic base and therefore, the bottleneck is how to obtain interspecific hybrids and introduce the germplasm directly from wild cotton to elite cultivars. Construction of genetic maps has provided insight into understanding the genome structure, interrelationships between organisms in relation to evolution, and discovery of genes that carry important agronomic traits in plants. In this study, we generated an interspecific hybrid between two wild diploid cottons, Gossypium davidsonii and Gossypium klotzschianum, and genotyped 188 F2:3 populations in order to develop a genetic map. We screened 12,560 SWU Simple Sequence Repeat (SSR) primers and obtained 1000 polymorphic markers which accounted for only 8%. A total of 928 polymorphic primers were successfully scored and only 728 were effectively linked across the 13 chromosomes, but with an asymmetrical distribution. The map length was 1480.23 cM, with an average length of 2.182 cM between adjacent markers. A high percentage of the markers on the map developed, and for the physical map of G. raimondii, exhibited highly significant collinearity, with two types of duplication. High level of segregation distortion was observed. A total of 27 key genes were identified with diverse roles in plant hormone signaling, development, and defense reactions. The achievement of developing the F2:3 population and its genetic map constructions may be a landmark in establishing a new tool for the genetic improvement of cultivars from wild plants in cotton. Our map had an increased recombination length compared to other maps developed from other D genome cotton species.

Overexpression of the Prunus sogdiana NBS-LRR Subgroup Gene PsoRPM2 Promotes Resistance to the Root-Knot Nematode Meloidogyne incognita in Tobacco

Root-knot nematodes (RKNs), particularly Meloidogyne incognita, are the most devastating soil-borne pathogens that significantly affect the production of Prunus spp. fruit. RKN infection is difficult to control and consequently causes massive yield losses each year. However, several germplasms of wild Prunus spp. have been shown to display resistance to M. incognita. Consequently, both the isolation of novel plant resistance (R) genes and the characterization of their resistance mechanisms are important strategies for future disease control. R proteins require the co-chaperone protein HSP90-SGT1-RAR1 to achieve correct folding, maturation, and stabilization. Here, we used homologous cloning to isolate the R gene PsoRPM2 from the RKN-resistant species Prunus sogdiana. PsoRPM2 was found to encode a TIR-NB-LRR-type protein and react with significantly elevated PsoRPM2 expression levels in response to RKN infection. Transient expression assays indicated PsoRPM2 to be located in both the cytoplasm and the nucleus. Four transgenic tobacco lines that heterologously expressed PsoRPM2 showed enhanced resistance to M. incognita. Yeast two-hybrid analysis and bimolecular fluorescence complementation analysis demonstrated that both PsoRAR1 and PsoRPM2 interacted with PsoHSP90-1 and PsoSGT1, but not with one another. These results indicate that the observed PsoRPM2-mediated RKN resistance requires both PsoHSP90-1 and PsoSGT1, further suggesting that PsoRAR1 plays a functionally redundant role in the HSP90-SGT1-RAR1 co-chaperone.

Construction of a genetic map using EST-SSR markers and QTL analysis of major agronomic characters in hexaploid sweet potato (Ipomoea batatas (L.) Lam)

The Sweet potato, Ipomoea batatas (L.) Lam, is difficult to study in genetics and genomics because it is a hexaploid. The sweet potato study not have been performed domestically or internationally. In this study was performed to construct genetic map and quantitative trait loci (QTL) analysis. A total of 245 EST-SSR markers were developed, and the map was constructed by using 210 of those markers. The total map length was 1508.1 cM, and the mean distance between markers was 7.2 cM. Fifteen characteristics were investigated for QTLs analysis. According to those, the Four QTLs were identified, and The LOD score was 3.0. Further studies need to develop molecular markers in terms of EST-SSR markers for doing to be capable of efficient breeding. The genetic map created here using EST-SSR markers will facilitate planned breeding of sweet potato cultivars with various desirable traits.

Genome size, cytogenetic data and transferability of EST-SSRs markers in wild and cultivated species of the genus Theobroma L. (Byttnerioideae, Malvaceae)

The genus Theobroma comprises several trees species native to the Amazon. Theobroma cacao L. plays a key economic role mainly in the chocolate industry. Both cultivated and wild forms are described within the genus. Variations in genome size and chromosome number have been used for prediction purposes including the frequency of interspecific hybridization or inference about evolutionary relationships. In this study, the nuclear DNA content, karyotype and genetic diversity using functional microsatellites (EST-SSR) of seven Theobroma species were characterized. The nuclear content of DNA for all analyzed Theobroma species was 1C = ~ 0.46 pg. These species presented 2n = 20 with small chromosomes and only one pair of terminal heterochromatic bands positively stained (CMA⁺/DAPI⁻ bands). The small size of Theobroma ssp. genomes was equivalent to other Byttnerioideae species, suggesting that the basal lineage of Malvaceae have smaller genomes and that there was an expansion of 2C values in the more specialized family clades. A set of 20 EST-SSR primers were characterized for related species of Theobroma, in which 12 loci were polymorphic. The polymorphism information content (PIC) ranged from 0.23 to 0.65, indicating a high level of information per locus. Combined results of flow cytometry, cytogenetic data and EST-SSRs markers will contribute to better describe the species and infer about the evolutionary relationships among Theobroma species. In addition, the importance of a core collection for conservation purposes is highlighted.

Allele Mining and Selective Patterns of Pi9 Gene in a Set of Rice Landraces from India

Allelic variants of the broad-spectrum blast resistance gene, Pi9 (nucleotide binding site-leucine-rich repeat region) have been analyzed in Indian rice landraces. They were selected from the list of 338 rice landraces phenotyped in the rice blast nursery at central Rainfed Upland Rice Research Station, Hazaribag. Six of them were further selected on the basis of their resistance and susceptible pattern for virulence analysis and selective pattern study of Pi9 gene. The sequence analysis and phylogenetic study illustrated that such sequences are vastly homologous and clustered into two groups. All the blast resistance Pi9 alleles were grouped into one cluster, whereas Pi9 alleles of susceptible landraces formed another cluster even though these landraces have a low level of DNA polymorphisms. A total number of 136 polymorphic sites comprising of transitions, transversions, and insertion and deletions (InDels) were identified in the 2.9 kb sequence of Pi9 alleles. Lower variation in the form of mutations (77) (Transition + Transversion), and InDels (59) were observed in the Pi9 alleles isolated from rice landraces studied. The results showed that the Pi9 alleles of the selected rice landraces were less variable, suggesting that the rice landraces would have been exposed to less number of pathotypes across the country. The positive Tajima's D (0.33580), P > 0.10 (not significant) was observed among the seven rice landraces, which suggests the balancing selection of Pi9 alleles. The value of synonymous substitution (-0.43337) was less than the non-synonymous substitution (0.78808). The greater non-synonymous substitution than the synonymous means that the coding region, mainly the leucine-rich repeat domain was under diversified selection. In this study, the Pi9 gene has been subjected to balancing selection with low nucleotide diversity which is different from the earlier reports, this may be because of the closeness of the rice landraces, cultivated in the same region, and under low pathotype pressure.

Characterization and development of EST-SSR markers in sweet potato (Ipomoea batatas (L.) Lam)

In this study, a cDNA library was constructed from the total RNA of sweet potato leaves. A total of 789 copies of the cDNA were cloned in Escherichia coli by employing the pGEM-T Easy vector. Sequencing was carried out by Solgent Co. (Korea). As many as 579 expressed sequence tag-simple sequence repeat (EST-SSR) markers were designed (73.38%) from the known cDNA nucleotide base sequences. The lengths of the developed EST-SSR markers ranged from 100 to 499 bp (average length 238 bp). Their motif sequence types were varied, with most being dinucleotides and pentanucleotides, and the most commonly found motifs were CAGAAT (29.0%) and TCT (2.8%). Based on these SSR-containing sequences, 619 pairs of high-quality SSR primers were designed using WebSat and Primer3web. The total number of primers designed was 144. Polymorphism was evident in 82 EST-SSR markers among 20 Korean sweet potato cultivars tested and in 90 EST-SSR markers in the two parents of a mapping population, Yeseumi and Annobeny. In this study, the hexaploid sweet potato (2n = 6x = 90) EST-SSR markers were developed in the absence of full-sequence data. Moreover, by acting as a molecular tag for particular traits, the EST-SSR marker can also simultaneously identify information about the corresponding gene. These EST-SSR markers will allow the molecular analysis of sweet potato to be done more efficiently. Thus, we can develop high-quality sweet potato while overcoming the challenges from climate change and other unfavorable conditions.

Genetic relationship and diversity among coconut (Cocos nucifera L.) accessions revealed through SCoT analysis

Coconut (Cocos nucifera L.) is one of the important palms grown both as a homestead and plantation crop in countries and most island territories of tropical regions. Different DNA-based marker systems have been utilized to assess the extent of genetic diversity in coconut. Advances in genomics research have resulted in the development of novel gene-targeted markers. In the present study, we have used a simple and novel marker system, start codon targeted polymorphism (SCoT), for its evaluation as a potential marker system in coconut. SCoT markers were utilized for assessment of genetic diversity in 23 coconut accessions (10 talls and 13 dwarfs), representing different geographical regions. Out of 25 SCoT primers screened, 15 primers were selected for this study based on their consistent amplification patterns. A total of 102 scorable bands were produced by the 15 primers, 88 % of which were polymorphic. The scored data were used to construct a similarity matrix. The similarity coefficient values ranged between 0.37 and 0.91. These coefficients were utilized to construct a dendrogram using the unweighted pair group of arithmetic means (UPGMA). The extent of genetic diversity observed based on SCoT analysis of coconut accessions was comparable to earlier findings using other marker systems. Tall and dwarf coconut accessions were clearly demarcated, and in general, coconut accessions from the same geographical region clustered together. The results indicate the potential of SCoT markers to be utilized as molecular markers to detect DNA polymorphism in coconut accessions.

De novo assembly and characterization of foot transcriptome and microsatellite marker development for Paphia textile

Paphia textile is an important, aquaculture bivalve clam species distributed mainly in China, Philippines, and Malaysia. Recent studies of P. textile have focused mainly on artificial breeding and nutrition analysis, and the transcriptome and genome of P. textile have rarely been reported. In this work, the transcriptome of P. textile foot tissue was sequenced on an Illumina HiSeq™ 2000 platform. A total of 20,219,795 reads were generated, resulting in 4.08 Gb of raw data. The raw reads were cleaned and assembled into 54,852 unigenes with an N50 length of 829 bp. Of these unigenes, 38.92% were successfully annotated based on their matches to sequences in seven public databases. Among the annotated unigenes, 14,571 were assigned Gene Ontology terms, 5448 were classified to Clusters of Orthologous Groups categories, and 6738 were mapped to 228 pathways in the Kyoto Encyclopedia of Genes and Genomes database. For functional marker development, 5605 candidate simple sequence repeats were identified in the transcriptome and 80 primer pairs were selected randomly and amplified in a wild population of P. textile. A total of 36 loci that exhibited obvious repeat length polymorphisms were detected. The transcriptomic data and microsatellite markers will provide valuable resources for future functional gene analyses, genetic map construction, and quantitative trait loci mapping in P. textile.

Genetic differentiation and hybrid identification using microsatellite markers in closely related wild species

Identifying the genetic basis of speciation is critical for understanding the evolutionary history of closely related wild species. Recently diverged species facilitate the study of speciation because many genetic and morphological characteristics are still shared by the organisms under study. The Petunia genus grows in South American grasslands and comprises both recently diverged wild species and commercial species. In this work, we analysed two closely related species: Petunia exserta, which has a narrow endemic range and grows exclusively in rocky shelters, and Petunia axillaris, which is widely distributed and comprises three allopatric subspecies. Petunia axillaris ssp. axillaris and P. exserta occur in sympatry, and putative hybrids between them have been identified. Here, we analysed 14 expressed sequence tag-simple sequence repeats (EST-SSRs) in 126 wild individuals and 13 putative morphological hybrids with the goals of identifying differentially encoded alleles to characterize their natural genetic diversity, establishing a genetic profile for each taxon and to verify the presence of hybridization signal. Overall, 143 alleles were identified and all taxa contained private alleles. Four major groups were identified in clustering analyses, which indicated that there are genetic distinctions among the groups. The markers evaluated here will be useful in evolutionary studies involving these species and may help categorize individuals by species, thus enabling the identification of hybrids between both their putative taxa. The individuals with intermediate morphology presented private alleles of their both putative parental species, although they showed a level of genetic mixing that was comparable with some of the individuals with typical P. exserta morphology. The EST-SSR markers scattered throughout the Petunia genome are very efficient tools for characterizing the genetic diversity in wild taxa of this genus and aid in identifying interspecific hybrids based on the presence of private alleles. These properties indicate that these markers will be helpful tools in evolutionary studies.

Functional marker development of miR1511-InDel and allelic diversity within the genus Glycine

BACKGROUND

Single-stranded non-protein coding small RNAs, 18-25 nucleotides in length, are ubiquitous throughout plants genomes and are involved in post-transcriptional gene regulation. Several types of DNA markers have been reported for the detection of genetic diversity or sequence variation in soybean, one of the most important legume crops in worldwide for seed protein and oil content. Recently, with the available of public genomic databases, there has been a shift from the labor-intensive development of PCR-based markers to sequence-based genotyping and the development of functional markers within genes, often coupled with the use of RNA information. But thus far miRNA-based markers have been only developed in rice and tobacco. Here we report the first functional molecular miRNA marker, miR1511-InDel, in soybean for a specific single copy locus used to assess genetic variation in domesticated soybean (Glycine max [L.] Merr) and its wild progenitor (Glycine soja Sieb. & Zucc.).

RESULTS

We genotyped a total of 1,669 accessions of domesticated soybean (G. max) and its wild progenitor G. soja which are native throughout the China and parts of Korea, Japan and Russia. The results indicate that the miR1511 locus is distributed in cultivated soybean and has three alleles in annual wild soybean. Based on this result, we proposed that miR-InDel marker technology can be used to assess genetic variation. The inclusion of geo-reference data with miR1511-InDel marker data corroborated that accessions from the Yellow River basin (Huanghuai) exhibited high genetic diversity which provides more molecular evidence for gene diversity in annual wild soybean and domestication of soybean.

CONCLUSIONS

These results provide evidence for the use of RNA marker, miRNA1511-InDel, as a soybean-specific functional maker for the study of genetic diversity, genotyping of germplasm and evolution studies. This is also the first report of functional marker developed from soybean miRNA located within the functional region of pre-miRNA1511.

De novo Assembly, Characterization of Immature Seed Transcriptome and Development of Genic-SSR Markers in Black Gram [Vigna mungo (L.) Hepper]

Black gram [V. mungo (L.) Hepper] is an important legume crop extensively grown in south and south-east Asia, where it is a major source of dietary protein for its predominantly vegetarian population. However, lack of genomic information and markers has become a limitation for genetic improvement of this crop. Here, we report the transcriptome sequencing of the immature seeds of black gram cv. TU94-2, by Illumina paired end sequencing technology to generate transcriptome sequences for gene discovery and genic-SSR marker development. A total of 17.2 million paired-end reads were generated and 48,291 transcript contigs (TCS) were assembled with an average length of 443 bp. Based on sequence similarity search, 33,766 TCS showed significant similarity to known proteins. Among these, only 29,564 TCS were annotated with gene ontology (GO) functional categories. A total number of 138 unique KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways were identified, of which majority of TCS are grouped into purine metabolism (678) followed by pyrimidine metabolism (263). A total of 48,291 TCS were searched for SSRs and 1,840 SSRs were identified in 1,572 TCS with an average frequency of one SSR per 11.9 kb. The tri-nucleotide repeats were most abundant (35%) followed by di-nucleotide repeats (32%). PCR primer pairs were successfully designed for 933 SSR loci. Sequences analyses indicate that about 64.4% and 35.6% of the SSR motifs were present in the coding sequences (CDS) and untranslated regions (UTRs) respectively. Tri-nucleotide repeats (57.3%) were preferentially present in the CDS. The rate of successful amplification and polymorphism were investigated using selected primers among 18 black gram accessions. Genic-SSR markers developed from the Illumina paired end sequencing of black gram immature seed transcriptome will provide a valuable resource for genetic diversity, evolution, linkage mapping, comparative genomics and marker-assisted selection in black gram.

Transciptome analysis reveals flavonoid biosynthesis regulation and simple sequence repeats in yam (Dioscorea alata L.) tubers

Background

Yam (Dioscorea alata L.) is an important tuber crop and purple pigmented elite cultivar has recently become popular because of associated health benefits. Identifying candidate genes responsible for flavonoid biosynthesis pathway (FBP) will facilitate understanding the molecular mechanism of controlling pigment formation in yam tubers. Here, we used Illumina sequencing to characterize the transcriptome of tubers from elite purple-flesh cultivar (DP) and conventional white-flesh cultivar (DW) of yam. In this process, we also designed high quality molecular markers to assist molecular breeding for tuber trait improvement.

Results

A total of 125,123 unigenes were identified from the DP and DW cDNA libraries, of which about 49.5% (60,020 unigenes) were annotated by BLASTX analysis using the publicly available protein database. These unigenes were further annotated functionally and subject to biochemical pathway analysis. 511 genes were identified to be more than 2-fold (FDR < 0.05) differentially expressed between the two yam cultivars, of which 288 genes were up-regulated and 223 genes were down-regulated in the DP tubers. Transcriptome analysis detected 61 unigenes encoding multiple well-known enzymes in the FBP. Furthermore, the unigenes encoding chalcone isomerase (CHS), flavanone 3-hydroxylase (F3H), flavonoid 3′-monooxygenase (F3’H), dihydroflavonol 4-reductase (DFR), leucoanthocyanidin dioxygenase (LDOX), and flavonol 3-O-glucosyltransferase (UF3GT) were found to be significantly up-regulated in the DP, implying that these genes were potentially associated with tuber color formation in this elite cultivar. The expression of these genes was further confirmed by qRT-PCR. Finally, 11,793 SSRs were successfully identified with these unigenes and 6,082 SSR markers were developed using Primer 3.

Conclusions

This study provides the first comprehensive transcriptomic dataset for yam tubers, which will significantly contribute to genomic research of this and other related species. Some key genes associated with purple-flesh trait were successfully identified, thus providing valuable information about molecular process of regulating pigment accumulation in elite yam tubers. In the future, this information might be directly used to genetically manipulate the conventional white-fleshed tuber cultivars to enable them to produce purple flesh. In addition, our SSR marker sets will facilitate identification of QTLs for various tuber traits in yam breeding programs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1547-8) contains supplementary material, which is available to authorized users.

Employing genome-wide SNP discovery and genotyping strategy to extrapolate the natural allelic diversity and domestication patterns in chickpea

The genome-wide discovery and high-throughput genotyping of SNPs in chickpea natural germplasm lines is indispensable to extrapolate their natural allelic diversity, domestication, and linkage disequilibrium (LD) patterns leading to the genetic enhancement of this vital legume crop. We discovered 44,844 high-quality SNPs by sequencing of 93 diverse cultivated desi, kabuli, and wild chickpea accessions using reference genome- and de novo-based GBS (genotyping-by-sequencing) assays that were physically mapped across eight chromosomes of desi and kabuli. Of these, 22,542 SNPs were structurally annotated in different coding and non-coding sequence components of genes. Genes with 3296 non-synonymous and 269 regulatory SNPs could functionally differentiate accessions based on their contrasting agronomic traits. A high experimental validation success rate (92%) and reproducibility (100%) along with strong sensitivity (93-96%) and specificity (99%) of GBS-based SNPs was observed. This infers the robustness of GBS as a high-throughput assay for rapid large-scale mining and genotyping of genome-wide SNPs in chickpea with sub-optimal use of resources. With 23,798 genome-wide SNPs, a relatively high intra-specific polymorphic potential (49.5%) and broader molecular diversity (13-89%)/functional allelic diversity (18-77%) was apparent among 93 chickpea accessions, suggesting their tremendous applicability in rapid selection of desirable diverse accessions/inter-specific hybrids in chickpea crossbred varietal improvement program. The genome-wide SNPs revealed complex admixed domestication pattern, extensive LD estimates (0.54-0.68) and extended LD decay (400-500 kb) in a structured population inclusive of 93 accessions. These findings reflect the utility of our identified SNPs for subsequent genome-wide association study (GWAS) and selective sweep-based domestication trait dissection analysis to identify potential genomic loci (gene-associated targets) specifically regulating important complex quantitative agronomic traits in chickpea. The numerous informative genome-wide SNPs, natural allelic diversity-led domestication pattern, and LD-based information generated in our study have got multidimensional applicability with respect to chickpea genomics-assisted breeding.

Genome-wide conserved non-coding microsatellite (CNMS) marker-based integrative genetical genomics for quantitative dissection of seed weight in chickpea

Phylogenetic footprinting identified 666 genome-wide paralogous and orthologous CNMS (conserved non-coding microsatellite) markers from 5'-untranslated and regulatory regions (URRs) of 603 protein-coding chickpea genes. The (CT)n and (GA)n CNMS carrying CTRMCAMV35S and GAGA8BKN3 regulatory elements, respectively, are abundant in the chickpea genome. The mapped genic CNMS markers with robust amplification efficiencies (94.7%) detected higher intraspecific polymorphic potential (37.6%) among genotypes, implying their immense utility in chickpea breeding and genetic analyses. Seventeen differentially expressed CNMS marker-associated genes showing strong preferential and seed tissue/developmental stage-specific expression in contrasting genotypes were selected to narrow down the gene targets underlying seed weight quantitative trait loci (QTLs)/eQTLs (expression QTLs) through integrative genetical genomics. The integration of transcript profiling with seed weight QTL/eQTL mapping, molecular haplotyping, and association analyses identified potential molecular tags (GAGA8BKN3 and RAV1AAT regulatory elements and alleles/haplotypes) in the LOB-domain-containing protein- and KANADI protein-encoding transcription factor genes controlling the cis-regulated expression for seed weight in the chickpea. This emphasizes the potential of CNMS marker-based integrative genetical genomics for the quantitative genetic dissection of complex seed weight in chickpea.

Fatty acid profile and unigene-derived simple sequence repeat markers in tung tree (Vernicia fordii)

Tung tree (Vernicia fordii) provides the sole source of tung oil widely used in industry. Lack of fatty acid composition and molecular markers hinders biochemical, genetic and breeding research. The objectives of this study were to determine fatty acid profiles and develop unigene-derived simple sequence repeat (SSR) markers in tung tree. Fatty acid profiles of 41 accessions showed that the ratio of α-eleostearic acid was increasing continuously with a parallel trend to the amount of tung oil accumulation while the ratios of other fatty acids were decreasing in different stages of the seeds and that α-eleostearic acid (18∶3) consisted of 77% of the total fatty acids in tung oil. Transcriptome sequencing identified 81,805 unigenes from tung cDNA library constructed using seed mRNA and discovered 6,366 SSRs in 5,404 unigenes. The di- and tri-nucleotide microsatellites accounted for 92% of the SSRs with AG/CT and AAG/CTT being the most abundant SSR motifs. Fifteen polymorphic genic-SSR markers were developed from 98 unigene loci tested in 41 cultivated tung accessions by agarose gel and capillary electrophoresis. Genbank database search identified 10 of them putatively coding for functional proteins. Quantitative PCR demonstrated that all 15 polymorphic SSR-associated unigenes were expressed in tung seeds and some of them were highly correlated with oil composition in the seeds. Dendrogram revealed that most of the 41 accessions were clustered according to the geographic region. These new polymorphic genic-SSR markers will facilitate future studies on genetic diversity, molecular fingerprinting, comparative genomics and genetic mapping in tung tree. The lipid profiles in the seeds of 41 tung accessions will be valuable for biochemical and breeding studies.

Development and characterization of polymorphic EST-SSR and genomic SSR markers for Tibetan annual wild barley

Tibetan annual wild barley is rich in genetic variation. This study was aimed at the exploitation of new SSRs for the genetic diversity and phylogenetic analysis of wild barley by data mining. We developed 49 novel EST-SSRs and confirmed 20 genomic SSRs for 80 Tibetan annual wild barley and 16 cultivated barley accessions. A total of 213 alleles were generated from 69 loci with an average of 3.14 alleles per locus. The trimeric repeats were the most abundant motifs (40.82%) among the EST-SSRs, while the majority of the genomic SSRs were di-nuleotide repeats. The polymorphic information content (PIC) ranged from 0.08 to 0.75 with a mean of 0.46. Besides this, the expected heterozygosity (He) ranged from 0.0854 to 0.7842 with an average of 0.5279. Overall, the polymorphism of genomic SSRs was higher than that of EST-SSRs. Furthermore, the number of alleles and the PIC of wild barley were both higher than that of cultivated barley, being 3.12 vs 2.59 and 0.44 vs 0.37. Indicating more polymorphism existed in the Tibetan wild barley than in cultivated barley. The 96 accessions were divided into eight subpopulations based on 69 SSR markers, and the cultivated genotypes can be clearly separated from wild barleys. A total of 47 SSR-containing EST unigenes showed significant similarities to the known genes. These EST-SSR markers have potential for application in germplasm appraisal, genetic diversity and population structure analysis, facilitating marker-assisted breeding and crop improvement in barley.

Characterization of EST-SSR markers in durum wheat EST library and functional analysis of SSR-containing EST fragments

The goal of this study is to identify characterization of expressed sequence tag (EST)-simple sequence repeats (SSR) markers from EST library of durum wheat and functional analysis of SSR-containing EST sequences for application in comparative genomics and breeding. 19,141 sequences were analyzed among which 18,937 ESTs were selected. Consistent with MISA results, 313 EST-SSRs were yielded. The final EST-SSRs were compared to the GenBank non-redundant database using BLASTX and classified based on these functions. Results indicated that the perfect EST-SSRs are the most frequent. The TTG/CTG imperfect EST-SSR had gamma-gliadin putative function that can be appropriate for durum wheat. Also, the mononucleotides and trinucleotides were the most frequent. Findings suggested that the identified EST-SSRs could be categorized into 83 types. Motifs TTG in trinucleotides and TC in dinucleotides had the highest frequency. TTG is the new motif in durum wheat identified in this study. We identified new EST-SSRs with more than trinucleotide and detected motifs that have potential to code amino acids. Arginine was the most frequent amino acid. Enzymes had the highest frequency among predicted functions. EST-SSRs have been identified in this study can be used for developing ESS-SSR-based detection tool for durum wheat in future studies and will be a useful resource for molecular breeding, genetics, genomics, and environmental stress studies. Motifs coding amino acids could be used as a new source of functional markers and biological study. In addition to, designed new PCR primer pairs are new resources for to identify useful alleles in transcription factors, storage proteins, and enzymes which incorporated them again into the cultivated material.

Development and characterization of simple sequence repeat (SSR) markers based on RNA-sequencing of Medicago sativa and in silico mapping onto the M. truncatula genome

Sufficient codominant genetic markers are needed for various genetic investigations in alfalfa since the species is an outcrossing autotetraploid. With the newly developed next generation sequencing technology, a large amount of transcribed sequences of alfalfa have been generated and are available for identifying SSR markers by data mining. A total of 54,278 alfalfa non-redundant unigenes were assembled through the Illumina HiSeqTM 2000 sequencing technology. Based on 3,903 unigene sequences, 4,493 SSRs were identified. Tri-nucleotide repeats (56.71%) were the most abundant motif class while AG/CT (21.7%), AGG/CCT (19.8%), AAC/GTT (10.3%), ATC/ATG (8.8%), and ACC/GGT (6.3%) were the subsequent top five nucleotide repeat motifs. Eight hundred and thirty- seven EST-SSR primer pairs were successfully designed. Of these, 527 (63%) primer pairs yielded clear and scored PCR products and 372 (70.6%) exhibited polymorphisms. High transferability was observed for ssp falcata at 99.2% (523) and 71.7% (378) in M. truncatula. In addition, 313 of 527 SSR marker sequences were in silico mapped onto the eight M. truncatula chromosomes. Thirty-six polymorphic SSR primer pairs were used in the genetic relatedness analysis of 30 Chinese alfalfa cultivated accessions generating a total of 199 scored alleles. The mean observed heterozygosity and polymorphic information content were 0.767 and 0.635, respectively. The codominant markers not only enriched the current resources of molecular markers in alfalfa, but also would facilitate targeted investigations in marker-trait association, QTL mapping, and genetic diversity analysis in alfalfa.

Transcriptome sequencing of two parental lines of cabbage (Brassica oleracea L. var. capitata L.) and construction of an EST-based genetic map

Background

Expressed sequence tag (EST)-based markers are preferred because they reflect transcribed portions of the genome. We report the development of simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers derived from transcriptome sequences in cabbage, and their utility for map construction.

Results

Transcriptome sequences were obtained from two cabbage parental lines, C1184 and C1234, which are susceptible and resistant to black rot disease, respectively, using the 454 platform. A total of 92,255 and 127,522 reads were generated and clustered into 34,688 and 40,947 unigenes, respectively. We identified 2,405 SSR motifs from the unigenes of the black rot-resistant parent C1234. Trinucleotide motifs were the most abundant (66.15%) among the repeat motifs. In addition, 1,167 SNPs were detected between the two parental lines. A total of 937 EST-based SSR and 97 SNP-based dCAPS markers were designed and used for detection of polymorphism between parents. Using an F₂ population, we built a genetic map comprising 265 loci, and consisting of 98 EST-based SSRs, 21 SNP-based dCAPS, 55 IBP markers derived from B. rapa genome sequence and 91 public SSRs, distributed on nine linkage groups spanning a total of 1,331.88 cM with an average distance of 5.03 cM between adjacent loci. The parental lines used in this study are elite breeding lines with little genetic diversity; therefore, the markers that mapped in our genetic map will have broad spectrum utility.

Conclusions

This genetic map provides additional genetic information to the existing B. oleracea map. Moreover, the new set of EST-based SSR and dCAPS markers developed herein is a valuable resource for genetic studies and will facilitate cabbage breeding. Additionally, this study demonstrates the usefulness of NGS transcriptomes for the development of genetic maps even with little genetic diversity in the mapping population.