SSR Locator: Tool for Simple Sequence Repeat Discovery Integrated with Primer Design and PCR Simulation

Studying Human Pathogenic Cryptococcus Gattii Lineages by Utilizing Simple Sequence Repeats to Create Diagnostic Markers and Analyzing Diversity

In this study, we compared the occurrence, relative abundance (RA), and density (RD) of simple sequence repeats (SSRs) among the lineages of human pathogenic Cryptococcus gattii using an in-silico approach to gain a deeper understanding of the structure and evolution of their genomes. C. gattii isolate MF34 showed the highest RA and RD of SSRs in both the genomic and transcriptomic sequences, followed by isolate WM276. In both the genomic (50%) and transcriptomic (65%) sequences, trinucleotide SSRs were the most common SSR class. A motif conservation study found that the isolates had stronger conservation (56.1%) of motifs, with isolate IND107 having the most (5.7%) unique motifs. We discovered the presence of SSRs in genes that are directly or indirectly associated with disease using gene enrichment analysis. Isolate-specific unique motifs identified in this study could be utilized as molecular probes for isolate identification. To improve genetic resources among C. gattii isolates, 6499 primers were developed. These genomic resources developed in this study could help with diversity analysis and the development of isolate-specific markers.

Comparative genome analysis reveals driving forces behind Monkeypox virus evolution and sheds light on the role of ATC trinucleotide motif

Monkeypox (MPOX), a zoonotic disease originating in Western and Central Africa in 1970, has seen a recent surge in outbreaks across 100+ countries. A comparative analysis of 404 Monkeypox virus (MPXV) genomes revealed notable changes in microsatellite abundance and density, especially within Clades I, IIa, and IIb. Each clade exhibited unique microsatellite motifs, with twenty-six conserved loci specific to MPXV, suggesting their potential as molecular markers in diagnostics. Additionally, nine genes in the MPXV genome featured ten variable hotspot microsatellite regions associated with surface protein synthesis and host control. Notably, gene OPG153, especially at the SSR locus '(ATC)n', exhibited the most pronounced variations among lineages over time and plays a role in virus pathogenesis within the host cell. These findings not only enhance our understanding of MPXV unique molecular profile but also offer valuable insights into potential pathogenic and evolutionary implications.

Distribution and conservation of simple sequence repeats in plant pathogenic species of Zymoseptoria and development of genomic resources for its orphaned species

To better understand the structure and evolution of the genomes of four plant pathogenic species of Zymoseptoria, we analyzed the occurrence, relative abundance (RA), and density (RD) of simple sequence repeats (SSRs) in their whole genome and transcriptome sequences. In this study, SSRs are defined as repeats of more than 12 bases in length. The genome and transcriptome sequences of Zymoseptoria ardabiliae show the highest RA (201.1 and 129.9) and RD (3229.4 and 1928.2) of SSRs, while those of Zymoseptoria pseudotritici show the lowest RA (167.2 and 118.5) and RD (2482.2 and 1687.0). The majority of SSRs in the genomic and transcriptome sequences of species were trinucleotide SSRs, while dinucleotide SSRs were the least common. The most common trinucleotide motifs in the transcriptomic sequences across all species were those that encoded the amino acid arginine. As per our motif conservation study, Zymoseptoria tritici (12.4%) possessed the most unique motifs, while Z. pseudotritici (3.9%) had the fewest. Overall, only 38.1% of the motifs were found to be conserved among the species. Gene enrichment studies reveal that three of the species, Z. ardabiliae, Zymoseptoria brevis, and Z. pseudotritici, have SSRs in their genes related to cellular metabolism, while the remaining Z. tritici harbors SSRs in genes related to DNA synthesis and gene expression. In an effort to improve the genetic resources for the orphan species of pathogenic Zymoseptoria, a total of 73,134 primers were created. The genomic resources developed in this study could help with analyses of genetic relatedness within the population and the development of species-specific markers.

EasySSR: a user-friendly web application with full command-line features for large-scale batch microsatellite mining and samples comparison

Microsatellites, also known as SSRs or STRs, are polymorphic DNA regions with tandem repetitions of a nucleotide motif of size 1-6 base pairs with a broad range of applications in many fields, such as comparative genomics, molecular biology, and forensics. However, the majority of researchers do not have computational training and struggle while running command-line tools or very limited web tools for their SSR research, spending a considerable amount of time learning how to execute the software and conducting the post-processing data tabulation in other tools or manually-time that could be used directly in data analysis. We present EasySSR, a user-friendly web tool with command-line full functionality, designed for practical use in batch identifying and comparing SSRs in sequences, draft, or complete genomes, not requiring previous bioinformatic skills to run. EasySSR requires only a FASTA and an optional GENBANK file of one or more genomes to identify and compare STRs. The tool can automatically analyze and compare SSRs in whole genomes, convert GenBank to PTT files, identify perfect and imperfect SSRs and coding and non-coding regions, compare their frequencies, abundancy, motifs, flanking sequences, and iterations, producing many outputs ready for download such as PTT files, interactive charts, and Excel tables, giving the user the data ready for further analysis in minutes. EasySSR was implemented as a web application, which can be executed from any browser and is available for free at https://computationalbiology.ufpa.br/easyssr/. Tutorials, usage notes, and download links to the source code can be found at https://github.com/engbiopct/EasySSR.

Comparative genomics reveals the presence of simple sequence repeats in genes related to virulence in plant pathogenic Pythium ultimum and Pythium vexans

In this study, we evaluated the occurrence, relative abundance (RA), and density (RD) of simple sequence repeats (SSRs) in the complete genome and transcriptomic sequences of the plant pathogenic species of Pythium to acquire a better knowledge of their genome structure and evolution. Among the species, P. ultimum had the highest RA and RD of SSRs in the genomic sequences, whereas P. vexans had the highest RA and RD in the transcriptomic sequences. The genomic and transcriptomic sequences of P. aphanidermatum showed the lowest RA and RD of SSRs. Trinucleotide SSRs were the most prevalent class in both genomic and transcriptomic sequences, while dinucleotide SSRs were the least prevalent. The G + C content of the transcriptomic sequences was found to be positively correlated with the number (r = 0.601) and RA (r = 0.710) of SSRs. A motif conservation study revealed the highest number of unique motifs in P. vexans (9.9%). Overall, a low conservation of motifs was observed among the species (25.9%). A gene enrichment study revealed P. vexans and P. ultimum carry SSRs in their genes that are directly connected to virulence, whereas the remaining two species, P. aphanidermatum and P. arrhenomanes, harbour SSRs in genes involved in transcription, translation, and ATP binding. In an effort to enhance the genomic resources, a total of 11,002 primers from the transcribed regions were designed for the pathogenic Pythium species. Furthermore, the unique motifs identified in this work could be employed as molecular probes for species identification.

QTL mapping for seed density per silique in Brassica napus

Seed density per silique (SDPS) and valid silique length (VSL) are two important yield-influencing traits in rapeseed. SDPS has a direct or indirect effect on rapeseed yield through its effect on seed per silique. In this study, a quantitative trait locus (QTL) for SDPS was detected on chromosome A09 using the QTL-seq approach and confirmed via linkage analysis in the mapping population obtained from 4263 × 3001 cross. Furthermore, one major QTL for SDPS (qSD.A9-1) was mapped to a 401.8 kb genomic interval between SSR markers Nys9A190 and Nys9A531. In the same genomic region, a QTL (qSL.A9) linked to VSL was also detected. The phenotypic variation of qSD.A9-1 and qSL.A9 was 53.1% and 47.6%, respectively. Results of the additive and dominant effects demonstrated that the expression of genes controlling SDPS and VSL were derived from a different parent in this population. Subsequently, we identified 56 genes that included 45 specific genes with exonic (splicing) variants. Further analysis identified specific genes containing mutations that may be related to seed density as well as silique length. These genes could be used for further studies to understand the details of these traits of rapeseed.

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.

