The rise, fall and renaissance of microsatellites in eukaryotic genomes

Development and characterization of polymorphic microsatellite markers for Periplaneta americana (Blattodea: Blattidae)

The American cockroach, Periplaneta americana (Blattodea: Blattidae), is a prevalent urban pest with significant public health implications. This study aimed to develop and validate novel microsatellite markers to understand the genetic diversity and population genetic structure of P. americana. In this study, a total of 397,898 microsatellite markers were developed based on 24.6 million genomic DNA sequences. Twenty microsatellite markers were selected and amplified with varying numbers of alleles ranging from 0 to 35. Seven out of 20 markers were characterized for their polymorphism and amplification efficiency. The polymorphic information content (PIC) values of these markers were high (0.669 to 0.950) implying their effectiveness. These markers also revealed 7 to 35 alleles per locus across tested samples, highlighting their utility in assessing the extensive genetic variation within Periplaneta americana populations. These results provide insightful information that may be applied to the genetic analysis of the American cockroach population using the developed species-specific microsatellite marker.

Construction of a Full-Length Transcriptome of Western Honeybee Midgut Tissue and Improved Genome Annotation

Honeybees are an indispensable pollinator in nature with pivotal ecological, economic, and scientific value. However, a full-length transcriptome for Apis mellifera, assembled with the advanced third-generation nanopore sequencing technology, has yet to be reported. Here, nanopore sequencing of the midgut tissues of uninoculated and Nosema ceranae-inoculated A. mellifera workers was conducted, and the full-length transcriptome was then constructed and annotated based on high-quality long reads. Next followed improvement of sequences and annotations of the current reference genome of A. mellifera. A total of 5,942,745 and 6,664,923 raw reads were produced from midguts of workers at 7 days post-inoculation (dpi) with N. ceranae and 10 dpi, while 7,100,161 and 6,506,665 raw reads were generated from the midguts of corresponding uninoculated workers. After strict quality control, 6,928,170, 6,353,066, 5,745,048, and 6,416,987 clean reads were obtained, with a length distribution ranging from 1 kb to 10 kb. Additionally, 16,824, 17,708, 15,744, and 18,246 full-length transcripts were respectively detected, including 28,019 nonredundant ones. Among these, 43,666, 30,945, 41,771, 26,442, and 24,532 full-length transcripts could be annotated to the Nr, KOG, eggNOG, GO, and KEGG databases, respectively. Additionally, 501 novel genes (20,326 novel transcripts) were identified for the first time, among which 401 (20,255), 193 (13,365), 414 (19,186), 228 (12,093), and 202 (11,703) were respectively annotated to each of the aforementioned five databases. The expression and sequences of three randomly selected novel transcripts were confirmed by RT-PCR and Sanger sequencing. The 5' UTR of 2082 genes, the 3' UTR of 2029 genes, and both the 5' and 3' UTRs of 730 genes were extended. Moreover, 17,345 SSRs, 14,789 complete ORFs, 1224 long non-coding RNAs (lncRNAs), and 650 transcription factors (TFs) from 37 families were detected. Findings from this work not only refine the annotation of the A. mellifera reference genome, but also provide a valuable resource and basis for relevant molecular and -omics studies.

Effects of trimer repeats on Psidium guajava L. gene expression and prospection of functional microsatellite markers

Most research on trinucleotide repeats (TRs) focuses on human diseases, with few on the impact of TR expansions on plant gene expression. This work investigates TRs' effect on global gene expression in Psidium guajava L., a plant species with widespread distribution and significant relevance in the food, pharmacology, and economics sectors. We analyzed TR-containing coding sequences in 1,107 transcripts from 2,256 genes across root, shoot, young leaf, old leaf, and flower bud tissues of the Brazilian guava cultivars Cortibel RM and Paluma. Structural analysis revealed TR sequences with small repeat numbers (5-9) starting with cytosine or guanine or containing these bases. Functional annotation indicated TR-containing genes' involvement in cellular structures and processes (especially cell membranes and signal recognition), stress response, and resistance. Gene expression analysis showed significant variation, with a subset of highly expressed genes in both cultivars. Differential expression highlighted numerous down-regulated genes in Cortibel RM tissues, but not in Paluma, suggesting interplay between tissues and cultivars. Among 72 differentially expressed genes with TRs, 24 form miRNAs, 13 encode transcription factors, and 11 are associated with transposable elements. In addition, a set of 20 SSR-annotated, transcribed, and differentially expressed genes with TRs was selected as phenotypic markers for Psidium guajava and, potentially for closely related species as well.

Ancient and Modern Genomes Reveal Microsatellites Maintain a Dynamic Equilibrium Through Deep Time

Microsatellites are widely used in population genetics, but their evolutionary dynamics remain poorly understood. It is unclear whether microsatellite loci drift in length over time. This is important because the mutation processes that underlie these important genetic markers are central to the evolutionary models that employ microsatellites. We identify more than 27 million microsatellites using a novel and unique dataset of modern and ancient Adélie penguin genomes along with data from 63 published chordate genomes. We investigate microsatellite evolutionary dynamics over 2 timescales: one based on Adélie penguin samples dating to ∼46.5 ka and the other dating to the diversification of chordates aged more than 500 Ma. We show that the process of microsatellite allele length evolution is at dynamic equilibrium; while there is length polymorphism among individuals, the length distribution for a given locus remains stable. Many microsatellites persist over very long timescales, particularly in exons and regulatory sequences. These often retain length variability, suggesting that they may play a role in maintaining phenotypic variation within populations.

GenoSSRFinder: a tool for rapid, precise, and targeted simple sequence repeat detection in genomic studies

The GenoSSRFinder is a new tool enables the research of Simple Sequence Repeats (SSRs) in DNA sequences and genomes much simpler and more precise in short time. The analysis is carried out by targeting a certain SSR in genome and gene sequences. This utility is quick, accurate, and does its function quite well. It quickly goes across the sequence, revealing all of the locations at which the selected SSR may be found. This tool will tell researchers where selected SSR begins and where it stops, how long it is, how often it repeats, and how long each repetition is. GenoSSRFinder gets the findings quickly, and they will be simple to comprehend. Therefore, when studying SSRs, researchers will have more time to use to thorough work as a result of this time savings. In addition, it provides a valuable information since it is highly precise. GenoSSRFinder is simple to use and produces high-quality findings. It is also accelerating SSRs gene research, which is a direct result of the new approach we use to analyse SSRs. Three case studies in this study demonstrated the usefulness of this program by immediately studying a particular SSR that was associated with genetic illness, biodiversity and criminal science in living organisms. This demonstration explains that GenoSSRFinder might be utilized in a wide variety of fields, such as the research of genetic illnesses, the biodiversity and genetic studies, or even in criminal investigations.

Transcriptome analysis of Haemaphysalis flava female using Illumina HiSeq 4000 sequencing: de novo assembly, functional annotation and discovery of SSR markers

BACKGROUND

Ticks are ectoparasites capable of directly damaging their hosts and transmitting vector-borne diseases. The ixodid tick Haemaphysalis flava has a broad distribution that extends from East to South Asia. This tick is a reservoir of severe fever with thrombocytopenia syndrome virus (SFTSV) that causes severe hemorrhagic disease, with cases reported from China, Japan and South Korea. Recently, the distribution of H. flava in South Korea was found to overlap with the occurrence of SFTSV.

METHODS

This study was undertaken to discover the molecular resources of H. flava female ticks using the Illumina HiSeq 4000 system, the Trinity de novo sequence assembler and annotation against public databases. The locally curated Protostome database (PANM-DB) was used to screen the putative adaptation-related transcripts classified to gene families, such as angiotensin-converting enzyme, aquaporin, adenylate cyclase, AMP-activated protein kinase, glutamate receptors, heat shock proteins, molecular chaperones, insulin receptor, mitogen-activated protein kinase and solute carrier family proteins. Also, the repeats and simple sequence repeats (SSRs) were screened from the unigenes using RepeatMasker (v4.0.6) and MISA (v1.0) software tools, followed by the designing of SSRs flanking primers using BatchPrimer 3 (v1.0) software.

