Transcriptome analysis reveals ginsenosides biosynthetic genes, microRNAs and simple sequence repeats in Panax ginseng C. A. Meyer

BACKGROUND

Panax ginseng C. A. Meyer is one of the most widely used medicinal plants. Complete genome information for this species remains unavailable due to its large genome size. At present, analysis of expressed sequence tags is still the most powerful tool for large-scale gene discovery. The global expressed sequence tags from P. ginseng tissues, especially those isolated from stems, leaves and flowers, are still limited, hindering in-depth study of P. ginseng.

RESULTS

Two 454 pyrosequencing runs generated a total of 2,423,076 reads from P. ginseng roots, stems, leaves and flowers. The high-quality reads from each of the tissues were independently assembled into separate and shared contigs. In the separately assembled database, 45,849, 6,172, 4,041 and 3,273 unigenes were only found in the roots, stems, leaves and flowers database, respectively. In the jointly assembled database, 178,145 unigenes were observed, including 86,609 contigs and 91,536 singletons. Among the 178,145 unigenes, 105,522 were identified for the first time, of which 65.6% were identified in the stem, leaf or flower cDNA libraries of P. ginseng. After annotation, we discovered 223 unigenes involved in ginsenoside backbone biosynthesis. Additionally, a total of 326 potential cytochrome P450 and 129 potential UDP-glycosyltransferase sequences were predicted based on the annotation results, some of which may encode enzymes responsible for ginsenoside backbone modification. A BLAST search of the obtained high-quality reads identified 14 potential microRNAs in P. ginseng, which were estimated to target 100 protein-coding genes, including transcription factors, transporters and DNA binding proteins, among others. In addition, a total of 13,044 simple sequence repeats were identified from the 178,145 unigenes.

CONCLUSIONS

This study provides global expressed sequence tags for P. ginseng, which will contribute significantly to further genome-wide research and analyses in this species. The novel unigenes identified here enlarge the available P. ginseng gene pool and will facilitate gene discovery. In addition, the identification of microRNAs and the prediction of targets from this study will provide information on gene transcriptional regulation in P. ginseng. Finally, the analysis of simple sequence repeats will provide genetic makers for molecular breeding and genetic applications in this species.

Identification of drought-responsive universal stress proteins in viridiplantae

Genes encoding proteins that contain the universal stress protein (USP) domain are known to provide bacteria, archaea, fungi, protozoa, and plants with the ability to respond to a plethora of environmental stresses. Specifically in plants, drought tolerance is a desirable phenotype. However, limited focused and organized functional genomic datasets exist on drought-responsive plant USP genes to facilitate their characterization. The overall objective of the investigation was to identify diverse plant universal stress proteins and Expressed Sequence Tags (ESTs) responsive to water-deficit stress. We hypothesize that cross-database mining of functional annotations in protein and gene transcript bioinformatics resources would help identify candidate drought-responsive universal stress proteins and transcripts from multiple plant species. Our bioinformatics approach retrieved, mined and integrated comprehensive functional annotation data on 511 protein and 1561 ESTs sequences from 161 viridiplantae taxa. A total of 32 drought-responsive ESTs from 7 plant genera Glycine, Hordeum, Manihot, Medicago, Oryza, Pinus and Triticum were identified. Two Arabidopsis USP genes At3g62550 and At3g53990 that encode ATP-binding motif were up-regulated in a drought microarray dataset. Further, a dataset of 80 simple sequence repeats (SSRs) linked to 20 singletons and 47 transcript assembles was constructed. Integrating the datasets on SSRs and drought-responsive ESTs identified three drought-responsive ESTs from bread wheat (BE604157), soybean (BM887317) and maritime pine (BX682209). The SSR sequence types were CAG, ATA and AT respectively. The datasets from cross-database mining provide organized resources for the characterization of USP genes as useful targets for engineering plant varieties tolerant to unfavorable environmental conditions.

Comparative transcriptome resources of eleven Primulina species, a group of 'stone plants' from a biodiversity hot spot

