Genome-wide distribution and organization of microsatellites in six species of birds

Comprehensive Comparative Analysis Sheds Light on the Patterns of Microsatellite Distribution across Birds Based on the Chromosome-Level Genomes

Microsatellites (SSRs) are widely distributed in the genomes of organisms and are an important genetic basis for genome evolution and phenotypic adaptation. Although the distribution patterns of microsatellites have been investigated in many phylogenetic lineages, they remain unclear within the morphologically and physiologically diverse avian clades. Here, based on high-quality chromosome-level genomes, we examined the microsatellite distribution patterns for 53 birds from 16 orders. The results demonstrated that each type of SSR had the same ratio between taxa. For example, the frequency of imperfect SSRs (I-SSRs) was 69.90-84.61%, while perfect SSRs (P-SSRs) were 14.86-28.13% and compound SSRs (C-SSRs) were 0.39-2.24%. Mononucleotide SSRs were dominant for perfect SSRs (32.66-76.48%) in most bird species (98.11%), and A(n) was the most abundant repeat motifs of P-SSRs in all birds (5.42-68.22%). Our study further confirmed that the abundance and diversity of microsatellites were less effected by evolutionary history but its length. The number of P-SSRs decreased with increasing repeat times, and longer P-SSRs motifs had a higher variability coefficient of the repeat copy number and lower diversity, indicating that longer motifs tended to have more stable preferences in avian genomes. We also found that P-SSRs were mainly distributed at the gene ends, and the functional annotation for these genes demonstrated that they were related to signal transduction and cellular process. In conclusion, our research provided avian SSR distribution patterns, which will help to explore the genetic basis for phenotypic diversity in birds.

Genome survey sequencing and mining of genome-wide microsatellite markers in yellow-billed babbler (Turdoides affinis)

Turdoides affinis is a species of group dwelling old world passerine of family Leiothrichidae. Unavailability of genome-wide sequence and species-specific molecular markers have hindered comprehensive understanding of cooperative breeding behaviour in T. affinis. Therefore, we generated genome-wide microsatellite markers through whole genome short read sequencing of T. affinis. A total of 68.8 gigabytes of paired-end raw data were sequenced containing 195,067,054 reads. Total sequenced reads spanned a coverage of 17X with genome size of 1.18 Gb. A large number of microsatellite markers (265,297) were mined in the T. affinis genome using Krait, and 50 most informative markers were identified and validated further. In-silico PCR results validated 47 markers. Of these 47 markers, five were randomly selected and validated in-vitro in twelve individuals of T. affinis. Genotyping data on these five loci estimated observed heterozygosity (H₀) and expected heterozygosity (H_e) ratios between 0.333 - 0.833 and 0.851-0.906, respectively. Effective allele size ranged from 6.698 to 10.667, inbreeding coefficient of the population ranged from 0.080 to 0.631 and null allele frequency was calculated at 0.055 to 0.303. Polymorphic information content of all the five loci varied between 0.850 and 0.906. Probabilities of exclusion and identity across 5 loci was estimated to be 0.95 and 0.0036, respectively. All the loci showed significant adherence to Hardy-Weinberg equilibrium. The microsatellite markers reported in this study will facilitate future population genetics studies on T. affinis and other congeneric species.

Microsatellite Characteristics of Silver Carp (Hypophthalmichthysmolitrix) Genome and Genetic Diversity Analysis in Four Cultured Populations

