(TG/CA)n repeats in human gene families: abundance and selective patterns of distribution according to function and gene length

ZSCAN4-binding motif-TGCACAC is conserved and enriched in CA/TG microsatellites in both mouse and human genomes

The Zinc finger and SCAN domain containing 4 (ZSCAN4) protein, expressed transiently in pluripotent stem cells, gametes, and early embryos, extends telomeres, enhances genome stability, and improves karyotypes in mouse embryonic stem (mES) cells. To gain insights into the mechanism of ZSCAN4 function, we identified genome-wide binding sites of endogenous ZSCAN4 protein using ChIP-seq technology in mouse and human ES cells, where the expression of endogenous ZSCAN4 was induced by treating cells with retinoic acids or by overexpressing DUX4. We revealed that both mouse and human ZSCAN4 bind to the TGCACAC motif located in CA/TG microsatellite repeats, which are known to form unstable left-handed duplexes called Z-DNA that can induce double-strand DNA breaks and mutations. These ZSCAN4 binding sites are mostly located in intergenic and intronic regions of the genomes. By generating ZSCAN4 knockout in human ES cells, we showed that ZSCAN4 does not seem to be involved in transcriptional regulation. We also found that ectopic expression of mouse ZSCAN4 enhances the suppression of chromatin at ZSCAN4-binding sites. These results together suggest that some of the ZSCAN4 functions are mediated by binding to the error-prone regions in mouse and human genomes.

Associations of BCL2 CA-Repeat Polymorphism and Breast Cancer Susceptibility in Isfahan Province of Iran

BCL2 apoptosis regulator (BCL2) is a cause of tumorigenesis whose CA-repeat promoter polymorphisms has inconsistent association with various types of cancers. The association of BCL2 polymorphism with breast cancer was investigated in the Isfahan province of Iran. PCRamplification of the CA-repeat was followed by polyacrylamide gel electrophoresis and direct sequencing for 120 breast cancer women and an equal number of corresponding healthy control individuals. Seven different alleles, ranging from 11 to 17 CA-repeats were observed. Short alleles with 11 to 14 repeats were protective (OR 0.363, P = 0.001), but large alleles with 15 to 17 repeats were threatening against breast cancer development (OR 2.780, P = 0.001). Accordingly, genotypes with large alleles showed a higher risk of breast cancer development (OR 3.400, P = 0.004). ERS1\ERBB2 positive breast cancer patients, but not PGRpositive ones, showed protection against breast cancer (OR 0.405, OR 0.346 respectively). In conclusion, women with at least one large allele of BCL2 were 3.4 times at higher risk of breast cancer development in the Isfahan province of Iran.

A comparative survey of microsatellites among wild and domestic cat provides valuable resources for marker development

Information on the level and distribution of genetic variation is important for conservation plan of captive population of an endangered species such as tiger and cheetah. We assayed the frequency of microsatellites in the genomic and genic sequences of wild cats (Panthera tigris, Acinonyx jubatus) and compared it with the domestic cat (Felis catus). Frequency, relative abundance and density of microsatellites were highest in the domestic cat when compared with wild cats. The frequency of microsatellites was positively correlated with the G+C content of genomic and genic sequences. The maximum frequency of microsatellites among all three sequence sets was of di-nucleotide repeats (genomic-88.1%; genic-70.4%), whereas the hexa-nucleotide repeat represents < 0.5%. Motif conservation study among the genomic and genic sequences revealed conservation of 81.3% and 51.0% motif within the members of family Felidae. A total 40,233 primers from genic sequences were designed in order to enrich the members of family Felidae with genomic resources. The designed primers could be useful in determining the molecular genetics of population structure and individualization of a particular cat.

Critical Mutation Rate has an Exponential Dependence on Population Size for Eukaryotic-length Genomes with Crossover

The critical mutation rate (CMR) determines the shift between survival-of-the-fittest and survival of individuals with greater mutational robustness ("flattest"). We identify an inverse relationship between CMR and sequence length in an in silico system with a two-peak fitness landscape; CMR decreases to no more than five orders of magnitude above estimates of eukaryotic per base mutation rate. We confirm the CMR reduces exponentially at low population sizes, irrespective of peak radius and distance, and increases with the number of genetic crossovers. We also identify an inverse relationship between CMR and the number of genes, confirming that, for a similar number of genes to that for the plant Arabidopsis thaliana (25,000), the CMR is close to its known wild-type mutation rate; mutation rates for additional organisms were also found to be within one order of magnitude of the CMR. This is the first time such a simulation model has been assigned input and produced output within range for a given biological organism. The decrease in CMR with population size previously observed is maintained; there is potential for the model to influence understanding of populations undergoing bottleneck, stress, and conservation strategy for populations near extinction.

