Analysis of Microsatellite Variation in Drosophila melanogaster with Population-Scale Genome Sequencing

Genome sequencing technologies promise to revolutionize our understanding of genetics, evolution, and disease by making it feasible to survey a broad spectrum of sequence variation on a population scale. However, this potential can only be realized to the extent that methods for extracting and interpreting distinct forms of variation can be established. The error profiles and read length limitations of early versions of next-generation sequencing technologies rendered them ineffective for some sequence variant types, particularly microsatellites and other tandem repeats, and fostered the general misconception that such variants are inherently inaccessible to these platforms. At the same time, tandem repeats have emerged as important sources of functional variation. Tandem repeats are often located in and around genes, and frequent mutations in their lengths exert quantitative effects on gene function and phenotype, rapidly degrading linkage disequilibrium between markers and traits. Sensitive identification of these variants in large-scale next-gen sequencing efforts will enable more comprehensive association studies capable of revealing previously invisible associations. We present a population-scale analysis of microsatellite repeats using whole-genome data from 158 inbred isolates from the Drosophila Genetics Reference Panel, a collection of over 200 extensively phenotypically characterized isolates from a single natural population, to uncover processes underlying repeat mutation and to enable associations with behavioral, morphological, and life-history traits. Analysis of repeat variation from next-generation sequence data will also enhance studies of genome stability and neurodegenerative diseases

The nuclear receptors of Biomphalaria glabrata and Lottia gigantea: Implications for developing new model organisms

© 2015 Kaur et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are creditedNuclear receptors (NRs) are transcription regulators involved in an array of diverse physiological functions including key roles in endocrine and metabolic function. The aim of this study was to identify nuclear receptors in the fully sequenced genome of the gastropod snail, Biomphalaria glabrata, intermediate host for Schistosoma mansoni and compare these to known vertebrate NRs, with a view to assessing the snail's potential as a invertebrate model organism for endocrine function, both as a prospective new test organism and to elucidate the fundamental genetic and mechanistic causes of disease. For comparative purposes, the genome of a second gastropod, the owl limpet, Lottia gigantea was also investigated for nuclear receptors. Thirty-nine and thirty-three putative NRs were identified from the B. glabrata and L. gigantea genomes respectively, based on the presence of a conserved DNA-binding domain and/or ligand-binding domain. Nuclear receptor transcript expression was confirmed and sequences were subjected to a comparative phylogenetic analysis, which demonstrated that these molluscs have representatives of all the major NR subfamilies (1-6). Many of the identified NRs are conserved between vertebrates and invertebrates, however differences exist, most notably, the absence of receptors of Group 3C, which includes some of the vertebrate endocrine hormone targets. The mollusc genomes also contain NR homologues that are present in insects and nematodes but not in vertebrates, such as Group 1J (HR48/DAF12/HR96). The identification of many shared receptors between humans and molluscs indicates the potential for molluscs as model organisms; however the absence of several steroid hormone receptors indicates snail endocrine systems are fundamentally different.The National Centre for the Replacement, Refinement and Reduction of Animals in Research, Grant Ref:G0900802 to CSJ, LRN, SJ & EJR [www.nc3rs.org.uk]

Characterization of simple sequence repeats (SSRs) from Phlebotomus papatasi (Diptera: Psychodidae) expressed sequence tags (ESTs)

<p>Abstract</p> <p>Background</p> <p><it>Phlebotomus papatasi </it>is a natural vector of <it>Leishmania major</it>, which causes cutaneous leishmaniasis in many countries. Simple sequence repeats (SSRs), or microsatellites, are common in eukaryotic genomes and are short, repeated nucleotide sequence elements arrayed in tandem and flanked by non-repetitive regions. The enrichment methods used previously for finding new microsatellite loci in sand flies remain laborious and time consuming; <it>in silico </it>mining, which includes retrieval and screening of microsatellites from large amounts of sequence data from sequence data bases using microsatellite search tools can yield many new candidate markers.</p> <p>Results</p> <p>Simple sequence repeats (SSRs) were characterized in <it>P. papatasi </it>expressed sequence tags (ESTs) derived from a public database, National Center for Biotechnology Information (NCBI). A total of 42,784 sequences were mined, and 1,499 SSRs were identified with a frequency of 3.5% and an average density of 15.55 kb per SSR. Dinucleotide motifs were the most common SSRs, accounting for 67% followed by tri-, tetra-, and penta-nucleotide repeats, accounting for 31.1%, 1.5%, and 0.1%, respectively. The length of microsatellites varied from 5 to 16 repeats. Dinucleotide types; AG and CT have the highest frequency. Dinucleotide SSR-ESTs are relatively biased toward an excess of (AX)n repeats and a low GC base content. Forty primer pairs were designed based on motif lengths for further experimental validation.</p> <p>Conclusion</p> <p>The first large-scale survey of SSRs derived from <it>P. papatasi </it>is presented; dinucleotide SSRs identified are more frequent than other types. EST data mining is an effective strategy to identify functional microsatellites in <it>P. papatasi</it>.</p