The genus Primulina is an emerging model system in studying the drivers and mechanisms of species diversification, for its high species richness and endemism, together with high degree of habitat specialization. In this study, we sequenced transcriptomes for eleven Primulina species across the phylogeny of the genus using the Illumina HiSeq 2000 platform. A total of 336 million clean reads were processed into 355 573 unigenes with a mean length of 1336 bp and an N50 value of 2191 bp after pooling and reassembling twelve individual pre-assembled unigene sets. Of these unigenes, 249 973 (70%) were successfully annotated and 256 601 (72%) were identified as coding sequences (CDSs). We identified a total of 38 279 simple sequence repeats (SSRs) and 367 123 single nucleotide polymorphisms (SNPs). Marker validation assay revealed that 354 (27.3%) of the 1296 SSR and 795 (39.6%) of the 2008 SNP loci showed successful genotyping performance and exhibited expected polymorphism profiles. We screened 834 putative single-copy nuclear genes and proved their high effectiveness in phylogeny construction and estimation of ancestral population parameters. We identified a total of 85 candidate orthologs under positive selection for 46 of the 66 species pairs. This study provided an efficient application of RNA-seq in development of genomic resources for a group of 'stone plants' from south China Karst regions, a biodiversity hot spot of the World. The assembled unigenes with annotations and the massive gene-associated molecular markers would help guide further molecular systematic, population genetic and ecological genomics studies in Primulina and its relatives.

Microsatellite landscape evolutionary dynamics across 450 million years of vertebrate genome evolution

The evolutionary dynamics of simple sequence repeats (SSRs or microsatellites) across the vertebrate tree of life remain largely undocumented and poorly understood. In this study, we analyzed patterns of genomic microsatellite abundance and evolution across 71 vertebrate genomes. The highest abundances of microsatellites exist in the genomes of ray-finned fishes, squamate reptiles, and mammals, while crocodilian, turtle, and avian genomes exhibit reduced microsatellite landscapes. We used comparative methods to infer evolutionary rates of change in microsatellite abundance across vertebrates and to highlight particular lineages that have experienced unusually high or low rates of change in genomic microsatellite abundance. Overall, most variation in microsatellite content, abundance, and evolutionary rate is observed among major lineages of reptiles, yet we found that several deeply divergent clades (i.e., squamate reptiles and mammals) contained relatively similar genomic microsatellite compositions. Archosauromorph reptiles (turtles, crocodilians, and birds) exhibit reduced genomic microsatellite content and the slowest rates of microsatellite evolution, in contrast to squamate reptile genomes that have among the highest rates of microsatellite evolution. Substantial branch-specific shifts in SSR content in primates, monotremes, rodents, snakes, and fish are also evident. Collectively, our results support multiple major shifts in microsatellite genomic landscapes among vertebrates.

De novo Assembly of Leaf Transcriptome in the Medicinal Plant Andrographis paniculata

Andrographis paniculata is an important medicinal plant containing various bioactive terpenoids and flavonoids. Despite its importance in herbal medicine, no ready-to-use transcript sequence information of this plant is made available in the public data base, this study mainly deals with the sequencing of RNA from A. paniculata leaf using Illumina HiSeq™ 2000 platform followed by the de novo transcriptome assembly. A total of 189.22 million high quality paired reads were generated and 1,70,724 transcripts were predicted in the primary assembly. Secondary assembly generated a transcriptome size of ~88 Mb with 83,800 clustered transcripts. Based on the similarity searches against plant non-redundant protein database, gene ontology, and eukaryotic orthologous groups, 49,363 transcripts were annotated constituting upto 58.91% of the identified unigenes. Annotation of transcripts-using kyoto encyclopedia of genes and genomes database-revealed 5606 transcripts plausibly involved in 140 pathways including biosynthesis of terpenoids and other secondary metabolites. Transcription factor analysis showed 6767 unique transcripts belonging to 97 different transcription factor families. A total number of 124 CYP450 transcripts belonging to seven divergent clans have been identified. Transcriptome revealed 146 different transcripts coding for enzymes involved in the biosynthesis of terpenoids of which 35 contained terpene synthase motifs. This study also revealed 32,341 simple sequence repeats (SSRs) in 23,168 transcripts. Assembled sequences of transcriptome of A. paniculata generated in this study are made available, for the first time, in the TSA database, which provides useful information for functional and comparative genomic analysis besides identification of key enzymes involved in the various pathways of secondary metabolism.

RNA-Seq of Guar (Cyamopsis tetragonoloba, L. Taub.) Leaves: De novo Transcriptome Assembly, Functional Annotation and Development of Genomic Resources

