Promoter microsatellites as modulators of human gene expression

Short tandem repeats bind transcription factors to tune eukaryotic gene expression

Short tandem repeats (STRs) are enriched in eukaryotic cis-regulatory elements and alter gene expression, yet how they regulate transcription remains unknown. We found that STRs modulate transcription factor (TF)-DNA affinities and apparent on-rates by about 70-fold by directly binding TF DNA-binding domains, with energetic impacts exceeding many consensus motif mutations. STRs maximize the number of weakly preferred microstates near target sites, thereby increasing TF density, with impacts well predicted by statistical mechanics. Confirming that STRs also affect TF binding in cells, neural networks trained only on in vivo occupancies predicted effects identical to those observed in vitro. Approximately 90% of TFs preferentially bound STRs that need not resemble known motifs, providing a cis-regulatory mechanism to target TFs to genomic sites.

Short tandem repeat polymorphism in the promoter region of cyclophilin 19B drives its transcriptional upregulation and contributes to drug resistance in the malaria parasite Plasmodium falciparum

Resistance of the human malaria parasites, Plasmodium falciparum, to artemisinins is now fully established in Southeast Asia and is gradually emerging in Sub-Saharan Africa. Although nonsynonymous SNPs in the pfk13 Kelch-repeat propeller (KREP) domain are clearly associated with artemisinin resistance, their functional relevance requires cooperation with other genetic factors/alterations of the P. falciparum genome, collectively referred to as genetic background. Here we provide experimental evidence that P. falciparum cyclophilin 19B (PfCYP19B) may represent one putative factor in this genetic background, contributing to artemisinin resistance via its increased expression. We show that overexpression of PfCYP19B in vitro drives limited but significant resistance to not only artemisinin but also piperaquine, an important partner drug in artemisinin-based combination therapies. We showed that PfCYP19B acts as a negative regulator of the integrated stress response (ISR) pathway by modulating levels of phosphorylated eIF2α (eIF2α-P). Curiously, artemisinin and piperaquine affect eIF2α-P in an inverse direction that in both cases can be modulated by PfCYP19B towards resistance. Here we also provide evidence that the upregulation of PfCYP19B in the drug-resistant parasites appears to be maintained by a short tandem repeat (SRT) sequence polymorphism in the gene's promoter region. These results support a model that artemisinin (and other drugs) resistance mechanisms are complex genetic traits being contributed to by altered expression of multiple genes driven by genetic polymorphism at their promoter regions.

Population-level genome-wide STR discovery and validation for population structure and genetic diversity assessment of Plasmodium species

