Identification of microsatellites in cattle unigenes

Development of 16 novel EST-SSR markers for species identification and cross-genus amplification in sambar, sika, and red deer

Deer genera around the globe are threatened by anthropogenic interference. The translocation of alien species and their subsequent genetic introgression into indigenous deer populations is particularly harmful to the species of greatest conservation concern. Products derived from deer, including venison and antler velvet, are also at risk of fraudulent labeling. The current molecular markers used to genetically identify deer species were developed from genome sequences and have limited applicability for cross-species amplification. The absence of efficacious diagnostic techniques for identifying deer species has hampered conservation and wildlife crime investigation efforts. Expressed sequence tag-simple sequence repeat (EST-SSR) markers are reliable tools for individual and species identification, especially in terms of cross-species genotyping. We conducted transcriptome sequencing of sambar (Rusa unicolor) antler velvet and acquired 11,190 EST-SSRs from 65,074 newly assembled unigenes. We identified a total of 55 unambiguous amplicons in sambar (n = 45), which were selected as markers to evaluate cross-species genotyping in sika deer (Cervus nippon, n = 30) and red deer (Cervus elaphus, n = 46), resulting in cross-species amplification rates of 94.5% and 89.1%, respectively. Based on polymorphic information content (>0.25) and genotyping fidelity, we selected 16 of these EST-SSRs for species identification. This marker set revealed significant genetic differentiation based on the fixation index and genetic distance values. Principal coordinate analysis and STRUCTURE analysis revealed distinct clusters of species and clearly identified red-sika hybrids. These markers showed applicability across different genera and proved suitable for identification and phylogenetic analyses across deer species.

Development of 49 novel microsatellite markers from Next-generation sequencing data and a robust method for parentage tests in the emu (Dromaius novaehollandiae)

The emu is a useful and new breed of poultry, but their genetic improvement has not advanced yet due to their very recent domestication. Pedigree information is difficult to record because of their complex reproduction system (polyandry). To identify parent-offspring relationships in the emu, parentage test based on polymorphic DNA markers have to be developed. In this study, we isolated more than 25,000 microsatellite (simple sequence repeat, SSR) regions from Next-generation sequencing data via the QDD pipeline and developed 49 SSR markers with polymorphism in the Japanese farmed emu. The dinucleotide motifs, (AC)_n, (AT)_n and (AG)_n, were the most frequently detected and were found on 10,167 (38.55%), 8,114 (30.76%) and 4,796 (18.18%) contigs, respectively. Forty-nine novel SSR markers were characterized in 20 individuals and showed N_A ranged from 2 to 12, with an average of 4.2. H_E/H_O ranged from 0.389/0.071 to 0.702/1.000 with an average of 0.601/0.515. PIC value ranged from 0.059 to 0.886 with an average of 0.528, and 17 of 49 markers showed a higher polymorphism than 0.500. Thirty-four individuals were genotyped using 12 markers, and CERVUS simulations based on genotype showed that parents of all offspring were identified with 0.9995-1.0 probability. Thus, 49 novel SSR markers and a robust method for parentage test for the Japanese emu were developed.

Reduced representation approach for identification of genome-wide SNPs and their annotation for economically important traits in Indian Tharparkar cattle

The present study was carried out in Tharparkar cattle for identification of genome-wide SNPs and microsatellites, and then annotate the identified high-quality SNPs to milk production, fertility, carcass, adaptability and immune response of economically important traits. A total of 146,011 SNPs were identified with respect to Bos taurus reference genome which are indicus specific, out of which 10,519 SNPs were found to be novel. Similarly, a total of 87,047 SNPs were identified with respect to Bos indicus reference genome. After final annotation of SNPs identified with respect to Bos indicus reference genome, 2871 SNPs were found to be associated in 383 candidate genes having to do with milk production, fertility, carcass, immune response and adaptability traits. Following that, 2571 microsatellites were identified. The information mined from the data might be of importance for the future breed improvement programs, conservation efforts and for enhancing the SNPs density of the existing bovine SNP chips.

Mining and comparative survey of EST-SSR markers among members of Euphorbiaceae family

Euphorbiaceae represents flowering plants family of tropical and sub-tropical region rich in secondary metabolites of economic importance. To understand and assess the genetic makeup among the members, this study was undertaken to characterize and compare SSR markers from publicly available ESTs and GSSs of nine selected species of the family. Mining of SSRs was performed by MISA, primer designing by Primer3, while functional annotation, gene ontology (GO) and enrichment analysis were performed by Blast2GO. A total 12,878 number of SSRs were detected from 101,701 number of EST sequences. SSR density ranged from 1 SSR/3.22 kb to 1 SSR/15.65 kb. A total of 1873 primer pairs were designed for the annotated SSR-Contigs. About 77.07% SSR-ESTs could be assigned a significant match to the protein database. 3037 unique SSR-FDM were assigned and IPR003657 (WRKY Domain) was found to be the most dominant FDM among the members. 1810 unique GO terms obtained were further subjected to enrichment analysis to obtain 513 statistically significant GO terms mapped to the SSR containing ESTs. Most frequent enriched GO terms were, GO:0003824 for molecular function, GO:0006350 for biological process and GO:0005886 for cellular component, justifying the richness of defensive secondary metabolites and phytomedicine within the family. The results from this study provides tangible insight to genetic make-up and distribution of SSRs. Functional annotation corresponded many genes of unknown functions which may be considered as novel genes or genes responsible for stress specific secondary metabolites. Further studies are required to understand stress specific genes accountable for leveraging the synthesis of secondary metabolites.

Genome-wide mapping and characterization of microsatellites in the swamp eel genome

