Molecular phylogeny reveals distinct evolutionary lineages of the banded krait, Bungarus fasciatus (Squamata, Elapidae) in Asia.

The banded krait, Bungarus fasciatus is a widespread elapid snake, likely to comprise several distinct species in different geographic regions of Asia. Therefore, based on molecular phylogenetics and comparative morphology data, we present an overview of the systematic composition of the species to delimit potential biogeographic boundaries. Our phylogenetic analyses, based on four mitochondrial genes, reveal the existence of at least three evolutionary lineages within B. fasciatus, corresponding to Indo-Myanmar, Sundaic and eastern Asian lineages. We are convinced that there are at least three taxonomic entities within the nomen B. fasciatus and restrict the distribution of B. fasciatus sensu stricto to the Indo-Myanmar region. We also provide additional natural history data of the taxon from eastern India. Finally, we advocate further studies to establish the degree of reproductive isolation among these diverging evolutionary lineages and to reassess the systematic status of this species complex especially the Sundaic and eastern Asian lineages.

Cytotoxin antibody-based colourimetric sensor for field-level differential detection of elapid among big four snake venom.

Development of a rapid, on-site detection tool for snakebite is highly sought after, owing to its clinically and forensically relevant medicolegal significance. Polyvalent antivenom therapy in the management of such envenomation cases is finite due to its poor venom neutralization capabilities as well as diagnostic ramifications manifested as untoward immunological reactions. For precise molecular diagnosis of elapid venoms of the big four snakes, we have developed a lateral flow kit using a monoclonal antibody (AB1; IgG1 - κ chain; Kd: 31 nM) generated against recombinant cytotoxin-7 (rCTX-7; 7.7 kDa) protein of the elapid venom. The monoclonal antibody specifically detected the venoms of Naja naja (p < 0.0001) and Bungarus caeruleus (p<0.0001), without showing any immunoreactivity against the viperidae snakes in big four venomous snakes. The kit developed attained the limit of quantitation of 170 pg/µL and 2.1 ng/µL in spiked buffer samples and 28.7 ng/µL and 110 ng/µL in spiked serum samples for detection of N. naja and B. caeruleus venoms, respectively. This kit holds enormous potential in identification of elapid venom of the big four snakes for effective prognosis of an envenomation; as per the existing medical guidelines.

A Wolf in Another Wolf's Clothing: Post-Genomic Regulation Dictates Venom Profiles of Medically-Important Cryptic Kraits in India.

The Common Krait (Bungarus caeruleus) shares a distribution range with many other 'phenotypically-similar' kraits across the Indian subcontinent. Despite several reports of fatal envenomings by other Bungarus species, commercial Indian antivenoms are only manufactured against B. caeruleus. It is, therefore, imperative to understand the distribution of genetically distinct lineages of kraits, the compositional differences in their venoms, and the consequent impact of venom variation on the (pre)clinical effectiveness of antivenom therapy. To address this knowledge gap, we conducted phylogenetic and comparative venomics investigations of kraits in Southern and Western India. Phylogenetic reconstructions using mitochondrial markers revealed a new species of krait, Romulus' krait (Bungarus romulusi sp. nov.), in Southern India. Additionally, we found that kraits with 17 mid-body dorsal scale rows in Western India do not represent a subspecies of the Sind Krait (B. sindanus walli) as previously believed, but are genetically very similar to B. sindanus in Pakistan. Furthermore, venom proteomics and comparative transcriptomics revealed completely contrasting venom profiles. While the venom gland transcriptomes of all three species were highly similar, venom proteomes and toxicity profiles differed significantly, suggesting the prominent role of post-genomic regulatory mechanisms in shaping the venoms of these cryptic kraits. In vitro venom recognition and in vivo neutralisation experiments revealed a strong negative impact of venom variability on the preclinical performance of commercial antivenoms. While the venom of B. caeruleus was neutralised as per the manufacturer's claim, performance against the venoms of B. sindanus and B. romulusi was poor, highlighting the need for regionally-effective antivenoms in India.

Mass spectrometric analysis to unravel the venom proteome composition of Indian snakes: opening new avenues in clinical research.