RESULTS

The transcriptome produced a total of 69,822 unigenes, of which 46,175 annotated to the homologous proteins in the PANM-DB. The unigenes were also mapped to the EuKaryotic Orthologous Groups (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) specializations. Promiscuous presence of protein kinase, zinc finger (C2H2-type), reverse transcriptase, and RNA recognition motif domains was observed in the unigenes. A total of 3480 SSRs were screened, of which 1907 and 1274 were found as tri- and dinucleotide repeats, respectively. A list of primer sequences flanking the SSR motifs was detailed for validation of polymorphism in H. flava and the related tick species.

CONCLUSIONS

The reference transcriptome information on H. flava female ticks will be useful for an enriched understanding of tick biology, its competency to act as a vector and the study of species diversity related to disease transmission.

Random amplified microsatellites (RAMS) analysis ascertains genetic variation of Alternaria alternata causing black spot disease on Carya illinoinensis in South Africa

Limited information regarding the occurrence of black spot disease of pecan (Carya illinoinensis), caused by A. alternata, in South Africa is known. The pecan industry is growing rapidly, so it is essential to understand the impact of the fungal pathogen to pecan health. In this study, the genetic variation of 364 A. alternata isolates was investigated by two RAMS primers (CCA5 and CGA5). In total, 6,525 alleles were produced, with a minimum of 3,182 alleles on the CGA5 primer and maximum of 3,343 alleles for CCA5 primer. Further analysis of the primers showed relatively low genetic diversity of A. alternata isolate populations, with mean values; (H = 0.12) and Shannon's information index (I = 0.20). The analysis of molecular variance (AMOVA) revealed significant differences between populations, with 88% of the genetic variation was found within populations (Nm = 3.59, PhiPT = 0.12), and were not significantly different (p > 0.001). While 12% variation was observed among populations (Nm = 2.89, PhiPT = 0.08) and the estimates were statistically significant (p < 0.001). STRUCTURE HARVESTER output showed that K value is K = 8, where ΔK cannot find the true number of populations because of less variation. The dendrogram cluster tree generated by Ward's analysis unveiled two main distinct clades and 10 sub-clades, revealing similar findings as those of PCoA analysis clusters. Therefore, it was evident that these analyses depicted no distinct relationship between the A. alternata isolates and their geographic locations or the prevalence of distribution among the populations.

Whole Genome Wide SSR Markers Identification Based on ddRADseq Data

The advent of advanced NGS technologies have led to the generation of enormous amount of sequence data which further aid in the discovery of the various type of markers such as SSRs, SNPs, InDels, etc. Among all these markers, microsatellite SSR markers can be mined from the ddRADseq data as certain properties of SSR markers make them ideal markers for study. These assist researchers and breeders in diversity analysis and producing new varieties with desired traits. To extract the markers, first, the ddRADseq data is assembled into consensus sequences using STACKS program which are further assembled for mining microsatellites using QDD along with MISA tool.

Tolerance to dietary linalool primarily involves co-expression of cytochrome P450s and cuticular proteins in Pagiophloeus tsushimanus (Coleoptera: Curculionidae) larvae using SMRT sequencing and RNA-seq

BACKGROUND

Pagiophloeus tsushimanus (Coleoptera: Curculionidae), an emerging forest pest exclusively infesting camphor trees, has recently caused severe ecological and economic damage in localized areas in China. Its population outbreak depends largely on the capacity to overcome the pressure of terpenoid-derived metabolites (e.g. linalool) from camphor trees. At present, the molecular basis of physiological adaptation of P. tsushimanus to dietary linalool is poorly understood, and there is no available reference genome or transcriptome.

RESULTS

Herein, we constructed the transcriptome profiling of P. tsushimanus larvae reared on linalool-infused diets using RNA sequencing and single-molecule real-time sequencing. A total of 20,325 high-quality full-length transcripts were identified as a reference transcriptome, of which 14,492 protein-coding transcripts including 130 transcription factors (TFs), and 5561 long non-coding RNAs (lncRNAs) were detected. Also, 30 alternative splicing events and 8049 simple sequence repeats were captured. Gene ontology enrichment of differential expressed transcripts revealed that overall up-regulation of both cytochrome P450s (CYP450s) and cuticular proteins (CPs), was the primary response characteristic against dietary linalool. Other physiological effects possibly caused by linalool exposure, such as increase in Reactive Oxygen Species (ROS) and hormetic stimulation, were compensated by a handful of induced genes encoding antioxidases, heat shock proteins (HSPs), juvenile hormone (JH) epoxide hydrolases, and digestive enzymes. Additionally, based on co-expression networks analysis, a diverse array of hub lncRNAs and TFs co-expressed with CYP450s and CPs were screened as the potential gene regulators. Temporal expression of candidate transcripts determined by quantitative real-time PCR also indicated a cooperative relationship between the inductions of CYP450s and CPs upon exposure to linalool.

CONCLUSIONS

Our present study provides an important transcriptome resource of P. tsushimanus, and lays a valuable foundation for understanding how this specialist pest copes with chemical challenges in its specific host environments.

SSR identification and phylogenetic analysis in four plant species based on complete chloroplast genome sequences

The effective utilization of traditional Chinese medicine (TCM) has been challenged by the difficulty to accurately distinguish between similar plant varieties. The stability and conservation of the chloroplast genome can aid in resolving genotypes. Previous studies using nuclear sequences and molecular markers have not effectively differentiated the species from related taxa, such as Machilus leptophylla, Hanceola exserta, Rubus bambusarum, and Rubus henryi. This study aimed to characterize the chloroplast genomes of these four plant species, and analyze their simple sequence repeats (SSRs) and phylogenetic positions. The results demonstrated the four chloroplast genomes consisted of 152.624 kb, 153.296 kb, 156.309 kb, and 158.953 kb in length, involving 124, 130, 129, and 131 genes, respectively. They also contained four specific regions with mononucleotide being the class with the most members. Moreover, these repeating types of SSR were various in individual class. Phylogenetic analysis showed that M. leptophylla was clustered with M. yunnanensis, and H. exserta was confirmed as belonging to the family Ocimeae. Additionally, R. bambusarum and R. henryi were grouped together but differed in their SSR features, indicating that they were not the same species. This research provides evidence for resolving species and contributes new genetic information for further studies.

Development of genome-specific SSR markers for the identification of introgressed segments of Sinapis alba in the Brassica juncea background

Sinapis alba L. (white mustard) is recognized for carrying host resistance against several biotic stresses including, Alternaria brassicae, which is responsible for blight disease in cultivated Brassica. However, another cultivated Brassica has a dearth for genetic resistance for these stresses due to its narrow genetic base. Therefore, we performed introgression of the genomic regions of S. alba into backcross progenies of B. juncea + S. alba somatic hybrids. These advanced generations with S. alba chromosomal segments are named B. juncea-S. alba introgression lines (ILs). In the present study, we developed the S. alba genome-specific microsatellites from the draft genome to track the S. alba genome introgressions and responsible regions for resistance to A. brassicae. For developing these SSR markers, the unique contigs of S. alba draft genome were identified through BLASTN with B. juncea, B. rapa, B. nigra, and B. oleracea reference genome assemblies, including mitochondrial and chloroplast genomes, and further used for marker development. Out of 403,423 contigs, we have identified 65,343 non-hit contigs of S. alba that yielded a total of 1231 genome-specific microsatellites, out of which 1107 were expected to produce a single allele upon amplification. Out of the total SSRs, 234 primer pairs were randomly picked from whole-genome and validated between B. juncea and S. alba genomes for their specificity. In the validation experiment, these markers gave a single amplicon into S. alba, while they did not amplify in B. juncea genome. Of these, 59 microsatellites were used to track S. alba introgressions in 80 BC2F3 lines. To the best of our knowledge, this is the first time that these two genetic resources are developed in the form of B. juncea-S. alba ILs and S. alba-specific markers. Therefore, both the resources unlock a new avenue of Brassica breeding for biotic and abiotic stresses along with quality traits.

Supplementary Information

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