Short tandem repeats (STRs) are highly informative genetic markers that have been used extensively in population genetics analysis. They are an important source of genetic diversity and can also have functional impact. Despite the availability of bioinformatic methods that permit large-scale genome-wide genotyping of STRs from whole genome sequencing data, they have not previously been applied to sequencing data from large collections of malaria parasite field samples. Here, we have genotyped STRs using HipSTR in more than 3,000 Plasmodium falciparum and 174 Plasmodium vivax published whole-genome sequence data from samples collected across the globe. High levels of noise and variability in the resultant callset necessitated the development of a novel method for quality control of STR genotype calls. A set of high-quality STR loci (6,768 from P. falciparum and 3,496 from P. vivax) were used to study Plasmodium genetic diversity, population structures and genomic signatures of selection and these were compared to genome-wide single nucleotide polymorphism (SNP) genotyping data. In addition, the genome-wide information about genetic variation and other characteristics of STRs in P. falciparum and P. vivax have been available in an interactive web-based R Shiny application PlasmoSTR (https://github.com/bahlolab/PlasmoSTR).

Revisiting tandem repeats in psychiatric disorders from perspectives of genetics, physiology, and brain evolution

Genome-wide association studies (GWASs) have revealed substantial genetic components comprised of single nucleotide polymorphisms (SNPs) in the heritable risk of psychiatric disorders. However, genetic risk factors not covered by GWAS also play pivotal roles in these illnesses. Tandem repeats, which are likely functional but frequently overlooked by GWAS, may account for an important proportion in the "missing heritability" of psychiatric disorders. Despite difficulties in characterizing and quantifying tandem repeats in the genome, studies have been carried out in an attempt to describe impact of tandem repeats on gene regulation and human phenotypes. In this review, we have introduced recent research progress regarding the genomic distribution and regulatory mechanisms of tandem repeats. We have also summarized the current knowledge of the genetic architecture and biological underpinnings of psychiatric disorders brought by studies of tandem repeats. These findings suggest that tandem repeats, in candidate psychiatric risk genes or in different levels of linkage disequilibrium (LD) with psychiatric GWAS SNPs and haplotypes, may modulate biological phenotypes related to psychiatric disorders (e.g., cognitive function and brain physiology) through regulating alternative splicing, promoter activity, enhancer activity and so on. In addition, many tandem repeats undergo tight natural selection in the human lineage, and likely exert crucial roles in human brain evolution. Taken together, the putative roles of tandem repeats in the pathogenesis of psychiatric disorders is strongly implicated, and using examples from previous literatures, we wish to call for further attention to tandem repeats in the post-GWAS era of psychiatric disorders.

STRs: Ancient Architectures of the Genome beyond the Sequence

Short tandem repeats (STRs) are commonly defined as short runs of repetitive nucleotides, consisting of tandemly repeating 2-6- bp motif units, which are ubiquitously distributed throughout genomes. Functional STRs are polymorphic in the population, and their variations influence gene expression, which subsequently may result in pathogenic phenotypes. To understand STR phenotypic effects and their functional roles, we describe four different mutational mechanisms including the unequal crossing-over model, gene conversion, retrotransposition mechanism and replication slippage. Due to the multi-allelic nature, small length, abundance, high variability, codominant inheritance, nearly neutral evolution, extensive genome coverage and simple assaying of STRs, these markers are widely used in various types of biological research, including population genetics studies, genome mapping, molecular epidemiology, paternity analysis and gene flow studies. In this review, we focus on the current knowledge regarding STR genomic distribution, function, mutation and applications.

Germline microsatellite genotypes differentiate children with medulloblastoma

BACKGROUND

The germline genetic events underpinning medulloblastoma (MB) initiation, and therefore the ability to determine who is at risk, are still unknown for the majority of cases. Microsatellites are short repeated sequences that make up ~3% of the genome. Repeat lengths vary among individuals and are often nonrandomly associated with disease, including several cancers such as breast, glioma, lung, and ovarian. Due to their effects on gene function, they have been called the "tuning knobs of the genome."

METHODS

We have developed a novel approach for identifying a microsatellite-based signature to differentiate MB patients from controls using germline DNA.

RESULTS

Analyzing germline whole exome sequencing data from a training set of 120 MB subjects and 425 controls, we identified 139 individual microsatellite loci whose genotypes differ significantly between the groups. Using a genetic algorithm, we identified a subset of 43 microsatellites that distinguish MB subjects from controls with a sensitivity and specificity of 92% and 88%, respectively. This microsatellite signature was validated in an independent dataset consisting of 102 subjects and 428 controls, with comparable sensitivity and specificity of 95% and 90%, respectively. Analysis of the allele genotypes of those 139 informative loci demonstrates that their association with MB is a consequence of individual microsatellites' genotypes rather than their hypermutability. Finally, an analysis of the genes harboring these microsatellite loci reveals cellular functions important for tumorigenesis.

CONCLUSION

This study demonstrates that MB-specific germline microsatellite variations mark those at risk for MB development and suggests mechanisms of predisposition.

Abundance of ethnically biased microsatellites in human gene regions

Microsatellites-a type of short tandem repeat (STR)-have been used for decades as putatively neutral markers to study the genetic structure of diverse human populations. However, recent studies have demonstrated that some microsatellites contribute to gene expression, cis heritability, and phenotype. As a corollary, some microsatellites may contribute to differential gene expression and RNA/protein structure stability in distinct human populations. To test this hypothesis, we investigate genotype frequencies, functional relevance, and adaptive potential of microsatellites in five super-populations (ethnicities) drawn from the 1000 Genomes Project. We discover 3,984 ethnically-biased microsatellite loci (EBML); for each EBML at least one ethnicity has genotype frequencies statistically different from the remaining four. South Asian, East Asian, European, and American EBML show significant overlap; on the contrary, the set of African EBML is mostly unique. We cross-reference the 3,984 EBML with 2,060 previously identified expression STRs (eSTRs); repeats known to affect gene expression (64 total) are over-represented. The most significant pathway enrichments are those associated with the matrisome: a broad collection of genes encoding the extracellular matrix and its associated proteins. At least 14 of the EBML have established links to human disease. Analysis of the 3,984 EBML with respect to known selective sweep regions in the genome shows that allelic variation in some of them is likely associated with adaptive evolution.

Conserved microsatellites may contribute to stem-loop structures in 5', 3' terminals of Ebolavirus genomes

Microsatellites (SSRs) are ubiquitous in coding and non-coding regions of the Ebolavirus genomes. We synthetically analyzed the microsatellites in whole-genome and terminal regions of 219 Ebolavirus genomes from five species. The Ebolavirus sequences were observed with small intraspecies variations and large interspecific variations, especially in the terminal non-coding regions. Only five conserved microsatellites were detected in the complete genomes, and four of them which well base-paired to help forming conserved stem-loop structures mainly appeared in the terminal non-coding regions. These results suggest that the conserved microsatellites may be evolutionary selected to form conserved secondary structures in 5′, 3′ terminals of Ebolavirus genomes. It may help to understand the biological significance of microsatellites in Ebolavirus and also other virus genomes.

•

Conserved microsatellites mainly occurred in 5′, 3′ terminal non-coding regions.

•

Conserved microsatellites may contribute to conserved stem-loop structures.

•

Conserved microsatellites might be preserved under greater evolutionary pressure.

CAGm: a repository of germline microsatellite variations in the 1000 genomes project

The human genome harbors an abundance of repetitive DNA; however, its function continues to be debated. Microsatellites-a class of short tandem repeat-are established as an important source of genetic variation. Array length variants are common among microsatellites and affect gene expression; but, efforts to understand the role and diversity of microsatellite variation has been hampered by several challenges. Without adequate depth, both long-read and short-read sequencing may not detect the variants present in a sample; additionally, large sample sizes are needed to reveal the degree of population-level polymorphism. To address these challenges we present the Comparative Analysis of Germline Microsatellites (CAGm): a database of germline microsatellites from 2529 individuals in the 1000 genomes project. A key novelty of CAGm is the ability to aggregate microsatellite variation by population, ethnicity (super population) and gender. The database provides advanced searching for microsatellites embedded in genes and functional elements. All data can be downloaded as Microsoft Excel spreadsheets. Two use-case scenarios are presented to demonstrate its utility: a mononucleotide (A) microsatellite at the BAT-26 locus and a dinucleotide (CA) microsatellite in the coding region of FGFRL1. CAGm is freely available at http://www.cagmdb.org/.

Regulatory Architecture of the Neuronal Cacng2/Tarpγ2 Gene Promoter: Multiple Repressive Domains, a Polymorphic Regulatory Short Tandem Repeat, and Bidirectional Organization with Co-regulated lncRNAs

CACNG2 (TARPγ2, Stargazin) is a multi-functional regulator of excitatory neurotransmission and has been implicated in the pathological processes of several brain diseases. Cacng2 function is dependent upon expression level, but currently, little is known about the molecular mechanisms that control expression of this gene. To address this deficit and investigate disease-related gene variants, we have cloned and characterized the rat Cacng2 promoter and have defined three major features: (i) multiple repressive domains that include an array of RE-1 silencing transcription factor (REST) elements, and a calcium regulatory element-binding factor (CaRF) element, (ii) a (poly-GA) short tandem repeat (STR), and (iii) bidirectional organization with expressed lncRNAs. Functional activity of the promoter was demonstrated in transfected neuronal cell lines (HT22 and PC12), but although selective removal of REST and CaRF domains was shown to enhance promoter-driven transcription, the enhanced Cacng2 promoter constructs were still about fivefold weaker than a comparable rat Synapsin-1 promoter sequence. Direct evidence of REST activity at the Cacng2 promoter was obtained through co-transfection with an established dominant-negative REST (DNR) construct. Investigation of the GA-repeat STR revealed polymorphism across both animal strains and species, and size variation was also observed in absence epilepsy disease model cohorts (Genetic Absence Epilepsy Rats, Strasbourg [GAERS] and non-epileptic control [NEC] rats). These data provide evidence of a genotype (STR)-phenotype correlation that may be unique with respect to proximal gene regulatory sequence in the demonstrated absence of other promoter, or 3' UTR variants in GAERS rats. However, although transcriptional regulatory activity of the STR was demonstrated in further transfection studies, we did not find a GAERS vs. NEC difference, indicating that this specific STR length variation may only be relevant in the context of other (Cacna1h and Kcnk9) gene variants in this disease model. Additional studies revealed further (bidirectional) complexity at the Cacng2 promoter, and we identified novel, co-regulated, antisense rat lncRNAs that are paired with Cacng2 mRNA. These studies have provided novel insights into the organization of a synaptic protein gene promoter, describing multiple repressive and modulatory domains that can mediate diverse regulatory inputs.

Identification of an (AC)n microsatellite in the Six1 gene promoter and its effect on production traits in Pietrain × Duroc × Landrace × Yorkshire pigs

The Sine oculis homeobox 1 (Six1) gene is important for skeletal muscle growth and fiber specification; therefore, it is considered as a promising candidate gene that may influence porcine growth and meat quality traits. Nevertheless, the association of Six1 with these processes and the mechanisms regulating its expression remain unclear. The objectives of this study were to identify variant sites of Six1 in different pig breeds, conduct association analysis to evaluate the relationship between polymorphisms of these variants and porcine production traits in Pietrain × Duroc × Landrace × Yorkshire commercial pigs, and explore the potential regulatory mechanisms of Six1 affecting production traits. A total of 12 variants were identified, including 10 single- nucleotide variations (SNVs), 1 insertion- deletion (Indel), and 1 (AC)n microsatellite. Association analysis demonstrated that the SNV, g.1595A>G, was significantly associated with meat color (redness, a*); individuals with the G allele had greater a* values (P < 0.05). Moreover, our results demonstrated that the (AC)n polymorphism in the Six1 promoter was significantly associated with weaning weight (P < 0.05), carcass weight (P < 0.05), and thoracic and lumbar back fat (P < 0.01).In addition, we found that the (AC)n variant was closely related with Six1 expression levels and demonstrated this polymorphism on promoter activity by in vitro experiments. Overall, this study provides novel evidence for elucidating the effects of Six1 on porcine production traits as promising candidate and describes two variants with these traits, which are potential reference markers for pig molecular breeding. In addition, our data on the relationship between porcine Six1 expression and the polymorphic (AC)n microsatellite in its promoter may facilitate similar studies in other species.

Of Men and Mice: Modeling the Fragile X Syndrome

The Fragile X Syndrome (FXS) is one of the most common forms of inherited intellectual disability in all human societies. Caused by the transcriptional silencing of a single gene, the fragile x mental retardation gene FMR1, FXS is characterized by a variety of symptoms, which range from mental disabilities to autism and epilepsy. More than 20 years ago, a first animal model was described, the Fmr1 knock-out mouse. Several other models have been developed since then, including conditional knock-out mice, knock-out rats, a zebrafish and a drosophila model. Using these model systems, various targets for potential pharmaceutical treatments have been identified and many treatments have been shown to be efficient in preclinical studies. However, all attempts to turn these findings into a therapy for patients have failed thus far. In this review, I will discuss underlying difficulties and address potential alternatives for our future research.

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.

Tandem repeats mediating genetic plasticity in health and disease

Accumulating evidence suggests that many classes of DNA repeats exhibit attributes that distinguish them from other genetic variants, including the fact that they are more liable to mutation; this enables them to mediate genetic plasticity. The expansion of tandem repeats, particularly of short tandem repeats, can cause a range of disorders (including Huntington disease, various ataxias, motor neuron disease, frontotemporal dementia, fragile X syndrome and other neurological disorders), and emerging data suggest that tandem repeat polymorphisms (TRPs) can also regulate gene expression in healthy individuals. TRPs in human genomes may also contribute to the missing heritability of polygenic disorders. A better understanding of tandem repeats and their associated repeatome, as well as their capacity for genetic plasticity via both germline and somatic mutations, is needed to transform our understanding of the role of TRPs in health and disease.

ZDHHC3 as a Risk and Mortality Marker for Breast Cancer in African American Women

African American woman are 43% more likely to die from breast cancer than white women and have increased the risk of tumor recurrence despite lower incidence. We investigate variations in microsatellite genomic regions-a type of repetitive DNA-and possible links to the breast cancer mortality gap. We screen 33 854 microsatellites in germline DNA of African American women with and without breast cancer: 4 are statistically significant. These are located in the 3' UTR (untranslated region) of gene ZDHHC3, an intron of transcribed pseudogene INTS4L1, an intron of ribosomal gene RNA5-8S5, and an intergenic region of chromosome 16. The marker in ZDHHC3 is interesting for 3 reasons: (a) the ZDHHC3 gene is located in region 3p21 which has already been linked to early invasive breast cancer, (b) the Kaplan-Meier estimator demonstrates that ZDHHC3 alterations are associated with poor breast cancer survival in all racial/ethnic groups combined, and (c) data from cBioPortal suggest that ZDHHC3 messenger RNA expression is significantly lower in African Americans compared with whites. These independent lines of evidence make ZDHHC3 a candidate for further investigation.

Functional Mechanisms of Microsatellite DNA in Eukaryotic Genomes

Microsatellite repeat DNA is best known for its length mutability, which is implicated in several neurological diseases and cancers, and often exploited as a genetic marker. Less well-known is the body of work exploring the widespread and surprisingly diverse functional roles of microsatellites. Recently, emerging evidence includes the finding that normal microsatellite polymorphism contributes substantially to the heritability of human gene expression on a genome-wide scale, calling attention to the task of elucidating the mechanisms involved. At present, these are underexplored, but several themes have emerged. I review evidence demonstrating roles for microsatellites in modulation of transcription factor binding, spacing between promoter elements, enhancers, cytosine methylation, alternative splicing, mRNA stability, selection of transcription start and termination sites, unusual structural conformations, nucleosome positioning and modification, higher order chromatin structure, noncoding RNA, and meiotic recombination hot spots.

Genetic variations at the human growth hormone receptor (GHR) gene locus are associated with idiopathic short stature

GH plays an essential role in the growing child by binding to the growth hormone receptor (GHR) on target cells and regulating multiple growth promoting and metabolic effects. Mutations in the GHR gene coding regions result in GH insensitivity (dwarfism) due to a dysfunctional receptor protein. However, children with idiopathic short stature (ISS) show growth impairment without GH or GHR defects. We hypothesized that decreased expression of the GHR gene may be involved. To test this, we investigated whether common genetic variants (microsatellites, SNPs) in regulatory regions of the GHR gene region were associated with the ISS phenotype. Genotyping of a GT‐repeat microsatellite in the GHR 5′UTR in a Montreal ISS cohort (n = 37 ISS, n = 105 controls) revealed that the incidence of the long/short (L/S) genotype was 3.3× higher in ISS children than controls (P = 0.04, OR = 3.85). In an Italian replication cohort (n = 143 ISS, n = 282 controls), the medium/short (M/S) genotype was 1.9× more frequent in the male ISS than controls (P = 0.017, OR = 2.26). In both ISS cohorts, logistic regression analysis of 27 SNPs showed an association of ISS with rs4292454, while haplotype analysis revealed specific risk haplotypes in the 3′ haploblocks. In contrast, there were no differences in GT genotype frequencies in a cohort of short stature (SS) adults versus controls (CARTaGENE: n = 168 SS, n = 207 controls) and the risk haplotype in the SS cohort was located in the most 5′ haploblock. These data suggest that the variants identified are potentially genetic markers specifically associated with the ISS phenotype.

Microsatellite polymorphisms associated with human behavioural and psychological phenotypes including a gene-environment interaction

Background

The genetic and environmental influences on human personality and behaviour are a complex matter of ongoing debate. Accumulating evidence indicates that short tandem repeats (STRs) in regulatory regions are good candidates to explain heritability not accessed by genome-wide association studies.

Methods

We tested for associations between the genotypes of four selected repeats and 18 traits relating to personality, behaviour, cognitive ability and mental health in a well-studied longitudinal birth cohort (n = 458-589) using one way analysis of variance. The repeats were a highly conserved poly-AC microsatellite in the upstream promoter region of the T-box brain 1 (TBR1) gene and three previously studied STRs in the activating enhancer-binding protein 2-beta (AP2-β) and androgen receptor (AR) genes. Where significance was found we used multiple regression to assess the influence of confounding factors.

Results

Carriers of the shorter, most common, allele of the AR gene’s GGN microsatellite polymorphism had fewer anxiety-related symptoms, which was consistent with previous studies, but in our study this was not significant following Bonferroni correction. No associations with two repeats in the AP2-β gene withstood this correction. A novel finding was that carriers of the minor allele of the TBR1 AC microsatellite were at higher risk of conduct problems in childhood at age 7-9 (p = 0.0007, which did pass Bonferroni correction). Including maternal smoking during pregnancy (MSDP) in models controlling for potentially confounding influences showed that an interaction between TBR1 genotype and MSDP was a significant predictor of conduct problems in childhood and adolescence (p < 0.001), and of self-reported criminal behaviour up to age 25 years (p ≤ 0.02). This interaction remained significant after controlling for possible confounders including maternal age at birth, socio-economic status and education, and offspring birth weight.

Conclusions

The potential functional importance of the TBR1 gene’s promoter microsatellite deserves further investigation. Our results suggest that it participates in a gene-environment interaction with MDSP and antisocial behaviour. However, previous evidence that mothers who smoke during pregnancy carry genes for antisocial behaviour suggests that epistasis may influence the interaction.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-017-0374-y) contains supplementary material, which is available to authorized users.