INTRODUCTION: The 'Big Four' venomous snakes - Daboia russelii, Naja naja, Bungarus caeruleus, and Echis carinatus - are primarily responsible for the majority of snake envenomation in India. Several other lesser-known venomous snake species also inflict severe envenomation in the country. AREAS COVERED: A comprehensive analysis of the venom proteome composition of the 'Big Four' and other medically important venomous snakes of India and the effect of regional variation in venom composition on immunorecognition and/or neutralization by commercial antivenom was undertaken by searching the literature (from 1985 to date) available in large public databases. Further, mass spectrometric identification of poorly immunogenic toxins of snake venom (against which commercial polyvalent antivenom contains a significantly lower proportion of antibodies) and its impact on antivenom therapy against snakebite are discussed. The application of mass spectrometry to identify protein (toxin) complexes as well as drug prototypes from Indian snake venoms and the clinical importance of such studies are also highlighted. EXPERT OPINION: Further detailed clinical and proteomic research is warranted to better understand the effects of regional snake venom composition on the clinical manifestation of envenomation and antivenom therapy and to improve the production of antibodies against poorly immunogenic venom components.

Cross-species transmission of the newly identified coronavirus 2019-nCoV.

The current outbreak of viral pneumonia in the city of Wuhan, China, was caused by a novel coronavirus designated 2019-nCoV by the World Health Organization, as determined by sequencing the viral RNA genome. Many initial patients were exposed to wildlife animals at the Huanan seafood wholesale market, where poultry, snake, bats, and other farm animals were also sold. To investigate possible virus reservoir, we have carried out comprehensive sequence analysis and comparison in conjunction with relative synonymous codon usage (RSCU) bias among different animal species based on the 2019-nCoV sequence. Results obtained from our analyses suggest that the 2019-nCoV may appear to be a recombinant virus between the bat coronavirus and an origin-unknown coronavirus. The recombination may occurred within the viral spike glycoprotein, which recognizes a cell surface receptor. Additionally, our findings suggest that 2019-nCoV has most similar genetic information with bat coronovirus and most similar codon usage bias with snake. Taken together, our results suggest that homologous recombination may occur and contribute to the 2019-nCoV cross-species transmission.

Expression and characterization of Flavikunin: A Kunitz-type serine protease inhibitor identified in the venom gland cDNA library of Bungarus flaviceps.

Trancriptomic analysis of the venom gland cDNA library of Bungarus flaviceps revealed Kunitz-type serine protease inhibitor as one of the major venom protein families with three groups A, B, C. One of the group B isoforms named Flavikunin, which lacked an extra cysteine residue involved in disulfide bond formation in ß-bungarotoxin, was synthesized, cloned, and overexpressed in Escherichia coli. To decipher the structure-function relationship, the P1 residue of Flavikunin, histidine, was mutated to alanine and arginine. Purified wild-type and mutant Flavikunins were screened against serine proteases-thrombin, factor Xa, trypsin, chymotrypsin, plasmin, and elastase. The wild-type and mutant Flavikunin (H∆R) inhibited plasmin with an IC 50 of 0.48 and 0.35 µM, respectively. The in-silico study showed that P1 residue of wild-type and mutant (H∆R) Flavikunin interacted with S1' and S1 site of plasmin, respectively. Thus, histidine at the P1 position was found to be involved in plasmin inhibition with mild anticoagulant activity.

Northward dispersal of sea kraits (Laticauda semifasciata) beyond their typical range.

Marine reptiles are declining globally, and recent climate change may be a contributing factor. The study of sea snakes collected beyond their typical distribution range provides valuable insight on how climate change affects marine reptile populations. Recently, we collected 12 Laticauda semifasciata (11 females, 1 male) from the waters around southern South Korea-an area located outside its typical distribution range (Japan, China including Taiwan, Philippines and Indonesia). We investigated the genetic origin of Korean specimens by analyzing mitochondrial cytochrome b gene (Cytb) sequences. Six individuals shared haplotypes with a group found in Taiwan-southern Ryukyu Islands, while the remaining six individuals shared haplotypes with a group encompassing the entire Ryukyu Archipelago. These results suggest L. semifasciata moved into Korean waters from the Taiwan-Ryukyu region via the Taiwan Warm Current and/or the Kuroshio Current, with extended survival facilitated by ocean warming. We highlight several contributing factors that increase the chances that L. semifasciata establishes new northern populations beyond the original distribution range.

Codon usage pattern and prediction of gene expression level in Bungarus species.