Genetic improvement in industrially important guar (Cyamopsis tetragonoloba, L. Taub.) crop has been hindered due to the lack of sufficient genomic or transcriptomic resources. In this study, RNA-Seq technology was employed to characterize the transcriptome of leaf tissues from two guar varieties, namely, M-83 and RGC-1066. Approximately 30 million high-quality pair-end reads of each variety generated by Illumina HiSeq platform were used for de novo assembly by Trinity program. A total of 62,146 non-redundant unigenes with an average length of 679 bp were obtained. The quality assessment of assembled unigenes revealed 87.50% of complete and 97.18% partial core eukaryotic genes (CEGs). Sequence similarity analyses and annotation of the unigenes against non-redundant protein (Nr) and Gene Ontology (GO) databases identified 175,882 GO annotations. A total of 11,308 guar unigenes were annotated with various enzyme codes (EC) and categorized in six categories with 55 subclasses. The annotation of biochemical pathways resulted in a total of 11,971 unigenes assigned with 145 KEGG maps and 1759 enzyme codes. The species distribution analysis of the unigenes showed highest similarity with Glycine max genes. A total of 5773 potential simple sequence repeats (SSRs) and 3594 high-quality single nucleotide polymorphisms (SNPs) were identified. Out of 20 randomly selected SSRs for wet laboratory validation, 13 showed consistent PCR amplification in both guar varieties. In silico studies identified 145 polymorphic SSR markers in two varieties. To the best of our knowledge, this is the first report on transcriptome analysis and SNPs identification in guar till date.

A second-generation genetic linkage map for bighead carp (Aristichthys nobilis) based on microsatellite markers

Bighead carp (Aristichthys nobilis) is an important aquaculture fish worldwide. Genetic linkage maps for the species were previously reported, but map resolution remained to be improved. In this study, a second-generation genetic linkage map was constructed for bighead carp through a pseudo-testcross strategy using interspecific hybrids between bighead carp and silver carp. Of the 754 microsatellites genotyped in two interspecific mapping families (with 77 progenies for each family), 659 markers were assigned to 24 linkage groups, which were equal to the chromosome numbers of the haploid genome. The consensus map spanned 1917.3 cM covering 92.8% of the estimated bighead carp genome with an average marker interval of 2.9 cM. The length of linkage groups ranged from 52.2 to 133.5 cM with an average of 79.9 cM. The number of markers per linkage group varied from 11 to 55 with an average of 27.5 per linkage group. Normality tests on interval distances of the map showed a non-normal marker distribution; however, significant correlation was found between the length of linkage group and the number of markers below the 0.01 significance level (two-tailed). The length of the female map was 1.12 times that of the male map, and the average recombination ratio of female to male was 1.10:1. Visual inspection showed that distorted markers gathered in some linkage groups and in certain regions of the male and female maps. This well-defined genetic linkage map will provide a basic framework for further genome mapping of quantitative traits, comparative mapping and marker-assisted breeding in bighead carp.

A Repeat-Associated Small RNA Controls the Major Virulence Factors of Helicobacter pylori

Many bacterial pathogens regulate their virulence genes via phase variation, whereby length-variable simple sequence repeats control the transcription or coding potential of those genes. Here, we have exploited this relationship between DNA structure and physiological function to discover a globally acting small RNA (sRNA) regulator of virulence in the gastric pathogen Helicobacter pylori. Our study reports the first sRNA whose expression is affected by a variable thymine (T) stretch in its promoter. We show the sRNA post-transcriptionally represses multiple major pathogenicity factors of H. pylori, including CagA and VacA, by base pairing to their mRNAs. We further demonstrate transcription of the sRNA is regulated by the nickel-responsive transcriptional regulator NikR (thus named NikS for nickel-regulated sRNA), thereby linking virulence factor regulation to nickel concentrations. Using in-vitro infection experiments, we demonstrate NikS affects host cell internalization and epithelial barrier disruption. Together, our results show NikS is a phase-variable, post-transcriptional global regulator of virulence properties in H. pylori.

Gene expression profile analysis of Ligon lintless-1 (Li1) mutant reveals important genes and pathways in cotton leaf and fiber development

Ligon lintless-1 (Li1) is a monogenic dominant mutant of Gossypium hirsutum (upland cotton) with a phenotype of impaired vegetative growth and short lint fibers. Despite years of research involving genetic mapping and gene expression profile analysis of Li1 mutant ovule tissues, the gene remains uncloned and the underlying pathway of cotton fiber elongation is still unclear. In this study, we report the whole genome-level deep-sequencing analysis of leaf tissues of the Li1 mutant. Differentially expressed genes in leaf tissues of mutant versus wild-type (WT) plants are identified, and the underlying pathways and potential genes that control leaf and fiber development are inferred. The results show that transcription factors AS2, YABBY5, and KANDI-like are significantly differentially expressed in mutant tissues compared with WT ones. Interestingly, several fiber development-related genes are found in the downregulated gene list of the mutant leaf transcriptome. These genes include heat shock protein family, cytoskeleton arrangement, cell wall synthesis, energy, H2O2 metabolism-related genes, and WRKY transcription factors. This finding suggests that the genes are involved in leaf morphology determination and fiber elongation. The expression data are also compared with the previously published microarray data of Li1 ovule tissues. Comparative analysis of the ovule transcriptomes of Li1 and WT reveals that a number of pathways important for fiber elongation are enriched in the downregulated gene list at different fiber development stages (0, 6, 9, 12, 15, 18dpa). Differentially expressed genes identified in both leaf and fiber samples are aligned with cotton whole genome sequences and combined with the genetic fine mapping results to identify a list of candidate genes for Li1.

