Strategies for microsatellite isolation: a review

Evolutionary and phylogenetic insights from the mitochondrial genomic analysis of Diceraeus melacanthus and D. furcatus (Hemiptera: Pentatomidae)

The mitochondrial genomes of D. melacanthus and D. furcatus were sequenced and used to investigate the phylogenetic relationships with 54 species of Pentatomidae. Their mitogenomes are 17,197 and 15,444 bp-long, respectively, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22/21 transfer RNA genes, with conserved gene arrangement. Leu, Lys, and Ser were the most common amino acids in their PCGs. PCGs evolutionary analysis indicated their mitogenomes are under purifying selection, and the most conserved genes are from the cytochrome complex, reinforcing their suitability as markers for molecular taxonomy. We identified 490 mtSSRs in 56 Pentatomidae species, with large variation and a positive correlation between mtSSR number and genome size. Three mtSSRs were identified in each Diceraeus species. Only the mtSSR in the nad6 (D. melacanthus) and nad4 (D. furcatus) appear to have application as molecular markers for species characterization. Phylogenetic analysis confirmed the monophyly of Pentatomidae. However, our analysis challenged the monophyly of Pentatominae and Podopinae. We also detected unexpected relationships among some tribes and genera, highlighting the complexity of the internal taxonomic structure of Pentatomidae. Both Diceraeus species were grouped in the same clade with the remaining Carpocorini analyzed.

Mfind: a tool for DNA barcode analysis in angiosperms and its relationship with microsatellites using a sliding window algorithm

MAIN CONCLUSION