Simple sequence repeats in Neurospora crassa: distribution, polymorphism and evolutionary inference

<p>Abstract</p> <p>Background</p> <p>Simple sequence repeats (SSRs) have been successfully used for various genetic and evolutionary studies in eukaryotic systems. The eukaryotic model organism <it>Neurospora crassa </it>is an excellent system to study evolution and biological function of SSRs.</p> <p>Results</p> <p>We identified and characterized 2749 SSRs of 963 SSR types in the genome of <it>N. crassa</it>. The distribution of tri-nucleotide (nt) SSRs, the most common SSRs in <it>N. crassa</it>, was significantly biased in exons. We further characterized the distribution of 19 abundant SSR types (AST), which account for 71% of total SSRs in the <it>N. crassa </it>genome, using a Poisson log-linear model. We also characterized the size variation of SSRs among natural accessions using Polymorphic Index Content (PIC) and ANOVA analyses and found that there are genome-wide, chromosome-dependent and local-specific variations. Using polymorphic SSRs, we have built linkage maps from three line-cross populations.</p> <p>Conclusion</p> <p>Taking our computational, statistical and experimental data together, we conclude that 1) the distributions of the SSRs in the sequenced N. crassa genome differ systematically between chromosomes as well as between SSR types, 2) the size variation of tri-nt SSRs in exons might be an important mechanism in generating functional variation of proteins in <it>N. crassa</it>, 3) there are different levels of evolutionary forces in variation of amino acid repeats, and 4) SSRs are stable molecular markers for genetic studies in <it>N. crassa</it>.</p

Analyses of carnivore microsatellites and their intimate association with tRNA-derived SINEs

BACKGROUND: The popularity of microsatellites has greatly increased in the last decade on account of their many applications. However, little is currently understood about the factors that influence their genesis and distribution among and within species genomes. In this work, we analyzed carnivore microsatellite clones from GenBank to study their association with interspersed repeats and elucidate the role of the latter in microsatellite genesis and distribution. RESULTS: We constructed a comprehensive carnivore microsatellite database comprising 1236 clones from GenBank. Thirty-three species of 11 out of 12 carnivore families were represented, although two distantly related species, the domestic dog and cat, were clearly overrepresented. Of these clones, 330 contained tRNA(Lys)-derived SINEs and 357 contained other interspersed repeats. Our rough estimates of tRNA SINE copies per haploid genome were much higher than published ones. Our results also revealed a distinct juxtaposition of AG and A-rich repeats and tRNA(Lys)-derived SINEs suggesting their coevolution. Both microsatellites arose repeatedly in two regions of the insterspersed repeat. Moreover, microsatellites associated with tRNA(Lys)-derived SINEs showed the highest complexity and less potential instability. CONCLUSION: Our results suggest that tRNA(Lys)-derived SINEs are a significant source for microsatellite generation in carnivores, especially for AG and A-rich repeat motifs. These observations indicate two modes of microsatellite generation: the expansion and variation of pre-existing tandem repeats and the conversion of sequences with high cryptic simplicity into a repeat array; mechanisms which are not specific to tRNA(Lys)-derived SINEs. Microsatellite and interspersed repeat coevolution could also explain different distribution of repeat types among and within species genomes. Finally, due to their higher complexity and lower potential informative content of microsatellites associated with tRNA(Lys)-derived SINEs, we recommend avoiding their use as genetic markers

PDK1 and HR46 Gene Homologs Tie Social Behavior to Ovary Signals

The genetic basis of division of labor in social insects is a central question in evolutionary and behavioral biology. The honey bee is a model for studying evolutionary behavioral genetics because of its well characterized age-correlated division of labor. After an initial period of within-nest tasks, 2–3 week-old worker bees begin foraging outside the nest. Individuals often specialize by biasing their foraging efforts toward collecting pollen or nectar. Efforts to explain the origins of foraging specialization suggest that division of labor between nectar and pollen foraging specialists is influenced by genes with effects on reproductive physiology. Quantitative trait loci (QTL) mapping of foraging behavior also reveals candidate genes for reproductive traits. Here, we address the linkage of reproductive anatomy to behavior, using backcross QTL analysis, behavioral and anatomical phenotyping, candidate gene expression studies, and backcross confirmation of gene-to-anatomical trait associations. Our data show for the first time that the activity of two positional candidate genes for behavior, PDK1 and HR46, have direct genetic relationships to ovary size, a central reproductive trait that correlates with the nectar and pollen foraging bias of workers. These findings implicate two genes that were not known previously to influence complex social behavior. Also, they outline how selection may have acted on gene networks that affect reproductive resource allocation and behavior to facilitate the evolution of social foraging in honey bees