The Complete Plastome Sequences of Eleven Capsicum Genotypes: Insights into DNA Variation and Molecular Evolution

Members of the genus Capsicum are of great economic importance, including both wild forms and cultivars of peppers and chilies. The high number of potentially informative characteristics that can be identified through next-generation sequencing technologies gave a huge boost to evolutionary and comparative genomic research in higher plants. Here, we determined the complete nucleotide sequences of the plastomes of eight Capsicum species (eleven genotypes), representing the three main taxonomic groups in the genus and estimated molecular diversity. Comparative analyses highlighted a wide spectrum of variation, ranging from point mutations to small/medium size insertions/deletions (InDels), with accD, ndhB, rpl20, ycf1, and ycf2 being the most variable genes. The global pattern of sequence variation is consistent with the phylogenetic signal. Maximum-likelihood tree estimation revealed that Capsicum chacoense is sister to the baccatum complex. Divergence and positive selection analyses unveiled that protein-coding genes were generally well conserved, but we identified 25 positive signatures distributed in six genes involved in different essential plastid functions, suggesting positive selection during evolution of Capsicum plastomes. Finally, the identified sequence variation allowed us to develop simple PCR-based markers useful in future work to discriminate species belonging to different Capsicum complexes.

Molecular Genotyping (SSR) and Agronomic Phenotyping for Utilization of Durum Wheat (Triticum durum Desf.) Ex Situ Collection from Southern Italy: A Combined Approach Including Pedigreed Varieties

In South Italy durum wheat (Triticum durum Desf.) has a long-time tradition of growing and breeding. Accessions collected and now preserved ex situ are a valuable genetic resource, but their effective use in agriculture and breeding programs remains very low. In this study, a small number (44) of simple sequence repeats (SSR) molecular markers were used to detect pattern of diversity for 136 accessions collected in South Italy over time, to identify the genepool of origin, and establish similarities with 28 Italian varieties with known pedigree grown in Italy over the same time-period. Phenotyping was conducted for 12 morphophysiological characters of agronomic interest. Based on discriminant analysis of principal components (DAPC) and STRUCTURE analysis six groups were identified, the assignment of varieties reflected the genetic basis and breeding strategies involved in their development. Some “old” varieties grown today are the result of evolution through natural hybridization and conservative pure line selection. A small number of molecular markers and little phenotyping coupled with powerful statistical analysis and comparison to pedigreed varieties can provide enough information on the genetic structure of durum wheat germplasm for a quick screening of the germplasm collection able to identify accessions for breeding or introduction in low input agriculture.

Association between simple sequence repeat-rich chromosome regions and intergenomic translocation breakpoints in natural populations of allopolyploid wild wheats

BACKGROUND AND AIMS

Repetitive DNA sequences are thought to be involved in the formation of chromosomal rearrangements. The aim of this study was to analyse the distribution of microsatellite clusters in Aegilops biuncialis and Aegilops geniculata, and its relationship with the intergenomic translocations in these allotetraploid species, wild genetic resources for wheat improvement.

METHODS

The chromosomal localization of (ACG)(n) and (GAA)(n) microsatellite sequences in Ae. biuncialis and Ae. geniculata and in their diploid progenitors Aegilops comosa and Aegilops umbellulata was investigated by sequential in situ hybridization with simple sequence repeat (SSR) probes and repeated DNA probes (pSc119·2, Afa family and pTa71) and by dual-colour genomic in situ hybridization (GISH). Thirty-two Ae. biuncialis and 19 Ae. geniculata accessions were screened by GISH for intergenomic translocations, which were further characterized by fluorescence in situ hybridization and GISH.

KEY RESULTS