Panmixia and active colonisation of the invasive palm Trachycarpus fortunei (Arecaceae) in Southern Switzerland and Northern Italy as inferred by microsatellites and SNP markers

Trachycarpus fortunei (Arecaceae: Coryphoideae) is an Asian palm that was introduced during the nineteenth century in southern Switzerland and northern Italy as an ornamental plant. In the recent decades, the palm has become an aggressive invasive species in the region. Before this study, the genetic structure and diversity of the naturalised populations were unknown. We aimed at understanding the dynamics of invasion and at comparing the results obtained with two types of markers. This genetic approach aimed at tracing back as far as possible the source of invasive populations comparing historical information found in literature and invasive genetic patterns. The genetic diversity was analysed using eight microsatellites (five were developed for that purpose) and 31′000 SNPs identified through GBS analyses. Genetic analyses were carried out for 200 naturalised individuals sampled from 21 populations in the Canton Ticino (Switzerland) and the provinces of Lombardy and Piedmont (Italy). The observed general panmixia indicates that the expansion of T. fortunei is active in its naturalised areas. The genetic pattern found for both SNPs and microsatellites appears to be related to the colonisation process, with a lack of geographic structure and bottleneck signatures occurring at the colonisation front, distantly from historical sites. This study gives a better understanding of the expansion of T. fortunei and adds new insights to its ecology.

Characterization of Microsatellites in the Akebia trifoliata Genome and Their Transferability and Development of a Whole Set of Effective, Polymorphic, and Physically Mapped Simple Sequence Repeat Markers

Akebia trifoliata is a perennial climbing woody liana plant with a high potential for commercial exploitation and theoretical research. Similarly, microsatellites (simple sequence repeats, SSRs) also have dual roles: as critical markers and as essential elements of the eukaryotic genome. To characterize the profile of SSRs and develop molecular markers, the high-quality assembled genome of A. trifoliata was used. Additionally, to determine the potential transferability of SSR loci, the genomes of Amborella trichopoda, Oryza sativa, Vitis vinifera, Arabidopsis thaliana, Papaver somniferum, and Aquilegia coerulea were also used. We identified 434,293 SSRs with abundant short repeats, such as 290,868 (66.98%) single-nucleotide repeats (SNRs) and 113,299 (26.09%) dinucleotide repeats (DNRs) in the A. trifoliata genome. 398,728 (91.81%) SSRs on 344,283 loci were physically mapped on the chromosomes, and a positive correlation (r = 0.98) was found between the number of SSRs and chromosomal length. Additionally, 342,916 (99.60%) potential SSR markers could be designed from the 344,283 physically mapped loci, while only 36,160 could be viewed as high-polymorphism-potential (HPP) markers, findings that were validated by PCR. Finally, SSR loci exhibited broad potential transferability, particularly DNRs such as the "AT/AT" and "AG/CT" loci, among all angiosperms, a finding that was not related to the genetic divergence distance. Practically, we developed a whole set of effective, polymorphic, and physically anchored molecular markers and found that, evolutionarily, DNRs could be responsible for microsatellite origin and protecting gene function.

Sex-specific splicing of Z- and W-borne nr5a1 alleles suggests sex determination is controlled by chromosome conformation

Pogona vitticeps has female heterogamety (ZZ/ZW), but the master sex-determining gene is unknown, as it is for all reptiles. We show that nr5a1 (Nuclear Receptor Subfamily 5 Group A Member 1), a gene that is essential in mammalian sex determination, has alleles on the Z and W chromosomes (Z-nr5a1 and W-nr5a1), which are both expressed and can recombine. Three transcript isoforms of Z-nr5a1 were detected in gonads of adult ZZ males, two of which encode a functional protein. However, ZW females produced 16 isoforms, most of which contained premature stop codons. The array of transcripts produced by the W-borne allele (W-nr5a1) is likely to produce truncated polypeptides that contain a structurally normal DNA-binding domain and could act as a competitive inhibitor to the full-length intact protein. We hypothesize that an altered configuration of the W chromosome affects the conformation of the primary transcript generating inhibitory W-borne isoforms that suppress testis determination. Under this hypothesis, the genetic sex determination (GSD) system of P. vitticeps is a W-borne dominant female-determining gene that may be controlled epigenetically.

Genetic Analysis of Plant Pathogens Natural Populations

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

Full-Length SMRT Transcriptome Sequencing and SSR Analysis of Bactrocera dorsalis (Hendel)

Simple Summary

In this study, a full-length transcriptome was analyzed with single-molecule real-time (SMRT) sequencing, which was first used to discover simple sequence repeat (SSR) genetic markers from B. dorsalis. Moreover, SSR markers from isoforms were screened for the identification of species diversity. These results could provide molecular biology methods for further population research.

Abstract

Bactrocera dorsalis (Hendel), as one of the most notorious and destructive invasive agricultural pests in the world, causes damage to over 250 different types of fruits and vegetables throughout tropical and subtropical areas. PacBio single-molecule real-time (SMRT) sequencing was used to generate the full-length transcriptome data of B. dorsalis. A total of 40,319,890 subreads (76.6 Gb, clean reads) were generated, including 535,241 circular consensus sequences (CCSs) and 386,916 full-length non-concatemer reads (FLNCs). Transcript cluster analysis of the FLNC reads revealed 22,780 high-quality reads (HQs). In total, 12,274 transcripts were functionally annotated based on four different databases. A total of 1978 SSR loci were distributed throughout 1714 HQ transcripts, of which 1926 were complete SSRs and 52 were complex SSRs. Among the total SSR loci, 2–3 nucleotide repeats were dominant, occupying 83.62%, of which di- and tri- nucleotide repeats were 39.38% and 44.24%, respectively. We detected 105 repeat motifs, of which AT/AT (50.19%), AC/GT (39.15%), CAA/TTG (32.46%), and ACA/TGT (10.86%) were the most common in di- and tri-nucleotide repeats. The repeat SSR motifs were 12–190 bp in length, and 1638 (88.02%) were shorter than 20 bp. According to the randomly selected microsatellite sequence, 80 pairs of primers were designed, and 174 individuals were randomly amplified by PCR using primers. The number of primers that had amplification products with clear bands and showed good polymorphism came to 41, indicating that this was a feasible way to explore SSR markers from the transcriptomic data of B. dorsalis. These results lay a foundation for developing highly polymorphic microsatellites for researching the functional genomics, population genetic structure, and genetic diversity of B. dorsalis.

Characterization and complexity of transcriptome in Gymnocypris przewalskii using single-molecule long-read sequencing and RNA-seq

The Tibetan Schizothoracinae fish Gymnocypris przewalskii has the ability to adapt to the extreme plateau environment, making it an ideal biological material for evolutionary biology research. However, the lack of well-annotated reference genomes has limited the study of the molecular genetics of G. przewalskii. To characterize its transcriptome features, we first used long-read sequencing technology in combination with RNA-seq for transcriptomic analysis. A total of 159,053 full-length (FL) transcripts were captured by Iso-Seq, having a mean length of 3,445 bp with N50 value of 4,348. Of all FL transcripts, 145,169 were well-annotated in the public database and 134,537 contained complete open reading frames. There were 4,149 pairs of alternative splicing events, of which three randomly selected were defined by RT–PCR and sequencing, and 13,293 long non-coding RNAs detected, based on all-vs.-all BLAST. A total of 118,185 perfect simple sequence repeats were identified from FL transcripts. The FL transcriptome might provide basis for further research of G. przewalskii.

PolyQ length co-evolution in neural proteins