Hypophthalmichthys molitrix is one of the four most important fish in China and has high breeding potential. However, simple sequence repeat (SSR) markers developed on H. molitrix genome level for genetic diversity analysis are limited. In this study, the distribution characteristics of SSRs in the assembled H. molitrix genome were analyzed, and new markers were developed to preliminarily evaluate the genetic diversity of the four breeding populations. A total of 368,572 SSRs were identified from the H. molitrix genome. The total length of SSRs was 6,492,076 bp, accounting for 0.77% of the total length of the genome sequence. The total frequency and total density were 437.73 loci/Mb and 7713.16 bp/Mb, respectively. Among the 2–6 different nucleotide repeat types, SSRs were dominated by di-nucleotide repeats (204,873, 55.59%), and AC/GT was the most abundant motif. The number of SSRs on each chromosome was positively correlated with the length. The 13 pairs of markers developed were used to analyze the genetic diversity of four cultivated populations in Hubei Province. The results showed that the genetic diversity of the four populations was low, and the ranges of alleles (Na), effective alleles (Ne), observed heterozygosity (Ho), and Shannon’s index information (I) were 3.538–4.462, 2.045–2.461, 0.392–0.450, and 0.879–0.954, respectively. Genetic variation occurs mainly among individuals within populations (95.35%). UPGMA tree and Bayesian analysis showed that four populations could be divided into two different branches. Therefore, the genome-wide SSRs were effectively in genetic diversity analysis on H. molitrix.

Genome-wide characterization of microsatellite DNA in fishes: survey and analysis of their abundance and frequency in genome-specific regions

BACKGROUND

Microsatellite repeats are ubiquitous in organism genomes and play an important role in the chromatin organization, regulation of gene activity, recombination and DNA replication. Although microsatellite distribution patterns have been studied in most phylogenetic lineages, they are unclear in fish species.

RESULTS

Here, we present the first systematic examination of microsatellite distribution in coding and non-coding regions of 14 fish genomes. Our study showed that the number and type of microsatellites displayed nonrandom distribution for both intragenic and intergenic regions, suggesting that they have potential roles in transcriptional or translational regulation and DNA replication slippage theories alone were insufficient to explain the distribution patterns. Our results showed that microsatellites are dominant in non-coding regions. The total number of microsatellites ranged from 78,378 to 1,012,084, and the relative density varied from 4925.76 bp/Mb to 25,401.97 bp/Mb. Overall, (A + T)-rich repeats were dominant. The dependence of repeat abundance on the length of the repeated unit (1-6 nt) showed a great similarity decrease, whereas more tri-nucleotide repeats were found in exonic regions than tetra-nucleotide repeats of most species. Moreover, the incidence of different repeated types appeared species- and genomic-specific. These results highlight potential mechanisms for maintaining microsatellite distribution, such as selective forces and mismatch repair systems.

CONCLUSIONS

Our data could be beneficial for the studies of genome evolution and microsatellite DNA evolutionary dynamics, and facilitate the exploration of microsatellites structural, function, composition mode and molecular markers development in these species.

Comparison of the Microsatellite Distribution Patterns in the Genomes of Euarchontoglires at the Taxonomic Level

Microsatellite or simple sequence repeat (SSR) instability within genes can induce genetic variation. The SSR signatures remain largely unknown in different clades within Euarchontoglires, one of the most successful mammalian radiations. Here, we conducted a genome-wide characterization of microsatellite distribution patterns at different taxonomic levels in 153 Euarchontoglires genomes. Our results showed that the abundance and density of the SSRs were significantly positively correlated with primate genome size, but no significant relationship with the genome size of rodents was found. Furthermore, a higher level of complexity for perfect SSR (P-SSR) attributes was observed in rodents than in primates. The most frequent type of P-SSR was the mononucleotide P-SSR in the genomes of primates, tree shrews, and colugos, while mononucleotide or dinucleotide motif types were dominant in the genomes of rodents and lagomorphs. Furthermore, (A)n was the most abundant motif in primate genomes, but (A)n, (AC)n, or (AG)n was the most abundant motif in rodent genomes which even varied within the same genus. The GC content and the repeat copy numbers of P-SSRs varied in different species when compared at different taxonomic levels, reflecting underlying differences in SSR mutation processes. Notably, the CDSs containing P-SSRs were categorized by functions and pathways using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes annotations, highlighting their roles in transcription regulation. Generally, this work will aid future studies of the functional roles of the taxonomic features of microsatellites during the evolution of mammals in Euarchontoglires.

