Genome-Wide Survey and Development of the First Microsatellite Markers Database (AnCorDB) in Anemone coronaria L

A Strategy of Assessing Gene Copy Number Differentiation Between Populations Using Ultra-Fast De Novo Assembly of Next-Generation Sequencing Data

Gene duplication and loss play pivotal roles in the evolutionary dynamics of genomes, contributing to species phenotypic diversity and adaptation. However, detecting copy number variations (CNVs) in homoploid populations and newly-diverged species using short reads from next-generation sequencing (NGS) with traditional methods can often be challenging due to uneven read coverage caused by variations in GC content and the presence of repetitive sequences. To address these challenges, we developed a novel pipeline, ST4gCNV, which leverages ultra-fast de novo assemblies of NGS data to detect gene-specific CNVs between populations. The pipeline effectively reduces the variance of read coverage due to technical factors such as GC bias, providing a reliable CNV detection with a minimum sequencing depth of 10. We successfully apply ST4gCNV to the resequencing analysis of homoploid species Nelumbo nucifera and Nelumbo lutea (lotus). We reveal significant CNV-driven differentiation between these species, particularly in genes related to petal colour diversity such as those involved in the anthocyanin pathway. By highlighting the extensive gene duplication and loss events in Nelumbo, our study demonstrates the utility of ST4gCNV in population genomics and underscores its potential of integrating genomic CNV analysis with traditional SNP-based resequencing analysis.

Development of Genome-Wide SSR Markers in Leymus chinensis with Genetic Diversity Analysis and DNA Fingerprints

Leymus chinensis, a major component of the plant community in the eastern Eurasian grasslands with a wide distribution, provides stability to grassland ecosystems and supports animal husbandry. This study aimed to bridge the gap between the molecular breeding and industrial application of L. chinensis by conducting a comprehensive simple sequence repeat (SSR) analysis. A total of 973,129 SSRs were identified in the L. chinensis whole genome, which was used to design 20 polymorphic pairs of SSR primers to further assess 105 L. chinensis accessions. On average, 33.55 alleles were detected per locus, with an average Shannon index of 2.939 and a polymorphic information content value of 0.910. Principal coordinate, maximum likelihood, and structure analyses consistently showed that all samples were coincidentally divided into four subclasses. In addition, Mantel test data indicated a weak correlation between genetic and geographical distances in L. chinensis, whose variability may be related to the pollination mode and natural selection pressures. Finally, we used the 20 pairs of selected markers to scan 105 accessions, constructing a fingerprint for them. These findings provide new foundations for identifying superior varieties, improving the management of genetic resources, and constructing a germplasm resource database for L. chinensis.

The roles of florivory and herbivory in maintaining intra-population flower colour variation in Anemone coronaria

Most flowering plants are colour monomorphic, while within-population flower colour variation is rare. Multiple selection agents on flower colour, each favouring a different colour morph, may drive such uncommon polymorphisms. We tested the role of biotic antagonistic interactions in maintaining flower colour variation in Anemone coronaria (Ranunculaceae), in colour-polymorphic populations comprised of red, purple, and white flowers. We estimated the extent of leaf herbivory and petal florivory in each flower colour morph in three populations over two flowering seasons. We categorized types of damage to four groups of herbivores and estimated the plant maternal fitness. We tested pollinator response to different levels (0-30%) of simulated florivory in experimental flower arrays. Leaf and petal damage did not differ between white- and purple-flowering plants. Red-flowering plants had higher leaf damage than white-flowering plants and higher petal damage than purple-flowering plants. Nevertheless, all colour morphs had similar fitness. Red flowers exhibited more petal scratches (attributed to glaphyrid beetles), but fewer petal bites (attributed to caterpillars or grasshoppers), than white and purple flowers. Experimentally induced florivory did not reduce visits by potential pollinators in any colour morph. Glaphyrid beetles are the major pollinators of red anemone flowers, suggesting that their service to red flowers as mutualists (pollinators) should be weighed against their disservice as antagonists (florivores). A balance between pollination service and petal scratch damage of red flowers, both mediated by Glaphyird beetles, may equalize fitness between the red and the purple/white colour morphs, contributing to colour polymorphism.