Intermolecular co-evolution optimizes physiological performance in functionally related proteins, ultimately increasing molecular co-adaptation and evolutionary fitness. Polyglutamine (polyQ) repeats, which are over-represented in nervous system-related proteins, are increasingly recognized as length-dependent regulators of protein function and interactions, and their length variation contributes to intraspecific phenotypic variability and interspecific divergence. However, it is unclear whether polyQ repeat lengths evolve independently in each protein or rather co-evolve across functionally related protein pairs and networks, as in an integrated regulatory system. To address this issue, we investigated here the length evolution and co-evolution of polyQ repeats in clusters of functionally related and physically interacting neural proteins in Primates. We observed function-/disease-related polyQ repeat enrichment and evolutionary hypervariability in specific neural protein clusters, particularly in the neurocognitive and neuropsychiatric domains. Notably, these analyses detected extensive patterns of intermolecular polyQ length co-evolution in pairs and clusters of functionally related, physically interacting proteins. Moreover, they revealed both direct and inverse polyQ length co-variation in protein pairs, together with complex patterns of coordinated repeat variation in entire polyQ protein sets. These findings uncover a whole system of co-evolving polyQ repeats in neural proteins with direct implications for understanding polyQ-dependent phenotypic variability, neurocognitive evolution and neuropsychiatric disease pathogenesis.

Full-length SMRT transcriptome sequencing and microsatellite characterization in Paulownia catalpifolia

Paulownia catalpifolia is an important, fast-growing timber species known for its high density, color and texture. However, few transcriptomic and genetic studies have been conducted in P. catalpifolia. In this study, single-molecule real-time sequencing technology was applied to obtain the full-length transcriptome of P. catalpifolia leaves treated with varying degrees of drought stress. The sequencing data were then used to search for microsatellites, or simple sequence repeats (SSRs). A total of 28.83 Gb data were generated, 25,969 high-quality (HQ) transcripts with an average length of 1624 bp were acquired after removing the redundant reads, and 25,602 HQ transcripts (98.59%) were annotated using public databases. Among the HQ transcripts, 16,722 intact coding sequences, 149 long non-coding RNAs and 179 alternative splicing events were predicted, respectively. A total of 7367 SSR loci were distributed throughout 6293 HQ transcripts, of which 763 complex SSRs and 6604 complete SSRs. The SSR appearance frequency was 28.37%, and the average distribution distance was 5.59 kb. Among the 6604 complete SSR loci, 1-3 nucleotide repeats were dominant, occupying 97.85% of the total SSR loci, of which mono-, di- and tri-nucleotide repeats were 44.68%, 33.86% and 19.31%, respectively. We detected 112 repeat motifs, of which A/T (42.64%), AG/CT (12.22%), GA/TC (9.63%), GAA/TTC (1.57%) and CCA/TGG (1.54%) were most common in mono-, di- and tri-nucleotide repeats, respectively. The length of the repeat SSR motifs was 10-88 bp, and 4997 (75.67%) were ≤ 20 bp. This study provides a novel full-length transcriptome reference for P. catalpifolia and will facilitate the identification of germplasm resources and breeding of new drought-resistant P. catalpifolia varieties.

Cytogenetic and genetic data support Crossodactylus aeneus Müller, 1924 as a new junior synonym of C. gaudichaudii Duméril and Bibron, 1841 (Amphibia, Anura)

The nominal anuran species Crossodactylus gaudichaudii Duméril and Bibron, 1841 and Crossodactylus aeneus Müller, 1924 are indistinguishable based on adult and larval morphology, being subject of taxonomic doubts. Here, we describe the karyotypes of C. gaudichaudii and C. aeneus, using classical and molecular cytogenetic markers. In addition, we used sequences of the H1 mitochondrial DNA to infer their phylogenetic relationships by Maximum Likelihood (ML) and Maximum Parsimony (MP) approaches and species delimitation test (by bPTP approach). The karyotypic data do not differentiate C. gaudichaudii and C. aeneus in any of the chromosome markers assessed. In both phylogenetic analyses, C. gaudichaudii and C. aeneus were recovered into a strongly supported clade. The species delimitation analysis recovered the specimens assigned to C. gaudichaudii and C. aeneus as a single taxonomic unit. Taken the cytogenetic and genetic results together with previous studies of internal and external morphology of tadpoles and biacoustic pattern, C. gaudichaudii and C. aeneus could not be differentiated, which supports the hypothesis that they correspond to the same taxonomic unit, with C. aeneus being a junior synonym of C. gaudichaudii.

Genome survey of Zanthoxylum bungeanum and development of genomic-SSR markers in congeneric species

Zanthoxylum bungeanum, a spice and medicinal plant, is cultivated in many parts of China and some countries in Southeast Asia; however, data on its genome are lacking. In the present study, we performed a whole-genome survey and developed novel genomic-SSR markers of Z. bungeanum. Clean data (∼197.16 Gb) were obtained and assembled into 11185221 scaffolds with an N50 of 183 bp. K-mer analysis revealed that Z. bungeanum has an estimated genome size of 3971.92 Mb, and the GC content, heterozygous rate, and repeat sequence rate are 37.21%, 1.73%, and 86.04%, respectively. These results indicate that the genome of Z. bungeanum is complex. Furthermore, 27153 simple sequence repeat (SSR) loci were identified from 57288 scaffolds with a minimum length > 1 kb. Mononucleotide repeats (19706) were the most abundant type, followed by dinucleotide repeats (5154). The most common motifs were A/T, followed by AT/AT; these SSRs accounted for 71.42% and 11.84% of all repeats, respectively. A total of 21243 non-repeating primer pairs were designed, and 100 were randomly selected and validated by PCR analysis using DNA from 10 Z. bungeanum individuals and 5 Zanthoxylum armatum individuals. Finally, 36 polymorphic SSR markers were developed with polymorphism information content (PIC) values ranging from 0.16 to 0.75. Cluster analysis revealed that Z. bungeanum and Z. armatum could be divided into two major clusters, suggesting that these newly developed SSR markers are useful for genetic diversity and germplasm resource identification in Z. bungeanum and Z. armatum.

Whole genome survey and microsatellite motif identification of Artemia franciscana

Artemia is an industrially important genus used in aquaculture as a nutritious diet for fish and as an aquatic model organism for toxicity tests. However, despite the significance of Artemia, genomic research remains incomplete and knowledge on its genomic characteristics is insufficient. In particular, Artemia franciscana of North America has been widely used in fisheries of other continents, resulting in invasion of native species. Therefore, studies on population genetics and molecular marker development as well as morphological analyses are required to investigate its population structure and to discriminate closely related species. Here, we used the Illumina Hi-Seq platform to estimate the genomic characteristics of A. franciscana through genome survey sequencing (GSS). Further, simple sequence repeat (SSR) loci were identified for microsatellite marker development. The predicted genome size was ∼867 Mb using K-mer (a sequence of k characters in a string) analysis (K = 17), and heterozygosity and duplication rates were 0.655 and 0.809%, respectively. A total of 421467 SSRs were identified from the genome survey assembly, most of which were dinucleotide motifs with a frequency of 77.22%. The present study will be a useful basis in genomic and genetic research for A. franciscana.

Guidelines for the Choice of Sequences for Molecular Plant Taxonomy

This chapter presents an overview of the major plant DNA sequences and molecular methods available for plant taxonomy. Guidelines are provided for the choice of sequences and methods to be used, based on the DNA compartment (nuclear, chloroplastic, mitochondrial), evolutionary mechanisms, and the level of taxonomic differentiation of the plants under survey.

Brainwide Genetic Sparse Cell Labeling to Illuminate the Morphology of Neurons and Glia with Cre-Dependent MORF Mice

Cajal recognized that the elaborate shape of neurons is fundamental to their function in the brain. However, there are no simple and generalizable genetic methods to study neuronal or glial cell morphology in the mammalian brain. Here, we describe four mouse lines conferring Cre-dependent sparse cell labeling based on mononucleotide repeat frameshift (MORF) as a stochastic translational switch. Notably, the optimized MORF3 mice, with a membrane-bound multivalent immunoreporter, confer Cre-dependent sparse and bright labeling of thousands of neurons, astrocytes, or microglia in each brain, revealing their intricate morphologies. MORF3 mice are compatible with imaging in tissue-cleared thick brain sections and with immuno-EM. An analysis of 151 MORF3-labeled developing retinal horizontal cells reveals novel morphological cell clusters and axonal maturation patterns. Our study demonstrates a conceptually novel, simple, generalizable, and scalable mouse genetic solution to sparsely label and illuminate the morphology of genetically defined neurons and glia in the mammalian brain.