ELFN1-AS1: a novel primate gene with possible microRNA function expressed predominantly in human tumors

Human gene LOC100505644 uncharacterized LOC100505644 [Homo sapiens] (Entrez Gene ID 100505644) is abundantly expressed in tumors but weakly expressed in few normal tissues. Till now the function of this gene remains unknown. Here we identified the chromosomal borders of the transcribed region and the major splice form of the LOC100505644-specific transcript. We characterised the major regulatory motifs of the gene and its splice sites. Analysis of the secondary structure of the major transcript variant revealed a hairpin-like structure characteristic for precursor microRNAs. Comparative genomic analysis of the locus showed that it originated in primates de novo. Taken together, our data indicate that human gene LOC100505644 encodes some non-protein coding RNA, likely a microRNA. It was assigned a gene symbol ELFN1-AS1 (ELFN1 antisense RNA 1 (non-protein coding)). This gene combines features of evolutionary novelty and predominant expression in tumors.

Comparative genomic analysis of dinucleotide repeats in Tritryps

The protozoans Trypanosoma cruzi, Trypanosoma brucei and Leishmania major (Tritryps), are evolutionarily ancient eukaryotes which cause worldwide human parasitosis. They present unique biological features. Indeed, canonical DNA/RNA cis-acting elements remain mostly elusive. Repetitive sequences, originally considered as selfish DNA, have been lately recognized as potentially important functional sequence elements in cell biology. In particular, the dinucleotide patterns have been related to genome compartmentalization, gene evolution and gene expression regulation. Thus, we perform a comparative analysis of the occurrence, length and location of dinucleotide repeats (DRs) in the Tritryp genomes and their putative associations with known biological processes. We observe that most types of DRs are more abundant than would be expected by chance. Complementary DRs usually display asymmetrical strand distribution, favoring TT and GT repeats in the coding strands. In addition, we find that GT repeats are among the longest DRs in the three genomes. We also show that specific DRs are non-uniformly distributed along the polycistronic unit, decreasing toward its boundaries. Distinctive non-uniform density patterns were also found in the intergenic regions, with predominance at the vicinity of the ORFs. These findings further support that DRs may control genome structure and gene expression.

Allelic polymorphism within the TAS1R3 promoter is associated with human taste sensitivity to sucrose

Human sweet taste perception is mediated by the heterodimeric G protein-coupled receptor encoded by the TAS1R2 and TAS1R3 genes. Variation in these genes has been characterized, but the functional consequences of such variation for sweet perception are unknown. We found that two C/T single-nucleotide polymorphisms (SNPs) located at positions -1572 (rs307355) and -1266 (rs35744813) upstream of the TAS1R3 coding sequence strongly correlate with human taste sensitivity to sucrose and explain 16% of population variability in perception. By using a luciferase reporter assay, we demonstrated that the T allele of each SNP results in reduced promoter activity in comparison to the C alleles, consistent with the phenotype observed in humans carrying T alleles. We also found that the distal region of the TAS1R3 promoter harbors a composite cis-acting element that has a strong silencing effect on promoter activity. We conclude that the rs307355 and rs35744813 SNPs affect gene transcription by altering the function of this regulatory element. A worldwide population survey reveals that the T alleles of rs307355 and rs35744813 occur at lowest frequencies in European populations. We propose that inherited differences in TAS1R3 transcription account for a substantial fraction of worldwide differences in human sweet taste perception.

Exploring the relationship between polymorphic (TG/CA)n repeats in intron 1 regions and gene expression