The mechanism of transactivation regulation due to polymorphic short tandem repeats (STRs) using IGF1 promoter as a model

Functional short tandem repeats (STR) are polymorphic in the population, and the number of repeats regulates the expression of nearby genes (known as expression STR, eSTR). STR in IGF1 promoter has been extensively studied for its association with IGF1 concentration in blood and various clinical traits and represents an important eSTR. We previously used an in-vitro luciferase reporter model to examine the interaction between STRs and SNPs in IGF1 promoter. Here, we further explored the mechanism how the number of repeats of the STR regulates gene transcription. An inverse correlation between the number of repeats and the extent of transactivation was found in a haplotype consisting of three promoter SNPs (C-STR-T-T). We showed that these adjacent SNPs located outside the STR were required for the STR to function as eSTR. The C allele of rs35767 provides a binding site for CCAAT/enhancer-binding-protein δ (C/EBPD), which is essential for the gradational transactivation property of eSTR and FOXA3 may also be involved. Therefore, we propose a mechanism in which the gradational transactivation by the eSTR is caused by the interaction of one or more transcriptional complexes located outside the STR, rather than by direct binding to a repeat motif of the STR.

Single Amino Acid Repeats in the Proteome World: Structural, Functional, and Evolutionary Insights

Microsatellites or simple sequence repeats (SSR) are abundant, highly diverse stretches of short DNA repeats present in all genomes. Tandem mono/tri/hexanucleotide repeats in the coding regions contribute to single amino acids repeats (SAARs) in the proteome. While SSRs in the coding region always result in amino acid repeats, a majority of SAARs arise due to a combination of various codons representing the same amino acid and not as a consequence of SSR events. Certain amino acids are abundant in repeat regions indicating a positive selection pressure behind the accumulation of SAARs. By analysing 22 proteomes including the human proteome, we explored the functional and structural relationship of amino acid repeats in an evolutionary context. Only ~15% of repeats are present in any known functional domain, while ~74% of repeats are present in the disordered regions, suggesting that SAARs add to the functionality of proteins by providing flexibility, stability and act as linker elements between domains. Comparison of SAAR containing proteins across species reveals that while shorter repeats are conserved among orthologs, proteins with longer repeats, >15 amino acids, are unique to the respective organism. Lysine repeats are well conserved among orthologs with respect to their length and number of occurrences in a protein. Other amino acids such as glutamic acid, proline, serine and alanine repeats are generally conserved among the orthologs with varying repeat lengths. These findings suggest that SAARs have accumulated in the proteome under positive selection pressure and that they provide flexibility for optimal folding of functional/structural domains of proteins. The insights gained from our observations can help in effective designing and engineering of proteins with novel features.