De novo transcriptome assembly and population genetic analyses of an important coastal shrub, Apocynum venetum L

BACKGROUND

Apocynum venetum L. is an important medicinal plant that is mainly distributed in the coastal areas and northwest of China. In addition to its high medical and economic value, its adaptation to saline-alkali and coastal saline lands makes A. venetum an ideal candidate for use in vegetation restoration. To date, the study of A. venetum has been limited in the northwest region of China, little attention has been paid to the genetic diversity and population structure of A. venetum populations in the coastal region. Here, we performed transcriptome sequencing of total RNA from A. venetum leaves and developed efficient expressed sequence tag-simple sequence repeat (EST-SSR) markers for analyzing the genetic diversity and population structure of A. venetum in the coastal region.

RESULTS

A total of 86,890 unigenes were generated after de novo assembly, and 68,751 of which were successfully annotated by searching against seven protein databases. Furthermore, 14,072 EST-SSR loci were detected and 10,243 primer pairs were successfully designed from these loci. One hundred primer pairs were randomly selected and synthesized, twelve primer pairs were identified as highly polymorphic and further used for population genetic analysis. Population genetic analyses showed that A. venetum exhibited low level of genetic diversity (mean alleles per locus, N_A = 3.3; mean expected heterozygosity, H_E = 0.342) and moderate level of genetic differentiation among the populations (genetic differentiation index, F_ST = 0.032-0.220) in the coastal region. Although the contemporary (mean m_c = 0.056) and historical (mean m_h = 0.106) migration rates among the six A. venetum populations were moderate, a decreasing trend over the last few generations was detected. Bayesian structure analysis clustered six populations into two major groups, and genetic bottlenecks were found to have occurred in two populations (QG, BH).

CONCLUSIONS

Using novel EST-SSR markers, we evaluated the genetic variation of A. venetum in the coastal region and determined conservation priorities based on these findings. The large dataset of unigenes and SSRs identified in our study, combining samples from a broader range, will support further research on the conservation and evolution of this important coastal plant and its related species.

Vector competence of Aedes albopictus populations for chikungunya virus is shaped by their demographic history

The mosquito Aedes albopictus is one of the most dangerous invasive species. Its worldwide spread has created health concerns as it is a major vector of arboviruses of public health significance such as chikungunya (CHIKV). Dynamics of different genetic backgrounds and admixture events may have impacted competence for CHIKV in adventive populations. Using microsatellites, we infer the genetic structure of populations across the expansion areas that we then associate with their competence for different CHIKV genotypes. Here we show that the demographic history of Ae. albopictus populations is a consequence of rapid complex patterns of historical lineage diversification and divergence that influenced their competence for CHIKV. The history of adventive populations is associated with CHIKV genotypes in a genotype-by-genotype interaction that impacts their vector competence. Thus, knowledge of the demographic history and vector competence of invasive mosquitoes is pivotal for assessing the risk of arbovirus outbreaks in newly colonized areas.

Being Merle: The Molecular Genetic Background of the Canine Merle Mutation

The intensity of the merle pattern is determined by the length of the poly(A) tail of a repeat element which has been inserted into the boundary of intron 10 and exon 11 of the PMEL17 locus in reverse orientation. This poly(A) tail behaves as a microsatellite, and due to replication slippage, longer and shorter alleles of it might be generated during cell divisions. The length of the poly(A) tail regulates the splicing mechanism. In the case of shorter tails, the removal of intron 10 takes place at the original splicing, resulting in a normal premelanosome protein (PMEL). Longer tails generate larger insertions, forcing splicing to a cryptic splice site, thereby coding for an abnormal PMEL protein, which is unable to form the normal fibrillar matrix of the eumelanosomes. Thus, eumelanin deposition ensuring the dark color formation is reduced. In summary, the longer the poly(A) tail, the lighter the coat color intensity of the melanocytes. These mutations can occur in the somatic cells and the resulting cell clones will shape the merle pattern of the coat. When they take place in the germ line, they occasionally produce offspring with unexpected color variations which are different from those of their parents.

Compound Dynamics and Combinatorial Patterns of Amino Acid Repeats Encode a System of Evolutionary and Developmental Markers

Homopolymeric amino acid repeats (AARs) like polyalanine (polyA) and polyglutamine (polyQ) in some developmental proteins (DPs) regulate certain aspects of organismal morphology and behavior, suggesting an evolutionary role for AARs as developmental "tuning knobs." It is still unclear, however, whether these are occasional protein-specific phenomena or hints at the existence of a whole AAR-based regulatory system in DPs. Using novel approaches to trace their functional and evolutionary history, we find quantitative evidence supporting a generalized, combinatorial role of AARs in developmental processes with evolutionary implications. We observe nonrandom AAR distributions and combinations in HOX and other DPs, as well as in their interactomes, defining elements of a proteome-wide combinatorial functional code whereby different AARs and their combinations appear preferentially in proteins involved in the development of specific organs/systems. Such functional associations can be either static or display detectable evolutionary dynamics. These findings suggest that progressive changes in AAR occurrence/combination, by altering embryonic development, may have contributed to taxonomic divergence, leaving detectable traces in the evolutionary history of proteomes. Consistent with this hypothesis, we find that the evolutionary trajectories of the 20 AARs in eukaryotic proteomes are highly interrelated and their individual or compound dynamics can sharply mark taxonomic boundaries, or display clock-like trends, carrying overall a strong phylogenetic signal. These findings provide quantitative evidence and an interpretive framework outlining a combinatorial system of AARs whose compound dynamics mark at the same time DP functions and evolutionary transitions.

Polymorphism of simple sequence repeats may quantitatively regulate gene transcription

The degree of polymorphism, i.e., DNA sequence divergence, of short AT-rich tandemly arranged simple sequence repeats at or near promoters and 5'- untranslated regions of mRNA may quantitatively regulate transcription of tissue-specific genes. Less polymorphic repeats allow greater gene expression. Preferential binding of hypophosphorylated H1 histone to these repeats may diminish binding of transcription factors. Preferential binding of hypophosphorylated high mobility group chromatin proteins would increase this binding. Shorter simple sequence repeats have undergone fewer point mutations than longer repeats, hence they are less polymorphic and more conserved. The role of transcribed simple sequence repeats in frog embryo germ layer determination is considered.

Finding and extending ancient simple sequence repeat-derived regions in the human genome

Background

Previously, 3% of the human genome has been annotated as simple sequence repeats (SSRs), similar to the proportion annotated as protein coding. The origin of much of the genome is not well annotated, however, and some of the unidentified regions are likely to be ancient SSR-derived regions not identified by current methods. The identification of these regions is complicated because SSRs appear to evolve through complex cycles of expansion and contraction, often interrupted by mutations that alter both the repeated motif and mutation rate. We applied an empirical, kmer-based, approach to identify genome regions that are likely derived from SSRs.

Results

The sequences flanking annotated SSRs are enriched for similar sequences and for SSRs with similar motifs, suggesting that the evolutionary remains of SSR activity abound in regions near obvious SSRs. Using our previously described P-clouds approach, we identified ‘SSR-clouds’, groups of similar kmers (or ‘oligos’) that are enriched near a training set of unbroken SSR loci, and then used the SSR-clouds to detect likely SSR-derived regions throughout the genome.

Conclusions

Our analysis indicates that the amount of likely SSR-derived sequence in the human genome is 6.77%, over twice as much as previous estimates, including millions of newly identified ancient SSR-derived loci. SSR-clouds identified poly-A sequences adjacent to transposable element termini in over 74% of the oldest class of Alu (roughly, AluJ), validating the sensitivity of the approach. Poly-A’s annotated by SSR-clouds also had a length distribution that was more consistent with their poly-A origins, with mean about 35 bp even in older Alus. This work demonstrates that the high sensitivity provided by SSR-Clouds improves the detection of SSR-derived regions and will enable deeper analysis of how decaying repeats contribute to genome structure.

New EST-SSR Markers for Individual Genotyping of Opium Poppy Cultivars (Papaver somniferum L.)