Genetic Analysis and Fine Mapping of ZmGHT1 Conferring Glufosinate Herbicide Tolerance in Maize (Zea mays L.)

Weed interference in the crop field is one of the major biotic stresses causing dramatic crop yield losses, and the development of herbicide-resistant crops is critical for weed control in the application of herbicide technologies. To identify herbicide-resistant germplasms, we screened 854 maize inbreed lines and 25,620 seedlings by spraying them with 1 g/L glufosinate. One plant (L336R), possibly derived from a natural variation of line L336, was identified to have the potential for glufosinate tolerance. Genetic analysis validated that the glufosinate tolerance of L336R is conferred by a single locus, which was tentatively designated as ZmGHT1. By constructing a bi-parental population derived from L336R, and a glufosinate sensitive line L312, ZmGHT1 was mapped between molecular markers M9 and M10. Interestingly, genomic comparation between the two sequenced reference genomes showed that large scale structural variations (SVs) occurred within the mapped region, resulting in 2.16 Mb in the inbreed line B73, and 11.5 kb in CML277, respectively. During the fine mapping process, we did not detect any additional recombinant, even by using more than 9500 F₂ and F₃ plants, suspecting that SVs should also have occurred between L336R and L312 in this region, which inhibited recombination. By evaluating the expression of the genes within the mapped interval and using functional annotation, we predict that the gene Zm00001eb361930, encoding an aminotransferase, is the most likely causative gene. After glufosinate treatment, lower levels of ammonia content and a higher activity of glutamine synthetase (GS) in L336R were detected compared with those of L336 and L312, suggesting that the target gene may participate in ammonia elimination involving GS activity. Collectively, our study can provide a material resource for maize herbicide resistant breeding, with the potential to reveal a new mechanism for herbicide resistance.

Simple sequence repeat insertion induced stability and potential 'gain of function' in the proteins of extremophilic bacteria

Here, we analysed the genomic evolution in extremophilic bacteria using long simple sequence repeats (SSRs). Frequencies of occurrence, relative abundance (RA) and relative density (RD) of long SSRs were analysed in the genomes of extremophilic bacteria. Thermus aquaticus had the most RA and RD of long SSRs in its coding sequences (110.6 and 1408.3), followed by Rhodoferax antarcticus (77.0 and 1187.4). A positive correlation was observed between G + C content and the RA-RD of long SSRs. Geobacillus kaustophilus, Geobacillus thermoleovorans, Halothermothrix orenii, R. antarcticus, and T. aquaticus preferred trinucleotide repeats within their genomes, whereas others preferred a higher number of tetranucleotide repeats. Gene enrichment showed the presence of these long SSRs in metabolic enzyme encoding genes related to stress tolerance. To analyse the functional implications of SSR insertions, three-dimensional protein structure modelling of SSR containing diguanylate cyclase (DGC) gene encoding protein was carried out. Removal of SSR sequence led to an inappropriate folding and instability of the modelled protein structure.

Genomic markers linked to Meloidogyne chitwoodi resistance introgressed from Solanum bulbocastanum to cultivated potato and their utility in marker-assisted selection

Meloidogyne chitwoodi is a major threat to potato production in the Pacific Northwest region of the United States. Infected tubers are rendered unmarketable; hence, growers' profitability is adversely affected. Breeding for nematode resistance is a long-term process and phenotyping the segregating populations for nematode resistance is the most time-consuming and laborious part of the process. Using DNA-based markers closely linked to the nematode resistance trait for marker-assisted selection (MAS) could enhance breeding efficiency and accuracy. In the present study, a pool of phenotyped progenies segregating for nematode resistance and susceptibility were fingerprinted using a 21 K single-nucleotide polymorphism (SNP) array. Eight candidate SNPs located on potato Chromosome 11, segregating with the nematode resistance trait, were identified and used as landmarks for discovery of other marker types such as simple sequence repeat (SSR) and insertion-deletion (INDEL) markers. Subsequently, a total of eight SNPs, 30 SSRs, and four INDELS located on scaffold 11 of Solanum bulbocastanum were used to design primers; markers were validated in a panel of resistant and susceptible clones. Two SNPs (SB_MC1Chr11-PotVar0066518 and SB_MC1Chr11-PotVar0064140), five SSRs (SB_MC1Chr11-SSR04, SB_MC1Chr11-SSR08, SB_MC1Chr11-SSR10, SB_MC1Chr11-SSR13, and SB_MC1Chr11-SSR20), and one INDEL (SB_MC1Chr11-INDEL4) markers showed polymorphism between the resistant and susceptible clones in the test panel and in other segregating progenies. These markers are robust, highly reproducible, and easy to use for MAS of nematode-resistant potato clones to enhance the breeding program.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01285-w.

QTL mapping for low temperature germination in rapeseed

Rapeseed, a major oil crop in the world, is easily affected by low-temperature stress. A low temperature delays seed germination and increases seedling mortality, adversely affecting rapeseed growth and production. In the present study, a tolerant cultivar (Huyou21) was crossed with a susceptible genotype (3429) to develop a mapping population consisting of 574 F₂ progenies and elucidate the genetic mechanisms of seed germination under low temperatures. Two quantitative trait loci (QTL) for low-temperature germination (LTG) were detected, one on chromosome A09 (named qLTGA9-1) and the other on chromosome C01 (named qLTGC1-1), using the QTL-seq approach and confirmed via linkage analysis in the mapping population. Further, qLTGA9-1 was mapped to a 341.86 kb interval between the SSR markers Nys9A212 and Nys9A215. In this region, 69 genes including six specific genes with moderate or high effect function variants were identified based on the Ningyou7 genome sequence. Meanwhile, qLTGC1-1 was mapped onto a 1.31 Mb interval between SSR markers Nys1C96 and Nys1C117. In this region, 133 genes including five specific genes with moderate effect function variants were identified. These specific genes within the two QTL could be used for further studies on cold tolerance and as targets in rapeseed breeding programs.

An empirical analysis of mtSSRs: could microsatellite distribution patterns explain the evolution of mitogenomes in plants?

Microsatellites (SSRs) are tandem repeat sequences in eukaryote genomes, including plant cytoplasmic genomes. The mitochondrial genome (mtDNA) has been shown to vary in size, number, and distribution of SSRs among different plant groups. Thus, SSRs contribute with genomic diversity in mtDNAs. However, the abundance, distribution, and evolutionary significance of SSRs in mtDNA from a wide range of algae and plants have not been explored. In this study, the mtDNAs of 204 plant and algal species were investigated related to the presence of SSRs. The number of SSRs was positively correlated with genome size. Its distribution is dependent on plant and algal groups analyzed, although the cluster analysis indicates the conservation of some common motifs in algal and terrestrial plants that reflect common ancestry of groups. Many SSRs in coding and non-coding regions can be useful for molecular markers. Moreover, mitochondrial SSRs are highly abundant, representing an important source for natural or induced genetic variation, i.e., for biotechnological approaches that can modulate mtDNA gene regulation. Thus, this comparative study increases the understanding of the plant and algal SSR evolution and brings perspectives for further studies.

Genetic Mapping of ms1s, a Recessive Gene for Male Sterility in Common Wheat

The utilization of heterosis is an important way to improve wheat yield, and the production of wheat hybrid seeds mainly relies on male-sterile lines. Male sterility in line 15 Fan 03 derived from a cross of 72,180 and Xiaoyan 6 is controlled by a single recessive gene. The gene was mapped to the distal region of chromosome 4BS in a genetic interval of 1.4 cM and physical distance of 6.57 Mb between SSR markers Ms4BS42 and Ms4BS199 using an F₂ population with 1205 individuals. Sterile individuals had a deletion of 4.57 Mb in the region presumed to carry the Ms1 locus. The allele for sterility was therefore named ms1s. Three CAPS markers were developed and verified from the region upstream of the deleted fragment and can be used for ms1s marker-assisted selection in wheat hybrid breeding. This work will enrich the utilization of male sterility genetic resources.