Genome-wide identification of SSR markers from coding regions for endangered Argania spinosa L. skeels and construction of SSR database: AsSSRdb

Microsatellites [simple sequence repeats (SSRs)] are one of the most widely used sources of genetic markers, particularly prevalent in plants. Despite their importance in various applications, a comprehensive genome-wide identification of coding sequence (CDS)-associated SSR markers in the Argania spinosa L. genome has yet to be conducted. In this study, 66 280 CDSs containing 5351 SSRs within 4535 A. spinosa L. CDSs were identified. Among these, tri-nucleotide motifs (58.96%) were the most common, followed by hexa-nucleotide (15.71%) and di-nucleotide motifs (13.32%). The predominant SSR motif in the tri-nucleotide category was AAG (24.4%), while AG (94.1%) was the most abundant among di-nucleotide repeats. Furthermore, the extracted CDSs containing SSRs were subjected to functional annotation; 3396 CDSs (74.88%) exhibited homology with known proteins, 3341 CDSs (73.7%) were assigned Gene Ontology terms, 1004 CDSs were annotated with Enzyme Commission numbers, and 832 (18.3%) were annotated with KEGG pathways. A total of 3475 primer pairs were designed, out of which 3264 were successfully validated in silico against the A. spinosa L. genome, with 99.6% representing high-resolution markers yielding no more than three products. Additionally, the SSR markers demonstrated a low rate of transferability through in-silico verification in two species within the Sapotaceae family. Furthermore, we developed an online database, the "Argania spinosa L. SSR database: https://as-fmmdb.shinyapps.io/asssrdb/" (AsSSRdb) to provide access to the CDS-associated SSRs identified in this study. Overall, this research provides valuable marker resources for DNA fingerprinting, genetic studies, and molecular breeding in argan and related species. Database URL: https://as-fmmdb.shinyapps.io/asssrdb/.

Genome-wide development of intra- and inter-specific transferable SSR markers and construction of a dynamic web resource for yam molecular breeding: Y2MD

Microsatellite markers are widely used in population genetics and breeding. Despite the economic significance of yams in developing countries, there is a paucity of microsatellite markers, and as of now, no comprehensive microsatellite marker database exists. In this study, we conducted genome-wide microsatellite marker development across four yam species, identified cross-species transferable markers, and designed an easy-to-use web portal for the yam researchers. The screening of Dioscorea alata, Dioscorea rotundata, Dioscorea dumetorum, and Dioscorea zingiberensis genomes resulted in 318,713, 322,501, 307,040, and 253,856 microsatellites, respectively. Mono-, di-, and tri-nucleotides were the most important types of repeats in the different species, and a total of 864,128 primer pairs were designed. Furthermore, we identified 1170 cross-species transferable microsatellite markers. Among them, 17 out of 18 randomly selected were experimentally validated with good discriminatory power, regardless of the species and ploidy levels. Ultimately, we created and deployed a dynamic Yam Microsatellite Markers Database (Y2MD) available at https://y2md.ucad.sn/. Y2MD is embedded with various useful tools such as JBrowse, Blast, insilicoPCR, and SSR Finder to facilitate the exploitation of microsatellite markers in yams. This study represents the first comprehensive microsatellite marker mining across several yam species and will contribute to advancing yam genetic research and marker-assisted breeding. The released user-friendly database constitutes a valuable platform for yam researchers.

Multi-genome comprehensive identification of SSR/SV and development of molecular markers database to serve Sorghum bicolor (L.) breeding

BACKGROUND