High-quality simple sequence repeat (SSR) markers are invaluable tools for revealing genetic variability which could be utilized for many purposes, such as breeding new varieties or the identifying current ones, among other applications. Based on the analysis of 3.7 million EST sequences and 15 genomic sequences from bacterial artificial chromosome (BAC) libraries, 200 trinucleotide genic (EST)-SSR and three genomic (gSSR) markers were tested, where 17 of them fulfilled all criteria for quality markers. Moreover, the reproducibility of these new markers was verified by two genetics laboratories, with a mean error rate per allele and per locus equal to 0.17%. These markers were tested on 38 accessions of Papaver somniferum and nine accessions of another five species of the Papaver and Argemone genera. In total, 118 alleles were detected for all accessions (median = 7; three to ten alleles per locus) and 88 alleles (median = 5; three to nine alleles per locus) within P. somniferum alone. Multivariate methods and identity analysis revealed high resolution capabilities of the new markers, where all but three pair accessions (41 out of 47) had a unique profile and opium poppy was distinguished from other species.

Conversion of DNA Sequences: From a Transposable Element to a Tandem Repeat or to a Gene

Eukaryotic genomes are rich in repetitive DNA sequences grouped in two classes regarding their genomic organization: tandem repeats and dispersed repeats. In tandem repeats, copies of a short DNA sequence are positioned one after another within the genome, while in dispersed repeats, these copies are randomly distributed. In this review we provide evidence that both tandem and dispersed repeats can have a similar organization, which leads us to suggest an update to their classification based on the sequence features, concretely regarding the presence or absence of retrotransposons/transposon specific domains. In addition, we analyze several studies that show that a repetitive element can be remodeled into repetitive non-coding or coding sequences, suggesting (1) an evolutionary relationship among DNA sequences, and (2) that the evolution of the genomes involved frequent repetitive sequence reshuffling, a process that we have designated as a "DNA remodeling mechanism". The alternative classification of the repetitive DNA sequences here proposed will provide a novel theoretical framework that recognizes the importance of DNA remodeling for the evolution and plasticity of eukaryotic genomes.

Hybridization, characterization and transferability of SSRs in the genus Morchella

Recently, Morchella importuna, M. sextelata, M. eximia, M. exuberans, Mel-13, and Mel-21 have been successfully cultivated in China and planting areas rapidly expanded because of their economic importance. Effective molecular markers are urgently needed for accurately identifying morel cultivars. Microsatellites are widely utilized for strain authentication in many fungal groups. To our knowledge, for the first time we characterized the distribution of microsatellites (simple sequence repeats, SSRs) in the M. importuna genome with 12902 SSRs and reported the first set of SSRs developed for Morchella species. Mono-nucleotides (66.2 %) were the most frequent motifs, followed by tri- (15.4 %), di- (12.1 %), tetra- (3.7 %), penta- (1.3 %) and hexa-nucleotides (1.3 %). We tested the cross-species amplification of 180 SSRs on 24 samples from the six species and high cross-species transferability of SSRs (87.7 %) was found. Among twenty-two microsatellites selected for genetic diversity analysis on 127 samples from the six species, fifteen to twenty polymorphic loci were identified in M. importuna, M. sextelata, M. eximia, M. exuberans, Mel-13 and Mel-21. Interspecific hybridization events were detected among morel species, indicating the potential application of morel crossbreeding. Ninety-one cultivated samples were characterized as new cultivars with different genotypes, but cultivar names used for these by farmers was confusing, with misnaming, synonyms and homonyms. Our results are not only helpful for cultivar identification and morel breeding programs in China, but also provide molecular tools for genetic studies in morels.

Genome-wide characterization of microsatellites and genetic diversity assessment of spinach in the Chinese germplasm collection

Spinach is a nutritional leafy green vegetable, and it also serves as a model species for studying sex chromosome evolution. Genetic marker development and genome structure analysis are important in breeding practice and theoretical evolution studies of spinach. In this study, the frequency and distribution of different microsatellites in the recently released draft spinach genome were characterized. A total of 261,002 perfect microsatellites were identified (estimated frequency: ~262.1 loci/Mbp). The most abundant microsatellites were tetranucleotide and trinucleotide, accounting for 33.2% and 27.7% of the total number of microsatellites, respectively. A total of 105 primer pairs were designed and screened, and 34 were polymorphic among the detected spinach cultivars. Combined with seven primer sets developed previously, 41 primer pairs were used to investigate genetic diversity among 43 spinach cultivars in China. The average polymorphism information content value of the 41 markers was 0.43, representing an intermediate level. The spinach cultivars had a low genetic diversity, and no detectable common factors were shared by each group in the UPGMA dendrogram. This study's findings facilitate further investigations on the organization of the microsatellites in spinach genome and provide clues for future breeding applications of spinach in China.

Isolation and characterization of microsatellite markers in Garcinia gummi-gutta by next-generation sequencing and cross-species amplification

Garcinia gummi-gutta (L.) Roxb. (Clusiaceae) is an endemic, semidomesticated, fruit-yielding tree species distributed in the Western Ghats of India and Sri Lanka. Various bioactive phytochemicals, such as garcinol, benzophenones and xanthones are isolated from G. gummi-gutta and have shown antibacterial, antiviral and antioxidant activities. We sequenced the total genomic DNA using Illumina Hiseq 2000 platform and examined 241,141,804 bp high quality data, assembled into 773,889 contigs. In these contigs, 27,313 simple-sequence repeats (SSRs) were identified, among which mononucleotide repeats were predominant (44.98%) followed by dinucleotide and trinucleotide repeats. Primers were designed for 9964 microsatellites among which 32 randomly selected SSR primer pairs were standardized for amplification. Polymerase chain reaction (PCR) amplification of genomic DNA in 30 G. gummi-gutta genotypes revealed polymorphic information content (PIC) across all 32 loci ranging from 0.867 to 0.951, with a mean value of 0.917. The observed and expected heterozygosity ranged from 0.00 to 0.63 and 0.896 to 0.974, respectively. Alleles per locus ranged from 12 to 27. This is the first report on the development of genomic SSR markers in G. gummi-gutta using next-generation sequencing technology. The genomic SSR markers developed in this study will be useful in identification, mapping, diversity and breeding studies.

Characterization of tri and tetra-nucleotide microsatellite loci for the freshwater snails Promenetus exacuous (Planorbidae) and Valvata tricarinata (Valvatidae) and their utility in population genetic studies

Objective

Promenetus exacuous and Valvata tricarinata are freshwater snail species with widespread distributions throughout North America. Information regarding their genetic diversity and population connectivity are currently lacking. We utilized next generation sequencing to develop the first microsatellites for each species to investigate genetic diversity within and differentiation among populations.

Results

Sixteen and seventeen microsatellite loci were developed for P. exacuous and V. tricarinata, respectively, and tested in a total of 43 P. exacuous and 48 V. tricarinata from two lakes approximately 183 km apart in New York State, USA. Fifteen P. exacuous loci were polymorphic in at least one lake and possessed 1–23 alleles and observed heterozygosities of 0.00–0.96 within individual lakes. Seventeen polymorphic V. tricarinata loci possessed 2–19 alleles and observed heterozygosities of 0.04–0.96 within lakes. Bayesian clustering using 12 loci for each species identified two distinct genetic populations, reflecting the two lakes. High assignment scores for individual snails to the lakes they were collected from supported strong population structure with minimal admixture at the scale of this study. These loci will be useful for investigating the genetic diversity and population structure of these species and indicate genetic differentiation may be common among their populations.