Comparison of microsatellite distribution patterns in twenty-nine beetle genomes

Simple sequence repeats (SSRs) represent an important source of genetic variation that provides a basis for adaptation to different environments in organisms. In this study, we examined the distribution patterns of SSRs in twenty-nine beetle genomes and carried out Gene Ontology (GO) analysis of CDSs embedded with perfect SSRs (P-SSRs). The results demonstrated that imperfect SSRs (I-SSRs) represented the most abundant SSR category in beetle genomes and in different genomic regions (CDS, exon, and intron regions). The numbers of P-SSRs, I-SSRs, compound SSRs, and variable number tandem repeats were positively correlated with beetle genome size, whereas neither the frequency nor the density of the SSRs was correlated with genome size. Moreover, our results demonstrated that common genomic features of P-SSRs within the same suborder or family of Coleoptera were rare. Mono-, di-, tri-, or tetranucleotide SSRs were the most abundant P-SSR categories in beetle genomes. The preferred predominant repeat motif among the mononucleotide P-SSRs was (A)n, but the most frequent repeat motifs for other length classes varied differentially among these genomes. Furthermore, the P-SSR type with the highest GC content differed in the beetle genomes and in different genomic regions. CV (coefficient of variability) analysis demonstrated that the repeat copy numbers of P-SSRs presented relatively higher variation in introns than in CDSs and exons. The GO terms of CDSs containing P-SSRs for molecular functions were mainly enriched in "binding" and "transcription". Our findings will be useful for studying the functional roles of microsatellite heterogeneity in beetle adaptation.

Comparative analyses of simple sequence repeats (SSRs) in 23 mosquito species genomes: Identification, characterization and distribution (Diptera: Culicidae)

Simple sequence repeats (SSRs) exist in both eukaryotic and prokaryotic genomes and are the most popular genetic markers, but the SSRs of mosquito genomes are still not well understood. In this study, we identified and analyzed the SSRs in 23 mosquito species using Drosophila melanogaster as reference at the whole-genome level. The results show that SSR numbers (33 076-560 175/genome) and genome sizes (574.57-1342.21 Mb) are significantly positively correlated (R² = 0.8992, P < 0.01), but the correlation in individual species varies in these mosquito species. In six types of SSR, mono- to trinucleotide SSRs are dominant with cumulative percentages of 95.14%-99.00% and densities of 195.65/Mb-787.51/Mb, whereas tetra- to hexanucleotide SSRs are rare with 1.12%-4.22% and 3.76/Mb-40.23/Mb. The (A/T)n, (AC/GT)n and (AGC/GCT)n are the most frequent motifs in mononucleotide, dinucleotide and trinucleotide SSRs, respectively, and the motif frequencies of tetra- to hexanucleotide SSRs appear to be species-specific. The 10-20 bp length of SSRs are dominant with the number of 110 561 ± 93 482 and the frequency of 87.25% ± 5.73% on average, and the number and frequency decline with the increase of length. Most SSRs (83.34% ± 7.72%) are located in intergenic regions, followed by intron regions (11.59% ± 5.59%), exon regions (3.74% ± 1.95%), and untranslated regions (1.32% ± 1.39%). The mono-, di- and trinucleotide SSRs are the main SSRs in both gene regions (98.55% ± 0.85%) and exon regions (99.27% ± 0.52%). An average of 42.52% of total genes contains SSRs, and the preference for SSR occurrence in different gene subcategories are species-specific. The study provides useful insights into the SSR diversity, characteristics and distribution in 23 mosquito species of genomes.

Characterization of Duck (Anas platyrhynchos) Short Tandem Repeat Variation by Population-Scale Genome Resequencing