Genome-wide analysis of Hsp70 and Hsp100 gene families in Ziziphus jujuba

The Ziziphus species are naturally tolerant to a range of abiotic stresses. Therefore, it is expected that they are an enriched source of genes conferring stress tolerance. Heat shock proteins (Hsps) play a significant role in plants in imparting tolerance against abiotic stress conditions. To get an insight into potential Hsp function in Ziziphus, we performed a genome-wide analysis and expression study of Hsp70 and Hsp100 gene families in Ziziphus jujuba. We identified 21 and 6 genes of the ZjHsp70 and ZjHsp100 families, respectively. Physiochemical properties, chromosomal location, gene structure, motifs, and protein domain organization were analysed for structural and functional characterization. We identified the contribution of tandem and segmental gene duplications in expansions of ZjHsp70s and ZjHsp100s in Z. jujuba. Promoter analysis suggested that ZjHsp70s and ZjHsp100s perform diverse functions related to abiotic stress. Furthermore, expression analyses revealed that most of the Z. jujuba Hsp genes are differentially expressed in response to heat, drought, and salinity stress. Our analyses suggested ZjHsp70-3, ZjHsp70-5, ZjHsp70-6, ZjHsp70-16, ZjHsp70-17, ZjHsp70-20, ZjHsp100-1, ZjHsp100-2, and ZjHsp100-3 are potential candidates for further functional analysis and with regard to breeding new more resilient strains. The present analysis laid the foundation for understanding the molecular mechanism of Hsps70 and Hsp100 gene families regulating abiotic stress tolerance in Z. jujuba.

Patterns of genetic structure and evidence of Egyptian Citrus rootstock based on informative SSR, LTR-IRAP and LTR-REMAP molecular markers

Background

Releasing the draft genome of sweet orange provides useful information on genetic structure and molecular marker association with heritable breeding traits in citrus species and their structures. Last decades, microsatellite and retrotransposons are well known as a significant diverse component of the structural evolution. They represented the most potent elements for assessing sustainable utilization of the complicated classification in citrus breeding. Our study was performed to verify the structure analysis and the parentage genetic diversity among the Egyptian citrus rootstocks and the related species.

Results

Here, the performance of 26 SSR and 14 LTR-IRAP in addition to 20 LTR-REMAP markers have been used to conduct the discriminating power and the status of the genetic structure analysis among twenty specimens of citrus genotypes. As a result, the three markers approach exhibited a remarkable variation among the tested genotypes. Overall, the three markers have different discrimination power; the co-dominant SSR markers can differentiate within the group level only in addition to the species level of sour orange, while the dominant markers LTR-IRAP had the ability to discriminate among the group level in addition to species level and the origin of acids. Similarly, LTR-REMAP is suitable for classifying the group level and species level for mandarins as well the origin of Egyptian acids; probably due to it is integration of SSR and LTR-IRAP techniques. Structure and PCoA results of LTR-REMAP marker in strong support for the group structure of citrus species have been divided into four sets: acids, grapefruit/pummelo, mandarin/orange, and sour orange.

Conclusion

Our findings of the genetic structure analysis support the monophyletic nature of the citrus species; are able to provide unambiguous identification and disposition of true species and related hybrids like lemon, lime, citron, sour orange, grapefruit, mandarin, sweet orange, pummelo, and fortunella; and resulted in their placement in individual or overlap groups. The outcomes of these results will offer helpful and potential information for breeding programs and conservation approaches as a key stage toward identifying the interspecific admixture and the inferred structure origins of Egyptian citrus rootstock and acid cultivars.

A wheat cysteine-rich receptor-like kinase confers broad-spectrum resistance against Septoria tritici blotch

The poverty of disease resistance gene reservoirs limits the breeding of crops for durable resistance against evolutionary dynamic pathogens. Zymoseptoria tritici which causes Septoria tritici blotch (STB), represents one of the most genetically diverse and devastating wheat pathogens worldwide. No fully virulent Z. tritici isolates against synthetic wheats carrying the major resistant gene Stb16q have been identified. Here, we use comparative genomics, mutagenesis and complementation to identify Stb16q, which confers broad-spectrum resistance against Z. tritici. The Stb16q gene encodes a plasma membrane cysteine-rich receptor-like kinase that was recently introduced into cultivated wheat and which considerably slows penetration and intercellular growth of the pathogen.

Genome-Wide Analysis of Heat Shock Transcription Factors in Ziziphus jujuba Identifies Potential Candidates for Crop Improvement Under Abiotic Stress

Plant heat shock transcription factors (Hsfs) play a significant role in adoption under abiotic stress conditions by modulating the expression of several stress-responsive genes. Analysis of the Hsf gene family will serve to understand the molecular mechanism which is involved in response to abiotic stress. The Ziziphus species grows in warm and dry regions and is inherently tolerant to abiotic stress conditions; thus, Ziziphus is a highly enriched source of genes conferring abiotic stress tolerance. Therefore, the present study provides a comprehensive genome-wide analysis of the Hsf gene family in Z. jujuba. Identified 21 non-redundant Hsf genes were grouped into three major classes (classes A, B, and C) based on the phylogenetic analysis. Promoter and gene ontology analysis suggested that ZjHsfs perform diverse functions in response to abiotic stress conditions. Two paralogous pairs resulting from tandem gene duplication events were identified. Also, physio-chemical properties of chromosomal locations, gene structure, motifs, and protein domain organization of Hsfs were analyzed. Real-time PCR expression analyses revealed that most of the Z. jujuba Hsf genes are differentially expressed in response to heat stress. The analysis suggested ZjHsf-2, ZjHsf-3, ZjHsf-5, ZjHsf-7, ZjHsf-8, ZjHsf-10, ZjHsf-12, ZjHsf-17, and ZjHsf-18 were the outstanding candidate genes for imparting heat stress tolerance and for future functional analysis. The present analysis laid the foundation for understanding the molecular mechanism of the Hsf gene family regulating Z. jujuba development and tolerance to abiotic stress conditions.

BSA-Seq Discovery and Functional Analysis of Candidate Hessian Fly (Mayetiola destructor) Avirulence Genes

The Hessian fly (HF, Mayetiola destructor) is a plant-galling parasite of wheat (Triticum spp.). Seven percent of its genome is composed of highly diversified signal-peptide-encoding genes that are transcribed in HF larval salivary glands. These observations suggest that they encode effector proteins that are injected into wheat cells to suppress basal wheat immunity and redirect wheat development towards gall formation. Genetic mapping has determined that mutations in four of these genes are associated with HF larval survival (virulence) on plants carrying four different resistance (R) genes. Here, this line of investigation was pursued further using bulked-segregant analysis combined with whole genome resequencing (BSA-seq). Virulence to wheat R genes H6, Hdic, and H5 was examined. Mutations associated with H6 virulence had been mapped previously. Therefore, we used H6 to test the capacity of BSA-seq to map virulence using a field-derived HF population. This was the first time a non-structured HF population had been used to map HF virulence. Hdic virulence had not been mapped previously. Using a structured laboratory population, BSA-seq associated Hdic virulence with mutations in two candidate effector-encoding genes. Using a laboratory population, H5 virulence was previously positioned in a region spanning the centromere of HF autosome 2. BSA-seq resolved H5 virulence to a 1.3 Mb fragment on the same chromosome but failed to identify candidate mutations. Map-based candidate effectors were then delivered to Nicotiana plant cells via the type III secretion system of Burkholderia glumae bacteria. These experiments demonstrated that the genes associated with virulence to wheat R genes H6 and H13 are capable of suppressing plant immunity. Results are consistent with the hypothesis that effector proteins underlie the ability of HFs to survive on wheat.

Genetic diversity analysis for narrow-leafed lupin (Lupinus angustifolius L.) by SSR markers