Resolving complex tandem repeats with long reads

MOTIVATION

Resolving tandemly repeated genomic sequences is a necessary step in improving our understanding of the human genome. Short tandem repeats (TRs), or microsatellites, are often used as molecular markers in genetics, and clinically, variation in microsatellites can lead to genetic disorders like Huntington's diseases. Accurately resolving repeats, and in particular TRs, remains a challenging task in genome alignment, assembly and variation calling. Though tools have been developed for detecting microsatellites in short-read sequencing data, these are limited in the size and types of events they can resolve. Single-molecule sequencing technologies may potentially resolve a broader spectrum of TRs given their increased length, but require new approaches given their significantly higher raw error profiles. However, due to inherent error profiles of the single-molecule technologies, these reads presents a unique challenge in terms of accurately identifying and estimating the TRs.

RESULTS

Here we present PacmonSTR, a reference-based probabilistic approach, to identify the TR region and estimate the number of these TR elements in long DNA reads. We present a multistep approach that requires as input, a reference region and the reference TR element. Initially, the TR region is identified from the long DNA reads via a 3-stage modified Smith-Waterman approach and then, expected number of TR elements is calculated using a pair-Hidden Markov Models-based method. Finally, TR-based genotype selection (or clustering: homozygous/heterozygous) is performed with Gaussian mixture models, using the Akaike information criteria, and coverage expectations.