The putative role of (TG/CA)n repeats in the regulation of transcription has recently been reported for several cancer- and disease-related genes, including the genes encoding the epidermal growth factor receptor (EGFR), hydroxysteroid (11-beta) dehydrogenase 2 (HSD11B2) and interferon-gamma (IFNG). These studies indicated a correlation between gene expression levels and the presence or length of (TG/CA)n repeats in their intron 1 regions. A genome-wide search for genes with similar features may provide evidence of whether these dinucleotide repeats represent a class of universal regulators of gene expression, which has recently begun to be investigated as a quantitative complex phenotype. Using a public database of simple repeats, we identified 330 genes containing potentially polymorphic long (TG/CA)n repeats (n >or= 12) in their intron 1 regions. One known physiological pathway, the calcium signalling pathway, was found to be enriched among the genes containing long repeats. In addition, certain biological processes, such as cation transport, signal transduction and ion transport, were found to be enriched in these genes. Genotyping of the long repeats showed that the majority of these dinucleotide repeats were polymorphic in the HapMap CEU (Caucasians from Utah, USA) samples of northern and western European ancestry. Evidence for a significant association between these repeats and gene expression was not observed in the genes selected based on their expression profiles in the HapMap CEU samples. Our current findings, therefore, do not support a role for these repeats as a class of universal gene expression regulators. A more comprehensive evaluation of the relationship between these repeats and gene expression, potentially in other tissues, may be necessary to illustrate their roles in gene regulation in the future.

Association of Microsatellite Polymorphisms of the Human 14q13.2 Region with Type 2 Diabetes Mellitus in Latvian and Finnish Populations

A polymorphic microsatellite in intron 6 of the human proteasome core particle PSMA6 gene (HSMS006), and four other microsatellites localized upstream on human chromosome 14q13.2 (HSMS801, HSMS702, HSMS701, HSMS602), were genotyped in 104 type 2 diabetic patients and 129 age-matched control subjects from Latvia and replicated in 91 type 2 diabetic patients and 88 age-matched healthy control subjects from the Botnia Study in Finland. In type 2 diabetic patients from both populations the HSMS006 (TG)22 allele was two times more frequent compared to the control group. In the Latvian population the (CAA)8 allele of the HSMS602 marker was less frequent in the diabetic group, as was the (AC)24 allele of microsatellite HSMS801. Allele frequencies of the HSMS701 and 702 repeats were similar in healthy controls and type 2 diabetic patients. In conclusion, our data suggest that variants in the PSMA6 gene on chromosome 14q13.2 are associated with type 2 diabetes.

Mining microsatellites in eukaryotic genomes

During recent decades, microsatellites have become the most popular source of genetic markers. More recently, the availability of enormous sequence data for a large number of eukaryotic genomes has accelerated research aimed at understanding the origin and functions of microsatellites and searching for new applications. This review presents recent developments of in silico mining of microsatellites to reveal various facets of the distribution and dynamics of microsatellites in eukaryotic genomes. Two aspects of microsatellite search strategies--using a suitable search tool and accessing a relevant microsatellite database--have been explored. Judicious microsatellite mining not only helps in addressing biological questions but also facilitates better exploitation of microsatellites for diverse applications.

Abundance of dinucleotide repeats and gene expression are inversely correlated: a role for gene function in addition to intron length

High and broad transcription of eukaryotic genes is facilitated by cost minimization, clustered localization in the genome, elevated G+C content, and low nucleosome formation potential. In this scenario, illumination of correlation between abundance of (TG/CA)(n>or=12) repeats, which are negative cis modulators of transcription, and transcriptional levels and other commonly occurring dinucleotide repeats, is required. Three independent microarray datasets were used to examine the correlation of (TG/CA)(n>or=12) and other dinucleotide repeats with gene expression. Compared with the expected equi-distribution pattern under neutral model, highly transcribed genes were poor in repeats, and conversely, weakly transcribed genes were rich in repeats. Furthermore, the inverse correlation between repeat abundance and transcriptional levels appears to be a global phenomenon encompassing all genes regardless of their breadth of transcription. This selective pattern of exclusion of (TG/CA)(n>or=12) and (AT)(n>or=12) repeats in highly transcribed genes is an additional factor along with cost minimization and elevated GC, and therefore, multiple factors govern high transcription of genes. We observed that even after controlling for the effects of GC and average intron lengths, the effect of repeats albeit somewhat weaker was persistent and definite. In the ribosomal protein coding genes, sequence analysis of orthologs suggests that negative selection for repeats perhaps occurred early in evolution. These observations suggest that negative selection of (TG/CA)(n>or=12) microsatellites in the evolution of the highly expressed genes was also controlled by gene function in addition to intron length.

Expoldb: expression linked polymorphism database with inbuilt tools for analysis of expression and simple repeats

BACKGROUND