Narrow-leafed lupin (Lupinus angustifolius L.) is used as grain legumes, fodder for livestock and green manure in the world and has a great potential to be developed as a new crop in China. In this study, we assessed the genetic diversity among a set of 109 newly introduced accessions of narrow-leafed lupin using 76 genomic SSR markers. Data analysis suggested that the average gene diversity index and average polymorphism information content (PIC) were 0.4758 and 0.4328, respectively. The mean allele number per loci (Na) was 6.3816. The population structure analysis identified two subgroups based on delta K (ΔK) values. This result is in accordance with that of a PCA. The AMOVA analysis showed that most of molecular variance were within population. These results will be useful to guide the genetic improvement of the narrow-leafed lupin crop in China.

Identification and characterization of SSR markers of Guadua chacoensis (Rojas) Londoño & P.M. Peterson and transferability to other bamboo species

The aim of this study was to develop simple sequence repeat (SSR) markers for genetic studies on G. chacoensis, as well as to evaluate their transferability to other bamboo species. Genomic DNA was isolated from G. chacoensis and its partial sequencing was used to find SSR loci. The obtained sequencing data were de novo assembled using the software CLC Genomics Workbench^® 8.0v. The SSR loci primers were identified and designed with the software SSRLocator. The selected markers were validated using 56 plants sampled in seven populations from southern Brazil. The markers with potential polymorphism were selected and fluorescently labeled for characterization by capillary electrophoresis. In total, 92 SSR loci were found in G. chacoensis contigs. Suitable primers were designed for 70 SSR loci, and the remaining 22 SSR loci did not have sequences for primer development. Out of 35 selected SSR markers, after PCR optimization, 10 with high polymorphism potential were characterized. These loci can be used in genetic analyses of G. chacoensis and all of them were successfully transferred to other bamboo species. Non-polymorphic loci require further tests with additional plants, from different populations, to identify possibilities of their use.

SSRome: an integrated database and pipelines for exploring microsatellites in all organisms

Over the past decade, many databases focusing on microsatellite mining on a genomic scale were released online with at least one of the following major deficiencies: (i) lacking the classification of microsatellites as genic or non-genic, (ii) not comparing microsatellite motifs at both genic and non-genic levels in order to identify unique motifs for each class or (iii) missing SSR marker development. In this study, we have developed ‘SSRome’ as a web-based, user-friendly, comprehensive and dynamic database with pipelines for exploring microsatellites in 6533 organisms. In the SSRome database, 158 million microsatellite motifs are identified across all taxa, in addition to all the mitochondrial and chloroplast genomes and expressed sequence tags available from NCBI. Moreover, 45.1 million microsatellite markers were developed and classified as genic or non-genic. All the stored motif and marker datasets can be downloaded freely. In addition, SSRome provides three user-friendly tools to identify, classify and compare motifs on either a genome- or transcriptome-wide scale. With the implementation of PHP, HTML and JavaScript, users can upload their data for analysis via a user-friendly GUI. SSRome represents a powerful database and mega-tool that will assist researchers in developing and dissecting microsatellite markers on a high-throughput scale.

De novo Transcriptome Assembly, Gene Annotation and SSR Marker Development in the Moon Seed Genus Menispermum (Menispermaceae)

The moonseed genus Menispermum L. (Menispermaceae) is disjunctly distributed in East Asia and eastern North America. Although Menispermum has important medicinal value, genetic and genomic information is scarce, with very few available molecular markers. In the current study, we used Illumina transcriptome sequencing and de novo assembly of the two Menispermum species to obtain in-depth genetic knowledge. From de novo assembly, 53,712 and 78,921 unigenes were generated for M. canadense and M. dauricum, with 37,527 (69.87%) and 55,211 (69.96%) showing significant similarities against the six functional databases, respectively. Moreover, 521 polymorphic EST-SSRs were identified. Of them, 23 polymorphic EST-SSR markers were selected to investigate the population genetic diversity within the genus. The newly developed EST-SSR markers also revealed high transferability among the three examined Menispermaceae species. Overall, we provide the very first transcriptomic analyses of this important medicinal genus. In addition, the novel microsatellite markers developed here will aid future studies on the population genetics and phylogeographic patterns of Menispermum at the intercontinental geographical scale.

Genome-wide mining and comparative analysis of fatty acid elongase gene family in Brassica napus and its progenitors

Very long chain fatty acids (VLCFAs) that are structural components of cell membrane lipid, cuticular waxes and seed oil, play crucial roles in plant growth, development and stress response. Fatty acid elongases (FAEs) comprising KCS and ELO, are key enzymes for VLCFA biosynthesis in plants. Although reference genomes of Brassica napus and its parental speices both have been sequenced, whole-genome analysis of FAE gene family in these Brassica speices is not reported. Here, 58, 33 and 30 KCS genes were identified in B. napus, B. rapa and B. oleracea genomes, respectively, whereas 14, 6 and 8 members were obtained for ELO genes. These KCS genes were unevenly located in 37 chromosomes and 3 scaffolds of 3 Brassica species, while these ELO genes were mapped to 19 chromosomes. The KCS and ELO proteins were divided into 8 and 4 subclasses, respectively. Gene structure and protein motifs remained highly conserved in each KCS or ELO subclass. Most promoters of KCS and ELO genes harbored various plant growth-, phytohormone-, and stress response-related cis-acting elements. 20 SSR loci existed in the KCS and ELO genes/promoters. The whole-genome duplication and segmental duplication mainly contributed to expansion of KCS and ELO genes in these genomes. Transcriptome analysis showed that KCS and ELO genes in 3 Brassica species were expressed in various tissues/organs with different levels, whereas 1 BnELO gene and 6 BnKCS genes might be pathogen-responsive genes. The qRT-PCR assay showed that BnKCS22 and BnELO04 responded to various phytohormone treatments and abiotic stresses. This work lays the foundation for further function identification of KCS and ELO genes in B. napus and its progenitors.

Natural variations in the promoter of Awn Length Inhibitor 1 (ALI-1) are associated with awn elongation and grain length in common wheat

Wheat awn plays a vital role in photosynthesis, grain production, and drought tolerance. However, the systematic identification or cloning of genes controlling wheat awn development is seldom reported. Here, we conducted a genome-wide association study (GWAS) with 364 wheat accessions and identified 26 loci involved in awn length development, including previously characterized B1, B2, Hd, and several rice homologs. The dominant awn suppressor B1 was fine mapped to a 125-kb physical interval, and a C₂ H₂ zinc finger protein Awn Length Inhibitor 1 (ALI-1) was confirmed to be the underlying gene of the B1 locus through the functional complimentary test with native awnless allele. ALI-1 expresses predominantly in the developing spike of awnless individuals, transcriptionally suppressing downstream genes. ALI-1 reduces cytokinin content and simultaneously restrains cytokinin signal transduction, leading to a stagnation of cell proliferation and reduction of cell numbers during awn development. Polymorphisms of four single nucleotide polymorphisms (SNPs) located in ALI-1 promoter region are diagnostic for the B1/b1 genotypes, and these SNPs are associated with awn length (AL), grain length (GL) and thousand-grain weight (TGW). More importantly, ali-1 was observed to increase grain length in wheat, which is a valuable attribute of awn on grain weight, aside from photosynthesis. Therefore, ALI-1 pleiotropically regulates awn and grain development, providing an alternative for grain yield improvement and addressing future climate changes.

A Novel Software and Method for the Efficient Development of Polymorphic SSR Loci Based on Transcriptome Data

