Identification and Characterization of Microsatellite Loci in Maqui (Aristotelia chilensis [Molina] Stunz) Using Next-Generation Sequencing (NGS)

Genome-wide survey and characterization of microsatellites in cashew and design of a web-based microsatellite database: CMDB

The cashew is an edible tree nut crop having a wide range of food and industrial applications. Despite great economic importance, the genome-wide characterization of microsatellites [simple sequence repeats (SSRs)] in cashew is lacking. In this study, we carried out the first comprehensive genome-wide microsatellites/SSRs characterization in cashew and developed polymorphic markers and a web-based microsatellite database. A total of 54526 SSRs were discovered in the cashew genome, with a mean frequency of 153 SSRs/Mb. Among the mined genome-wide SSRs (2-6 bp size motifs), the dinucleotide repeat motifs were dominant (68.98%) followed by the trinucleotides (24.56%). The Class I type of SSRs (≥20 bp) were 45.10%, while Class II repeat motifs (≥12-<20 bp) were 54.89% of the total genomic SSRs discovered here. Further, the AT-rich SSRs occurred more frequently in the cashew genome (84%) compared to the GC-rich SSRs. The validation of the in silico-mined genome-wide SSRs by PCR screening in cashew genotypes resulted in the development of 59 polymorphic SSR markers, and the polymorphism information content (PIC) of the polymorphic SSR markers ranged from 0.19 to 0.84. Further, a web-based database, "Cashew Microsatellite Database (CMDB)," was constructed to provide access to the genome-wide SSRs mined in this study as well as transcriptome-based SSRs from our previous study to the research community through a user-friendly searchable interface. Besides, CMDB provides information on experimentally validated SSRs. CMDB permits the retrieval of SSR markers information with the customized search options. Altogether, the genome-wide SSRs characterization, the polymorphic markers and CMDB database developed in this study would serve as valuable marker resources for DNA fingerprinting, germplasm characterization, genetic studies, and molecular breeding in cashew and related Anacardium species.

Small cardamom genome: development and utilization of microsatellite markers from a draft genome sequence of Elettaria cardamomum Maton

Small cardamom (Elettaria cardamomum Maton), the queen of spices, is the third most expensive spice in the world after saffron and vanilla, valued highly for its aroma and taste. This perennial herbaceous plant is a native of coastal parts of Southern India and displays a significant amount of morphological diversity. Its genetic potential has not been exploited due to lack of genomic resources limiting our understanding of the genome and important metabolic pathways which give it the economic advantage in the spice industry. Here, we report upon the de novo assembled, draft whole genome sequence of cardamom variety, Njallani Green Gold. We used a hybrid assembly strategy using the reads from the Oxford Nanopore, Illumina and 10x Genomics GemCode sequencing chemistries. The assembled genome length was 1.06 Gb (gigabases) which is close to the estimated genome size of cardamom. More than 75% of the genome was captured in 8000 scaffolds with a N50 of 0.15 Mb. The genome appears to have a high repeat content and 68055 gene models were predicted. The genome is close to Musa species and displays an expansion and contraction in different gene families. The draft assembly was used for in silico mining of simple sequence repeats (SSRs). A total of 2,50,571 SSRs were identified of which 2,18,270 were perfect SSRs and 32,301 were compound SSRs. Among the perfect SSRs, trinucleotides were most abundant (1,25,329) and hexanucleotide repeats appear least (2,380). From the 2,50,571 SSRs mined, 2,27,808 primer pairs were designed based on flanking sequence information. Wet lab validation was performed for 246 SSR loci and based on their amplification profiles, 60 SSR markers were used for diversity analysis of a set of 60 diverse cardamom accessions. The average number of alleles detected per locus were 14.57 with a minimum of 4 and maximum of 30 alleles. Population structure analysis revealed the presence of high degree of admixtures which could primarily be due to cross-pollination prevalent in this species. The SSR markers identified would help in the development of gene or trait-linked markers which can be subsequently used for marker-assisted breeding for crop improvement in cardamom. The information on utilization of the SSR loci for generation of markers has been developed into a public database, 'cardamomSSRdb' that is freely available for use by the cardamom community.

Genome survey sequencing and mining of genome-wide microsatellite markers in yellow-billed babbler (Turdoides affinis)