Genome wide survey and analysis of small repetitive sequences in caulimoviruses

Microsatellites are known to exhibit ubiquitous presence across all kingdoms of life including viruses. Members of the Caulimoviridae family severely affect growth of vegetable and fruit plants and reduce economic yield in diverse cropping systems worldwide. Here, we analyzed the nature and distribution of both simple and complex microsatellites present in complete genome of 44 species of Caulimoviridae. Our results showed, in all analyzed genomes, genome size and GC content had a weak influence on number, relative abundance and relative density of microsatellites, respectively. For each genome, mono- and dinucleotide repeats were found to be highly predominant and are overrepresented in genome of majority of caulimoviruses. AT/TA and GAA/AAG/AGA was the most abundant di- and trinucleotide repeat motif, respectively. Repeats larger than trinucleotide were rarely found in these genomes. Comparative study of occurrence, abundance and density of microsatellite among available RNA and DNA viral genomes indicated that simple repeats were least abundant in genomes of caulimoviruses. Polymorphic repeats even though rare were observed in the large intergenic region of the genome, indicating strand slippage and/or unequal recombination processes do occur in caulimoviruses. To our knowledge, this is the first analysis of microsatellites occurring in any dsDNA viral genome. Characterization of such variations in repeat sequences would be important in deciphering the origin, mutational processes, and role of repeat sequences in viral genomes.