Traditional methods for developing polymorphic microsatellite loci without reference sequences are time-consuming and labor-intensive, and the polymorphisms of simple sequence repeat (SSR) loci developed from expressed sequence tag (EST) databases are generally poor. To address this issue, in this study, we developed a new software (PSSRdt) and established an effective method for directly obtaining polymorphism details of SSR loci by analyzing diverse transcriptome data. The new method includes three steps, raw data processing, PSSRdt application, and loci extraction and verification. To test the practicality of the method, we successfully obtained 1940 potential polymorphic SSRs from the transcript dataset combined with 44 pea aphid transcriptomes. Fifty-two SSR loci obtained by the new method were selected for validating the polymorphic characteristics by genotyping in pea aphid individuals. The results showed that over 92% of SSR loci were polymorphic and 73.1% of loci were highly polymorphic. Our new software and method provide an innovative approach to microsatellite development based on RNA-seq data, and open a new path for the rapid mining of numerous loci with polymorphism to add to the body of research on microsatellites.

Evaluating the genetic structure of wild and commercial red cusk-eel (Genypterus chilensis) populations through the development of novel microsatellite markers from a reference transcriptome

The red cusk-eel (Genypterus chilensis) is a native Chilean species with a high-value market, with the potential to diversify Chilean aquaculture. The objective of this study was to develop a set of microsatellite markers, estimate genetic parameters, determine population differentiation, and identify the population structure of wild and commercial populations of G. chilensis. We discovered 6427 microsatellites markers from RNA-seq data, of which 54.9%, 20.2% and 16.8% were di-, tri-, and tetranucleotides, respectively. We used 12 of these markers to genotype two sets of broodstock, one group from commercial fish, and one group from wild fish from the Coquimbo Region of G. chilensis. We estimate the genetic parameters of the markers, selecting ten polymorphic markers (PIC > 0.5). We observed differences in the inbreeding coefficient among populations, with high values of inbreeding in one broodstock set and lower values in the other groups. The evaluation of population differentiation using Fst showed small (0.0195) to large (0.1888) genetic differentiation between the groups. The structure analysis showed that commercial and wild groups were formed by three clusters, without relevant evidence of admixture process, suggesting that groups evaluated in this study are formed of at least three subpopulations of G. chilensis, which could be explained by the low or lack of migration suggested for this species. This is the first study that identifies a high number of molecular markers in G. chilensis, providing relevant information of the genetic structure of commercial and wild population of this species.

Population Genetic Analysis of Phytophthora parasitica From Tobacco in Chongqing, Southwestern China

Tobacco black shank, caused by Phytophthora parasitica, is one of the most notorious tobacco diseases and causes huge economic losses worldwide. Understanding the genetic variation of P. parasitica populations is essential to the development of disease control measures. In this research, 210 simple sequence repeat (SSR) markers for P. parasitica were identified, 10 of which were polymorphic among nine reference strains. We further performed population genetic analysis of 245 P. parasitica isolates randomly collected from tobacco fields in Chongqing for mating type, molecular variation at 14 SSR loci (four of which were identified previously), and sensitivity to the fungicide metalaxyl. The results showed that the A2 mating type was dominant and no A1 mating type isolate was discovered. SSR genotyping distinguished 245 P. parasitica isolates into 46 genotypes, four of which were dominant in the population. Low genotypic diversity and excess heterozygosity were common in nearly all of the populations from Chongqing. Population analysis showed that no differentiation existed among different populations. All isolates tested were highly sensitive to metalaxyl. Taken together, our results showed that the P. parasitica populations from tobacco fields in Chongqing belonged to a clonal lineage and were highly sensitive to metalaxyl.

IDSSR: An Efficient Pipeline for Identifying Polymorphic Microsatellites from a Single Genome Sequence

Simple sequence repeats (SSRs) are known as microsatellites, and consist of tandem 1–6-base motifs. They have become one of the most popular molecular markers, and are widely used in molecular ecology, conservation biology, molecular breeding, and many other fields. Previously reported methods identify monomorphic and polymorphic SSRs and determine the polymorphic SSRs via experimental validation, which is potentially time-consuming and costly. Herein, we present a new strategy named insertion/deletion (INDEL) SSR (IDSSR) to identify polymorphic SSRs by integrating SSRs with nucleotide insertions/deletions (INDEL) solely based on a single genome sequence and the sequenced pair-end reads. These INDEL indexes and polymorphic SSRs were identified, as well as the number of repeats, repeat motifs, chromosome location, annealing temperature, and primer sequences, enabling future experimental approaches to determine the correctness and polymorphism. Experimental validation with the giant panda demonstrated that our method has high reliability and stability. The efficient SSR pipeline would help researchers obtain high-quality genetic markers for plants and animals of interest, save labor, and reduce costly marker-screening experiments. IDSSR is freely available at https://github.com/Allsummerking/IDSSR.

A high-quality genome of Eragrostis curvula grass provides insights into Poaceae evolution and supports new strategies to enhance forage quality

The Poaceae constitute a taxon of flowering plants (grasses) that cover almost all Earth’s inhabitable range and comprises some of the genera most commonly used for human and animal nutrition. Many of these crops have been sequenced, like rice, Brachypodium, maize and, more recently, wheat. Some important members are still considered orphan crops, lacking a sequenced genome, but having important traits that make them attractive for sequencing. Among these traits is apomixis, clonal reproduction by seeds, present in some members of the Poaceae like Eragrostis curvula. A de novo, high-quality genome assembly and annotation for E. curvula have been obtained by sequencing 602 Mb of a diploid genotype using a strategy that combined long-read length sequencing with chromosome conformation capture. The scaffold N50 for this assembly was 43.41 Mb and the annotation yielded 56,469 genes. The availability of this genome assembly has allowed us to identify regions associated with forage quality and to develop strategies to sequence and assemble the complex tetraploid genotypes which harbor the apomixis control region(s). Understanding and subsequently manipulating the genetic drivers underlying apomixis could revolutionize agriculture.

ImtRDB: a database and software for mitochondrial imperfect interspersed repeats annotation

BACKGROUND

Mitochondria is a powerhouse of all eukaryotic cells that have its own circular DNA (mtDNA) encoding various RNAs and proteins. Somatic perturbations of mtDNA are accumulating with age thus it is of great importance to uncover the main sources of mtDNA instability. Recent analyses demonstrated that somatic mtDNA deletions depend on imperfect repeats of various nature between distant mtDNA segments. However, till now there are no comprehensive databases annotating all types of imperfect repeats in numerous species with sequenced complete mitochondrial genome as well as there are no algorithms capable to call all types of imperfect repeats in circular mtDNA.

RESULTS

We implemented naïve algorithm of pattern recognition by analogy to standard dot-plot construction procedures allowing us to find both perfect and imperfect repeats of four main types: direct, inverted, mirror and complementary. Our algorithm is adapted to specific characteristics of mtDNA such as circularity and an excess of short repeats - it calls imperfect repeats starting from the length of 10 b.p. We constructed interactive web available database ImtRDB depositing perfect and imperfect repeats positions in mtDNAs of more than 3500 Vertebrate species. Additional tools, such as visualization of repeats within a genome, comparison of repeat densities among different genomes and a possibility to download all results make this database useful for many biologists. Our first analyses of the database demonstrated that mtDNA imperfect repeats (i) are usually short; (ii) associated with unfolded DNA structures; (iii) four types of repeats positively correlate with each other forming two equivalent pairs: direct and mirror versus inverted and complementary, with identical nucleotide content and similar distribution between species; (iv) abundance of repeats is negatively associated with GC content; (v) dinucleotides GC versus CG are overrepresented on light chain of mtDNA covered by repeats.

CONCLUSIONS

ImtRDB is available at http://bioinfodbs.kantiana.ru/ImtRDB/ . It is accompanied by the software calling all types of interspersed repeats with different level of degeneracy in circular DNA. This database and software can become a very useful tool in various areas of mitochondrial and chloroplast DNA research.

Discovery and characterization of SSR markers in Eugenia uniflora L. (Myrtaceae) using low coverage genome sequencing