Turdoides affinis is a species of group dwelling old world passerine of family Leiothrichidae. Unavailability of genome-wide sequence and species-specific molecular markers have hindered comprehensive understanding of cooperative breeding behaviour in T. affinis. Therefore, we generated genome-wide microsatellite markers through whole genome short read sequencing of T. affinis. A total of 68.8 gigabytes of paired-end raw data were sequenced containing 195,067,054 reads. Total sequenced reads spanned a coverage of 17X with genome size of 1.18 Gb. A large number of microsatellite markers (265,297) were mined in the T. affinis genome using Krait, and 50 most informative markers were identified and validated further. In-silico PCR results validated 47 markers. Of these 47 markers, five were randomly selected and validated in-vitro in twelve individuals of T. affinis. Genotyping data on these five loci estimated observed heterozygosity (H₀) and expected heterozygosity (H_e) ratios between 0.333 - 0.833 and 0.851-0.906, respectively. Effective allele size ranged from 6.698 to 10.667, inbreeding coefficient of the population ranged from 0.080 to 0.631 and null allele frequency was calculated at 0.055 to 0.303. Polymorphic information content of all the five loci varied between 0.850 and 0.906. Probabilities of exclusion and identity across 5 loci was estimated to be 0.95 and 0.0036, respectively. All the loci showed significant adherence to Hardy-Weinberg equilibrium. The microsatellite markers reported in this study will facilitate future population genetics studies on T. affinis and other congeneric species.

A Pipeline for the Development of Microsatellite Markers using Next Generation Sequencing Data

Background: Also known as Simple Sequence Repetitions (SSRs), microsatellites are profoundly informative molecular markers and powerful tools in genetics and ecology studies on plants. Objective: This research presents a workflow for developing microsatellite markers using genome skimming. Methods: The pipeline was proposed in several stages that must be performed sequentially: obtaining DNA sequences, identifying microsatellite regions, designing primers, and selecting candidate microsatellite regions to develop the markers. Our pipeline efficiency was analyzed using Illumina sequencing data from the non-model tree species Pterodon emarginatus Vog. Results: The pipeline revealed 4,382 microsatellite regions and drew 7,411 pairs of primers for P. emarginatus. However, a much larger number of microsatellite regions with the potential to develop markers were discovered from our pipeline. We selected 50 microsatellite regions with high potential for developing markers and organized 29 microsatellite regions in sets for multiplex PCR. Conclusion: The proposed pipeline is a powerful tool for fast and efficient development of microsatellite markers on a large scale in several species, especially nonmodel plant species.

Salicylic Acid Improves Antioxidant Defense System and Photosynthetic Performance in Aristotelia chilensis Plants Subjected to Moderate Drought Stress

Salicylic acid (SA) has been shown to ameliorate drought stress. However, physiological and biochemical mechanisms involved in drought stress tolerance induced by SA in plants have not been well understood. Thus, this study aimed to study the role of SA application on enzymatic and non-enzymatic antioxidants, photosynthetic performance, and plant growth in A. chilensis plants subjected to moderate drought stress. One-year-old A. chilensis plants were subjected to 100% and 60% of field capacity. When plants reached moderate drought stress (average of stem water potential of -1.0 MPa, considered as moderate drought stress), a single SA application was performed on plants. Then, physiological and biochemical features were determined at different times during 14 days. Our study showed that SA application increased 13.5% plant growth and recovered 41.9% A_N and 40.7% g_s in drought-stressed plants on day 3 compared to drought-stressed plants without SA application. Interestingly, SOD and APX activities were increased 85% and 60%, respectively, in drought-stressed SA-treated plants on day 3. Likewise, SA improved 30% total phenolic content and 60% antioxidant capacity in drought-stressed A. chilensis plants. Our study provides insight into the SA mechanism to tolerate moderate drought stress in A. chilensis plants.

Phenotypic variation of fruit and ecophysiological traits among maqui (Aristotelia chilensis [Molina] Stuntz) provenances established in a common garden