Single pericentromeric (ACG)(n) signals were localized on most U and on some M genome chromosomes, whereas strong pericentromeric and several intercalary and telomeric (GAA)(n) sites were observed on the Aegilops chromosomes. Three Ae. biuncialis accessions carried 7U(b)-7M(b) reciprocal translocations and one had a 7U(b)-1M(b) rearrangement, while two Ae. geniculata accessions carried 7U(g)-1M(g) or 5U(g)-5M(g) translocations. Conspicuous (ACG)(n) and/or (GAA)(n) clusters were located near the translocation breakpoints in eight of the ten translocated chromosomes analysed, SSR bands and breakpoints being statistically located at the same chromosomal site in six of them.

CONCLUSIONS

Intergenomic translocation breakpoints are frequently mapped to SSR-rich chromosomal regions in the allopolyploid species examined, suggesting that microsatellite repeated DNA sequences might facilitate the formation of those chromosomal rearrangements. The (ACG)(n) and (GAA)(n) SSR motifs serve as additional chromosome markers for the karyotypic analysis of UM genome Aegilops species.

MSDB: A Comprehensive Database of Simple Sequence Repeats

Microsatellites, also known as Simple Sequence Repeats (SSRs), are short tandem repeats of 1-6 nt motifs present in all genomes, particularly eukaryotes. Besides their usefulness as genome markers, SSRs have been shown to perform important regulatory functions, and variations in their length at coding regions are linked to several disorders in humans. Microsatellites show a taxon-specific enrichment in eukaryotic genomes, and some may be functional. MSDB (Microsatellite Database) is a collection of >650 million SSRs from 6,893 species including Bacteria, Archaea, Fungi, Plants, and Animals. This database is by far the most exhaustive resource to access and analyze SSR data of multiple species. In addition to exploring data in a customizable tabular format, users can view and compare the data of multiple species simultaneously using our interactive plotting system. MSDB is developed using the Django framework and MySQL. It is freely available at http://tdb.ccmb.res.in/msdb.

Discovery of miRNAs and Development of Heat-Responsive miRNA-SSR Markers for Characterization of Wheat Germplasm for Terminal Heat Tolerance Breeding

A large proportion of the Asian population fulfills their energy requirements from wheat (Triticum aestivum L.). Wheat quality and yield are critically affected by the terminal heat stress across the globe. It affects approximately 40% of the wheat-cultivating regions of the world. Therefore, there is a critical need to develop improved terminal heat-tolerant wheat varieties. Marker-assisted breeding with genic simple sequence repeats (SSR) markers have been used for developing terminal heat-tolerant wheat varieties; however, only few studies involved the use of microRNA (miRNA)-based SSR markers (miRNA-SSRs) in wheat, which were found as key players in various abiotic stresses. In the present study, we identified 104 heat-stress-responsive miRNAs reported in various crops. Out of these, 70 miRNA-SSR markers have been validated on a set of 20 terminal heat-tolerant and heat-susceptible wheat genotypes. Among these, only 19 miRNA-SSR markers were found to be polymorphic, which were further used to study the genetic diversity and population structure. The polymorphic miRNA-SSRs amplified 61 SSR loci with an average of 2.9 alleles per locus. The polymorphic information content (PIC) value of polymorphic miRNA-SSRs ranged from 0.10 to 0.87 with a mean value of 0.48. The dendrogram constructed using unweighted neighbor-joining method and population structure analysis clustered these 20 wheat genotypes into 3 clusters. The target genes of these miRNAs are involved either directly or indirectly in providing tolerance to heat stress. Furthermore, two polymorphic markers miR159c and miR165b were declared as very promising diagnostic markers, since these markers showed specific alleles and discriminated terminal heat-tolerant genotypes from the susceptible genotypes. Thus, these identified miRNA-SSR markers will prove useful in the characterization of wheat germplasm through the study of genetic diversity and population structural analysis and in wheat molecular breeding programs aimed at terminal heat tolerance of wheat varieties.

Twentieth-century changes in the genetic composition of Swedish field pea metapopulations

Landrace crops are formed by local adaptation, genetic drift and gene flow through seed exchange. In reverse, the study of genetic structure between landrace populations can reveal the effects of these forces over time. We present here the analysis of genetic diversity in 40 Swedish field pea (Pisum sativum L.) populations, either available as historical seed samples from the late nineteenth century or as extant gene bank accessions assembled in the late twentieth century. The historical material shows constant high levels of within-population diversity, whereas the extant accessions show varying, and overall lower, levels of within-population diversity. Structure and principal component analysis cluster most accessions, both extant and historical, in groups after geographical origin. County-wise analyses of the accessions show that the genetic diversity of the historical accessions is largely overlapping. In contrast, most extant accessions show signs of genetic drift. They harbor a subset of the alleles found in the historical accessions and are more differentiated from each other. These results reflect how, historically present metapopulations have been preserved during the twentieth century, although as genetically isolated populations.