Codon bias study in an organism gains significance in understanding the molecular mechanism as well as the functional conservation of gene expression during the course of evolution. The prime focus in this study is to compare the codon usage patterns among the four species belonging to the genus Bungarus (B. multicinctus, B. fasciatus, B. candidus and B. flaviceps) using several codon bias parameters. Our results suggested that relatively low codon bias exists in the coding sequences of the selected species. The compositional constraints together with gene expression level might influence the patterns of codon usage among the genes of Bungarus species. Both natural selection and mutation pressure affect the codon usage pattern in Bungarus species as evident from correspondence analysis. Neutrality plot indicates that natural selection played a major role while mutation pressure played a minor role in codon usage pattern of the genes in Bungarus species.

The Primary Structure of ß(I)-Chain of Hemoglobin from Snake Sindhi Krait (Bungarus sindanus sindanus).

The amino acid sequence of ß(I)-globin chain from Sindhi Krait (Bungarus sindanus sindanus) was determined to study the molecular evolution among snakes. The hemoglobin was isolated from the red blood cells and was analyzed by ion-exchange chromatography (IEX). The crude globin was subjected to reversed phased-high performance liquid chromatography (RP-HPLC) using C4 column. The N-terminal sequences of intact globin chains and tryptic peptides were determined by Edman degradation in a pulsed liquid gas phase sequencer using an online Phenylthiohydantoin analyzer. Sindhi Krait is expected to express three hemoglobin components that are composed of ß(II), ß(I), α(D) and α(A)-globin chains, as apparent by IEX, RP-HPLC and N-terminal sequence analyses. Sequence alignment and phylogenetic analyses of ß(I) globin chain from Sindhi Krait showed closest relationship with ß(I) globin chain from Rattlesnake, Water snake and Indigo snake. Interestingly, comparison of primary sequence of ß(I) globin chain of Sindhi Krait with human ß chain revealed 63 % similarity along with the retention of all heme contact points. Variations among the two sequences were prominent at αß contact points and in regions directly not important for function.

Venom gland transcriptomes of two elapid snakes (Bungarus multicinctus and Naja atra) and evolution of toxin genes.

BACKGROUND: Kraits (genus Bungarus) and cobras (genus Naja) are two representative toxic genera of elapids in the old world. Although they are closely related genera and both of their venoms are very toxic, the compositions of their venoms are very different. To unveil their detailed venoms and their evolutionary patterns, we constructed venom gland cDNA libraries and genomic bacterial artificial chromosome (BAC) libraries for Bungarus multicinctus and Naja atra, respectively. We sequenced about 1500 cDNA clones for each of the venom cDNA libraries and screened BAC libraries of the two snakes by blot analysis using four kinds of toxin probes; i.e., three-finger toxin (3FTx), phospholipase A2 (PLA2), kunitz-type protease inhibitor (Kunitz), and natriuretic peptide (NP). RESULTS: In total, 1092 valid expressed sequences tags (ESTs) for B. multicinctus and 1166 ESTs for N. atra were generated. About 70% of these ESTs can be annotated as snake toxin transcripts. 3FTx (64.5%) and ß bungarotoxin (25.1%) comprise the main toxin classes in B. multicinctus, while 3FTx (95.8%) is the dominant toxin in N. atra. We also observed several less abundant venom families in B. multicinctus and N. atra, such as PLA2, C-type lectins, and Kunitz. Peculiarly a cluster of NP precursors with tandem NPs was detected in B. multicinctus. A total of 71 positive toxin BAC clones in B. multicinctus and N. atra were identified using four kinds of toxin probes (3FTx, PLA2, Kunitz, and NP), among which 39 3FTx-positive BACs were sequenced to reveal gene structures of 3FTx toxin genes. CONCLUSIONS: Based on the toxin ESTs and 3FTx gene sequences, the major components of B. multicinctus venom transcriptome are neurotoxins, including long chain alpha neurotoxins (α-ntx) and the recently originated ß bungarotoxin, whereas the N. atra venom transcriptome mainly contains 3FTxs with cytotoxicity and neurotoxicity (short chain α-ntx). The data also revealed that tandem duplications contributed the most to the expansion of toxin multigene families. Analysis of nonsynonymous to synonymous nucleotide substitution rate ratios (dN/dS) indicates that not only multigene toxin families but also other less abundant toxins might have been under rapid diversifying evolution.

Transcriptomic analysis of the venom gland of the red-headed krait (Bungarus flaviceps) using expressed sequence tags.