The domestication of forest species has traditionally relied on productivity issues. However, today there are concerns about the potential responses of natural populations and new cultivars to extreme climatic conditions derived from climate change and how to incorporate this knowledge into the domestication programs. Aristotelia chilensis (Molina) Stuntz ('Maqui') is a widely distributed native species in Chile. Its berry is considered a "super fruit" with an increasing interest in the food industry. This study investigated the phenotypic variation of growth, fruit, and ecophysiological traits of 20 A. chilensis clones originated from six provenances along the latitudinal gradient and established in a common-garden experiment in the Mediterranean zone of central Chile (center part of the species distribution). Differences among provenances were observed for most of the  traits under study, especially between the northern and southernmost provenances (i.e., San Fernando versus Entre Lagos). Northern provenances showed higher development of vegetative tissue and fruit yield but lower intrinsic water use efficiency (WUE_int) compared with southern ones. Clonal variation within provenances was found significant for the ripening index, WUE_int, and fruit number and weight but not significant for traits related to the crown and leaf morphology. A genetic differentiation due to latitudinal cline was not evident in this study, but differences among provenances suggest local adaptation for some traits. The genotypic variation in productive traits must be considered in the outgoing domestication of the species and future selection programs.

J, A. S, Bhat KV, Gaikwad AB. Genome wide identification and characterization of microsatellite markers in black pepper (Piper nigrum): A valuable resource for boosting genomics applications

Black pepper is one of the most valued and widely used spices in the world and dominates multi-billion dollar global spices trade. India is amongst the major producers, consumers and exporters of black pepper. In spite of its commercial and cultural importance, black pepper has received meagre attention in terms of generation of genomic resources. Availability of markers distributed throughout the genome would facilitate and accelerate genetic studies, QTL identification, genetic enhancement and crop improvement in black pepper. In this perspective, the sequence information from the recently sequenced black pepper (Piper nigrum) genome has been used for identification and characterisation of Simple Sequence Repeats (SSRs). Total 69,126 SSRs were identified from assembled genomic sequence of P. nigrum. The SSR frequency was 158 per MB making it, one SSR for every 6.3 kb in the assembled genome. Among the different types of microsatellite repeat motifs, dinucleotides were the most abundant (48.6%), followed by trinucleotide (23.7%) and compound repeats (20.62%). A set of 85 SSRs were used for validation, of which 74 produced amplification products of expected size. Genetic diversity of 30 black pepper accessions using 50 SSRs revealed four distinct clusters. Further, the cross species transferability of the SSRs was checked in nine other Piper species. Out of 50 SSRs used, 19 and 31 SSRs were amplified in nine and seven species, respectively. Thus the identified SSRs may have application in other species of the genus Piper where genome sequence is not available yet. Present study reports the first NGS based genomic SSRs in black pepper and thus constitute a valuable resource for a whole fleet of applications in genetics and plant breeding studies such as genetic map construction, QTL identification, map-based gene cloning, marker-assisted selection and evolutionary studies in Piper nigrum and related species.

Development of microsatellite markers for the threatened species Coleocephalocereus purpureus (Cactaceae) using next-generation sequencing

Ten microsatellite loci were developed and validated for the endangered cactus species Coleocephalocereus purpureus. The markers were obtained from sequences generated by whole genome shotgun sequencing approaches. A testing group of 36 specimens of the main population were genotyped and all described markers presented suitable outcomes to population genetic studies, showing polymorphic status for C. purpureus testing group with clean and reproducible amplification. No evidence for scoring errors, null alleles or linkage disequilibrium was detected. Number of alleles per locus ranged from 3 to 6 and expected heterozygosity ranged from 0.78 to 0.99. These new microsatellite loci are suitable to be used in future diversity and structure population studies of C. purpureus.

Draft genome sequence data of maqui (Aristotelia chilensis) and identification of SSR markers

Maqui (Aristotelia chilensis [Molina] Stunz) is a small dioecious tree, belonging to the Elaeocarpaceae family. Maqui fruit has high levels of antioxidant activity, which are due to elevated anthocyanin and polyphenol content. Here we describe a draft genome sequence data of maqui (A. chilensis). The genomic sequence datasets were obtained using Illumina NextSeq platform. Nucleotide sequences of raw reads and the assembled draft genome are available at NCBI's Sequence Read Archive as BioProject PRJNA544858. Also, a total of 210067 microsatellite or simple sequence repeat (SSR) markers were identified.

Optimization of Microwave-Assisted Extraction for the Recovery of Bioactive Compounds from the Chilean Superfruit (Aristotelia chilensis (Mol.) Stuntz)