Eugenia uniflora L. (Myrtaceae) is a tree species widely distributed in South America suffering the effects of the exploitation of natural populations. In this study, we employed low coverage sequencing of the E. uniflora genome for mining of SSR markers. The de novo assembly generated 2,601 contigs with an average length of 1139 bp and spans 3.15 Mb. A total of 76 dimer, 33 trimer and two compound SSR loci were identified. Twelve selected SSR loci were employed to genotype 30 individuals from two natural populations. A total of 73 alleles were detected (mean A= 6.1) were observed, the mean effective number of alleles was Ae = 3.91, mean Ho was 0.23 and mean HE was 0.70). The mean Wright's within population fixation index was FIS = 0.66 and significant deviation of HWE was observed in all loci, except one. The FST between populations equaled 0.27. The levels of genetic diversity and structure estimated with these 12 SSR markers are in accordance with data from genetics studies performed on other tree species of the Pampa biome, presenting moderate to high polymorphism and may be employed in studies of species conservation measures and breeding programs.

Improved genetic resolution for linkage mapping of resistance to potato wart in monoparental dihaploids with potential diagnostic value in tetraploid potato varieties

KEY MESSAGE

We achieved improved mapping resolution of the major wart resistance locus Xla-TNL containing also Sen1 in a dihaploid population using SNP data and developed additional markers with diagnostic value in tetraploid varieties. We analyzed a segregating monoparental dihaploid potato population comprising 215 genotypes derived from a tetraploid variety that is highly resistant to Synchytrium endobioticum pathotypes 18 and 6. The clear bimodal segregation for both pathotypes indicated that a major dominant resistance factor in a simplex allele configuration was present in the tetraploid donor genotype. Compared to that in previous analyses of the same tetraploid donor in conventional crosses with susceptible tetraploid genotypes, a segregation pattern with a reduced genetic complexity of resistance in dihaploids was observed here. Using the 12.8 k SolCAP SNP array, we mapped a resistance locus to the Xla-TNL region containing also Sen1 on potato chromosome 11. The improved mapping resolution provided by the monoparental dihaploids allowed for the localization of the genes responsible for the resistance to both pathotypes in an interval spanning less than 800 kbp on the reference genome. Furthermore, we identified eight molecular markers segregating without recombination to pathotype 18 and pathotype 6 resistance. Also, two developed markers display improved diagnostic properties in an independent panel of tetraploid varieties. Overall, our data provide the highest resolution mapping of wart resistance genes at the Xla-TNL locus thus far.

Genome-wide regulatory gene-derived SSRs reveal genetic differentiation and population structure in fiber flax genotypes

We designed a set of 580 simple sequence repeat markers; 506 from transcription factor-coding genes, and 74 from long non-coding RNAs and designated them as regulatory gene-derived simple sequence repeat (ReG-SSR) markers. From this set, we could anchor 559 ReG-SSR markers on 15 flax chromosomes with an average marker distance of 0.56 Mb. Thirty-one polymorphic ReG-SSR primers, amplifying SSR loci length of at least 20 bp were chosen from 134 screened primers. This primer set was used to characterize a diversity panel of 93 flax accessions. The panel included 33 accessions from India, including released varieties, dual-purpose lines and landraces, and 60 fiber flax accessions from the global core collection. Thirty-one ReG-SSR markers generated 76 alleles, with an average of 2.5 alleles per primer and a mean allele frequency of 0.77. These markers recorded 0.32 average gene diversity, 0.26 polymorphism information content and 1.35% null alleles. All the 31 ReG-SSR loci were found selectively neutral and showed no evidence of population reduction. A model-based clustering analysis separated the flax accessions into two sub-populations-Indian and global, with some accessions showing admixtures. The distinct clustering pattern of the Indian accessions compared to the global accessions, conforms to the principal coordinate analysis, genetic dissimilarity-based unweighted neighbor-joining tree and analysis of molecular variance. Fourteen flax accessions with 99.3% allelic richness were found optimum to adopt in breeding programs. In summary, the genome-wide ReG-SSR markers will serve as a functional marker resource for genetic and phenotypic relationship studies, marker-assisted selections, and provide a basis for selection of accessions from the Indian and global gene pool in fiber flax breeding programs.

Characterization of Plastidial and Nuclear SSR Markers for Understanding Invasion Histories and Genetic Diversity of Schinus molle L

Invasive plant species are expected to display high dispersal capacity but low levels of genetic diversity due to the founder effect occurring at each invasion episode. Understanding the history of invasions and the levels of genetic diversity of such species is an important task for planning management and monitoring strategy for these events. Peruvian Peppertree (Schinus molle L.) is a pioneer tree species native from South America which was introduced in North America, Europe and Africa, becoming a threat to these non-native habitats. In this study, we report the discovery and characterization of 17 plastidial (ptSSR) and seven nuclear (nSSR) markers for S. molle based on low-coverage whole-genome sequencing data acquired through next-generation sequencing. The markers were tested in 56 individuals from two natural populations sampled in the Brazilian Caatinga and Pampa biomes. All loci are moderately to highly polymorphic and revealed to be suitable for genetic monitoring of new invasions, for understanding the history of old invasions, as well as for genetic studies of native populations in their natural occurrence range and of orchards established with commercial purposes.

Red Clover (Trifolium pratense) and Zigzag Clover (T. medium) - A Picture of Genomic Similarities and Differences

The genus clover (Trifolium sp.) is one of the most economically important genera in the Fabaceae family. More than 10 species are grown as manure plants or forage legumes. Red clover's (T. pratense) genome size is one of the smallest in the Trifolium genus, while many clovers with potential breeding value have much larger genomes. Zigzag clover (T. medium) is closely related to the sequenced red clover; however, its genome is approximately 7.5x larger. Currently, almost nothing is known about the architecture of this large genome and differences between these two clover species. We sequenced the T. medium genome (2n = 8x = 64) with ∼23× coverage and managed to partially assemble 492.7 Mbp of its genomic sequence. A thorough comparison between red clover and zigzag clover sequencing reads resulted in the successful validation of 7 T. pratense- and 45 T. medium-specific repetitive elements. The newly discovered repeats led to the set-up of the first partial T. medium karyotype. Newly discovered red clover and zigzag clover tandem repeats were summarized. The structure of centromere-specific satellite repeat resembling that of T. repens was inferred in T. pratense. Two repeats, TrM300 and TrM378, showed a specific localization into centromeres of a half of all zigzag clover chromosomes; TrM300 on eight chromosomes and TrM378 on 24 chromosomes. A comparison with the red clover draft sequence was also used to mine more than 105,000 simple sequence repeats (SSRs) and 1,170,000 single nucleotide variants (SNVs). The presented data obtained from the sequencing of zigzag clover represent the first glimpse on the genomic sequence of this species. Centromeric repeats indicated its allopolyploid origin and naturally occurring homogenization of the centromeric repeat motif was somehow prevented. Using various repeats, highly uniform 64 chromosomes were separated into eight types of chromosomes. Zigzag clover genome underwent substantial chromosome rearrangements and cannot be counted as a true octoploid. The resulting data, especially the large number of predicted SSRs and SNVs, may have great potential for further research of the legume family and for rapid advancements in clover breeding.

Mutations in CsPID encoding a Ser/Thr protein kinase are responsible for round leaf shape in cucumber (Cucumis sativus L.)