BACKGROUND: The Red-headed krait (Bungarus flaviceps, Squamata: Serpentes: Elapidae) is a medically important venomous snake that inhabits South-East Asia. Although the venoms of most species of the snake genus Bungarus have been well characterized, a detailed compositional analysis of B. flaviceps is currently lacking. RESULTS: Here, we have sequenced 845 expressed sequence tags (ESTs) from the venom gland of a B. flaviceps. Of the transcripts, 74.8% were putative toxins; 20.6% were cellular; and 4.6% were unknown. The main venom protein families identified were three-finger toxins (3FTxs), Kunitz-type serine protease inhibitors (including chain B of beta-bungarotoxin), phospholipase A2 (including chain A of beta-bungarotoxin), natriuretic peptide (NP), CRISPs, and C-type lectin. CONCLUSION: The 3FTxs were found to be the major component of the venom (39%). We found eight groups of unique 3FTxs and most of them were different from the well-characterized 3FTxs. We found three groups of Kunitz-type serine protease inhibitors (SPIs); one group was comparable to the classical SPIs and the other two groups to chain B of beta-bungarotoxins (with or without the extra cysteine) based on sequence identity. The latter group may be functional equivalents of dendrotoxins in Bungarus venoms. The natriuretic peptide (NP) found is the first NP for any Asian elapid, and distantly related to Australian elapid NPs. Our study identifies several unique toxins in B. flaviceps venom, which may help in understanding the evolution of venom toxins and the pathophysiological symptoms induced after envenomation.

[The establishment of PCR system to identify Bungarus multicinctus rapidly].

The purpose of the present study is to establish a rapid and effective PCR method for the identification of B. multicinctus. Based on sequence alignment of B. multicinctus and its adulterants, we found that Cyt b gene is a good molecular genetic marker for the authentication of B. multicinctus. On the basis of the sequence data, a pair of highly specialized primers was designed. The templates were extracted by the DNA purification system. Key factors such as annealing temperature, concentration of Taq enzyme and cycle numbers were analyzed and optimized. The modified PCR program consisted of an initial denaturation step at 95 degrees C for 5 min, followed by 30 cycles of 95 degrees C for 30 s and 55 degrees C for 45 s and a final extension at 72 degrees C for 5 min. Thirteen samples of B. multicinctus were identified accurately from their 20 adulterants in 4 hours. The results indicated it is a highly accurate, rapid and applicable method for the authentication of B. multicinctus.

Genetic organization of Bungarus multicinctus protease inhibitor-like proteins.

The structural organization of the genes encoding Bungarus multicinctus protease inhibitor-like proteins (PILPs), PILP-1, PILP-2 and PILP-3, are reported in this study. Unlike PILP-2 and PILP-3, recombinant PILP-1 exhibited inhibitory activity on trypsin. PILP genes and B chain genes shared identical organization with three exons interrupted by two introns in similar positions. On the contrary, intron 1 of these genes had a similar size, a notable variation with the size of intron 2 was observed. It was found that two regions at the second intron of B1 chain and B2 chain genes were absent in that of PILP genes. Noticeably, intronic insertion in the second intron of B chain genes appeared in the promoter region of PILP-1 gene, but not in that of PILP-2 and PILP-3 genes. Comparative analyses of PILP genes and B chain genes showed that the protein-coding regions of the exons are more diverse than introns, except for in the signal peptide domain. These results suggest that PILP genes and B chain genes originate from a common ancestor, and that accelerated evolution may diversify PILP and B chain genes as that proposed for snake venom phospholipase A(2), neurotoxin and cardiotoxin genes.

Expression pattern of three-finger toxin and phospholipase A2 genes in the venom glands of two sea snakes, Lapemis curtus and Acalyptophis peronii: comparison of evolution of these toxins in land snakes, sea kraits and sea snakes.