Parallel Evolution of Auditory Genes for Echolocation in Bats and Toothed Whales

The ability of bats and toothed whales to echolocate is a remarkable case of convergent evolution. Previous genetic studies have documented parallel evolution of nucleotide sequences in Prestin and KCNQ4, both of which are associated with voltage motility during the cochlear amplification of signals. Echolocation involves complex mechanisms. The most important factors include cochlear amplification, nerve transmission, and signal re-coding. Herein, we screen three genes that play different roles in this auditory system. Cadherin 23 (Cdh23) and its ligand, protocadherin 15 (Pcdh15), are essential for bundling motility in the sensory hair. Otoferlin (Otof) responds to nerve signal transmission in the auditory inner hair cell. Signals of parallel evolution occur in all three genes in the three groups of echolocators—two groups of bats (Yangochiroptera and Rhinolophoidea) plus the dolphin. Significant signals of positive selection also occur in Cdh23 in the Rhinolophoidea and dolphin, and Pcdh15 in Yangochiroptera. In addition, adult echolocating bats have higher levels of Otof expression in the auditory cortex than do their embryos and non-echolocation bats. Cdh23 and Pcdh15 encode the upper and lower parts of tip-links, and both genes show signals of convergent evolution and positive selection in echolocators, implying that they may co-evolve to optimize cochlear amplification. Convergent evolution and expression patterns of Otof suggest the potential role of nerve and brain in echolocation. Our synthesis of gene sequence and gene expression analyses reveals that positive selection, parallel evolution, and perhaps co-evolution and gene expression affect multiple hearing genes that play different roles in audition, including voltage and bundle motility in cochlear amplification, nerve transmission, and brain function

Population genomics applications for conservation: the case of the tropical dry forest dweller Peromyscus melanophrys

Recent advances in genomic sequencing have opened new horizons in the study of population genetics and evolution in non-model organisms. However, very few population genomic studies have been performed on wild mammals to understand how the landscape affects the genetic structure of populations, useful information for the conservation of biodiversity. Here, we applied a genomic approach to evaluate the relationship between habitat features and genetic patterns at spatial and temporal scales in an endangered ecosystem, the Tropical Dry Forest (TDF). We studied populations of the Plateau deer mouse Peromyscus melanophrys to analyse its genomic diversity and structure in a TDF protected area in the Huautla Mountain Range (HMR), Mexico based on 8,209 SNPs obtained through Genotyping-by-Sequencing. At a spatial scale, we found a significant signature of isolation-by-distance, few significant differences in genetic diversity indices among study sites, and no significant differences between habitats with different levels of human perturbation. At a temporal scale, while genetic diversity levels fluctuated significantly over time, neither seasonality nor disturbance levels had a significant effect. Also, outlier analysis revealed loci potentially under selection. Our results suggest that the population genetics of P. melanophrys may be little impacted by anthropogenic disturbances, or by natural spatial and temporal habitat heterogeneity in our study area. The genome-wide approach adopted here provides data of value for conservation planning, and a baseline to be used as a reference for future studies on the effects of habitat fragmentation and seasonality in the HMR and in TDF

Intronic Cis-Regulatory Modules Mediate Tissue-Specific and Microbial Control of angptl4/fiaf Transcription

The intestinal microbiota enhances dietary energy harvest leading to increased fat storage in adipose tissues. This effect is caused in part by the microbial suppression of intestinal epithelial expression of a circulating inhibitor of lipoprotein lipase called Angiopoietin-like 4 (Angptl4/Fiaf). To define the cis-regulatory mechanisms underlying intestine-specific and microbial control of Angptl4 transcription, we utilized the zebrafish system in which host regulatory DNA can be rapidly analyzed in a live, transparent, and gnotobiotic vertebrate. We found that zebrafish angptl4 is transcribed in multiple tissues including the liver, pancreatic islet, and intestinal epithelium, which is similar to its mammalian homologs. Zebrafish angptl4 is also specifically suppressed in the intestinal epithelium upon colonization with a microbiota. In vivo transgenic reporter assays identified discrete tissue-specific regulatory modules within angptl4 intron 3 sufficient to drive expression in the liver, pancreatic islet β-cells, or intestinal enterocytes. Comparative sequence analyses and heterologous functional assays of angptl4 intron 3 sequences from 12 teleost fish species revealed differential evolution of the islet and intestinal regulatory modules. High-resolution functional mapping and site-directed mutagenesis defined the minimal set of regulatory sequences required for intestinal activity. Strikingly, the microbiota suppressed the transcriptional activity of the intestine-specific regulatory module similar to the endogenous angptl4 gene. These results suggest that the microbiota might regulate host intestinal Angptl4 protein expression and peripheral fat storage by suppressing the activity of an intestine-specific transcriptional enhancer. This study provides a useful paradigm for understanding how microbial signals interact with tissue-specific regulatory networks to control the activity and evolution of host gene transcription