Two round-leaf mutants, rl-1 and rl-2, were identified from EMS-induced mutagenesis. High throughput sequencing and map-based cloning suggested CsPID encoding a Ser/Thr protein kinase as the most possible candidate for rl-1. Rl-2 was allelic to Rl-1. Leaf shape is an important plant architecture trait that is affected by plant hormones, especially auxin. In Arabidopsis, PINOID (PID), a regulator for the auxin polar transporter PIN (PIN-FORMED) affects leaf shape formation, but this function of PID in crop plants has not been well studied. From an EMS mutagenesis population, we identified two round-leaf (rl) mutants, C356 and C949. Segregation analysis suggested that both mutations were controlled by single recessive genes, rl-1 and rl-2, respectively. With map-based cloning, we show that CsPID as the candidate gene of rl-1; a non-synonymous SNP in the second exon of CsPID resulted in an amino acid substitution and the round leaf phenotype. As compared in the wild type plant, CsPID had significantly lower expression in the root, leaf and female flowers in C356, which may result in the less developed roots, round leaves and abnormal female flowers, respectively in the rl-1 mutant. Among the three copies of PID genes, CsPID, CsPID2 and CSPID2L (CsPID2-like) in the cucumber genome, CsPID was the only one with significantly differential expression in adult leaves between WT and C356 suggesting CsPID plays a main role in leaf shape formation. The rl-2 mutation in C949 was also cloned, which was due to another SNP in a nearby location of rl-1 in the same CsPID gene. The two round leaf mutants and the work presented herein provide a good foundation for understanding the molecular mechanisms of CsPID in cucumber leaf development.

Characterization of Heterobasidion occidentale transcriptomes reveals candidate genes and DNA polymorphisms for virulence variations

Characterization of genes involved in differentiation of pathogen species and isolates with variations of virulence traits provides valuable information to control tree diseases for meeting the challenges of sustainable forest health and phytosanitary trade issues. Lack of genetic knowledge and genomic resources hinders novel gene discovery, molecular mechanism studies and development of diagnostic tools in the management of forest pathogens. Here, we report on transcriptome profiling of Heterobasidion occidentale isolates with contrasting virulence levels. Comparative transcriptomic analysis identified orthologous groups exclusive to H. occidentale and its isolates, revealing biological processes involved in the differentiation of isolates. Further bioinformatics analyses identified an H. occidentale secretome, CYPome and other candidate effectors, from which genes with species- and isolate-specific expression were characterized. A large proportion of differentially expressed genes were revealed to have putative activities as cell wall modification enzymes and transcription factors, suggesting their potential roles in virulence and fungal pathogenesis. Next, large numbers of simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) were detected, including more than 14 000 interisolate non-synonymous SNPs. These polymorphic loci and species/isolate-specific genes may contribute to virulence variations and provide ideal DNA markers for development of diagnostic tools and investigation of genetic diversity.

Chloroplast genome analyses and genomic resource development for epilithic sister genera Oresitrophe and Mukdenia (Saxifragaceae), using genome skimming data

Background

Epilithic sister genera Oresitrophe and Mukdenia (Saxifragaceae) have an epilithic habitat (rocky slopes) and a parapatric distribution in East Asia, which makes them an ideal model for a more comprehensive understanding of the demographic and divergence history and the influence of climate changes in East Asia. However, the genetic background and resources for these two genera are scarce.

Results

The complete chloroplast (cp) genomes of two Oresitrophe rupifraga and one Mukdenia rossii individuals were reconstructed and comparative analyses were conducted to examine the evolutionary pattern of chloroplast genomes in Saxifragaceae. The cp genomes ranged from 156,738 bp to 156,960 bp in length and had a typical quadripartite structure with a conserved genome arrangement. Comparative analysis revealed the intron of rpl2 has been lost in Heuchera parviflora, Tiarella polyphylla, M. rossii and O. rupifraga but presents in the reference genome of Penthorum chinense. Seven cp hotspot regions (trnH-psbA, trnR-atpA, atpI-rps2, rps2-rpoC2, petN-psbM, rps4-trnT and rpl33-rps18) were identified between Oresitrophe and Mukdenia, while four hotspots (trnQ-psbK, trnR-atpA, trnS-psbZ and rpl33-rps18) were identified within Oresitrophe. In addition, 24 polymorphic cpSSR loci were found between Oresitrophe and Mukdenia. Most importantly, we successfully developed 126 intergeneric polymorphic gSSR markers between Oresitrophe and Mukdenia, as well as 452 intrageneric ones within Oresitrophe. Twelve randomly selected intergeneric gSSRs have shown that these two genera exhibit a significant genetic structure.

Conclusions

In this study, we conducted genome skimming for Oresitrophe rupifraga and Mukdenia rossii. Using these data, we were able to not only assemble their complete chloroplast genomes, but also develop abundant genetic resources (cp hotspots, cpSSRs, polymorphic gSSRs). The genomic patterns and genetic resources presented here will contribute to further studies on population genetics, phylogeny and conservation biology in Saxifragaceae.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4633-x) contains supplementary material, which is available to authorized users.

Mining and Development of Novel SSR Markers Using Next Generation Sequencing (NGS) Data in Plants

Microsatellites, or simple sequence repeats (SSRs), are one of the most informative and multi-purpose genetic markers exploited in plant functional genomics. However, the discovery of SSRs and development using traditional methods are laborious, time-consuming, and costly. Recently, the availability of high-throughput sequencing technologies has enabled researchers to identify a substantial number of microsatellites at less cost and effort than traditional approaches. Illumina is a noteworthy transcriptome sequencing technology that is currently used in SSR marker development. Although 454 pyrosequencing datasets can be used for SSR development, this type of sequencing is no longer supported. This review aims to present an overview of the next generation sequencing, with a focus on the efficient use of de novo transcriptome sequencing (RNA-Seq) and related tools for mining and development of microsatellites in plants.

Quantitative Trait Locus Mapping and Candidate Gene Analysis for Verticillium Wilt Resistance Using Gossypium barbadense Chromosomal Segment Introgressed Line

Verticillium wilt (VW) is a soil-borne fungal disease that is caused by Verticillium dahliae Kleb and seriously damages cotton production annually in China. To date, many efforts have been made to improve the resistance of upland cotton against VW, but little progress has been achieved because of a lack of resistant upland cotton to VW. G. barbadense is known to carry high resistance to VW; however, it is difficult to transfer the resistance trait from G. barbadense to upland cotton because of linkage drag and distortion in the interspecific hybrid. In this study, a chromosomal segment introgression line (CSIL), SuVR043, containing a single and homozygous chromosome segment of G. barbadense cv. H7124 D04 (Chr 22), was created and used to construct an F₂ population for mapping of VW resistance quantitative trait loci (QTLs) in the greenhouse. Two major resistance QTLs against nondefoliating V. dahliae isolate Bp2, called qVW-Bp2-1 and qVW-Bp2-2, which were flanked by the markers cgr6409-ZHX37 and ZHX57-ZHX70 and explained an average of 16.38 and 22.36% of the observed phenotypic variation, respectively, were detected in three independent replicate experiments. The genetic distances from cgr6409 to ZHX37 and from ZHX57 to ZHX70 were 2.4 and 0.8 cM, respectively. By analyzing the genome sequence of the qVW-Bp2-1 and qVW-Bp2-2 regions, we determined that the accurate physical distances from cgr6409 to ZHX37 and from ZHX57 to ZHX70 in the G. barbadense genome are 254 and 140 kb, and that those spans 36 and 20 putative genes, respectively. The results of the expression analysis showed significant differences in the expression profiles of GbCYP450, GbTMEM214, and GbRLK among G. barbadense cv. H7124, CSIL SuVR043 and G. hirsutum acc. Sumian 8 at different times after inoculation with V. dahliae isolate Bp2. Virus-induced gene silencing (VIGS) analysis showed that silencing of GbCYP450 and GbTMEM214 decreased H7124 and CSIL SuVR043 resistance to VW. These results form a solid foundation for fine mapping and cloning of resistance genes in the substituted segment and will provide valuable assistance in future efforts to breed for VW resistance.

De novo transcriptome sequencing and assembly from apomictic and sexual Eragrostis curvula genotypes