Sequencing, de novo assembly and comparative analysis of Raphanus sativus transcriptome

Raphanus sativus is an important Brassicaceae plant and also an edible vegetable with great economic value. However, currently there is not enough transcriptome information of R. sativus tissues, which impedes further functional genomics research on R. sativus. In this study, RNA-seq technology was employed to characterize the transcriptome of leaf tissues. Approximately 70 million clean pair-end reads were obtained and used for de novo assembly by Trinity program, which generated 68,086 unigenes with an average length of 576 bp. All the unigenes were annotated against GO and KEGG databases. In the meanwhile, we merged leaf sequencing data with existing root sequencing data and obtained better de novo assembly of R. sativus using Oases program. Accordingly, potential simple sequence repeats (SSRs), transcription factors (TFs) and enzyme codes were identified in R. sativus. Additionally, we detected a total of 3563 significantly differentially expressed genes (DEGs, P = 0.05) and tissue-specific biological processes between leaf and root tissues. Furthermore, a TFs-based regulation network was constructed using Cytoscape software. Taken together, these results not only provide a comprehensive genomic resource of R. sativus but also shed light on functional genomic and proteomic research on R. sativus in the future.

Simple Repeats as Building Blocks for Genetic Computers

Processing of RNA involves heterogeneous nuclear ribonucleoproteins. The simple sequence repeats (SSRs) they bind can also adopt alternative DNA structures, like Z DNA, triplexes, G quadruplexes, and I motifs. Those SSRs capable of switching conformation under physiological conditions (called flipons) are genetic elements that can encode alternative RNA processing by their effects on RNA processivity, most likely as DNA:RNA hybrids. Flipons are elements of a binary, instructive genetic code directing how genomic sequences are compiled into transcripts. The combinatorial nature of this code provides a rich set of options for creating genetic computers able to reproduce themselves and use a heritable and evolvable code to optimize survival. The underlying computational logic potentiates a diverse set of genetic programs that modify cis-mediated heritability and disease risk.

De Novo Sequencing and Hybrid Assembly of the Biofuel Crop Jatropha curcas L.: Identification of Quantitative Trait Loci for Geminivirus Resistance

Jatropha curcas is an important perennial, drought tolerant plant that has been identified as a potential biodiesel crop. We report here the hybrid de novo genome assembly of J. curcas generated using Illumina and PacBio sequencing technologies, and identification of quantitative loci for Jatropha Mosaic Virus (JMV) resistance. In this study, we generated scaffolds of 265.7 Mbp in length, which correspond to 84.8% of the gene space, using Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis. Additionally, 96.4% of predicted protein-coding genes were captured in RNA sequencing data, which reconfirms the accuracy of the assembled genome. The genome was utilized to identify 12,103 dinucleotide simple sequence repeat (SSR) markers, which were exploited in genetic diversity analysis to identify genetically distinct lines. A total of 207 polymorphic SSR markers were employed to construct a genetic linkage map for JMV resistance, using an interspecific F₂ mapping population involving susceptible J. curcas and resistant Jatropha integerrima as parents. Quantitative trait locus (QTL) analysis led to the identification of three minor QTLs for JMV resistance, and the same has been validated in an alternate F₂ mapping population. These validated QTLs were utilized in marker-assisted breeding for JMV resistance. Comparative genomics of oil-producing genes across selected oil producing species revealed 27 conserved genes and 2986 orthologous protein clusters in Jatropha. This reference genome assembly gives an insight into the understanding of the complex genetic structure of Jatropha, and serves as source for the development of agronomically improved virus-resistant and oil-producing lines.

Microsatellite Diversity, Complexity, and Host Range of Mycobacteriophage Genomes of the Siphoviridae Family