We described genome-wide screening and characterization of microsatellites in the swamp eel genome. A total of 99,293 microsatellite loci were identified in the genome with an overall density of 179 microsatellites per megabase of genomic sequences. The dinucleotide microsatellites were the most abundant type representing 71% of the total microsatellite loci and the AC-rich motifs were the most recurrent in all repeat types. Microsatellite frequency decreased as numbers of repeat units increased, which was more obvious in long than short microsatellite motifs. Most of microsatellites were located in non-coding regions, whereas only approximately 1% of the microsatellites were detected in coding regions. Trinucleotide repeats were most abundant microsatellites in the coding regions, which represented amino acid repeats in proteins. There was a chromosome-biased distribution of microsatellites in non-coding regions, with the highest density of 203.95/Mb on chromosome 8 and the least on chromosome 7 (164.06/Mb). The most abundant dinucleotides (AC)n was mainly located on chromosome 8. Notably, genomic mapping showed that there was a chromosome-biased association of genomic distributions between microsatellites and transposon elements. Thus, the novel dataset of microsatellites in swamp eel provides a valuable resource for further studies on QTL-based selection breeding, genetic resource conservation and evolutionary genetics.

Utilization of in silico EST-SSR markers for diversity studies in castor (Ricinus communis L.)

Castor (Ricinus communis L.) a chief non-edible oilseed crop has numerous industrial applications. Systematic genetic diversity analysis utilizing DNA based markers has been quick and reliable method that ensures selection of diverse parents for exploitation of higher levels of heterosis in breeding programs. From NCBI database, 63,852 EST sequences of castor were mined. One thousand one hundred and five (1105) EST-SSRs and 1652 repeat motifs sequences were identified from 20,495 non-redundant unigene sequences. Repeat motifs consisted of 29.7 % mono nucleotide repeats, 24.8 % di nucleotide repeats, 27.27 % tri nucleotide repeats and 3.94 % tetra nucleotide repeats. Twenty eight primer pairs were chosen from SSR-containing ESTs to determine genetic diversity among 27 castor accessions. Twelve EST-SSRs showed polymorphism. Number of alleles detected were 2-3 with an average of 2.33 per locus. 150-400 bp was the size of an allele. Dendrogram analysis grouped the 27 accessions into two separate clusters. Genetic similarity coefficient of dendrogram ranged from 0.24 to 0.83. The polymorphic information content value of 0.28-0.49 revealed medium level of diversity in castor. Results of present study indicated that EST-SSRs to be efficient markers for genetic diversity studies. Knowledge on level of diversity existing in castor genotypes would be useful for breeders to plan efficient hybrid breeding programme.

De Novo Transcriptome Assembly of the Chinese Swamp Buffalo by RNA Sequencing and SSR Marker Discovery

The Chinese swamp buffalo (Bubalis bubalis) is vital to the lives of small farmers and has tremendous economic importance. However, a lack of genomic information has hampered research on augmenting marker assisted breeding programs in this species. Thus, a high-throughput transcriptomic sequencing of B. bubalis was conducted to generate transcriptomic sequence dataset for gene discovery and molecular marker development. Illumina paired-end sequencing generated a total of 54,109,173 raw reads. After trimming, de novo assembly was performed, which yielded 86,017 unigenes, with an average length of 972.41 bp, an N50 of 1,505 bp, and an average GC content of 49.92%. A total of 62,337 unigenes were successfully annotated. Among the annotated unigenes, 27,025 (43.35%) and 23,232 (37.27%) unigenes showed significant similarity to known proteins in NCBI non-redundant protein and Swiss-Prot databases (E-value < 1.0E-5), respectively. Of these annotated unigenes, 14,439 and 15,813 unigenes were assigned to the Gene Ontology (GO) categories and EuKaryotic Ortholog Group (KOG) cluster, respectively. In addition, a total of 14,167 unigenes were assigned to 331 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Furthermore, 17,401 simple sequence repeats (SSRs) were identified as potential molecular markers. One hundred and fifteen primer pairs were randomly selected for amplification to detect polymorphisms. The results revealed that 110 primer pairs (95.65%) yielded PCR amplicons and 69 primer pairs (60.00%) presented polymorphisms in 35 individual buffaloes. A phylogenetic analysis showed that the five swamp buffalo populations were clustered together, whereas two river buffalo breeds clustered separately. In the present study, the Illumina RNA-seq technology was utilized to perform transcriptome analysis and SSR marker discovery in the swamp buffalo without using a reference genome. Our findings will enrich the current SSR markers resources and help spearhead molecular genetic research studies on the swamp buffalo.

Developing new SSR markers from ESTs of pea (Pisum sativum L.)

The development of expressed sequence tags (ESTs) from pea has provided a useful source for mining novel simple sequence repeat (SSR) markers. In the present research, in order to find EST-derived SSR markers, 18 552 pea ESTs from the National Center for Biotechnology Information (NCBI) database were downloaded and assembled into 10 086 unigenes. A total of 586 microsatellites in 530 unigenes were identified, indicating that merely 5.25% of sequences contained SSRs. The most abundant SSRs within pea were tri-nucleotide repeat motifs, and among all the tri-nucleotide repeats, the motif GAA was the most abundant type. In total, 49 SSRs were used for primer design. EST-SSR loci were subsequently screened on 10 widely adapted varieties in China. Of these, nine loci showed polymorphic profiles that revealed two to three alleles per locus. The polymorphism information content value ranged from 0.18 to 0.58 with an average of 0.41. Furthermore, transferable analysis revealed that some of these loci showed transferability to faba bean. Because of their polymorphism and transferability, these nine novel EST-SSRs will be valuable tools for marker-assisted breeding and comparative mapping of pea in the future.