A long-standing goal in plant breeding has been the ability to confer apomixis to agriculturally relevant species, which would require a deeper comprehension of the molecular basis of apomictic regulatory mechanisms. Eragrostis curvula (Schrad.) Nees is a perennial grass that includes both sexual and apomictic cytotypes. The availability of a reference transcriptome for this species would constitute a very important tool toward the identification of genes controlling key steps of the apomictic pathway. Here, we used Roche/454 sequencing technologies to generate reads from inflorescences of E. curvula apomictic and sexual genotypes that were de novo assembled into a reference transcriptome. Near 90% of the 49568 assembled isotigs showed sequence similarity to sequences deposited in the public databases. A gene ontology analysis categorized 27448 isotigs into at least one of the three main GO categories. We identified 11475 SSRs, and several of them were assayed in E curvula germoplasm using SSR-based primers, providing a valuable set of molecular markers that could allow direct allele selection. The differential contribution to each library of the spliced forms of several transcripts revealed the existence of several isotigs produced via alternative splicing of single genes. The reference transcriptome presented and validated in this work will be useful for the identification of a wide range of gene(s) related to agronomic traits of E. curvula, including those controlling key steps of the apomictic pathway in this species, allowing the extrapolation of the findings to other plant species.

Quantitative trait locus mapping under irrigated and drought treatments based on a novel genetic linkage map in mungbean (Vigna radiata L.)

A novel genetic linkage map was constructed using SSR markers and stable QTLs were identified for six drought tolerance related-traits using single-environment analysis under irrigation and drought treatments. Mungbean (Vigna radiata L.) is one of the most important leguminous food crops. However, mungbean production is seriously constrained by drought. Isolation of drought-responsive genetic elements and marker-assisted selection breeding will benefit from the detection of quantitative trait locus (QTLs) for traits related to drought tolerance. In this study, we developed a full-coverage genetic linkage map based on simple sequence repeat (SSR) markers using a recombinant inbred line (RIL) population derived from an intra-specific cross between two drought-resistant varieties. This novel map was anchored with 313 markers. The total map length was 1010.18 cM across 11 linkage groups, covering the entire genome of mungbean with a saturation of one marker every 3.23 cM. We subsequently detected 58 QTLs for plant height (PH), maximum leaf area (MLA), biomass (BM), relative water content, days to first flowering, and seed yield (Yield) and 5 for the drought tolerance index of 3 traits in irrigated and drought environments at 2 locations. Thirty-eight of these QTLs were consistently detected two or more times at similar linkage positions. Notably, qPH5A and qMLA2A were consistently identified in marker intervals from GMES5773 to MUS128 in LG05 and from Mchr11-34 to the HAAS_VR_1812 region in LG02 in four environments, contributing 6.40-20.06% and 6.97-7.94% of the observed phenotypic variation, respectively. None of these QTLs shared loci with previously identified drought-related loci from mungbean. The results of these analyses might facilitate the isolation of drought-related genes and help to clarify the mechanism of drought tolerance in mungbean.

A draft genome sequence of the rose black spot fungus Diplocarpon rosae reveals a high degree of genome duplication

BACKGROUND

Black spot is one of the most severe and damaging diseases of garden roses. We present the draft genome sequence of its causative agent Diplocarpon rosae as a working tool to generate molecular markers and to analyze functional and structural characteristics of this fungus.

RESULTS

The isolate DortE4 was sequenced with 191x coverage of different read types which were assembled into 2457 scaffolds. By evidence supported genome annotation with the MAKER pipeline 14,004 gene models were predicted and transcriptomic data indicated that 88.5% of them are expressed during the early stages of infection. Analyses of k-mer distributions resulted in unexpectedly large genome size estimations between 72.5 and 91.4 Mb, which cannot be attributed to its repeat structure and content of transposable elements alone, factors explaining such differences in other fungal genomes. In contrast, different lines of evidences demonstrate that a huge proportion (approximately 80%) of genes are duplicated, which might indicate a whole genome duplication event. By PCR-RFLP analysis of six paralogous gene pairs of BUSCO orthologs, which are expected to be single copy genes, we could show experimentally that the duplication is not due to technical error and that not all isolates tested possess all of the paralogs.

CONCLUSIONS

The presented genome sequence is still a fragmented draft but contains almost the complete gene space. Therefore, it provides a useful working tool to study the interaction of D. rosae with the host and the influence of a genome duplication outside of the model yeast in the background of a phytopathogen.

A genome-wide transcriptome map of pistachio (Pistacia vera L.) provides novel insights into salinity-related genes and marker discovery

Background

Pistachio (Pistacia vera L.) is one of the most important commercial nut crops worldwide. It is a salt-tolerant and long-lived tree, with the largest cultivation area in Iran. Climate change and subsequent increased soil salt content have adversely affected the pistachio yield in recent years. However, the lack of genomic/global transcriptomic sequences on P. vera impedes comprehensive researches at the molecular level. Hence, whole transcriptome sequencing is required to gain insight into functional genes and pathways in response to salt stress.

Results

RNA sequencing of a pooled sample representing 24 different tissues of two pistachio cultivars with contrasting salinity tolerance under control and salt treatment by Illumina Hiseq 2000 platform resulted in 368,953,262 clean 100 bp paired-ends reads (90 Gb). Following creating several assemblies and assessing their quality from multiple perspectives, we found that using the annotation-based metrics together with the length-based parameters allows an improved assessment of the transcriptome assembly quality, compared to the solely use of the length-based parameters. The generated assembly by Trinity was adopted for functional annotation and subsequent analyses. In total, 29,119 contigs annotated against all of five public databases, including NR, UniProt, TAIR10, KOG and InterProScan. Among 279 KEGG pathways supported by our assembly, we further examined the pathways involved in the plant hormone biosynthesis and signaling as well as those to be contributed to secondary metabolite biosynthesis due to their importance under salinity stress. In total, 11,337 SSRs were also identified, which the most abundant being dinucleotide repeats. Besides, 13,097 transcripts as candidate stress-responsive genes were identified. Expression of some of these genes experimentally validated through quantitative real-time PCR (qRT-PCR) that further confirmed the accuracy of the assembly. From this analysis, the contrasting expression pattern of NCED3 and SOS1 genes were observed between salt-sensitive and salt-tolerant cultivars.

Conclusion

This study, as the first report on the whole transcriptome survey of P. vera, provides important resources and paves the way for functional and comparative genomic studies on this major tree to discover the salinity tolerance-related markers and stress response mechanisms for breeding of new pistachio cultivars with more salinity tolerance.

Electronic supplementary material

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

Genome-Wide Association Study Identifies NBS-LRR-Encoding Genes Related with Anthracnose and Common Bacterial Blight in the Common Bean

Nucleotide-binding site and leucine-rich repeat (NBS-LRR) genes represent the largest and most important disease resistance genes in plants. The genome sequence of the common bean (Phaseolus vulgaris L.) provides valuable data for determining the genomic organization of NBS-LRR genes. However, data on the NBS-LRR genes in the common bean are limited. In total, 178 NBS-LRR-type genes and 145 partial genes (with or without a NBS) located on 11 common bean chromosomes were identified from genome sequences database. Furthermore, 30 NBS-LRR genes were classified into Toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) types, and 148 NBS-LRR genes were classified into coiled-coil (CC)-NBS-LRR (CNL) types. Moreover, the phylogenetic tree supported the division of these PvNBS genes into two obvious groups, TNL types and CNL types. We also built expression profiles of NBS genes in response to anthracnose and common bacterial blight using qRT-PCR. Finally, we detected nine disease resistance loci for anthracnose (ANT) and seven for common bacterial blight (CBB) using the developed NBS-SSR markers. Among these loci, NSSR24, NSSR73, and NSSR265 may be located at new regions for ANT resistance, while NSSR65 and NSSR260 may be located at new regions for CBB resistance. Furthermore, we validated NSSR24, NSSR65, NSSR73, NSSR260, and NSSR265 using a new natural population. Our results provide useful information regarding the function of the NBS-LRR proteins and will accelerate the functional genomics and evolutionary studies of NBS-LRR genes in food legumes. NBS-SSR markers represent a wide-reaching resource for molecular breeding in the common bean and other food legumes. Collectively, our results should be of broad interest to bean scientists and breeders.