BACKGROUND: Snake venom composition varies widely both among closely related species and within the same species, based on ecological variables. In terrestrial snakes, such variation has been proposed to be due to snakes' diet. Land snakes target various prey species including insects (arthropods), lizards (reptiles), frogs and toads (amphibians), birds (aves), and rodents (mammals), whereas sea snakes target a single vertebrate class (fishes) and often specialize on specific types of fish. It is therefore interesting to examine the evolution of toxins in sea snake venoms compared to that of land snakes. RESULTS: Here we describe the expression of toxin genes in the venom glands of two sea snakes, Lapemis curtus (Spine-bellied Sea Snake) and Acalyptophis peronii (Horned Sea Snake), two members of a large adaptive radiation which occupy very different ecological niches. We constructed cDNA libraries from their venom glands and sequenced 214 and 192 clones, respectively. Our data show that despite their explosive evolutionary radiation, there is very little variability in the three-finger toxin (3FTx) as well as the phospholipase A2 (PLA2) enzymes, the two main constituents of Lapemis curtus and Acalyptophis peronii venom. To understand the evolutionary trends among land snakes, sea snakes and sea kraits, pairwise genetic distances (intraspecific and interspecific) of 3FTx and PLA2 sequences were calculated. Results show that these proteins appear to be highly conserved in sea snakes in contrast to land snakes or sea kraits, despite their extremely divergent and adaptive ecological radiation. CONCLUSION: Based on these results, we suggest that streamlining in habitat and diet in sea snakes has possibly kept their toxin genes conserved, suggesting the idea that prey composition and diet breadth may contribute to the diversity and evolution of venom components.

Cloning, expression and characterization of two C-type lectins from the venom gland of Bungarus multicinctus.

C-type lectins found in many animals are non-enzymatic proteins and able to bind with mono- and oligosaccharides in a Ca(2+)-dependent fashion. Here, we report the cloning of two C-type lectins named BML-1 and BML-2 from the venom gland of Bungarus multicinctus, and expression of their mature peptides with 135 and 137 amino acids as inclusion bodies. Recombinant BML-1 and BML-2 proteins with 135 amino acids formed monomers, and those with 137 amino acids formed homodimers and monomers and both of them displayed certain hemagglutinating activity to rabbit erythrocytes. The results of Western blotting and immuno-affinity chromatography demonstrated that C-type lectins in B. multicinctus formed dimers in physiological conditions, and their molecular weight is lower than previous predictions. This is the first report of the cloning of the BML-2 gene from the venom gland of B. multicinctus, as well as an investigation of its confirmation and biological functions.

[High specific PCR identification of Bungarus multicinctus and its adulterants].

OBJECTIVE: To develop a convenient and effective method for the identification of Bungarus multicinctus. METHOD: Based on the sequence of Cyt b gene fragment of B. multicinctus and its adulterants, a pair of highly specific primer (HJL- and HJH-) were designed for distinguishing B. ulticinctus from other species of snake. To establish specific PCR reaction condition, the primers were employed to amplify the DNA templates extracted from B. multicinctus and 6 other species of snake, under different annealing temperature. Using this method, B. multicinctus was identified from 18 samples bought from many drugstores. RESULT: A 230 bp DNA fragment was amplified from B. multicinctus in PCR with annealed temperature at 67 degrees C, whereas no DNA fragment was amplified from other snake samples under the same reaction condition, B. multicinctus could be clearly distinguished from others by PCR reaction with the highly specific primers. In the present study, 18 sample, bought from different drugstores, were also identified by the highly specific PCR with the primers. The results indicated that 14 samples were B. multicinctus and the other 4 were adulterant, which was consistent with the conclusion of authentication based on morphological. CONCLUSION: The primers designed in the present study were highly specific for B. multicinctus.

[Molecular cloning and expression of an isotoxin gene, alpha-bungarotoxin, from Bungarus multicinctus].

Abstract: Snake venom contains a number of small proteins,enzymes and other components,which displays a broad spectrum of biological activities. With the ability of specifically binding on acetylcholine acceptor, alpha-bungarotoxins are not only useful molecular probes in investigating the mechanism of neural signal transmission, but also potential pharmic preparations for neural disease treatment. In current research,cDNAs of Bungarus multicinutus venom gland were synthesized using SMART cDNA amplification kit and then, alpha-bungarotoxin genes were cloned and sequenced. Total of 20 clones were sequenced representing 14 isotoxin mRNAs of alpha-bungarotoxins. Among those clones, a novel isotoxin gene was subcloned into two expression plasmids, alpha-BgTX/pQE30a and alpha-BgTX/pGEX-4T-1, and transformed into E. coli. After inducing with IPTG, fused protein of GST-alpha-BgTX was successfully expressed at level of 30% gross proteins of bacteria. More than 25% of fused protein was in the soluble fraction and the rest in inclusion body.

Novel neurotoxins from Taiwan banded krait (Bungarus multicinctus) venom: purification, characterization and gene organization.