Maqui berry (Aristotelia chilensis) is being explored in the food industry, and is considered to be one of the healthiest berries due to its bioactive components and high commercial value. Microwave-assisted extraction (MAE) was developed for the determination of total phenolic compounds and anthocyanins from maqui. A Box–Behnken experimental design was employed in conjunction with a response surface methodology to optimize the conditions based on 27 different experiments. The extractions were carried out with four factors (i.e., methanol percentage, pH, temperature, and solvent volume:sample mass ratio), and two responses—total phenolics and anthocyanins. Temperature and methanol percentage were found to be the most influential parameters for total phenolic compounds and anthocyanins, respectively. The optimum MAE conditions were: 65% MeOH in water at pH 2, temperature of 100 °C, and a ratio of 10:0.5 for total phenolics; and 60% MeOH in water at pH 2, temperature of 50 °C, and a ratio of 14:0.5 for anthocyanins. Kinetics assays were carried out and an optimum time of only 2 min was identified for the extractions. Repeatability and intermediate precision were also evaluated, and coefficients of variation below 5% were obtained. The new methods were successfully applied to a foodstuff made with maqui.

Transcriptome characterization and screening of molecular markers in ecologically important Himalayan species (Rhododendron arboreum)

Rhododendron arboreum is an ecologically prominent species, which also lends commercial and medicinal benefits in the form of palatable juices and useful herbal drugs. Local abundance and survival of the species under a highly fluctuating climate make it an ideal model for genetic structure and functional analysis. However, a lack of genomic data has hampered additional research. In the present study, cDNA libraries from floral and foliar tissues of the species were sequenced to provide a foundation for understanding the functional aspects of the genome and to construct an enriched repository that will promote genomics studies in the genera. Illumina’s platform facilitated the generation of ∼100 million high-quality paired-end reads. De novo assembly, clustering, and filtering out of shorter transcripts predicted 113 167 non-redundant transcripts with an average length of 1164.6 bases. Of these, 71 961 transcripts were categorized based on functional annotations in the Gene Ontology database, whereby 5710 were grouped into 141 pathways and 23 746 encoded for different transcription factors. Transcriptome screening further identified 35 419 microsatellite regions, of which, 43 polymorphic loci were characterized on 30 genotypes. Seven hundred and nineteen transcripts had 811 high-quality single-nucleotide polymorphic variants with a minimum coverage of 10, a total score of 20, and SNP% of 50.

Age-related mechanism and its relationship with secondary metabolism and abscisic acid in Aristotelia chilensis plants subjected to drought stress

Drought stress is the most important stress factor for plants, being the main cause of agricultural crop loss in the world. Plants have developed complex mechanisms for preventing water loss and oxidative stress such as synthesis of abscisic acid (ABA) and non-enzymatic antioxidant compounds such as anthocyanins, which might help plants to cope with abiotic stress as antioxidants and for scavenging reactive oxygen species. A. chilensis (Mol.) is a pioneer species, colonizing and growing on stressed and disturbed environments. In this research, an integrated analysis of secondary metabolism in Aristotelia chilensis was done to relate ABA effects on anthocyanins biosynthesis, by comparing between young and fully-expanded leaves under drought stress. Plants were subjected to drought stress for 20 days, and physiological, biochemical, and molecular analyses were performed. The relative growth rate and plant water status were reduced in stressed plants, with young leaves significantly more affected than fully-expanded leaves beginning from the 5th day of drought stress. A. chilensis plants increased their ABA and total anthocyanin content and showed upregulation of gene expression when they were subjected to severe drought (day 20), with these effects being higher in fully-expanded leaves. Multivariate analysis indicated a significant positive correlation between transcript levels for NCED1 (9-cis-epoxycarotenoid dioxygenase) and UFGT (UDP glucose: flavonoid-3-O-glucosyltransferase) with ABA and total anthocyanin, respectively. Thus, this research provides a more comprehensive analysis of the mechanisms that allow plants to cope with drought stress. This is highlighted by the differences between young and fully-expanded leaves, showing different sensibility to stress due to their ability to synthesize anthocyanins. In addition, this ability to synthesize different and high amounts of anthocyanins could be related to higher NCED1 and MYB expression and ABA levels, enhancing drought stress tolerance.