As an important food and cash crop, identification of DNA molecular markers is of great significance for molecular marker-assisted breeding of Sorghum (Sorghum bicolor (L.) moench). Although some sorghum-related mutation databases have been published, the special SSR and SV databases still need to be constructed and updated.

RESULTS

In this study, the quality of 18 different sorghum genomes was evaluated, and two genomes were assembled at chromosome level. Through the identification and comparative analysis of SSR loci in these genomes, the distribution characteristics of SSR in the above sorghum genomes were initially revealed. At the same time, five representative reference genomes were selected to identify the structural variation of sorghum. Finally, a convenient SSR/SV database of sorghum was constructed by integrating the above results ( http://www.sorghum.top:8079/ ; http://43.154.129.150:8079/ ; http://47.106.184.91:8079/ ). Users can query the information of related sites and primer pairs.

CONCLUSIONS

Anyway, our research provides convenience for sorghum researchers and will play an active role in sorghum molecular marker-assisted breeding.

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.

Diversity analyses in two ornamental and large-genome Ranunculaceae species based on a low-cost Klenow NGS-based protocol

Persian buttercup (Ranunculus asiaticus L.) and poppy anemone (Anemone coronaria L.) are ornamental, outcrossing, perennial species belonging to the Ranunculaceae family, characterized by large and highly repetitive genomes. We applied K-seq protocol in both species to generate high-throughput sequencing data and produce a large number of genetic polymorphisms. The technique entails the application of Klenow polymerase-based PCR using short primers designed by analyzing k-mer sets in the genome sequence. To date the genome sequence of both species has not been released, thus we designed primer sets based on the reference the genome sequence of the related species Aquilegia oxysepala var. kansuensis (Brühl). A whole of 11,542 SNPs were selected for assessing genetic diversity of eighteen commercial varieties of R. asiaticus, while 1,752 SNPs for assessing genetic diversity in six cultivars of A. coronaria. UPGMA dendrograms were constructed and in R. asiaticus integrated in with PCA analysis. This study reports the first molecular fingerprinting within Persian buttercup, while the results obtained in poppy anemone were compared with a previously published SSR-based fingerprinting, proving K-seq to be an efficient protocol for the genotyping of complex genetic backgrounds.

Complete Chloroplast Genome Determination of Ranunculus sceleratus from Republic of Korea (Ranunculaceae) and Comparative Chloroplast Genomes of the Members of the Ranunculus Genus

Ranunculus sceleratus (family: Ranunculaceae) is a medicinally and economically important plant; however, gaps in taxonomic and species identification limit its practical applicability. This study aimed to sequence the chloroplast genome of R. sceleratus from Republic of Korea. Chloroplast sequences were compared and analyzed among Ranunculus species. The chloroplast genome was assembled from Illumina HiSeq 2500 sequencing raw data. The genome was 156,329 bp and had a typical quadripartite structure comprising a small single-copy region, a large single-copy region, and two inverted repeats. Fifty-three simple sequence repeats were identified in the four quadrant structural regions. The region between the ndhC and trnV-UAC genes could be useful as a genetic marker to distinguish between R. sceleratus populations from Republic of Korea and China. The Ranunculus species formed a single lineage. To differentiate between Ranunculus species, we identified 16 hotspot regions and confirmed their potential using specific barcodes based on phylogenetic tree and BLAST-based analyses. The ndhE, ndhF, rpl23, atpF, rps4, and rpoA genes had a high posterior probability of codon sites in positive selection, while the amino acid site varied between Ranunculus species and other genera. Comparison of the Ranunculus genomes provides useful information regarding species identification and evolution that could guide future phylogenetic analyses.

Genome-Wide Development of Polymorphic Microsatellite Markers and Genetic Diversity Analysis for the Halophyte Suaeda aralocaspica (Amaranthaceae)

Suaeda aralocaspica, which is an annual halophyte, grows in saline deserts in Central Asia with potential use in saline soil reclamation and salt tolerance breeding. Studying its genetic diversity is critical for effective conservation and breeding programs. In this study, we aimed to develop a set of polymorphic microsatellite markers to analyze the genetic diversity of S. aralocaspica. We identified 177,805 SSRs from the S. aralocaspica genome, with an average length of 19.49 bp, which were present at a density of 393.37 SSR/Mb. Trinucleotide repeats dominated (75.74%) different types of motifs, and the main motif was CAA/TTG (44.25%). We successfully developed 38 SSR markers that exhibited substantial polymorphism, displaying an average of 6.18 alleles with accompanying average polymorphism information content (PIC) value of 0.516. The markers were used to evaluate the genetic diversity of 52 individuals collected from three populations of S. aralocaspica in Xinjiang, China. The results showed that the genetic diversity was moderate to high, with a mean expected heterozygosity (He) of 0.614, a mean Shannon's information index (I) of 1.23, and a mean genetic differentiation index (Fst) of 0.263. The SSR markers developed in this study provide a valuable resource for future genetic studies and breeding programs of S. aralocaspica, and even other species in Suaeda.

Breeding and Growth Performance of 'Ningzhi 4', a New Blackberry Cultivar with High Yield Potential and Good Quality in China

The thornless blackberry cultivar 'Ningzhi 4' was developed by the Institute of Botany, Jiangsu Province and the Chinese Academy of Sciences (Nanjing Botanical Garden Mem. Sun Yat-Sen). The new blackberry cultivar was selected from the 'Kiowa' (female parent) and 'Hull Thornless' (male parent) F1 hybrid. 'Ningzhi 4' had excellent plant characteristics, including thornlessness, semi-erect to erect canes, vigorous growth and good disease resistance. 'Ningzhi 4' had large fruit and high yield. In addition, the parents of the superior hybrid plant were further identified by SSR markers, which provided the basis for the fingerprint of the new blackberry cultivar 'Ningzhi 4'. This is a commercial cultivar to be grown for fruit production for either shipping or local sales. It also has value as a home-garden plant. This unique type of blackberry fruit was a traditional summer fruit. This new cultivar has thornless semi-erect to erect canes and produces high-quality berries with large size, good firmness, excellent flavor, and potential for shipping and postharvest storage. The new blackberry cultivar 'Ningzhi 4' is adapted to all areas of southern China and is expected to replace or complement 'Kiowa', 'Hull Thornless', 'Chester Thornless' and 'Triple Crown'. A local cultivar patent has been approved by the Jiangsu Variety Approval Committee as 'Rubus spp. Ningzhi 4' in 2020 (S-SV-RS-014-2020). In the future, 'Ningzhi 4' could be promoted as an advantageous thornless blackberry cultivar in the main production regions of China.

Identifying potential ligand molecules EGFR mediated TNBC targeting the kinase domain-identification of customized drugs through in silico methods

Background and Purpose

Triple-negative breast cancer (TNBC) is an aggressive subtype of breast cancer in which three hormone receptors are negative. This work aimed at identifying customized potential molecules inhibiting epidermal growth factor receptor (EGFR) by exploring variants using the pharmacogenomics approaches.

Experimental approach

The pharmacogenomics approach has been followed to identify the genetic variants across the 1000 genomes continental population. Model proteins for the populations have been designed by including genetic variants in the reported positions. The 3D structures of the mutated proteins have been generated through homology modeling. The kinase domain present in the parent and the model protein molecules has been investigated. The docking study has been performed with the protein molecules against the kinase inhibitors evaluated by the molecular dynamic simulation studies. Molecular evolution has been performed to generate the potential derivatives of these kinase inhibitors suitable for the conserved region of the kinase domain. This study considered variants within the kinase domain as the sensitive region and remaining residues as the conserved region.

Findings/Results

The results reveal that few kinase inhibitors interact with the sensitive region. Among the derivatives of these kinase inhibitors molecules, the potential kinase inhibitor that interacts with the different population models has been identified.

Conclusions and implications

This study encompasses the importance of genetic variants in drug action as well as in the design of customized drugs. This research gives way to designing customized potential molecules inhibiting EGFR by exploring variants using the pharmacogenomics approaches.

Genome-Wide Survey and Analysis of Microsatellites in Waterlily, and Potential for Polymorphic Marker Development

Waterlily (Nymphaeaceae), a diploid dicotyledon, is an ornamental aquatic plant. In 2020, the complete draft genome for the blue-petal waterlily (Nymphaea colorata) was made available in GenBank. To date, the genome-wide mining of microsatellites or simple sequence repeats (SSRs) in waterlily is still absent. In the present study, we investigated the characteristics of genome-wide microsatellites for N. colorata and developed polymorphic SSR markers across tropical and hardy waterlilies. A total of 238,816 SSRs were identified in 14 N. colorata chromosomes with an average density of 662.60 SSRs per Mb, and the largest number of SSRs were present on chromosome 1 (n = 30,426, 705.94 SSRs per Mb). The dinucleotide was the most common type, and AT-rich repeats prevail in the N. colorata genome. The SSR occurrence frequencies decreased as the number of motif repeats increased. Among 2442 protein-coding region SSRs, trinucleotides, accounting for 63.84%, were the most abundant. Gene ontology terms for signal transduction (e.g., GO: 0045859 and GO: 0019887) and the lipoic acid metabolism (ko00785,) were overrepresented in GO and KEGG enrichment analysis, respectively. In addition, 107,152 primer pairs were identified, and 13 novel polymorphism SSR markers were employed to distinguish among nine waterlily cultivars, of which Ny-5.2 and Ny-10.1 were the most informative SSR loci. This study contributes the first detailed characterization of SSRs in N. colorata genomes and delivers 13 novel polymorphism markers, which are useful for the molecular breeding strategies, genetic diversity and population structure analysis of waterlily.

First genetic maps development and QTL mining in Ranunculus asiaticus L. through ddRADseq

Persian Buttercup (Ranunculus asiaticus L.; 2x=2n=16; estimated genome size: 7.6Gb) is an ornamental and perennial crop native of Asia Minor and Mediterranean basin, marketed both as cut flower or potted plant. Currently new varieties are developed by selecting plants carrying desirable traits in segregating progenies obtained by controlled mating, which are propagated through rhizomes or micro-propagated in vitro. In order to escalate selection efficiency and respond to market requests, more knowledge of buttercup genetics would facilitate the identification of markers associated with loci and genes controlling key ornamental traits, opening the way for molecular assisted breeding programs. Reduced-representation sequencing (RRS) represents a powerful tool for plant genotyping, especially in case of large genomes such as the one of buttercup, and have been applied for the development of high-density genetic maps in several species. We report on the development of the first molecular-genetic maps in R. asiaticus based on of a two-way pseudo-testcross strategy. A double digest restriction-site associated DNA (ddRAD) approach was applied for genotyping two F₁ mapping populations, whose female parents were a genotype of a so called 'ponpon' and of a 'double flower' varieties, while the common male parental ('Cipro') was a genotype producing a simple flower. The ddRAD generated a total of ~2Gb demultiplexed reads, resulting in an average of 8,3M reads per line. The sstacks pipeline was applied for the construction of a mock reference genome based on sequencing data, and SNP markers segregating in only one of the parents were retained for map construction by treating the F₁ population as a backcross. The four parental maps (two of the female parents and two of the common male parent) were aligned with 106 common markers and 8 linkage groups were identified, corresponding to the haploid chromosome number of the species. An average of 586 markers were associated with each parental map, with a marker density ranging from 1 marker/cM to 4.4 markers/cM. The developed maps were used for QTL analysis for flower color, leading to the identification of major QTLs for purple pigmentation. These results contribute to dissect on the genetics of Persian buttercup, enabling the development of new approaches for future varietal development.