The incidence, distribution, and variation of simple sequence repeats (SSRs) in viruses is instrumental in understanding the functional and evolutionary aspects of repeat sequences. Full-length genome sequences retrieved from NCBI were used for extraction and analysis of repeat sequences using IMEx software. We have also developed two MATLAB-based tools for extraction of gene locations from GenBank in tabular format and simulation of this data with SSR incidence data. Present study encompassing 147 Mycobacteriophage genomes revealed 25,284 SSRs and 1,127 compound SSRs (cSSRs) through IMEx. Mono- to hexa-nucleotide motifs were present. The SSR count per genome ranged from 78 (M100) to 342 (M58) while cSSRs incidence ranged from 1 (M138) to 17 (M28, M73). Though cSSRs were present in all the genomes, their frequency and SSR to cSSR conversion percentage varied from 1.08 (M138 with 93 SSRs) to 8.33 (M116 with 96 SSRs). In terms of localization, the SSRs were predominantly localized to coding regions (∼78%). Interestingly, genomes of around 50 kb contained a similar number of SSRs/cSSRs to that in a 110 kb genome, suggesting functional relevance for SSRs which was substantiated by variation in motif constitution between species with different host range. The three species with broad host range (M97, M100, M116) have around 90% of their mono-nucleotide repeat motifs composed of G or C and only M16 has both A and T mononucleotide motifs. Around 20% of the di-nucleotide repeat motifs in the genomes exhibiting a broad host range were CT/TC, which were either absent or represented to a much lesser extent in the other genomes.

Molecular structure and phylogenetic analysis of complete chloroplast genomes of medicinal species Paeonia lactiflora from Zhejiang Province

Paeonia lactiflora is a geo-authentic and superior medicinal plant in Zhejiang province. Here, we report the complete chloroplast genome sequence of P. lactiflora. The total genome size of P. lactiflora is 153,405 bp in length, including a small single-copy (SSC) region of 16,969 bp, a large single-copy (LSC) region of 84,340 bp separated by a pair of inverted repeats (IRs) of 26,048 bp. The overall annotated gene number is 109, containing 76 protein-coding genes, 29 tRNAs and 4 rRNAs. The entire GC content of P. lactiflora is 38.43%, with the highest GC content of 42.99% in IR region. A total of 52 simple sequence repeats are identified in the cp genome of P. lactiflora. Phylogenetic analysis indicated a sister relationship between P. lactiflora and P. veitchii, and supported a unique evolutionary status of Family Paeoniaceae. This work provides a valuable genetic resource to develop robust markers and investigate the population genetics diversities for this famous medicinal species.

Transcriptomic Analysis of the Endangered Neritid Species Clithon retropictus: De Novo Assembly, Functional Annotation, and Marker Discovery

An aquatic gastropod belonging to the family Neritidae, Clithon retropictus is listed as an endangered class II species in South Korea. The lack of information on its genomic background limits the ability to obtain functional data resources and inhibits informed conservation planning for this species. In the present study, the transcriptomic sequencing and de novo assembly of C. retropictus generated a total of 241,696,750 high-quality reads. These assembled to 282,838 unigenes with mean and N50 lengths of 736.9 and 1201 base pairs, respectively. Of these, 125,616 unigenes were subjected to annotation analysis with known proteins in Protostome DB, COG, GO, and KEGG protein databases (BLASTX; E ≤ 0.00001) and with known nucleotides in the Unigene database (BLASTN; E ≤ 0.00001). The GO analysis indicated that cellular process, cell, and catalytic activity are the predominant GO terms in the biological process, cellular component, and molecular function categories, respectively. In addition, 2093 unigenes were distributed in 107 different KEGG pathways. Furthermore, 49,280 simple sequence repeats were identified in the unigenes (>1 kilobase sequences). This is the first report on the identification of transcriptomic and microsatellite resources for C. retropictus, which opens up the possibility of exploring traits related to the adaptation and acclimatization of this species.

PSR: polymorphic SSR retrieval

BACKGROUND

With the advent of high-throughput sequencing technologies large-scale identification of microsatellites became affordable and was especially directed to non-model species. By contrast, few efforts have been published toward the automatic identification of polymorphic microsatellites by exploiting sequence redundancy. Few tools for genotyping microsatellite repeats have been implemented so far that are able to manage huge amount of sequence data and handle the SAM/BAM file format. Most of them have been developed for and tested on human or model organisms with high quality reference genomes.

RESULTS

In this note we describe polymorphic SSR retrieval (PSR), a read counter and simple sequence repeat (SSR) length polymorphism detection tool. It is written in Perl and was developed to identify length polymorphisms in perfect microsatellites exploiting next generation sequencing (NGS) data. PSR has been developed bearing in mind plant non-model species for which de novo transcriptome assembly is generally the first sequence resource available to be used for SSR-mining. PSR is divided into two modules: the read-counting module (PSR_read_retrieval) identifies all the reads that cover the full-length of perfect microsatellites; the comparative module (PSR_poly_finder) detects both heterozygous and homozygous alleles at each microsatellite locus across all genotypes under investigation. Two threshold values to call a length polymorphism and reduce the number of false positives can be defined by the user: the minimum number of reads overlapping the repetitive stretch and the minimum read depth. The first parameter determines if the microsatellite-containing sequence must be processed or not, while the second one is decisive for the identification of minor alleles. PSR was tested on two different case studies. The first study aims at the identification of polymorphic SSRs in a set of de novo assembled transcripts defined by RNA-sequencing of two different plant genotypes. The second research activity aims to investigate sequence variations within a collection of newly sequenced chloroplast genomes. In both the cases PSR results are in agreement with those obtained by capillary gel separation.