Microsatellite tandem repeats are abundant in human promoters and are associated with regulatory elements

Tandem repeats are genomic elements that are prone to changes in repeat number and are thus often polymorphic. These sequences are found at a high density at the start of human genes, in the gene’s promoter. Increasing empirical evidence suggests that length variation in these tandem repeats can affect gene regulation. One class of tandem repeats, known as microsatellites, rapidly alter in repeat number. Some of the genetic variation induced by microsatellites is known to result in phenotypic variation. Recently, our group developed a novel method for measuring the evolutionary conservation of microsatellites, and with it we discovered that human microsatellites near transcription start sites are often highly conserved. In this study, we examined the properties of microsatellites found in promoters. We found a high density of microsatellites at the start of genes. We showed that microsatellites are statistically associated with promoters using a wavelet analysis, which allowed us to test for associations on multiple scales and to control for other promoter related elements. Because promoter microsatellites tend to be G/C rich, we hypothesized that G/C rich regulatory elements may drive the association between microsatellites and promoters. Our results indicate that CpG islands, G-quadruplexes (G4) and untranslated regulatory regions have highly significant associations with microsatellites, but controlling for these elements in the analysis does not remove the association between microsatellites and promoters. Due to their intrinsic lability and their overlap with predicted functional elements, these results suggest that many promoter microsatellites have the potential to affect human phenotypes by generating mutations in regulatory elements, which may ultimately result in disease. We discuss the potential functions of human promoter microsatellites in this context.