Two novel neurotoxins BM10-1 and BM10-2 were isolated from Bungarus multicinctus (Taiwan banded krait) venom using the combinations of chromatography on a SP-Sephadex C-25 column and a reverse phase HPLC column. BM10-1 contained 66 amino acid residues including 10 Cys residues, while BM10-2 consisted of 65 amino acid residues with 8 Cys residues. The secondary structure of both BM10-1 and BM10-2 was dominated with beta-sheet, but their gross conformation differed as evidenced by CD spectra and acrylamide quenching studies. BM10-1 inhibited carbachol-induced muscle contraction in a reversible manner and the dose for achieving 50% inhibition was approximately fourfold that of alpha-bungarotoxin. BM10-2 exhibited an irreversible but weak inhibition on carbachol-induced muscle contraction. Sequence alignment of neurotoxins with BM10-1 and BM10-2 suggested that the manner in the manifestation of their activity could be partly elucidated by the residues at toxin second loop. The genomic DNAs encoding BM10-1 and BM10-1-like protein (BM10-1L) were amplified by PCR. The two genes shared virtually identical structural organization and high degree of sequence identity with B. multicinctus neurotoxin genes. Compared to intron sequences of these genes, the protein-coding regions were highly variable. The difference between BM10-1 gene and BM10-1L gene notably arose from the third exon. These results suggest the evolution of B. multicinctus neurotoxins via the path of gene duplication.

Identification of alpha-bungarotoxin (A31) as the major postsynaptic neurotoxin, and complete nucleotide identity of a genomic DNA of Bungarus candidus from Java with exons of the Bungarus multicinctus alpha-bungarotoxin (A31) gene.

The Malayan krait (Bungarus candidus) is one of the most medically significant snake species in Southeast Asia. No specific antivenom exists to treat envenoming by this species. Death within 30 min after its bite has been reported from Java, suggesting the presence of highly lethal postsynaptic neurotoxins in the venom of these snakes. We purified and identified the major postsynaptic toxin in the venom of B. candidus from Java. The toxin was indistinguishable from alpha-bungarotoxin (A31), a toxin originally isolated from Bungarus multicinctus, in its mass (7983.75 Da), LD50 (0.23 microg/g in mice i.p.), affinity to nicotinic acetylcholine receptors, and by its 40 N-terminal amino acid residues as determined by Edman degradation. Identity with alpha-bungarotoxin was confirmed by cloning and sequencing a genomic DNA from B. candidus which encodes the 74 amino acid sequence of alpha-bungarotoxin (A31) and part of its signal peptide, revealing complete identity to the alpha-bungarotoxin (A31) gene in exon and 98.9% identity in intron sequences. The entire mitochondrial cytochrome b gene of the krait species B. candidus from Java and B. multicinctus from Taiwan was sequenced for comparison, suggesting that these snakes are phylogenetically closely related. alpha-Bungarotoxin appears to be widely present and conserved in Southeast and East Asian black-and-white kraits across populations and taxa.

Genome-wide analysis of Bkm sequences (GATA repeats): predominant association with sex chromosomes and potential role in higher order chromatin organization and function.

MOTIVATION: Bkm (Banded krait minor) satellite DNA sequences (GATA repeats) have been shown to be associated with the sex determining chromosomes of various eukaryotes and have been implicated in the evolution and differentiation of sex chromosomes in snakes. The objective of the study is to analyze the GATA repeats of human genome specifically, the Y-chromosome, and other model organisms to understand the possible function and potential role in higher order chromatin organization. RESULTS: Our extensive analysis of GATA repeats in the prokaryotic and eukaryotic genomes, which have been completely sequenced so far, has revealed that GATA repeats are absent in prokaryotes and have been gradually accumulated in higher organisms during the course of evolution. In human, the Y-chromosome has the highest GATA repeat density, which predominantly exists in the Yq centromeric region. Generally, occurrence of repeats in the genomes decreases steadily as the length of the repeat increases. In contrast, we report, that the occurrence of GATA repeats increases as the length of the repeat increases from six tandem repeats onwards and peaks at (GATA)(10-12). This has not been observed with any other simple repeat. Distribution of (GATA)(10-12) along the chromosome and their close proximity to Matrix Associated Regions (GATA-MAR) suggests that it may be demarking chromatin domains for a coordinated expression of genes residing in these domains.