CONCLUSION

PSR has been specifically developed from the need to automate the gene-based and genome-wide identification of polymorphic microsatellites from NGS data. It overcomes the limits related to the existing and time-consuming efforts based on tools developed in the pre-NGS era.

Fungal Infection Increases the Rate of Somatic Mutation in Scots Pine (Pinus sylvestris L.)

Somatic mutations are transmitted during mitosis in developing somatic tissue. Somatic cells bearing the mutations can develop into reproductive (germ) cells and the somatic mutations are then passed on to the next generation of plants. Somatic mutations are a source of variation essential to evolve new defense strategies and adapt to the environment. Stem rust disease in Scots pine has a negative effect on wood quality, and thus adversely affects the economy. It is caused by the 2 most destructive fungal species in Scandinavia: Peridermium pini and Cronartium flaccidum. We studied nuclear genome stability in Scots pine under biotic stress (fungus-infected, 22 trees) compared to a control population (plantation, 20 trees). Stability was assessed as accumulation of new somatic mutations in 10 microsatellite loci selected for genotyping. Microsatellites are widely used as molecular markers in population genetics studies of plants, and are particularly used for detection of somatic mutations as their rate of mutation is of a much higher magnitude when compared with other DNA markers. We report double the rate of somatic mutation per locus in the fungus-infected trees (4.8×10(-3) mutations per locus), as compared to the controls (2.0×10(-3) mutations per locus) when individual samples were analyzed at 10 different microsatellite markers. Pearson's chi-squared test indicated a significant effect of the fungal infection which increased the number of mutations in the fungus-infected trees (χ(2) = 12.9883, df = 1, P = 0.0003134).

De Novo Transcriptome Sequencing Analysis of cDNA Library and Large-Scale Unigene Assembly in Japanese Red Pine (Pinus densiflora)

Japanese red pine (Pinus densiflora) is extensively cultivated in Japan, Korea, China, and Russia and is harvested for timber, pulpwood, garden, and paper markets. However, genetic information and molecular markers were very scarce for this species. In this study, over 51 million sequencing clean reads from P. densiflora mRNA were produced using Illumina paired-end sequencing technology. It yielded 83,913 unigenes with a mean length of 751 bp, of which 54,530 (64.98%) unigenes showed similarity to sequences in the NCBI database. Among which the best matches in the NCBI Nr database were Picea sitchensis (41.60%), Amborella trichopoda (9.83%), and Pinus taeda (4.15%). A total of 1953 putative microsatellites were identified in 1784 unigenes using MISA (MicroSAtellite) software, of which the tri-nucleotide repeats were most abundant (50.18%) and 629 EST-SSR (expressed sequence tag- simple sequence repeats) primer pairs were successfully designed. Among 20 EST-SSR primer pairs randomly chosen, 17 markers yielded amplification products of the expected size in P. densiflora. Our results will provide a valuable resource for gene-function analysis, germplasm identification, molecular marker-assisted breeding and resistance-related gene(s) mapping for pine for P. densiflora.

Phase variable DNA repeats in Neisseria gonorrhoeae influence transcription, translation, and protein sequence variation

There are many types of repeated DNA sequences in the genomes of the species of the genus Neisseria, from homopolymeric tracts to tandem repeats of hundreds of bases. Some of these have roles in the phase-variable expression of genes. When a repeat mediates phase variation, reversible switching between tract lengths occurs, which in the species of the genus Neisseria most often causes the gene to switch between on and off states through frame shifting of the open reading frame. Changes in repeat tract lengths may also influence the strength of transcription from a promoter. For phenotypes that can be readily observed, such as expression of the surface-expressed Opa proteins or pili, verification that repeats are mediating phase variation is relatively straightforward. For other genes, particularly those where the function has not been identified, gathering evidence of repeat tract changes can be more difficult. Here we present analysis of the repetitive sequences that could mediate phase variation in the Neisseria gonorrhoeae strain NCCP11945 genome sequence and compare these results with other gonococcal genome sequences. Evidence is presented for an updated phase-variable gene repertoire in this species, including a class of phase variation that causes amino acid changes at the C-terminus of the protein, not previously described in N. gonorrhoeae.