Quantitative variation in gene expression has been proposed to underlie phenotypic variation among human individuals. A facilitating step towards understanding the basis for gene expression variability is associating genome wide transcription patterns with potential cis modifiers of gene expression.

DESCRIPTION

EXPOLDB, a novel Database, is a new effort addressing this need by providing information on gene expression levels variability across individuals, as well as the presence and features of potentially polymorphic (TG/CA)n repeats. EXPOLDB thus enables associating transcription levels with the presence and length of (TG/CA)n repeats. One of the unique features of this database is the display of expression data for 5 pairs of monozygotic twins, which allows identification of genes whose variability in expression, are influenced by non-genetic factors including environment. In addition to queries by gene name, EXPOLDB allows for queries by a pathway name. Users can also upload their list of HGNC (HUGO (The Human Genome Organisation) Gene Nomenclature Committee) symbols for interrogating expression patterns. The online application 'SimRep' can be used to find simple repeats in a given nucleotide sequence. To help illustrate primary applications, case examples of Housekeeping genes and the RUNX gene family, as well as one example of glycolytic pathway genes are provided.

CONCLUSION

The uniqueness of EXPOLDB is in facilitating the association of genome wide transcription variations with the presence and type of polymorphic repeats while offering the feature for identifying genes whose expression variability are influenced by non genetic factors including environment. In addition, the database allows comprehensive querying including functional information on biochemical pathways of the human genes. EXPOLDB can be accessed at http://expoldb.igib.res.in/expol.

Specific Interactions of Divalent Metal Ions with a DNA Duplex Containing the d(CA)n/(GT)nTandem Repeat

Divalent metal ions are essential for maintaining functional states of the DNA molecule. Their participation in DNA structure is modulated by the base sequence and varies depending on the nature of the ion. The present investigation addresses the interaction of Ca2+ ions with a tandem repeat of two CA dinucleotides, (CA)2/(TG)2. The binding of Ca2+ to the repeat is monitored by nuclear magnetic resonance (NMR) spectroscopy using chemical shift mapping. Parallel experiments monitor binding of Mg2+ ions to the repeat as well as binding of each ion to a DNA duplex in which the (CA)2/(TG)2 repeat is eliminated. The results reveal that the direction and the magnitude of chemical shift changes induced by Ca2+ ions in the NMR spectra of the repeat are different from those induced by Mg2+ ions. The differences between the two cations are significantly diminished by the elimination of the (CA)2/(TG)2 repeat. These findings suggest a specific interaction of Ca2+ ions with the (CA)2/(TG)2 motif. The specificity of the interaction resides in the two A-T base pairs of the repeat, and it involves the major groove of the first A-T base pair and both grooves of the second A-T base pair.

Microsatellite profile in hormonal receptor genes associated with breast cancer

Given that breast cancer is depending on multiple hormonal influences, the nuclear receptors, estrogen receptor alpha, estrogen receptor beta and androgen receptor, are candidates for cancer susceptibility markers. We conducted an association study in a case-control population (139 cases and 145 controls) by genotyping three potentially functional microsatellites (TA)n, (CA)n and (CAG)n in the ERa, ERb and AR genes respectively. For (CAG)n polymorphism, a significant difference was observed using a cut-off 15 repeats CAG between genotypes short-short/short-long/long-long in cases and control subjects (p = 0.009) and also between the distribution of short/long allele in the two groups of individuals (p = 0.001). Genotypes comprising one or two short (CAG)n sequences had higher risk of breast cancer compared to genotypes with two long allele (odds ratio = 1,93; confidence interval = 1.05-3.55; p = 0.03). No significant difference was observed in allele frequency or in short/long allele percentage for (CA)n or (TA)n polymorphism (cut-off 22 CA and 19 TA repeats), neither in genotype frequencies (short-short, short-long or long-long). When the three microsatellite genotype were taken in analysis, the profile short CA-long TA-short CAG could clearly discriminate between cases and controls (p = 0.006). Also, this combined genotype profile has greater predictive values for breast cancer than (CAG)n genotype alone (predictive positive value 57,1% versus 53,7% and predictive negative value 53% versus 23% respectively). Our results sustain a polygenic model of breast cancer with gene-gene interactions; combined effects of three low-risk polymorphisms conferred significant genetic predisposition. Genotyping hormonal receptor genes ERa, ERb and AR could be a useful genetic marker for defining disease risk.