The Complete Chloroplast Genome Sequence of Tree of Heaven (Ailanthus altissima (Mill.) (Sapindales: Simaroubaceae), an Important Pantropical Tree

Ailanthus altissima (Mill.) Swingle (Simaroubaceae) is a deciduous tree widely distributed throughout temperate regions in China, hence suitable for genetic diversity and evolutionary studies. Previous studies in A. altissima have mainly focused on its biological activities, genetic diversity and genetic structure. However, until now there is no published report regarding genome of this plant species or Simaroubaceae family. Therefore, in this paper, we first characterized A. altissima complete chloroplast genome sequence. The tree of heaven chloroplast genome was found to be a circular molecule 160,815 base pairs (bp) in size and possess a quadripartite structure. The A. altissima chloroplast genome contains 113 unique genes of which 79 and 30 are protein coding and transfer RNA (tRNA) genes respectively and also 4 ribosomal RNA genes (rRNA) with overall GC content of 37.6%. Microsatellite marker detection identified A/T mononucleotides as majority SSRs in all the seven analyzed genomes. Repeat analyses of seven Sapindales revealed a total of 49 repeats in A. altissima, Rhus chinensis, Dodonaea viscosa, Leitneria floridana, while Azadirachta indica, Boswellia sacra, and Citrus aurantiifolia had a total of 48 repeats. The phylogenetic analysis using protein coding genes revealed that A. altissima is a sister to Leitneria floridana and also suggested that Simaroubaceae is a sister to Rutaceae family. The genome information reported here could be further applied for evolution and invasion, population genetics, and molecular studies in this plant species and family.

A genome-wide investigation of microsatellite mismatches and the association with body mass among bird species

Mutation rate is usually found to covary with many life history traits of animals such as body mass, which has been readily explained by the higher number of mutation opportunities per unit time. Although the precise reason for the pattern is not yet clear, to determine the universality of this pattern, we tested whether life history traits impact another form of genetic mutation, the motif mismatches in microsatellites. Employing published genome sequences from 65 avian species, we explored the motif mismatches patterns of microsatellites in birds on a genomic level and assessed the relationship between motif mismatches and body mass in a phylogenetic context. We found that small-bodied species have a higher average mismatches and we suggested that higher heterozygosity in imperfect microsatellites lead to the increase of motif mismatches. Our results obtained from this study imply that a negative body mass trend in mutation rate may be a general pattern of avian molecular evolution.

High-Throughput Sequencing Strategy for Microsatellite Genotyping Using Neotropical Fish as a Model

Genetic diversity and population studies are essential for conservation and wildlife management programs. However, monitoring requires the analysis of multiple loci from many samples. These processes can be laborious and expensive. The choice of microsatellites and PCR calibration for genotyping are particularly daunting. Here we optimized a low-cost genotyping method using multiple microsatellite loci for simultaneous genotyping of up to 384 samples using next-generation sequencing (NGS). We designed primers with adapters to the combinatorial barcoding amplicon library and sequenced samples by MiSeq. Next, we adapted a bioinformatics pipeline for genotyping microsatellites based on read-length and sequence content. Using primer pairs for eight microsatellite loci from the fish Prochilodus costatus, we amplified, sequenced, and analyzed the DNA of 96, 288, or 384 individuals for allele detection. The most cost-effective methodology was a pseudo-multiplex reaction using a low-throughput kit of 1 M reads (Nano) for 384 DNA samples. We observed an average of 325 reads per individual per locus when genotyping eight loci. Assuming a minimum requirement of 10 reads per loci, two to four times more loci could be tested in each run, depending on the quality of the PCR reaction of each locus. In conclusion, we present a novel method for microsatellite genotyping using Illumina combinatorial barcoding that dispenses exhaustive PCR calibrations, since non-specific amplicons can be eliminated by bioinformatics analyses. This methodology rapidly provides genotyping data and is therefore a promising development for large-scale conservation-genetics studies.

The effects of transcription and recombination on mutational dynamics of short tandem repeats

Short tandem repeats (STR) are ubiquitous components of the genomic architecture of most living organisms. Recent work has highlighted the widespread functional significance of such repeats, particularly around gene regulation, but the mutational processes underlying the evolution of these highly abundant and highly variable sequences are not fully understood. Traditional models assume that strand misalignment during replication is the predominant mechanism, but empirical data suggest the involvement of other processes including recombination and transcription. Despite this evidence, the relative influences of these processes have not previously been tested experimentally on a genome-wide scale. Using deep sequencing, we identify mutations at >200 microsatellites, across 700 generations in replicated populations of two otherwise identical sexual and asexual Saccharomyces cerevisiae strains. Using generalized linear models, we investigate correlates of STR mutability including the nature of the mutation, STR composition and contextual factors including recombination, transcription and replication origins. Sexual capability was not a significant predictor of microsatellite mutability, but, intriguingly, we identify transcription as a significant positive predictor. We also find that STR density is substantially increased in regions neighboring, but not within, recombination hotspots.

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.

Systems biology approach in plant abiotic stresses

Plant abiotic stresses are the major constraint on plant growth and development, causing enormous crop losses across the world. Plants have unique features to defend themselves against these challenging adverse stress conditions. They modulate their phenotypes upon changes in physiological, biochemical, molecular and genetic information, thus making them tolerant against abiotic stresses. It is of paramount importance to determine the stress-tolerant traits of a diverse range of genotypes of plant species and integrate those traits for crop improvement. Stress-tolerant traits can be identified by conducting genome-wide analysis of stress-tolerant genotypes through the highly advanced structural and functional genomics approach. Specifically, whole-genome sequencing, development of molecular markers, genome-wide association studies and comparative analysis of interaction networks between tolerant and susceptible crop varieties grown under stress conditions can greatly facilitate discovery of novel agronomic traits that protect plants against abiotic stresses.

Genomic data reveal a loss of diversity in two species of tuco-tucos (genus Ctenomys) following a volcanic eruption

Marked reductions in population size can trigger corresponding declines in genetic variation. Understanding the precise genetic consequences of such reductions, however, is often challenging due to the absence of robust pre- and post-reduction datasets. Here, we use heterochronous genomic data from samples obtained before and immediately after the 2011 eruption of the Puyehue-Cordón Caulle volcanic complex in Patagonia to explore the genetic impacts of this event on two parapatric species of rodents, the colonial tuco-tuco (Ctenomys sociabilis) and the Patagonian tuco-tuco (C. haigi). Previous analyses using microsatellites revealed no post-eruption changes in genetic variation in C. haigi, but an unexpected increase in variation in C. sociabilis. To explore this outcome further, we used targeted gene capture to sequence over 2,000 putatively neutral regions for both species. Our data revealed that, contrary to the microsatellite analyses, the eruption was associated with a small but significant decrease in genetic variation in both species. We suggest that genome-level analyses provide greater power than traditional molecular markers to detect the genetic consequences of population size changes, particularly changes that are recent, short-term, or modest in size. Consequently, genomic analyses promise to generate important new insights into the effects of specific environmental events on demography and genetic variation.

An accurate and efficient method for large-scale SSR genotyping and applications

Accurate and efficient genotyping of simple sequence repeats (SSRs) constitutes the basis of SSRs as an effective genetic marker with various applications. However, the existing methods for SSR genotyping suffer from low sensitivity, low accuracy, low efficiency and high cost. In order to fully exploit the potential of SSRs as genetic marker, we developed a novel method for SSR genotyping, named as AmpSeq-SSR, which combines multiplexing polymerase chain reaction (PCR), targeted deep sequencing and comprehensive analysis. AmpSeq-SSR is able to genotype potentially more than a million SSRs at once using the current sequencing techniques. In the current study, we simultaneously genotyped 3105 SSRs in eight rice varieties, which were further validated experimentally. The results showed that the accuracies of AmpSeq-SSR were nearly 100 and 94% with a single base resolution for homozygous and heterozygous samples, respectively. To demonstrate the power of AmpSeq-SSR, we adopted it in two applications. The first was to construct discriminative fingerprints of the rice varieties using 3105 SSRs, which offer much greater discriminative power than the 48 SSRs commonly used for rice. The second was to map Xa21, a gene that confers persistent resistance to rice bacterial blight. We demonstrated that genome-scale fingerprints of an organism can be efficiently constructed and candidate genes, such as Xa21 in rice, can be accurately and efficiently mapped using an innovative strategy consisting of multiplexing PCR, targeted sequencing and computational analysis. While the work we present focused on rice, AmpSeq-SSR can be readily extended to animals and micro-organisms.

RNA biology of disease-associated microsatellite repeat expansions

Microsatellites, or simple tandem repeat sequences, occur naturally in the human genome and have important roles in genome evolution and function. However, the expansion of microsatellites is associated with over two dozen neurological diseases. A common denominator among the majority of these disorders is the expression of expanded tandem repeat-containing RNA, referred to as xtrRNA in this review, which can mediate molecular disease pathology in multiple ways. This review focuses on the potential impact that simple tandem repeat expansions can have on the biology and metabolism of RNA that contain them and underscores important gaps in understanding. Merging the molecular biology of repeat expansion disorders with the current understanding of RNA biology, including splicing, transcription, transport, turnover and translation, will help clarify mechanisms of disease and improve therapeutic development.

Abundance and Characterization of Perfect Microsatellites on the Cattle Y Chromosome

Microsatellites or simple sequence repeats (SSRs) are found in most organisms and play an important role in genomic organization and function. To characterize the abundance of SSRs (1-6 base-pairs [bp]) on the cattle Y chromsome, the relative frequency and density of perfect or uninterrupted SSRs based on the published Y chromosome sequence were examined. A total of 17,273 perfect SSRs were found, with total length of 324.78 kb, indicating that approximately 0.75% of the cattle Y chromosome sequence (43.30 Mb) comprises perfect SSRs, with an average length of 18.80 bp. The relative frequency and density were 398.92 loci/Mb and 7500.62 bp/Mb, respectively. The proportions of the six classes of perfect SSRs were highly variable on the cattle Y chromosome. Mononucleotide repeats had a total number of 8073 (46.74%) and an average length of 15.45 bp, and were the most abundant SSRs class, while the percentages of di-, tetra-, tri-, penta-, and hexa-nucleotide repeats were 22.86%, 11.98%, 11.58%, 6.65%, and 0.19%, respectively. Different classes of SSRs varied in their repeat number, with the highest being 42 for dinucleotides. Results reveal that repeat categories A, AC, AT, AAC, AGC, GTTT, CTTT, ATTT, and AACTG predominate on the Y chromosome. This study provides insight into the organization of cattle Y chromosome repetitive DNA, as well as information useful for developing more polymorphic cattle Y-chromosome-specific SSRs.

Characterization of the whole chloroplast genome of Chikusichloa mutica and its comparison with other rice tribe (Oryzeae) species

Chloroplast genomes are a significant genomic resource in plant species and have been used in many research areas. The complete genomic information from wild crop species could supply a valuable genetic reservoir for breeding. Chikusichloa mutica is one of the most important wild distant relatives of cultivated rice. In this study, we sequenced and characterized its complete chloroplast (cp) genome and compared it with other species in the same tribe. The whole cp genome sequence is 136,603 bp in size and exhibits a typical quadripartite structure with large and small single-copy regions (LSC, 82,327 bp; SSC, 12,598 bp) separated by a pair of 20,839-bp inverted repeats (IRA, B). A total of 110 unique genes are annotated, including 76 protein-coding genes, 4 ribosomal RNA genes and 30 tRNA genes. The genome structure, gene order, GC content, and other features are similar to those of other angiosperm cp genomes. When comparing the cp genomes between Oryzinae and Zizaniinae subtribes, the main differences were found between the junction regions and distribution of simple sequence repeats (SSRs). In comparing the two Chikusichloa species, the genomes were only 40 bp different in length and 108 polymorphic sites, including 83 single nucleotide substitutions (SNPs) and 25 insertion-deletions (Indels), were found between the whole cp genomes. The complete cp genome of C. mutica will be an important genetic tool for future breeding programs and understanding the evolution of wild rice relatives.

Leaf Transcriptome Sequencing for Identifying Genic-SSR Markers and SNP Heterozygosity in Crossbred Mango Variety 'Amrapali' (Mangifera indica L.)

Mango (Mangifera indica L.) is called “king of fruits” due to its sweetness, richness of taste, diversity, large production volume and a variety of end usage. Despite its huge economic importance genomic resources in mango are scarce and genetics of useful horticultural traits are poorly understood. Here we generated deep coverage leaf RNA sequence data for mango parental varieties ‘Neelam’, ‘Dashehari’ and their hybrid ‘Amrapali’ using next generation sequencing technologies. De-novo sequence assembly generated 27,528, 20,771 and 35,182 transcripts for the three genotypes, respectively. The transcripts were further assembled into a non-redundant set of 70,057 unigenes that were used for SSR and SNP identification and annotation. Total 5,465 SSR loci were identified in 4,912 unigenes with 288 type I SSR (n ≥ 20 bp). One hundred type I SSR markers were randomly selected of which 43 yielded PCR amplicons of expected size in the first round of validation and were designated as validated genic-SSR markers. Further, 22,306 SNPs were identified by aligning high quality sequence reads of the three mango varieties to the reference unigene set, revealing significantly enhanced SNP heterozygosity in the hybrid Amrapali. The present study on leaf RNA sequencing of mango varieties and their hybrid provides useful genomic resource for genetic improvement of mango.

Differentiation Among Blumeria graminis f. sp. tritici Isolates Originating from Wild Versus Domesticated Triticum Species in Israel

Israel and its vicinity constitute a center of diversity of domesticated wheat species (Triticum aestivum and T. durum) and their sympatrically growing wild relatives, including wild emmer wheat (T. dicoccoides). We investigated differentiation within the forma specialis of their obligate powdery mildew pathogen, Blumeria graminis f. sp. tritici. A total of 61 B. graminis f. sp. tritici isolates were collected from the three host species in four geographic regions of Israel. Genetic relatedness of the isolates was characterized using both virulence patterns on 38 wheat lines (including 21 resistance gene differentials) and presumptively neutral molecular markers (simple-sequence repeats and single-nucleotide polymorphisms). All isolates were virulent on at least some genotypes of all three wheat species tested. All assays divided the B. graminis f. sp. tritici collection into two distinct groups, those from domesticated hosts and those from wild emmer wheat. One-way migration was detected from the domestic wheat B. graminis f. sp. tritici population to the wild emmer B. graminis f. sp. tritici population at a rate of five to six migrants per generation. This gene flow may help explain the overlap between the distinct domestic and wild B. graminis f. sp. tritici groups. Overall, B. graminis f. sp. tritici is significantly differentiated into wild-emmer and domesticated-wheat populations, although the results do not support the existence of a separate f. sp. dicocci.

Identification and Characterization of Microsatellite Loci in Maqui (Aristotelia chilensis [Molina] Stunz) Using Next-Generation Sequencing (NGS)

Maqui (Aristotelia chilensis [Molina] Stunz) is a small dioecious tree native to South America with edible fruit characterized by very high antioxidant capacity and anthocyanin content. To preserve maqui as a genetic resource it is essential to study its genetic diversity. However, the complete genome is unknown and only a few gene sequences are available in databases. Simple sequence repeats (SSR) markers, which are neutral, co-dominant, reproducible and highly variable, are desirable to support genetic studies in maqui populations. By means of identification and characterization of microsatellite loci from a maqui genotype, using 454 sequencing technology, we develop a set of SSR for this species. Obtaining a total of 165,043 shotgun genome sequences, with an average read length of 387 bases, we covered 64 Mb of the maqui genome. Reads were assembled into 4,832 contigs, while 98,546 reads remained as singletons, generating a total of 103,378 consensus genomic sequences. A total of 24,494 SSR maqui markers were identified. Of them, 15,950 SSR maqui markers were classified as perfects. The most common SSR motifs were dinucleotide (31%), followed by tetranucleotide (26%) and trinucleotide motifs (24%). The motif AG/CT (28.4%) was the most abundant, while the motif AC (89 bp) was the largest. Eleven polymorphic SSRs were selected and used to analyze a population of 40 maqui genotypes. Polymorphism information content (PIC) ranged from 0.117 to 0.82, with an average of 0.58. Non-significant groups were observed in the maqui population, showing a panmictic genetic structure. In addition, we also predicted 11150 putative genes and 3 microRNAs (miRNAs) in maqui sequences. This results, including partial sequences of genes, some miRNAs and SSR markers from high throughput next generation sequencing (NGS) of maqui genomic DNA, constitute the first platform to undertake genetic and molecular studies of this important species.