Mfind is a tool to analyze the impact of microsatellite presence on DNA barcode specificity. We found a significant correlation between barcode entropy and microsatellite count in angiosperm. Genetic barcodes and microsatellites are some of the identification methods in taxonomy and biodiversity research. It is important to establish a relationship between microsatellite quantification and genetic information in barcodes. In order to clarify the association between the genetic information in barcodes (expressed as Shannon's Measure of Information, SMI) and microsatellites count, a total of 330,809 DNA barcodes from the BOLD database (Barcode of Life Data System) were analyzed. A parallel sliding-window algorithm was developed to compute the Shannon entropy of the barcodes, and this was compared with the quantification of microsatellites like (AT)n, (AC)n, and (AG)n. The microsatellite search method utilized an algorithm developed in the Java programming language, which systematically examined the genetic barcodes from an angiosperm database. For this purpose, a computational tool named Mfind was developed, and its search methodology is detailed. This comprehensive study revealed a broad overview of microsatellites within barcodes, unveiling an inverse correlation between the sumz of microsatellites count and barcodes information. The utilization of the Mfind tool demonstrated that the presence of microsatellites impacts the barcode information when considering entropy as a metric. This effect might be attributed to the concise length of DNA barcodes and the repetitive nature of microsatellites, resulting in a direct influence on the entropy of the barcodes.

The complete Chloroplast genome of Stachys geobombycis and comparative analysis with related Stachys species

Herb genomics, at the forefront of traditional Chinese medicine research, combines genomics with traditional practices, facilitating the scientific validation of ancient remedies. This integration enhances public understanding of traditional Chinese medicine's efficacy and broadens its scope in modern healthcare. Stachys species encompass annual or perennial herbs or small shrubs, exhibiting simple petiolate or sessile leaves. Despite their wide-ranging applications across various fields, molecular data have been lacking, hindering the precise identification and taxonomic elucidation of Stachys species. To address this gap, we assembled the complete chloroplast (CP) genome of Stachys geobombycis and conducted reannotation and comparative analysis of seven additional species within the Stachys genus. The findings demonstrate that the CP genomes of these species exhibit quadripartite structures, with lengths ranging from 14,523 to 150,599 bp. Overall, the genome structure remains relatively conserved, hosting 131 annotated genes, including 87 protein coding genes, 36 tRNA genes, and 8 rRNA genes. Additionally, 78 to 98 SSRs and long repeat sequences were detected , and notably, 6 highly variable regions were identified as potential molecular markers in the CP genome through sequence alignment. Phylogenetic analysis based on Bayesian inference and maximum likelihood methods strongly supported the phylogenetic position of the genus Stachys as a member of Stachydeae tribe. Overall, this comprehensive bioinformatics study of Stachys CP genomes lays the groundwork for phylogenetic classification, plant identification, genetic engineering, evolutionary studies, and breeding research concerning medicinal plants within the Stachys genus.

Plastid genome data provide new insights into the dynamic evolution of the tribe Ampelopsideae (Vitaceae)

BACKGROUND

Ampelopsideae J. Wen & Z.L. Nie is a small-sized tribe of Vitaceae Juss., including ca. 47 species from four genera showing a disjunct distribution worldwide across all the continents except Antarctica. There are numerous species from the tribe that are commonly used as medicinal plants with immune-modulating, antimicrobial, and anti-hypertensive properties. The tribe is usually recognized into three clades, i.e., Ampelopsis Michx., Nekemias Raf., and the Southern Hemisphere clade. However, the relationships of the three clades differ greatly between the nuclear and the plastid topologies. There has been limited exploration of the chloroplast phylogenetic relationships within Ampelopsideae, and studies on the chloroplast genome structure of this tribe are only available for a few individuals. In this study, we aimed to investigate the evolutionary characteristics of plastid genomes of the tribe, including their genome structure and evolutionary insights.

RESULTS

We sequenced, assembled, and annotated plastid genomes of 36 species from the tribe and related taxa in the family. Three main clades were recognized within Ampelopsideae, corresponding to Ampelopsis, Nekemias, and the Southern Hemisphere lineage, respectively, and all with 100% bootstrap supports. The genome sequences and content of the tribe are highly conserved. However, comparative analyses suggested that the plastomes of Nekemias demonstrate a contraction in the large single copy region and an expansion in the inverted repeat region, and possess a high number of forward and palindromic repeat sequences distinct from both Ampelopsis and the Southern Hemisphere taxa.

CONCLUSIONS

Our results highlighted plastome variations in genome length, expansion or contraction of the inverted repeat region, codon usage bias, and repeat sequences, are corresponding to the three lineages of the tribe, which probably faced with different environmental selection pressures and evolutionary history. This study provides valuable insights into understanding the evolutionary patterns of plastid genomes within the Ampelopsideae of Vitaceae.

Managing genomic diversity in conservation programs of Chinese domestic chickens

BACKGROUND

Effective conservation and utilization of farm animals are fundamental for realizing sustainable increases in food production. In situ and ex situ conservation are the two main strategies that are currently used to protect the genetic integrity of Chinese domestic chicken breeds. However, genomic diversity and population structure have not been compared in these conserved populations.

RESULTS

Three hundred and sixty-one individuals from three Chinese domestic chicken breeds were collected from populations conserved in situ and ex situ and genotyped using genotyping-by-sequencing (GBS). First, we used different parameters based on heterozygosity, genomic inbreeding, and linkage disequilibrium to estimate the genomic diversity of these populations, and applied principal component analysis (PCA), neighbor-joining tree, and ADMIXTURE to analyze population structure. We found that the small ex situ conserved populations, which have been maintained in controlled environments, retained less genetic diversity than the in situ conserved populations. In addition, genetic differentiation was detected between the in situ and ex situ conserved populations of the same breed. Next, we analyzed signatures of selection using three statistical methods (fixation index (FST), nucleotide diversity (Pi), and cross-population extended haplotype homozygosity (XP-EHH) to study the genetic footprints that underlie the differentiation between in situ and ex situ conserved populations. We concluded that, in these small populations, differentiation might be caused by genetic drift or by mutations from the original populations. The differentiation observed in the population of Beijing You chicken probably reflects adaptation to environmental changes in temperature and humidity that the animals faced when they were moved from their place of origin to the new site for ex situ conservation.

CONCLUSIONS

Conservation programs of three Chinese domestic chicken breeds have maintained their genomic diversity to a sustainable degree. The small ex situ conserved populations, which are maintained in controlled environments, retain less genetic diversity than populations conserved in situ. In addition, the transfer of populations from their place of origin to another site for conservation purposes results in genetic differentiation, which may be caused by genetic drift or adaptation. This study provides a basis for further optimization of in situ and ex situ conservation programs for domestic chicken breeds in China.

A Multiplex PCR System of Novel Microsatellite Loci for Population Genetic Application in Walnuts

Multiplex polymerase chain reaction (PCR) of microsatellite loci allows for simultaneous amplification of two or more pairs of primers in a single PCR reaction; hence, it is cost and time effective. However, very few attempts have been reported in non-model species. In this study, by combining a genome-based de novo development and cross-species application approach, a multiplex PCR system comprising 5 PCR reactions of 33 microsatellites consisting of 26 novel genomic and 7 literature-sourced loci was tested for polymorphisms, cross-species transferability, and the ability to assess genetic diversity and population structure of three walnut species (Juglans spp.). We found that the genome-based approach is more efficient than other methods. An allelic ladder was developed for each locus to enhance consistent genotyping among laboratories. The population genetic analysis results showed that all 33 loci were successfully transferred across the three species, showing high polymorphism and a strong genetic structure. Hence, the multiplex PCR system is highly applicable in walnut species. Furthermore, we propose an efficient pipeline to characterize and genotype polymorphic microsatellite loci. The novel toolbox developed here will aid future ecology and evolution studies in walnut and could serve as a model for other plant species.

Comparative Analysis of Physical and Chemical Mutagenesis in Chrysanthemum cv. 'Candid': Assessing Genetic Variation and Breeding Potential

In this study, we developed a mutagenesis protocol specifically designed for chrysanthemum cv. "Candid" in order to introduce genetic variation. By subjecting chrysanthemum shoots to different doses of physical and chemical mutagens, we successfully generated a total of 24 mutants, each with unique genetic compositions. We observed that the mortality rate was lowest when the shoots were exposed to 10 Gy gamma irradiation and 1.00% EMS. To assess the diversity and relatedness among the mutants, we employed RAPD and SSR markers. The combination of these markers allowed us to construct a dendrogram that effectively categorized the mutant population into distinct clusters based on the specific mutagen treatments. Interestingly, the mutants induced by 10 Gy gamma irradiation exhibited greater genetic diversity in terms of flower colors. On the other hand, mutants created with 1.00% EMS displayed a higher level of variation and yielded more viable mutants. To determine the optimal markers for studying genetic diversity, we analyzed the polymorphic information content (PIC) of different markers. Among the tested markers, OPA-07 (RAPD) and JH47 (SSR) showed the highest PIC values, indicating their effectiveness in capturing genetic variability within the mutant population. Conversely, the PIC values of OPD-07 and JH20 demonstrated the lowest among the markers tested. Our results revealed a percentage of polymorphism ranging from 81.81% to 100% for RAPD markers and 66.66% to 100% for SSR markers. These findings indicate that physical mutation induced by 10 Gy gamma irradiation can be clearly distinguished from chemical mutation induced by EMS at concentrations of 1% and 0.75% in chrysanthemum cv. "Candid.″ Overall, this study provides valuable insights into the genetic composition of the generated mutants and highlights their potential for enhancing chrysanthemum-breeding programs. The identified markers, particularly, OPA-07 and JH47, can serve as valuable tools for future studies aimed at exploring and exploiting the genetic diversity within the chrysanthemum population.

The novel developed microsatellite markers revealed potential hybridization among Cymbidium species and the interspecies sub-division of C. goeringii and C. ensifolium

BACKGROUND

Orchids (Cymbidium spp.) exhibit significant variations in floral morphology, pollinator relations, and ecological habitats. Due to their exceptional economic and ornamental value, Cymbidium spp. have been commercially cultivated for centuries. SSR markers are extensively used genetic tools for biology identification and population genetics analysis.

RESULT

In this study, nine polymorphic EST-SSR loci were isolated from Cymbidium goeringii using RNA-Seq technology. All nine SSR loci showed transferability in seven other congeneric species, including 51 cultivars. The novel SSR markers detected inter-species gene flow among the Cymbidium species and intra-species sub-division of C. goeringii and C. ensifolium, as revealed by neighborhood-joining and Structure clustering analyses.

CONCLUSION

In this study, we developed nine microsatellites using RNA-Seq technology. These SSR markers aided in detecting potential gene flow among Cymbidium species and identified the intra-species sub-division of C. goeringii and C. ensifolium.

A first insight into the genomic background of Ilex pubescens (Aquifoliaceae) by flow cytometry and genome survey sequencing

BACKGROUND

Ilex pubescens is an important traditional Chinese medicinal plant with many naturally occurring compounds and multiple pharmacological effects. However, the lack of reference genomic information has led to tardiness in molecular biology research and breeding programs of this plant.

RESULTS

To obtain knowledge on the genomic information of I. pubescens, a genome survey was performed for the first time by next generation sequencing (NGS) together with genome size estimation using flow cytometry. The whole genome survey of I. pubescens generated 46.472 Gb of sequence data with approximately 82.2 × coverage. K-mer analysis indicated that I. pubescens has a small genome of approximately 553 Mb with 1.93% heterozygosity rate and 39.1% repeat rate. Meanwhile, the genome size was estimated to be 722 Mb using flow cytometry, which was possibly more precise for assessment of genome size than k-mer analysis. A total of 45.842 Gb clean reads were assembled into 808,938 scaffolds with a relatively short N50 of 760 bp. The average guanine and cytosine (GC) content was 37.52%. In total, 197,429 microsatellite motifs were detected with a frequency of 2.8 kb, among which mononucleotide motifs were the most abundant (up to 62.47% of the total microsatellite motifs), followed by dinucleotide and trinucleotide motifs.

CONCLUSION

In summary, the genome of I. pubescens is small but complex with a high level of heterozygosity. Even though not successfully applied for estimation of genome size due to its complex genome, the survey sequences will help to design whole genome sequencing strategies and provide genetic information support for resource protection, genetic diversity analysis, genetic improvement and artificial breeding of I. pubescens.

Mapping the root systems of individual trees in a natural community using genotyping-by-sequencing

The architecture of root systems is an important driver of plant fitness, competition and ecosystem processes. However, the methodological difficulty of mapping roots hampers the study of these processes. Existing approaches to match individual plants to belowground samples are low throughput and species specific. Here, we developed a scalable sequencing-based method to map the root systems of individual trees across multiple species. We successfully applied it to a tropical dry forest community in the Brazilian Caatinga containing 14 species. We sequenced all 42 individual shrubs and trees in a 14 × 14 m plot using double-digest restriction site-associated sequencing (ddRADseq). We identified species-specific markers and individual-specific haplotypes from the data. We matched these markers to the ddRADseq data from 100 mixed root samples from across the centre (10 × 10 m) of the plot at four different depths using a newly developed R package. We identified individual root samples for all species and all but one individual. There was a strong significant correlation between belowground and aboveground size measurements, and we also detected significant species-level root-depth preference for two species. The method is more scalable and less labour intensive than the current techniques and is broadly applicable to ecology, forestry and agricultural biology.

Genome survey sequencing-based SSR marker development and their validation in Dendrocalamus longispathus

Bamboo is an important genetic resource in India, supporting rural livelihood and industries. Unfortunately, most Indian bamboo taxa are devoid of basic genomic or marker information required to comprehend the genetic processes for further conservation and management. In this study, we perform genome survey sequencing for development of de novo genomic SSRs in Dendrocalamus longispathus, a socioeconomically important bamboo species of northeast India. Using Illumina platform, 69.49 million raw reads were generated and assembled into 1,145,321 contig with GC content 43% and N50 1228 bp. In total, 46,984 microsatellite repeats were mined-out wherein di-nucleotide repeats were most abundant (54.71%) followed by mono- (31.91%) and tri-repeats (9.85%). Overall, AT-rich repeats were predominant in the genome, but GC-rich motifs were more frequent in tri-repeats. Afterwards, 21,596 SSR loci were successfully tagged with the primer pairs, and a subset of 50 were validated through polymerase chain reaction amplification. Of these, 36 SSR loci were successfully amplified, and 16 demonstrated polymorphism. Using 13 polymorphic SSRs, a moderate level of gene diversity (He = 0.480; Ar = 3.52) was recorded in the analysed populations of D. longispathus. Despite the high gene flow (Nm = 4.928) and low genetic differentiation (F_ST = 0.119), severe inbreeding (F_IS = 0.407) was detected. Further, genetic clustering and STRUCTURE analysis revealed that the entire genetic variability is captured under two major gene pools. Conclusively, we present a comprehensive set of novel SSR markers in D. longispathus as well as other taxa of tropical woody bamboos.

Pattern and variation in simple sequence repeat (SSR) at different genomic regions and its implications to maize evolution and breeding

BACKGROUND

Repetitive DNA sequences accounts for over 80% of maize genome. Although simple sequence repeats (SSRs) account for only 0.03% of the genome, they have been widely used in maize genetic research and breeding as highly informative codominant DNA markers. The genome-wide distribution and polymorphism of SSRs are not well studied due to the lack of high-quality genome DNA sequence data.

RESULTS

In this study, using data from high-quality de novo-sequenced maize genomes of five representative maize inbred lines, we revealed that SSRs were more densely present in telomeric region than centromeric region, and were more abundant in genic sequences than intergenic sequences. On genic sequences, tri- and hexanucleotide motifs were more abundant in CDS sequence and some mono- and dinucleotide motifs were more abundant in UTR sequences. Median length and chromosomal density of SSRs were both narrowly range-bound, with median length of 14-18 bp and genome-wide average density of 3355.77 bp/Mbp. LTR-RTs of < 0.4 Mya had higher SSR density (4498-4992 bp/Mbp). The genome-specific and motif-specific SSR polymorphism were studied. Their potential breeding applications were discussed.

CONCLUSIONS

We found that the median length of SSR sequences of different SSR motifs was nearly constant. SSR density in genic regions was much higher than intergenic regions. In addition, SSR density at LTR-RTs of different evolutionary ages varied in a narrow range. The SSRs and their LTR-RT carriers evolved at an equal rate. All these observations indicated that SSR length and density were under control of yet unknown evolutionary forces. The chromosome region-specific and motif-specific SSR polymorphisms we observed supported the notion that SSR polymorphism was invaluable genome resource for developing highly informative genome and gene markers in maize genetic research and molecular breeding.

Development of twenty-one novel microsatellite loci for Gila topminnow, Poeciliopsis occidentalis occidentalis

BACKGROUND

Gila topminnow (Poeciliopsis occidentalis occidentalis) was once highly abundant throughout the Lower Colorado River Basin of the southwestern United States. However, this Sonoran Desert endemic suffered extreme population declines over the past century because of habitat degradation and nonnative species introductions. Much of the prior conservation genetic work conducted on the species relied upon a small number of microsatellite loci; many exhibiting low variability in extant populations. Consequently, there was a need for additional microsatellite loci to provide high-resolution delimitation of populations for conservation purposes.

METHODS AND RESULTS

Paired-end Illumina sequencing was utilized to screen the Gila topminnow genome for novel microsatellite loci. We identified 21 novel loci that exhibited no deviations from expectations of genetic equilibrium, and cross-amplified in Yaqui topminnow (P. o. sonoriensis). These loci were amplified from 401 samples representing eight populations of Gila topminnow and Yaqui topminnow. Although diversity was low for all populations (observed heterozygosity = 0.12 to 0.45), these novel markers provided ample power to identify population of origin for each individual in Bayesian assignment tests.

CONCLUSIONS

This novel set of microsatellite loci provide a useful genetic tool to assess population genetic parameters of the endangered Gila topminnow and delineate populations for identifying conservation priorities. The cross-amplification of these loci in Yaqui topminnow shows promise for application to other Poeciliopsis species of Mexico and Central America.

Microsatellite Genome-Wide Database Development for the Commercial Blackhead Seabream (Acanthopagrus schlegelii)

Simple sequence repeats (SSRs), the markers with the highest polymorphism and co-dominance degrees, offer a crucial genetic research resource. Limited SSR markers in blackhead seabream have been reported. The availability of the blackhead seabream genome assembly provided the opportunity to carry out genome-wide identification for all microsatellite markers, and bioinformatic analyses open the way for developing a microsatellite genome-wide database in blackhead seabream. In this study, a total of 412,381 SSRs were identified in the 688.08 Mb genome by Krait software. Whole-genome sequences (10×) of 42 samples were aligned against the reference genome and genotyped using the HipSTR tools by comparing and counting repeat number variation across the SSR loci. A total of 156,086 SSRs with a 2–4 bp repeat were genotyped by HipSTR tools, which accounted for 55.78% of the 2–4 bp SSRs in the reference genome. High accuracy of genotyping was observed by comparing HipSTR tools and PCR amplification. A set of 109,131 loci with a number of alleles ≥ 3 and with a number of genotyped individuals ≥ 6 were reserved to constitute the polymorphic SSR database. Fifty-one polymorphic SSR loci were identified through PCR amplification. This strategy to develop polymorphic SSR markers not only obtained a large set of polymorphic SSRs but also eliminated the need for laborious experimental screening. SSR markers developed in this study may facilitate blackhead seabream research, which lays a certain foundation for further gene tagging and genetic linkage analysis, such as marker-assisted selection, genetic mapping, as well as comparative genomic analysis.

Genetic Diversity and Fine-Scale Genetic Structure of Spodoptera litura Fabricius (Lepidoptera: Noctuidae) in Southern China Based on Microsatellite Markers

Population genetic structure is strongly affected by dispersal events, especially for migratory species. The investigation of population structure is therefore conducive to increasing our understanding of species dispersal. Spodoptera litura (Fabricius) (Lepidoptera: Noctuidae) is an important tobacco pest in China causing serious damage to multiple crops. In this study, we explore its dispersal dynamics by clarifying the fine-scale population genetics using 545 S. litura samples collected from tobacco plantations at 24 locations (mainly in Baise, Hechi, and Hezhou, Southern China). We analyzed the genetic diversity, genetic structure, and gene flow of these populations using seven microsatellite loci. Our results revealed high genetic diversity and low population genetic structure among S. litura. The genetic distance was uncorrelated with geographical distance, indicating the complete randomness of dispersal among the local populations. Our results suggest that the movement scope of contemporary S. litura might be much higher than the local-level spatial scale, which will provide a theoretical basis for pest management.

Chloroplast genomes of four Carex species: Long repetitive sequences trigger dramatic changes in chloroplast genome structure

The chloroplast genomes of angiosperms usually have a stable circular quadripartite structure that exhibits high consistency in genome size and gene order. As one of the most diverse genera of angiosperms, Carex is of great value for the study of evolutionary relationships and speciation within its genus, but the study of the structure of its chloroplast genome is limited due to its highly expanded and restructured genome with a large number of repeats. In this study, we provided a more detailed account of the chloroplast genomes of Carex using a hybrid assembly of second- and third-generation sequencing and examined structural variation within this genus. The study revealed that chloroplast genomes of four Carex species are significantly longer than that of most angiosperms and are characterized by high sequence rearrangement rates, low GC content and gene density, and increased repetitive sequences. The location of chloroplast genome structural variation in the species of Carex studied is closely related to the positions of long repeat sequences; this genus provides a typical example of chloroplast structural variation and expansion caused by long repeats. Phylogenetic relationships constructed based on the chloroplast protein-coding genes support the latest taxonomic system of Carex, while revealing that structural variation in the chloroplast genome of Carex may have some phylogenetic significance. Moreover, this study demonstrated a hybrid assembly approach based on long and short reads to analyze complex chloroplast genome assembly and also provided an important reference for the analysis of structural rearrangements of chloroplast genomes in other taxa.

Development of Diagnostic Markers and Applied for Genetic Diversity Study and Population Structure of Bipolaris sorokiniana Associated with Leaf Blight Complex of Wheat

Bipolaris sorokiniana, a key pathogenic fungus in the wheat leaf blight complex, was the subject of research that resulted in the development of fifty-five polymorphic microsatellite markers. These markers were then used to examine genetic diversity and population structure in Indian geographical regions. The simple sequence repeat (SSR) like trinucleotides, dinucleotides, and tetranucleotides accounted for 43.37% (1256), 23.86% (691), and 16.54% (479) of the 2896 microsatellite repeats, respectively. There were 109 alleles produced by these loci overall, averaging 2.36 alleles per microsatellite marker. The average polymorphism information content value was 0.3451, with values ranging from 0.1319 to 0.5932. The loci's Shannon diversity varied from 0.2712 to 1.2415. These 36 isolates were divided into two main groups using population structure analysis and unweighted neighbour joining. The groupings were not based on where the isolates came from geographically. Only 7% of the overall variation was found to be between populations, according to an analysis of molecular variance. The high amount of gene flow estimate (N_M = 3.261 per generation) among populations demonstrated low genetic differentiation in the entire populations (F_ST = 0.071). The findings indicate that genetic diversity is often minimal. In order to examine the genetic diversity and population structure of the B. sorokiniana populations, the recently produced microsatellite markers will be helpful. This study's findings may serve as a foundation for developing improved management plans for the leaf blight complex and spot blotch of wheat diseases in India.

Molecular Characterization of Tinospora cordifolia (Willd.) Miers Using Novel g-SSR Markers and Their Comparison with EST-SSR and SCoT Markers for Genetic Diversity Study

In the present study, novel genomic-SSR (g-SSR) markers generated in our laboratory were used to characterize Tinospora cordifolia and related species. The g-SSR marker was also compared with EST-SSR and SCoT markers used earlier in our laboratory to assess the genetic diversity of T. cordifolia. A total of 26 accessions of T. cordifolia and 1 accession each of Tinospora rumphii and Tinospora sinensis were characterized using 65 novel g-SSR markers. A total of 125 alleles were detected with 49 polymorphic g-SSR markers. The number of alleles per locus varied from 1–4 with a mean value of 2.55 alleles per locus. Mean PIC, gene diversity, and heterozygosity were estimated to be 0.33, 0.41, and 0.65, respectively. The two species, namely T. rumphii and T. sinensis, showed cross-species transferability of g-SSRs developed in T. cordifolia. The success rate of cross-species transferability in T. rumphii was 95.3% and 93.8% in T. sinensis, proving the usefulness of this marker in genetic diversity studies of related species. The Tinospora accessions were also used for molecular characterization using SCoT and EST-SSR markers and compared for genetic diversity and cross-species transferability. The PIC, gene diversity, heterozygosity, and principal coordinate analysis showed that g-SSR is the better maker for a genetic diversity study of T. cordifolia. Additionally, high cross-species transferability of g-SSRs was found (95.3% and 93.8%) compared to EST-SSRs (68.8% and 67.7%) in T. rumphii and T. sinensis, respectively.

Pipeline for developing polymorphic microsatellites in species without reference genomes

Microsatellites, also known as simple sequence repeats (SSRs), are the preferred type of marker for many genetic applications. In conjunction with the ongoing development of next-generation sequencing, several bioinformatic tools have been developed for identifying SSRs from genomic or transcriptomic sequences. Although these tools are handy for generating polymorphic SSRs, their application almost always depends on an existing reference genome or self-assembly of the reference genome. With this in mind, we propose a pipeline for developing polymorphic SSRs that may be applied to species without reference genomes. Using a species without a reference genome (black Amur bream; Megalobrama terminalis Richardson, 1846) as a model, our pipeline was able to effectively discover polymorphic SSRs. Under different R parameters of a reference-free single nucleotide polymorphisms (SNPs) caller (ebwt2InDel), a total of 258, 208, 102, and 11 polymorphic SSRs were mined. To quantify the accuracy of the polymorphic SSRs detected using our pipeline, we analyzed 25 SSRs with PCR experiments. All primers were successfully amplified, and most SSRs (23 SSRs, 92%) were polymorphic. From the 36 individual black Amur bream, we acquired an average of 3.36 alleles per locus, ranging from one to 11. This demonstrates the effectiveness of our pipeline in identifying polymorphic SSRs and designing primers for SSR genotyping. Ultimately, our pipeline can effectively mine polymorphic SSRs for species without reference genomes, complementing SSR mining approaches based on reference genomes and helping to resolve biological issues that accompany these methods.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03313-0.

Mining and analysis of microsatellites in human coronavirus genomes using the in-house built Java pipeline

Microsatellites or simple sequence repeats are motifs of 1 to 6 nucleotides in length present in both coding and non-coding regions of DNA. These are found widely distributed in the whole genome of prokaryotes, eukaryotes, bacteria, and viruses and are used as molecular markers in studying DNA variations, gene regulation, genetic diversity and evolutionary studies, etc. However, in vitro microsatellite identification proves to be time-consuming and expensive. Therefore, the present research has been focused on using an in-house built java pipeline to identify, analyse, design primers and find related statistics of perfect and compound microsatellites in the seven complete genome sequences of coronavirus, including the genome of coronavirus disease 2019, where the host is Homo sapiens. Based on search criteria among seven genomic sequences, it was revealed that the total number of perfect simple sequence repeats (SSRs) found to be in the range of 76 to 118 and compound SSRs from 01 to10, thus reflecting the low conversion of perfect simple sequence to compound repeats. Furthermore, the incidence of SSRs was insignificant but positively correlated with genome size (R² = 0.45, p > 0.05), with simple sequence repeats relative abundance (R² = 0.18, p > 0.05) and relative density (R² = 0.23, p > 0.05). Dinucleotide repeats were the most abundant in the coding region of the genome, followed by tri, mono, and tetra. This comparative study would help us understand the evolutionary relationship, genetic diversity, and hypervariability in minimal time and cost.

Whole-Genome Survey and Microsatellite Marker Detection of Antarctic Crocodile Icefish, Chionobathyscus dewitti

Simple Summary

Crocodile icefish inhabit the deep sea around the Southern Ocean and belong to the family Channichthyidae. The species lacks hemoglobin and has evolved an antifreeze protein, unlike other teleosts. In this study, the whole-genome survey and microsatellite motifs were analyzed, which provide relevant information on genetic diversity, population genetics, and the genomic study of crocodile icefish.

Abstract

The crocodile icefish, Chionobathyscus dewitti, belonging to the family Channichthyidae, is an endemic species of the Southern Ocean. The study of its biological features and genetics is challenging as the fish inhabits the deep sea around Antarctic waters. The icefish, the sole cryopelagic species, shows unique physiological and genetic features, unlike other teleosts. It lacks hemoglobin and has evolved antifreeze proteins. Here, we report the genome sequencing data of crocodile icefish produced using the Illumina Novaseq 6000 platform. The estimated genome size was 0.88 Gb with a K-value of 19, and the unique sequence, heterozygosity, error, and duplication rates were 57.4%, 0.421%, 0.317%, and 0.738%, respectively. A genome assembly of 880.69 Mb, with an N50 scaffold length of 2401 bp, was conducted. We identified 2,252,265 microsatellite motifs from the genome assembly data, and dinucleotide repeats (1,920,127; 85.25%) had the highest rate. We selected 84 primer pairs from the genome survey assembly and randomly selected 30 primer pairs for validation. As a result, 15 primer pairs were validated as microsatellite markers.

Large-scale and small-scale population genetic structure of the medically important gastropod species Bulinus truncatus (Gastropoda, Heterobranchia)

BACKGROUND

Gastropod snails remain strongly understudied, despite their important role in transmitting parasitic diseases. Knowledge of their distribution and population dynamics increases our understanding of the processes driving disease transmission. We report the first study to use high-throughput sequencing (HTS) to elucidate the population genetic structure of the hermaphroditic snail Bulinus truncatus (Gastropoda, Heterobranchia) on a regional (17-150 km) and inter-regional (1000-5400 km) scale. This snail species acts as an intermediate host of Schistosoma haematobium and Schistosoma bovis, which cause human and animal schistosomiasis respectively.

METHODS

Bulinus truncatus snails were collected in Senegal, Cameroon, Egypt and France and identified through DNA barcoding. A single-end genotyping-by-sequencing (GBS) library, comprising 87 snail specimens from the respective countries, was built and sequenced on an Illumina HiSeq 2000 platform. Reads were mapped against S. bovis and S. haematobium reference genomes to identify schistosome infections, and single nucleotide polymorphisms (SNPs) were scored using the Stacks pipeline. These SNPs were used to estimate genetic diversity, assess population structure and construct phylogenetic trees of B. truncatus.

RESULTS

A total of 10,750 SNPs were scored and used in downstream analyses. The phylogenetic analysis identified five clades, each consisting of snails from a single country but with two distinct clades within Senegal. Genetic diversity was low in all populations, reflecting high selfing rates, but varied between locations due to habitat variability. Significant genetic differentiation and isolation by distance patterns were observed at both spatial scales, indicating that gene flow is not strong enough to counteract the effects of population bottlenecks, high selfing rates and genetic drift. Remarkably, the population genetic differentiation on a regional scale (i.e. within Senegal) was as large as that between populations on an inter-regional scale. The blind GBS technique was able to pick up parasite DNA in snail tissue, demonstrating the potential of HTS techniques to further elucidate the role of snail species in parasite transmission.

CONCLUSIONS

HTS techniques offer a valuable toolbox to further investigate the population genetic patterns of intermediate schistosome host snails and the role of snail species in parasite transmission.

Genetic structure of the small yellow croaker (Larimichthys polyactis) across the Yellow Sea and the East China Sea by microsatellite DNA variation: implications for the division of management units

The small yellow croaker, Larimichthys polyactis, is a commercial fish of the order Perciformes that mainly inhabit estuaries and coastal waters.In recent years, the resources and catch of L. polyactis have undergone huge fluctuations. To detect genetic variations caused by the fluctuation of resources, genetic diversity of L. polyactis in the coastal waters of China were analyzed in this study using microsatellite DNA marker. The results revealed high genetic diversity of this species. The STRUCTURE, DAPC and F ST results all indicated that there was no genetic structure consistent with the distribution pattern. Overall, our main findings are in agreement with previous studies, indicating that L. polyactis showed high genetic diversity and low genetic differentiation. Our results for high genetic connectivity among L. polyactis localities provide insights into the development of management strategies, that is, to manage this species as a single management unit.

Turtles in Malaysia: A Review of Conservation Status and a Call for Research

Simple Summary

Turtles are threatened all over the world. Malaysia has 24 species of turtles. This review focuses on current conservation status and some requirements for sustainability. We propose integrating concepts of ecology and molecular biology to provide almost comprehensive turtle reviews in Malaysia.

Abstract

Approximately 356 species of turtles inhabit saltwater and freshwater habitats globally, except in Antarctica. Twenty-four species of turtles have been reported in Malaysia, four of which are sea turtles. The state of Terengganu harbored the highest number of turtles, with 17 different reported species. Based on the IUCN Red List, 29% of turtle species in Malaysia are critically endangered. In comparison, another 25% are classified as endangered. Likewise, CITES reported that 67% of Malaysia’s turtles are threatened, while 25% are classified as critically endangered. This review discusses the checklists, molecular genetics work, conservation status, recent trends, and recommendations for future research. Factors contributing to their population declines and current endangered status are also discussed.

Comprehensive Analyses of Simple Sequence Repeat (SSR) in Bamboo Genomes and Development of SSR Markers with Peroxidase Genes

Simple sequence repeats (SSRs) are one of the most important molecular markers, which are widespread in plants. Bamboos are important forest resources worldwide. Here, the comprehensive identification and comparative analysis of SSRs were performed in three woody and two herbaceous bamboo species. Altogether 567,175 perfect SSRs and 71,141 compound SSRs were identified from 5737.8 Mb genome sequences of five bamboo species. Di-nucleotide SSRs were the most predominant type, with an average of ~50,152.2 per species. Most SSRs were located in intergenic regions, while those located in genic regions were relatively less. Moreover, the results of annotation distribution indicated that terms with P450, peroxidase and ATP-binding cassette transporter related to lignin biosynthesis might play important roles in woody and herbaceous bamboos under the mediation of SSRs. Furthermore, the peroxidase gene family consisted of a large number of genes containing SSRs was selected for the evolutionary relationship analysis and SSR markers development. Fifteen SSR markers derived from peroxidase family genes of Phyllostachys edulis were identified as polymorphic in 34 accessions belonging to seven genera in Bambusoideae. These results provided a comprehensive insight of SSR markers into bamboo genomes, which would facilitate bamboo research related to comparative genomics, evolution and marker-assisted selection.

Development and characterisation of SSR markers in the potato rot nematode Ditylenchus destructor

The potato rot nematode, Ditylenchus destructor, causes serious disease limiting the production of many crops. This disease usually decreases sweet potato yield by 20-50%, and in heavily infested fields the crop may be completely lost. Although the nematode has economic importance in China, its transmission route and genetic diversity are unknown. In this study, a collection of 1761 contigs of the D. destructor genome was mined for simple sequence repeat (SSR) markers, which resulted in the identification of 9745 SSRs. A total of 150 pairs of SSR primers were further developed and used for validation of the amplification rate and assessment of the polymorphism. Nine SSR markers were finally identified and analysed using 96 individual specimens of D. destructor sampled from four provinces in China. These loci were found to be moderately polymorphic with 2-8 alleles per locus. The observed and expected heterozygosity across the four populations ranged from 0.000 to 0.833 and from 0.000 to 0.666, respectively. This is the first report of the development and characterisation of genomic SSR markers in D. destructor. Our study demonstrated the obvious gene differentiation among different populations of D. destructor in China. This suggests that D. destructor in China may have been introduced from multiple origins. Much more work is needed on this species to identify patterns of spread, and the microsatellite loci we develop here should be useful in many regions for modelling range expansion, studying the evolution of resistance, and increasing the effectiveness of pest management strategies.

Parthenogenesis without costs in a grasshopper with hybrid origins

The rarity of parthenogenetic species is typically attributed to the reduced genetic variability that accompanies the absence of sex, yet natural parthenogens can be surprisingly successful. Ecological success is often proposed to derive from hybridization through enhanced genetic diversity from repetitive origins or enhanced phenotypic breadth from heterosis. Here, we tested and rejected both hypotheses in a classic parthenogen, the diploid grasshopper Warramaba virgo. Genetic data revealed a single hybrid mating origin at least 0.25 million years ago, and comparative analyses of 14 physiological and life history traits showed no evidence for altered fitness relative to its sexual progenitors. Our findings imply that the rarity of parthenogenesis is due to constraints on origin rather than to rapid extinction.

High genetic differentiation of Indo-Pacific humpback dolphins (Sousa chinensis) along the Asian Coast of the Pacific Ocean

The Indo-Pacific humpback dolphin (Sousa chinensis) is a vulnerable marine mammal species that inhabits shallow, coastal waters from Southeast China, southward throughout Southeast Asia, and westward around the Bay of Bengal to eastern India. Polymorphic microsatellites are useful for elucidating ecological and population genetics-related questions. Here, 18 new polymorphic microsatellites were developed from S. chinensis genomic DNA by Illumina MiSeq sequencing. Population genetic analyses were conducted on 42 S. chinensis individuals from three geographic locations, including the Xiamen Bay of China, the Western Gulf of Thailand, and Andaman Sea. Our microsatellite data revealed a strong and significant population structure among the three sampling regions (overall F ST = 0.371, p = .001). Pairwise mutual information index also demonstrated high levels of genetic differentiation between different region pairs (values range from 0.272 to 0.339, p < .001). Moreover, Structure analysis inferred three genetic clusters, with the high assignment probabilities of 95.92%, 99.47%, and 99.68%, respectively. Principal coordinate analysis plots of individuals divided entire genotypes into three clusters, indicating high level of genetic differentiation. Our results indicated the strong genetic structure in S. chinensis populations is a result of geographic distances. Other factors such as environmental variables, anthropogenic interference, and social behavior may also have contributed to population differentiation.

A pipeline for effectively developing highly polymorphic simple sequence repeats markers based on multi-sample genomic data

Simple sequence repeats (SSRs) are widely used genetic markers in ecology, evolution, and conservation even in the genomics era, while a general limitation to their application is the difficulty of developing polymorphic SSR markers. Next-generation sequencing (NGS) offers the opportunity for the rapid development of SSRs; however, previous studies developing SSRs using genomic data from only one individual need redundant experiments to test the polymorphisms of SSRs. In this study, we designed a pipeline for the rapid development of polymorphic SSR markers from multi-sample genomic data. We used bioinformatic software to genotype multiple individuals using resequencing data, detected highly polymorphic SSRs prior to experimental validation, significantly improved the efficiency and reduced the experimental effort. The pipeline was successfully applied to a globally threatened species, the brown eared-pheasant (Crossoptilon mantchuricum), which showed very low genomic diversity. The 20 newly developed SSR markers were highly polymorphic, the average number of alleles was much higher than the genomic average. We also evaluated the effect of the number of individuals and sequencing depth on the SSR mining results, and we found that 10 individuals and ~10X sequencing data were enough to obtain a sufficient number of polymorphic SSRs, even for species with low genetic diversity. Furthermore, the genome assembly of NGS data from the optimal number of individuals and sequencing depth can be used as an alternative reference genome if a high-quality genome is not available. Our pipeline provided a paradigm for the application of NGS technology to mining and developing molecular markers for ecological and evolutionary studies.

Assessment of Genetic Diversity and Relatedness in an Andean Potato Collection from Argentina by High-Density Genotyping

Native potatoes are the most diverse among cultivated potato species and thus constitute a valuable source for identifying genes for potato improvement. Nevertheless, high-density mapping, needed to reveal allelic diversity, has not been performed for native Argentinian potatoes. We present a study of the genetic variability and population structure of 96 Andigena potatoes from Northwestern Argentina performed using a subset of 5035 SNPs with no missing data and full reproducibility. These high-density markers are distributed across the genome and present a good coverage of genomic regions. A Bayesian approach revealed the presence of: (I) a major group comprised of most of the Andean accessions; (II) a smaller group containing the out-group cv. Spunta and the sequenced genotype DM; and (III) a third group containing colored flesh potatoes. This grouping was also consistent when maximum likelihood trees were constructed and further confirmed by a principal coordinate analysis. A group of 19 accessions stored as Andean varieties clustered consistently with group Tuberosum accessions. This was in agreement with previous studies and we hypothesize that they may be reintroductions of European-bred long day-adapted potatoes. The present study constitutes a valuable source for allele mining of genes of interest and thus provides a tool for association mapping studies.

Genetic Analysis of Plant Pathogens Natural Populations

Population genetics allow to address fundamental questions about the biology of plant pathogens. By testing specific hypotheses, population genetics provide insights into the population genetic variability of pathogens across different geographical areas, time, and associated plant hosts, as well as on the structure and differentiation of populations, and on the possibility that a population is introduced and from where it has originated. In this chapter, basic concepts of population genetics are introduced, as well as the five evolutionary factors affecting populations, that is, mutations, recombination, variation in population size, gene flow, and natural selection. A step-by-step workflow, from sampling to data analysis, on how to perform a genetic analysis of natural populations of plant pathogens is discussed. Increased knowledge of the population biology of pathogens is pivotal to improve management strategies of diseases in agricultural and forest ecosystems.

Whole Genome Sequencing and Annotation of Naematelia aurantialba (Basidiomycota, Edible-Medicinal Fungi)

Naematelia aurantialba is a rare edible fungus with both nutritional and medicinal values and especially rich in bioactive polysaccharides. However, due to the lack of genomic information, researches on the mining of active compounds, artificial breeding and cultivation, genetics, and molecular biology are limited. To facilitate the medicinal and food applications of N. aurantialba, we sequenced and analyzed the whole genome of N. aurantialba for the first time. The 21-Mb genome contained 15 contigs, and a total of 5860 protein-coding genes were predicted. The genome sequence shows that 296 genes are related to polysaccharide synthesis, including 15 genes related to nucleoside-activated sugar synthesis and 11 genes related to glucan synthesis. The genome also contains genes and gene clusters for the synthesis of other active substances, including terpenoids, unsaturated fatty acids, and bioactive proteins. In addition, it was also found that N. aurantialba was more closely related to Naematelia encephala than to Tremella fuciformis. In short, this study provides a reference for molecular cognition of N. aurantialba and related researches.

The genetic variation and relationship among the natural hybrids of Mangifera casturi Kosterm

Mangifera casturi Kosterm., a mango plant from Kalimantan Selatan, Indonesia, has limited genetic information, severely limiting the research on its genetic variation and phylogeny. We collected M. casturi's genomic information using next-generation sequencing, developed microsatellite markers and performed Sanger sequencing for DNA barcoding analysis. These markers were used to confirm parental origin and genetic diversity of M. casturi hybrids. The clean reads of the Kasturi accession were assembled de novo, producing 259 872 scaffolds (N50 = 1 445 bp). Fourteen polymorphic microsatellite markers were developed from 11 040 microsatellite motif-containing sequences. In total, 58 alleles were produced with a mean of 4.14 alleles per locus. Microsatellite marker analysis revealed broad genetic variation in M. casturi. Phylogenetic analysis was performed using internal transcribed spacers (ITS), matK, rbcL, and trnH-psbA. The phylogenetic tree of chloroplast markers placed Kasturi, Cuban, Pelipisan, Pinari, and Hambawang in one group, with M. indica as the female ancestor. Meanwhile, the phylogenetic tree of ITS markers indicated several Mangifera species as ancestors of M. casturi. Thus, M. casturi very likely originated from the cross-hybridization of multiple ancestors. Furthermore, crossing the F1 hybrids of M. indica and M. quadrifida with other Mangifera spp. may have generated much genetic variation. The genetic information for M. casturi will be a resource for breeding improvement, and conservation studies.

Development of a new set of genic SSR markers in the genus Gentiana: in silico mining, characterization and validation

Gentiana is an important genus of around 360 medicinally important species, majority of which are not well characterized. Despite its importance, very few genomic resources are available for Gentiana L. Till date, the number of informative and robust simple sequence repeat (SSR)-based markers is limited and very few efforts have been made for their development. A set of robust, freely accessible and informative SSR markers for Gentiana is a pre-requisite for any molecular systematic as well as improvement studies in this group of pharmacologically valuable plants. In view of the importance of these plants, Expressed Sequence Tag (EST) sequences of 18 Gentiana species were surveyed for the development of a large set of non-redundant SSR markers. A total of 5808 perfect SSR with an average length of 17 bp and relative abundance of 214 loci/Mb were identified in the analysed 47,487 EST sequences using Krait software. Mapping of the ESTs resulted in gene ontology annotations of 49.14% of the sequences. Based on these perfect SSRs, 2902 primer pairs were designed, and 60 markers were randomly selected and validated on a set of Gentiana kurroo Royle accessions. Among the screened markers, 39 (65%) were found to be cross-species transferable. This is the first report of the largest set of functional, novel genic SSR markers in Gentiana, which will be a valuable resource for future characterization, genotype identification, conservation and genomic studies in the various species of this group of important medicinal plants.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02969-4.

Characterization of novel genotyping-by-sequencing (GBS)-based simple sequence repeats (SSRs) and their application for population genomics of Capoeta aculeata (Valenciennes, 1844)

BACKGROUND

The species Capoeta aculeata (Valenciennes, 1844) is one of the most important freshwater species endemic to Iran. However, the investigation of the population genetic structure of this species is limited by the low number of molecular markers currently described.

METHODS AND RESULTS

In this study, we implemented next generation sequencing technology to identify polymorphic microsatellite markers and investigate the population genetic structure of C. aculeata sampled from three geographical sites in Iran. We characterized and developed 36 novel polymorphic microsatellite markers and these loci were examined in 120 individuals from three populations occurring in the Zagros basin. The average number of alleles per locus varied from 1.7 to 16 (average = 7.89). The results showed that, the polymorphism information content (PIC) of these simple sequence repeat (SSR) loci varied from 0.254 to 0.888. The observed heterozygosity (H_O) per locus ranged from 0.170 to 0.881, while the expected heterozygosity (H_E) per locus was from 0.170 to 0.881. Among these SSR loci, 20 loci deviated significantly from the Hardy-Weinberg equilibrium after Bonferroni correction (p < 0.05).

CONCLUSIONS

These microsatellite markers could provide a valuable tool for future population and conservation genetics studies of C. aculeate and other closely related species.

Plastome evolution in the Caesalpinia group (Leguminosae) and its application in phylogenomics and populations genetics

The chloroplast genomes of Caesalpinia group species are structurally conserved, but sequence level variation is useful for both phylogenomic and population genetic analyses. Variation in chloroplast genomes (plastomes) has been an important source of information in plant biology. The Caesalpinia group has been used as a model in studies correlating ecological and genomic variables, yet its intergeneric and infrageneric relationships are not fully solved, despite densely sampled phylogenies including nuclear and plastid loci by Sanger sequencing. Here, we present the de novo assembly and characterization of plastomes from 13 species from the Caesalpinia group belonging to eight genera. A comparative analysis was carried out with 13 other plastomes previously available, totalizing 26 plastomes and representing 15 of the 26 known Caesalpinia group genera. All plastomes showed a conserved quadripartite structure and gene repertoire, except for the loss of four ndh genes in Erythrostemon gilliesii. Thirty polymorphic regions were identified for inter- or intrageneric analyses. The 26 aligned plastomes were used for phylogenetic reconstruction, revealing a well-resolved topology, and dividing the Caesalpinia group into two fully supported clades. Sixteen microsatellite (cpSSR) loci were selected from Cenostigma microphyllum for primer development and at least two were cross-amplified in different Leguminosae subfamilies by in vitro or in silico approaches. Four loci were used to assess the genetic diversity of C. microphyllum in the Brazilian Caatinga. Our results demonstrate the structural conservation of plastomes in the Caesalpinia group, offering insights into its systematics and evolution, and provides new genomic tools for future phylogenetic, population genetics, and phylogeographic studies.

Development, characterization, functional annotation and validation of genomic and genic-SSR markers using de novo next generation sequencing in Melia dubia Cav

Melia dubia Cav. (Meliaceae), a fast-growing tropical tree finds use in plywood, pulp and high-value solid wood products. To increase its productivity, we must essentially capture genetic diversity and identify genotypes with superior wood properties. This study aimed to develop novel microsatellite markers from genomic data and validate the markers in M. dubia. Direct Seq-to-SSR approach was adopted and using an in-house Perl script, 426,390 SSR markers identified. For validation, selected 151 markers, of which 50 were genomic markers chosen randomly, and 101 were genic markers identified through BLAST2GO. Amplification was observed in all loci, and 81.4% generated high-quality, reproducible amplicons of the expected size. Out of 50 genomic markers, we used ten highly polymorphic markers to assess genetic diversity among 75 genotypes from three populations. One hundred fourteen alleles were recorded, with a moderate level of diversity and a positive fixation index. Twenty-nine genic markers representing 13 enzymes showing polymorphism for wood stiffness were selected for diversity assessment of 24 genotypes (12 genotypes each with high and low-stress wave velocity). The product size ranged from 87 to 279, covering the majority of the genome. Cluster and structure analysis segregated ~ 80% of the genotypes based on the trait. This is the first report of the development of genic markers from a genomic survey and has proved efficient in differentiating genotypes based on the trait. The markers developed in this study will be useful for genetic mapping, diversity estimation, marker-assisted selection for desired traits and breeding for wood traits in M. dubia.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02858-w.

Development of polymorphic microsatellites for genetic studies of white scar oyster (Crassostrea belcheri) using paired-end shotgun sequencing

White scar oyster Crassostrea belcheri is a commercially important bivalve species in Thailand. Appropriate genetic markers are needed for effective management to elevate its production efficiency. Type II microsatellites of C. belcheri were identified and characterized using an Illumina paired-end shotgun sequencing. A total of 14,743,710 reads were generated for which 198,849 reads containing microsatellites and 217,998 microsatellite loci were found. Twenty out of 60 microsatellite loci (33.33%) were polymorphic and these microsatellites were further tested against DNA bulks (N = 10 each) originating from 7 different geographic locations in Thai waters. Results indicated that newly developed microsatellites can be used for genetic diversity analysis of C. belcheri. Genotyping of C. belcheri collected from Surat Thani (Gulf of Thailand; N = 50) were performed. The number of alleles per locus ranged from 2 to 12 (average = 4.95). Observed and expected heterozygosities ranged from 0.0000 to 0.9400 (average = 0.3419) and 0.1139 to 0.8190 (average = 0.5844), respectively. Genome information and 20 newly isolated microsatellites will facilitate further studies in population genetics, stock management, and genetic improvement of C. belcheri in Thailand.

Recent demographic histories of temperate deciduous trees inferred from microsatellite markers

Background

Accurate inference of demographic histories for temperate tree species can aid our understanding of current climate change as a driver of evolution. Microsatellites are more suitable for inferring recent historical events due to their high mutation rates. However, most programs analyzing microsatellite data assume a strict stepwise mutation model (SMM), which could cause false detection of population shrinkage when microsatellite mutation does not follow SMM.

Results

This study aims to reconstruct the recent demographic histories of five cool-temperate tree species in Eastern Asia, Quercus mongolica, Q. liaotungensis, Juglans cathayensis, J. mandshurica and J. ailantifolia, by using 19 microsatellite markers with two methods considering generalized stepwise mutation model (GSM) (MIGRAINE and VarEff). Both programs revealed that all the five species experienced expansions after the Last Glacial Maximum (LGM). Within butternuts, J. cathayensis experienced a more serious bottleneck than the other species, and within oaks, Q. mongolica showed a moderate increase in population size and remained stable after the expansion. In addition, the point estimates of the multistep mutation proportion in the GSM model (p_GSM) for all five species were between 0.50 and 0.65, indicating that when inferring population demographic history of the cool-temperate forest species using microsatellite markers, it is better to assume a GSM rather than a SMM.

Conclusions

This study provides the first direct evidence that five cool-temperate tree species in East Asia have experienced expansions after the LGM with microsatellite data. Considering the mutation model of microsatellite has a vital influence on demographic inference, combining multiple programs such as MIGRAINE and VarEff can effectively reduce errors caused by inappropriate model selection and prior setting.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01805-w.

Whole genome survey analysis and microsatellite motif identification of Sebastiscus marmoratus

The marbled rockfish Sebastiscus marmoratus is an ecologically and economically important marine fish species distributed along the northwestern Pacific coast from Japan to the Philippines. Here, next-generation sequencing was used to generate a whole genome survey dataset to provide fundamental information of its genome and develop genome-wide microsatellite markers for S. marmoratus. The genome size of S. marmoratus was estimated as approximate 800 Mb by using K-mer analyses, and its heterozygosity ratio and repeat sequence ratio were 0.17% and 39.65%, respectively. The preliminary assembled genome was nearly 609 Mb with GC content of 41.3%, and the data were used to develop microsatellite markers. A total of 191,592 microsatellite motifs were identified. The most frequent repeat motif was dinucleotide with a frequency of 76.10%, followed by 19.63% trinucleotide, 3.91% tetranucleotide, and 0.36% pentanucleotide motifs. The AC, GAG, and ATAG repeats were the most abundant motifs of dinucleotide, trinucleotide, and tetranucleotide motifs, respectively. In summary, a wide range of candidate microsatellite markers were identified and characterized in the present study using genome survey analysis. High-quality whole genome sequence based on the “Illumina+PacBio+Hi-C” strategy is warranted for further comparative genomics and evolutionary biology studies in this species.

Development and characterization of Novel EST-based single-copy genic microsatellite DNA markers in white spruce and black spruce

Due mainly to large genome size and prevalence of repetitive sequences in the nuclear genome of spruce (Picea Mill.), it is very difficult to develop single-copy genomic microsatellite markers. We have developed and characterized 25 polymorphic, single-copy genic microsatellites from white spruce (Picea glauca (Moench) Voss) EST sequences and determined their informativeness in white spruce and black spruce (Picea mariana (Mill.) B.S.P.) and inheritance in black spruce. White spruce EST sequences from NCBI dbEST were searched for the presence of microsatellite repeats. Forty-seven sequences containing dinucleotide, trinucleotide, tetranucleotide and compound repeats were selected to develop primers. Twenty-five of the designed primer pairs yielded scorable amplicons, with single-locus patterns, and were characterized in 20 individuals each of white spruce and black spruce. All 25 microsatellites were polymorphic in white spruce and 24 in black spruce. The number of alleles at a locus ranged from two to 18, with a mean of 8.8 in white spruce, and from one to 17, with a mean of 7.6 in black spruce. The expected heterozygosity/polymorphic information content ranged from 0.10 to 0.92, with a mean of 0.67 in white spruce, and from 0 to 0.93, with a mean of 0.59 in black spruce. Microsatellites with dinucleotide and compound repeats were more informative than those with trinucleotide and tetranucleotide repeats. Eighteen microsatellite markers polymorphic between the parents of a black spruce controlled cross inherited in a single-locus Mendelian fashion. The microsatellite markers developed can be applied for various genetics, genomics, breeding, and conservation studies and applications.

Insight into morphological and molecular variations across Iranian spinach landraces

Spinach is salt-tolerant leafy vegetable stemmed from central Asia. It has been well adapted to different climates of Iran. We aimed to study genetic diversity between several landraces as a prerequisite for crop improvement programs using molecular markers including microsatellites (SSRs) and morphological traits. Genetic diversity was studied among 22 spinach landraces using morphological and molecular tools. We developed 17 genic and genomic SSR markers based on the information acquired from NCBI resources. Morphological evaluation indicated high variability for economic traits including leaf color, leaf thickness, leaf wrinkle and seed type, across Iranian landraces. The molecular results exhibited that 12 out of 17 primer pairs successfully amplified genomic DNA with explicit bands. The results verified that genic markers were superior to genomic markers to detect polymorphism and genetic diversity. In this regard, PIC for genomic and genic SSRs was in the range of 0.4616-0.6621 and 0.5188-0.7394, respectively. Polymorphic genic SSRs were identified to be directly and indirectly involved in biotic/abiotic stresses. High degree of polymorphism, which was detected across the landraces by genic SSRs, could assist us to select 11 landraces for the second experiment. The second experiment was designed to evaluate the response of selected landraces to salinity stress. The results confirmed genetic variability among the landraces in terms of salinity tolerance. A highly diverse germplasm of Iranian spinach based on molecular and morphological characteristics along with the tolerance to oxidative stress provides an ample opportunity for plant breeders to select superior genotypes.

Population Genetic Structure of A Marine Pelagic Egg Producer and Popular Marine Aquarium Species, the Mandarinfish Synchiropus splendidus

The mandarinfish Synchiropus splendidus is extensively collected in Southeast Asia (mainly in the Philippines) and highly favoured for the marine aquarium trade. Males are more popular than females for their large first dorsal fins and the fishery is not managed. To examine possible population replenishment dynamics arising as a result of selective fishing, the effects of sex-selective fishing on sex ratios and population connectivity were considered. This study determined the sex ratios and analyzed the population genetic structure from mandarinfish collected at six locations: one from Palau, where the species is not exploited, and five from Bohol in the Philippines, where the species has long been heavily fished. The findings reported very low male to female ratios (0.12 to 0.30) from four of the five locations in Bohol, with relatively more males to females in the specimens collected from Palau (2.3). The analyses from allozymes (43 alleles from 10 loci) and microsatellites (118 alleles from 5 loci) revealed that genetic connectivity was high among the five locations in the Philippines as well as with the specimens collected from the more-distant Palau. The genetic homogeneity observed across the geographical range considered is inconsistent with the hypothesized limited dispersal ability of the species and could be explained by recent species range expansion associated with sea level rise in the region. The results suggest that the present genetic structure, at least in the geographic region considered, may not be determined by current patterns of gene flow, but may, instead, be driven by recent sea-level changes associated with periods of glaciation. Caution is suggested to ensure that heavily localized fishing does not produce excessively biased adult sex ratios.

citSATdb: Genome-Wide Simple Sequence Repeat (SSR) Marker Database of Citrus Species for Germplasm Characterization and Crop Improvement

Microsatellites or simple sequence repeats (SSRs) are popular co-dominant markers that play an important role in crop improvement. To enhance genomic resources in general horticulture, we identified SSRs in the genomes of eight citrus species and characterized their frequency and distribution in different genomic regions. Citrus is the world's most widely cultivated fruit crop. We have implemented a microsatellite database, citSATdb, having the highest number (~1,296,500) of putative SSR markers from the genus Citrus, represented by eight species. The database is based on a three-tier approach using MySQL, PHP, and Apache. The markers can be searched using multiple search parameters including chromosome/scaffold number(s), motif types, repeat nucleotides (1-6), SSR length, patterns of repeat motifs and chromosome/scaffold location. The cross-species transferability of selected markers can be checked using e-PCR. Further, the markers can be visualized using the Jbrowse feature. These markers can be used for distinctness, uniformity, and stability (DUS) tests of variety identification, marker-assisted selection (MAS), gene discovery, QTL mapping, and germplasm characterization. citSATdb represents a comprehensive source of markers for developing/implementing new approaches for molecular breeding, required to enhance Citrus productivity. The potential polymorphic SSR markers identified by cross-species transferability could be used for genetic diversity and population distinction in other species.

Mining the red deer genome (CerEla1.0) to develop X-and Y-chromosome-linked STR markers

Microsatellites are widely applied in population and forensic genetics, wildlife studies and parentage testing in animal breeding, among others, and recently, high-throughput sequencing technologies have greatly facilitated the identification of microsatellite markers. In this study the genomic data of Cervus elaphus (CerEla1.0) was exploited, in order to identify microsatellite loci along the red deer genome and for designing the cognate primers. The bioinformatics pipeline identified 982,433 microsatellite motifs genome-wide, assorted along the chromosomes, from which 45,711 loci mapped to the X- and 1096 to the Y-chromosome. Primers were successfully designed for 170,873 loci, and validated with an independently developed autosomal tetranucleotide STR set. Ten X- and five Y-chromosome-linked microsatellites were selected and tested by two multiplex PCR setups on genomic DNA samples of 123 red deer stags. The average number of alleles per locus was 3.3, and the average gene diversity value of the markers was 0.270. The overall observed and expected heterozygosities were 0.755 and 0.832, respectively. Polymorphic Information Content (PIC) ranged between 0.469 and 0.909 per locus with a mean value of 0.813. Using the X- and Y-chromosome linked markers 19 different Y-chromosome and 72 X-chromosome lines were identified. Both the X- and the Y-haplotypes split to two distinct clades each. The Y-chromosome clades correlated strongly with the geographic origin of the haplotypes of the samples. Segregation and admixture of subpopulations were demonstrated by the use of the combination of nine autosomal and 16 sex chromosomal STRs concerning southwestern and northeastern Hungary. In conclusion, the approach demonstrated here is a very efficient method for developing microsatellite markers for species with available genomic sequence data, as well as for their use in individual identifications and in population genetics studies.

The genetic legacy of fragmentation and overexploitation in the threatened medicinal African pepper-bark tree, Warburgia salutaris

The pepper-bark tree (Warburgia salutaris) is one of the most highly valued medicinal plant species worldwide. Native to southern Africa, this species has been extensively harvested for the bark, which is widely used in traditional health practices. Illegal harvesting coupled with habitat degradation has contributed to fragmentation of populations and a severe decline in its distribution. Even though the species is included in the IUCN Red List as Endangered, genetic data that would help conservation efforts and future re-introductions are absent. We therefore developed new molecular markers to understand patterns of genetic diversity, structure, and gene flow of W. salutaris in one of its most important areas of occurrence (Mozambique). In this study, we have shown that, despite fragmentation and overexploitation, this species maintains a relatively high level of genetic diversity supporting the existence of random mating. Two genetic groups were found corresponding to the northern and southern locations. Our study suggests that, if local extinctions occurred in Mozambique, the pepper-bark tree persisted in sufficient numbers to retain a large proportion of genetic diversity. Management plans should concentrate on maintaining this high level of genetic variability through both in and ex-situ conservation actions.