Short tandem repeats (STRs) are usually associated with genetic diseases and gene regulatory functions, and are also important genetic markers for analysis of evolutionary, genetic diversity and forensic. However, for the majority of STRs in the duck genome, their population genetic properties and functional impacts remain poorly defined. Recent advent of next generation sequencing (NGS) has offered an opportunity for profiling large numbers of polymorphic STRs. Here, we reported a population-scale analysis of STR variation using genome resequencing in mallard and Pekin duck. Our analysis provided the first genome-wide duck STR reference including 198,022 STR loci with motif size of 2–6 base pairs. We observed a relatively uneven distribution of STRs in different genomic regions, which indicates that the occurrence of STRs in duck genome is not random, but undergoes a directional selection pressure. Using genome resequencing data of 23 mallard and 26 Pekin ducks, we successfully identified 89,891 polymorphic STR loci. Intensive analysis of this dataset suggested that shorter repeat motif, longer reference tract length, higher purity, and residing outside of a coding region are all associated with an increase in STR variability. STR genotypes were utilized for population genetic analysis, and the results showed that population structure and divergence patterns among population groups can be efficiently captured. In addition, comparison between Pekin duck and mallard identified 3,122 STRs with extremely divergent allele frequency, which overlapped with a set of genes related to nervous system, energy metabolism and behavior. The evolutionary analysis revealed that the genes containing divergent STRs may play important roles in phenotypic changes during duck domestication. The variation analysis of STRs in population scale provides valuable resource for future study of genetic diversity and genome evolution in duck.

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.

Population diversity and relatedness in Sugarbirds (Promeropidae: Promerops spp.)

Sugarbirds are a family of two socially-monogamous passerine species endemic to southern Africa. Cape and Gurney’s Sugarbird (Promerops cafer and P. gurneyi) differ in abundance, dispersion across their range and in the degree of sexual dimorphism in tail length, factors that affect breeding systems and potentially genetic diversity. According to recent data, P. gurneyi are in decline and revision of the species’ IUCN conservation status to a threatened category may be warranted. It is therefore necessary to understand genetic diversity and risk of inbreeding in this species. We used six polymorphic microsatellite markers and one mitochondrial gene (ND2) to compare genetic diversity in P. cafer from Helderberg Nature Reserve and P. gurneyi from Golden Gate Highlands National Park, sites at the core of each species distribution. We describe novel universal avian primers which amplify the entire ND2 coding sequence across a broad range of bird orders. We observed high mitochondrial and microsatellite diversity in both sugarbird populations, with no detectable inbreeding and large effective population sizes.

Genome-wide mining of perfect microsatellites and tetranucleotide orthologous microsatellites estimates in six primate species

Advancement in genome sequencing and in silico mining tools have provided new opportunities for comparative primate genomics of microsatellites. The SSRs (simple sequence repeats) numbers were not correlated with the genome size (Pearson, r=0.310, p=0.550), and were positively correlated with the total length of SSRs (Pearson, r=0.992, p=0.00). A total of 224,289 tetranucleotide orthologous microsatellites families and 367 single-copy orthologous SSRs loci were found in six primate species by homologous alignment. The inner mutation types of single-copy orthologous SSRs loci included the copy number variance, point mutation, and chromosomal translocation. The accumulated repeat times and average length of tetranucleotide orthologous microsatellites in Rhinopithecus roxellana, Papio anubis and Macaca mulatta were longer than Homo sapiens and Pan troglodytes, which showed the tetranucleotide orthologous SSRs loci had more repeat times and longer average length on the branches with earlier divergence time, one exception may be Microcebus murinus as a primitive monkey with a smallest morphology in Malagasy. Our conclusion indicated that single-copy tetranucleotide orthologous SSRs sequences accumulated individual mutation more slowly through time in H. sapiens and P. troglodytes than in R. roxellanae, P. anubis and M. mulatta. However, such divergence wouldn't arise uniformly in all branches of the primate tree. A comparison of genomic sequence assemblages would offer remarkable insights about comparisons and contrasts, and the evolutionary processes of the microsatellites involved in human and nonhuman primate species.