Use of an ordered cosmid library to deduce the genomic organization of Mycobacterium leprae

In an attempt to unify the genetic and biological research on Mycobacterium leprae, the aetiological agent of leprosy, a cosmid library was constructed and then ordered by a combination of fingerprinting and hybridization techniques. The genome of M. leprae is represented by four contigs of overlapping clones which, together, account for nearly 2.B Mb of DNA. Several arguments suggest that the gaps between the contigs are small in size and that virtually complete coverage of the chromosome has been obtained. All of the cloned M. leprae genes have been positioned on the contig maps together with the 29 copies of the dispersed repetitive element, RLEP. These have been classified into four groups on the basis of differences in their organization. Several key housekeeping genes were identified and mapped by hybridization with heterologous probes, and the current genome map of this uncultivable pathogen comprises 72 loci.     

The Use of Bayesian Networks to Assess the Quality of Evidence from Research Synthesis: 1.

Background

The grades of recommendation, assessment, development and evaluation (GRADE) approach is widely implemented in systematic reviews, health technology assessment and guideline development organisations throughout the world. A key advantage to this approach is that it aids transparency regarding judgments on the quality of evidence. However, the intricacies of making judgments about research methodology and evidence make the GRADE system complex and challenging to apply without training.

Methods

We have developed a semi-automated quality assessment tool (SAQAT) l based on GRADE. This is informed by responses by reviewers to checklist questions regarding characteristics that may lead to unreliability. These responses are then entered into the Bayesian network to ascertain the probabilities of risk of bias, inconsistency, indirectness, imprecision and publication bias conditional on review characteristics. The model then combines these probabilities to provide a probability for each of the GRADE overall quality categories. We tested the model using a range of plausible scenarios that guideline developers or review authors could encounter.

Results

Overall, the model reproduced GRADE judgements for a range of scenarios. Potential advantages over standard assessment are use of explicit and consistent weightings for different review characteristics, forcing consideration of important but sometimes neglected characteristics and principled downgrading where small but important probabilities of downgrading are accrued across domains.

Conclusions

Bayesian networks have considerable potential for use as tools to assess the validity of research evidence. The key strength of such networks lies in the provision of a statistically coherent method for combining probabilities across a complex framework based on both belief and evidence. In addition to providing tools for less experienced users to implement reliability assessment, the potential for sensitivity analyses and automation may be beneficial for application and the methodological development of reliability tools.     

Use of short tandem repeat sequences to study Mycobacterium leprae in leprosy patients in Malawi and India

Background

Inadequate understanding of the transmission of Mycobacterium leprae makes it difficult to predict the impact of leprosy control interventions. Genotypic tests that allow tracking of individual bacterial strains would strengthen epidemiological studies and contribute to our understanding of the disease.

Methodology/Principal Findings

Genotyping assays based on variation in the copy number of short tandem repeat sequences were applied to biopsies collected in population-based epidemiological studies of leprosy in northern Malawi, and from members of multi-case households in Hyderabad, India. In the Malawi series, considerable genotypic variability was observed between patients, and also within patients, when isolates were collected at different times or from different tissues. Less within-patient variability was observed when isolates were collected from similar tissues at the same time. Less genotypic variability was noted amongst the closely related Indian patients than in the Malawi series.

Conclusions/Significance

Lineages of M. leprae undergo changes in their pattern of short tandem repeat sequences over time. Genetic divergence is particularly likely between bacilli inhabiting different (e.g., skin and nerve) tissues. Such variability makes short tandem repeat sequences unsuitable as a general tool for population-based strain typing of M. leprae, or for distinguishing relapse from reinfection. Careful use of these markers may provide insights into the development of disease within individuals and for tracking of short transmission chains.     

The Phylogeographic History of the New World Screwworm Fly,Inferred by Approximate Bayesian Computation Analysis

Insect pest phylogeography might be shaped both by biogeographic events and by human influence. Here, we conducted an approximate Bayesian computation (ABC) analysis to investigate the phylogeography of the New World screwworm fly, Cochliomyia hominivorax, with the aim of understanding its population history and its order and time of divergence. Our ABC analysis supports that populations spread from North to South in the Americas, in at least two different moments. The first split occurred between the North/Central American and South American populations in the end of the Last Glacial Maximum (15,300-19,000 YBP). The second split occurred between the North and South Amazonian populations in the transition between the Pleistocene and the Holocene eras (9,100-11,000 YBP). The species also experienced population expansion. Phylogenetic analysis likewise suggests this north to south colonization and Maxent models suggest an increase in the number of suitable areas in South America from the past to present. We found that the phylogeographic patterns observed in C. hominivorax cannot be explained only by climatic oscillations and can be connected to host population histories. Interestingly we found these patterns are very coincident with general patterns of ancient human movements in the Americas, suggesting that humans might have played a crucial role in shaping the distribution and population structure of this insect pest. This work presents the first hypothesis test regarding the processes that shaped the current phylogeographic structure of C. hominivorax and represents an alternate perspective on investigating the problem of insect pests.     

Detection of Mycobacterium leprae DNA by polymerase chain reaction in the blood and nasal secretion of Brazilian household contacts

DNA samples from blood and nasal swabs of 125 healthy household contacts was submitted to amplification by polymerase chain reaction (PCR) using a Mycobacterium leprae-specific sequence as a target for the detection of subclinical infection with M. leprae. All samples were submitted to hybridization analysis in order to exclude any false positive or negative results. Two positive samples were confirmed from blood out of 119 (1.7%) and two positive samples from nasal secretion out of 120 (1.7%). The analysis of the families with positive individuals showed that 2.5% (n = 3) of the contacts were relatives of multibacilary patients while 0.8% of the cases (n = 1) had a paucibacilary as an index case. All positive contacts were followed up and after one year none of them presented clinical signs of the disease. In spite of the PCR sensitivity to detect the presence of the M. leprae in a subclinical stage, this molecular approach did not seem to be a valuable tool to screen household contacts, since we determined a spurious association of the PCR positivity and further development of leprosy.     

Comparative Sequence Analysis of Mycobacterium leprae and the New Leprosy-Causing Mycobacterium lepromatosis

Mycobacterium lepromatosis is a newly discovered leprosy-causing organism. Preliminary phylogenetic analysis of its 16S rRNA gene and a few other gene segments revealed significant divergence from Mycobacterium leprae, a well-known cause of leprosy, that justifies the status of M. lepromatosis as a new species. In this study we analyzed the sequences of 20 genes and pseudogenes (22,814 nucleotides). Overall, the level of matching of these sequences with M. leprae sequences was 90.9%, which substantiated the species-level difference; the levels of matching for the 16S rRNA genes and 14 protein-encoding genes were 98.0% and 93.1%, respectively, but the level of matching for five pseudogenes was only 79.1%. Five conserved protein-encoding genes were selected to construct phylogenetic trees and to calculate the numbers of synonymous substitutions (dS values) and nonsynonymous substitutions (dN values) in the two species. Robust phylogenetic trees constructed using concatenated alignment of these genes placed M. lepromatosis and M. leprae in a tight cluster with long terminal branches, implying that the divergence occurred long ago. The dS and dN values were also much higher than those for other closest pairs of mycobacteria. The dS values were 14 to 28% of the dS values for M. leprae and Mycobacterium tuberculosis, a more divergent pair of species. These results thus indicate that M. lepromatosis and M. leprae diverged ∼10 million years ago. The M. lepromatosis pseudogenes analyzed that were also pseudogenes in M. leprae showed nearly neutral evolution, and their relative ages were similar to those of M. leprae pseudogenes, suggesting that they were pseudogenes before divergence. Taken together, the results described above indicate that M. lepromatosis and M. leprae diverged from a common ancestor after the massive gene inactivation event described previously for M. leprae.Leprosy, one of the oldest human diseases, remains a significant public health problem in many developing countries (8). Mycobacterium leprae was the only known cause of leprosy until recently, when a new mycobacterium, Mycobacterium lepromatosis, was found to be the cause of diffuse lepromatous leprosy (DLL), a unique form of leprosy endemic in Mexico and the Caribbean (17). The discovery of this new species may provide an explanation for the clinical and geographic variability of leprosy.The initial phylogenetic analysis of M. lepromatosis was carried out using the sequences of the 16S rRNA gene and segments of groEL, rpoB, and other genes (total, 4.99 kb) (17). This study revealed significant sequence differences between M. lepromatosis and all known Mycobacterium species and placed M. lepromatosis closest to M. leprae. However, the sequence variation justified assigning a new species for the new organism instead of classifying it as a variant of M. leprae. All M. leprae strains collected worldwide have been found to be clonal and to differ by only single-nucleotide polymorphism or variable numbers of tandem repeats (24). Also, the genomes of two M. leprae strains, strain TN from India (GenBank accession numbers {"type":"entrez-nucleotide-range","attrs":{"text":"AL583917 to AL583926","start_term":"AL583917","end_term":"AL583926","start_term_id":"13092412","end_term_id":"13093786"}}AL583917 to AL583926) (4) and strain Br4923 from Brazil (GenBank accession number {"type":"entrez-nucleotide","attrs":{"text":"FM211192","term_id":"219932450","term_text":"FM211192"}}FM211192) (N. Honore et al., unpublished data) share 99.98% identity.Like M. leprae, M. lepromatosis has not been cultivated on artificial media. In addition, our previous study also showed other similarities between these organisms, such as degeneration of mmaA3 into a pseudogene, the presence of unique AT-rich inserted sequences in the 16S rRNA gene, identical six-base tandem repeats in rpoT, similar G+C contents, and great evolutionary distance from other mycobacteria (17).The M. leprae genome (3.3 Mb) is much smaller than the Mycobacterium tuberculosis genome (4.4 Mb) (3, 4). More intriguingly, the M. leprae genome has undergone reductive evolution; ∼40% of the genes are inactivated (4), and ∼50% of the genes of the last common ancestor of M. leprae and M. tuberculosis have been lost (13). On the other hand, the M. leprae genome has been far more stable than the M. tuberculosis genome, and the worldwide clonality of the M. leprae strains paralleled the global spread of M. leprae strains that occurred via human activity and migration during the last ∼100,000 years (24). Recently, by comparing the genomes of M. leprae and M. tuberculosis and by analyzing the ages of the M. leprae pseudogenes, Gomez-Valero et al. (13) estimated that a massive gene inactivation event took place in the M. leprae genome in the last 20 million years.The discovery of M. lepromatosis and its differences from M. leprae make it relevant for further study for diagnosis, treatment, and prevention of DLL. Likewise, the many similarities between these two organisms prompted questions about their evolutionary histories and about how M. lepromatosis became endemic mainly in Mexico, while M. leprae occurs worldwide. In this study, we extended and refined our previous phylogenetic study by determining and analyzing the sequences of 20 genes and pseudogenes of M. lepromatosis. Our findings solidified the phylogeny of this new organism and provided new insights into the history of pseudogenes.     

Understanding the recent colonization history of a plant pathogenic fungus using population genetic tools and Approximate Bayesian Computation

Understanding the processes by which new diseases are introduced in previously healthy areas is of major interest in elaborating prevention and management policies, as well as in understanding the dynamics of pathogen diversity at large spatial scale. In this study, we aimed to decipher the dispersal processes that have led to the emergence of the plant pathogenic fungus Microcyclus ulei, which is responsible for the South American Leaf Blight (SALB). This fungus has devastated rubber tree plantations across Latin America since the beginning of the twentieth century. As only imprecise historical information is available, the study of population evolutionary history based on population genetics appeared most appropriate. The distribution of genetic diversity in a continental sampling of four countries (Brazil, Ecuador, Guatemala and French Guiana) was studied using a set of 16 microsatellite markers developed specifically for this purpose. A very strong genetic structure was found (F_st=0.70), demonstrating that there has been no regular gene flow between Latin American M. ulei populations. Strong bottlenecks probably occurred at the foundation of each population. The most likely scenario of colonization identified by the Approximate Bayesian Computation (ABC) method implemented in 𝒟ℐ𝒴𝒜ℬ𝒞 suggested two independent sources from the Amazonian endemic area. The Brazilian, Ecuadorian and Guatemalan populations might stem from serial introductions through human-mediated movement of infected plant material from an unsampled source population, whereas the French Guiana population seems to have arisen from an independent colonization event through spore dispersal.     

On the Use of Bayesian Methods for Evaluating Compartmental Neural Models

Computational modeling is being used increasingly in neuroscience. In deriving such models, inference issues such as model selection, model complexity, and model comparison must be addressed constantly. In this article we present briefly the Bayesian approach to inference. Under a simple set of commonsense axioms, there exists essentially a unique way of reasoning under uncertainty by assigning a degree of confidence to any hypothesis or model, given the available data and prior information. Such degrees of confidence must obey all the rules governing probabilities and can be updated accordingly as more data becomes available. While the Bayesian methodology can be applied to any type of model, as an example we outline its use for an important, and increasingly standard, class of models in computational neuroscience—compartmental models of single neurons. Inference issues are particularly relevant for these models: their parameter spaces are typically very large, neurophysiological and neuroanatomical data are still sparse, and probabilistic aspects are often ignored. As a tutorial, we demonstrate the Bayesian approach on a class of one-compartment models with varying numbers of conductances. We then apply Bayesian methods on a compartmental model of a real neuron to determine the optimal amount of noise to add to the model to give it a level of spike time variability comparable to that found in the real cell.     

Relation of Infection to Tissue Temperature in Mice Infected with Mycobacterium marinum and Mycobacterium leprae

Intravenous and footpad infections with Mycobacterium marinum and footpad infections with M. leprae were compared in the following mouse strains: A/He, BALB/C, CBA, C3H, C57BL, C57L, DBA, 101, and CFW. The results varied a great deal according to mouse strain used. Intravenous injection of high doses of M. marinum resulted in deaths after 28 days of 100% of strain A/He, and none of strain 101; 27 days after injection, the feet and noses of all strain CBA mice, but few of the C57BL, 101, or CFW mice, were involved. Injection of a small dose of M. marinum into the footpad produced visible disease in 5 days in all of the C57BL and 101 mice, but in not more than 60% of the A/He, DBA, and CFW mice; the average amount of swelling at 17 days varied from 4.40 mm in strain C57L to 0.92 in strain 101. After footpad injection of M. leprae, the average plateau harvests varied from 1.3 x 10(7) acid-fast bacteria in strain CBA to 6.5 x 10(5) in strain C57L. The infections in CBA mice extended from the site of inoculation throughout the foot. The temperature was measured rectally, in the footpad, and in the tail. Analysis of all the results revealed little correlation among the three types of infection. There was a strong negative correlation between the tail temperature and the death rate after intravenous injection of M. marinum, and a strong positive correlation between footpad temperature and plateau harvest of M. leprae.     

Use of Stochastic Models To Assess the Effect of Environmental Factors on Microbial Growth

We present a novel application of a stochastic ecological model to the study and analysis of microbial growth dynamics as influenced by environmental conditions in an extensive experimental data set. The model proved to be useful in bridging the gap between theoretical ideas in ecology and an applied problem in microbiology. The data consisted of recorded growth curves of Escherichia coli grown in triplicate in a base medium with all 32 possible combinations of five supplements: glucose, NH₄Cl, HCl, EDTA, and NaCl. The potential complexity of 2⁵ experimental treatments and their effects was reduced to 2² as just the metal chelator EDTA, the presumed osmotic pressure imposed by NaCl, and the interaction between these two factors were enough to explain the variability seen in the data. The statistical analysis showed that the positive and negative effects of the five chemical supplements and their combinations were directly translated into an increase or decrease in time required to attain stationary phase and the population size at which the stationary phase started. The stochastic ecological model proved to be useful, as it effectively explained and summarized the uncertainty seen in the recorded growth curves. Our findings have broad implications for both basic and applied research and illustrate how stochastic mathematical modeling coupled with rigorous statistical methods can be of great assistance in understanding basic processes in microbial ecology.     

Genetic control of murine T cell proliferative responses to Mycobacterium leprae. V. Evidence for cross-reactivity between host antigens and Mycobacterium leprae

T cell proliferative responses to Mycobacterium leprae were measured by immunization of mice at the base of the tail with Ag and challenging lymphocytes from draining lymph nodes in culture with M. leprae. C57BL/10J and B10.BR mice were identified as low responder mice and the congenic strains B10.M, B10.Q, and B10.AKM as high responders whereas F1 (high x low) hybrid mice were found to be low responders. The cellular basis of low responsiveness did not appear to result from a defect in Ag-presenting cells or the activation of suppressor T cells by M. leprae. The influence of the environment in which T cells developed on responsiveness to M. leprae was analyzed in chimeric mice prepared by irradiating F1(C57BL/10J x B10.M) mice and reconstituting with bone marrow from C57BL/10J, B10.M, or F1 donors. Six weeks later, chimeric mice were immunized with M. leprae, lymph node cells were subsequently prepared, and H-2 phenotyped and challenged in culture with M. leprae Ag. T cell proliferative responses were found to be low in all cases, similar to those observed using lymph node cells from F1 hybrid mice. These results suggested that high responder B10.M lymphocytes developing in the irradiated F1 mice became tolerized to antigenic determinants found on M. leprae. This implied cross-reactive epitopes existed between some mouse strains and M. leprae. Low responsiveness to M. leprae in low responder and F1 hybrid mice may result from tolerance to H-2-encoded Ag that show cross-reactivity with M. leprae.     

Comparison of the UDP-N-acetylmuramate:L-alanine ligase enzymes from Mycobacterium tuberculosis and Mycobacterium leprae

In the peptidoglycan of Mycobacterium leprae, L-alanine of the side chain is replaced by glycine. When expressed in Escherichia coli, MurC (UDP-N-acetyl-muramate:L-alanine ligase) of M. leprae showed K(m) and V(max) for L-alanine and glycine similar to those of Mycobacterium tuberculosis MurC, suggesting that another explanation should be sought for the presence of glycine.     

Quantitative comparison of the mycolic and fatty acid compositions of Mycobacterium leprae and Mycobacterium gordonae

The mycolic and fatty acids of three samples each of Mycobacterium leprae and Mycobacterium gordonae were compared. Acids released by whole-organism alkaline hydrolysis were converted to 4-nitrobenzyl esters and mycolic acids were further derivatized to t-butyldimethylsilyl ethers. Thin-layer chromatography of the derivatized long-chain extracts showed that all three M. leprae preparations contained so-called alpha-mycolates and ketomycolates but that the M. gordonae samples had a methoxymycolate in addition to the above types. Silica gel normal-phase high-performance liquid chromatography of the total mycolic acid derivatives confirmed the lack of detectable amounts of methoxymycolates in M. leprae and reverse-phase chromatography of the individual mycolate types demonstrated the homogeneity of the chain lengths of the mycolic acids in each species. Non-hydroxylated fatty acid 4-nitrobenzyl esters were transformed to methyl esters and examined by gas chromatography. Tuberculostearic (10-methyloctadecanoic) acid was a major component of the lipids of all three M. leprae preparations but it was absent in one M. gordonae strain and a very minor component in the other representatives of this latter species. On the basis of fatty and mycolic acid compositions, therefore, a previously suggested close relationship between M. leprae and M. gordonae was not supported.     

Quantitative Reconstruction of Weaning Ages in Archaeological Human Populations Using Bone Collagen Nitrogen Isotope Ratios and Approximate Bayesian Computation

BackgroundNitrogen isotope analysis of bone collagen has been used to reconstruct the breastfeeding practices of archaeological human populations. However, weaning ages have been estimated subjectively because of a lack of both information on subadult bone collagen turnover rates and appropriate analytical models.MethodologyTemporal changes in human subadult bone collagen turnover rates were estimated from data on tissue-level bone metabolism reported in previous studies. A model for reconstructing precise weaning ages was then developed using a framework of approximate Bayesian computation and incorporating the estimated turnover rates. The model is presented as a new open source R package, WARN (Weaning Age Reconstruction with Nitrogen isotope analysis), which computes the age at the start and end of weaning, ¹⁵N-enrichment through maternal to infant tissue, and value of collagen synthesized entirely from weaning foods with their posterior probabilities. The model was applied to 39 previously reported Holocene skeletal populations from around the world, and the results were compared with weaning ages observed in ethnographic studies.ConclusionsThere were no significant differences in the age at the end of weaning between the archaeological (2.80±1.32 years) and ethnographic populations. By comparing archaeological populations, it appears that weaning ages did not differ with the type of subsistence practiced (i.e., hunting–gathering or not). Most of -enrichment (2.44±0.90‰) was consistent with biologically valid values. The nitrogen isotope ratios of subadults after the weaning process were lower than those of adults in most of the archaeological populations (−0.48±0.61‰), and this depletion was greater in non-hunter–gatherer populations. Our results suggest that the breastfeeding period in humans had already been shortened by the early Holocene compared with those in extant great apes.     

Local comparison of the genomes of Mycobacterium tuberculosis and Mycobacterium leprae using the polymerase chain reaction

Abstract A gene fusion between the Saccharomyces cerevisiae actin gene promoter and the cDNA of the Fusarium solani f. sp. pisi pelA gene has been constructed. This expression cassette has been introduced into the industrial wine yeast strain T₇₃. The resulting recombinant strain is able to secrete active PELA enzyme into the culture medium. In preliminary microvinification experiments the wine produced by this pectinolytic strain is indistinguishable from wine produced using the non-transformed strain on the basis of the chemical analyses. Large scale fermentations need to be carried out in order to assess the effects on filtrability.     

Characterization and taxonomic implications of the rRNA genes of Mycobacterium leprae.

The number of rRNA genes of Mycobacterium leprae was determined by restriction analysis of M. leprae total chromosomal DNA. A single set of rRNA genes was found. This set was subcloned from a cosmid library of M. leprae DNA into pUC13 and was characterized by restriction analysis and hybridization with Escherichia coli rRNA genes. The 16S, 23S, and 5S genes of M. leprae were clustered on a 5.3-kilobase DNA fragment. On one hand, restriction analysis of the set of rRNA genes showed the uniqueness of M. leprae among mycobacteria, but on the other hand, it suggested that M. leprae strains of several origins are very much alike. Quantitative hybridization studies between M. leprae rDNA and total DNA of various bacteria demonstrated a close relatedness between M. leprae and corynebacteria, nocardia, and mycobacteria, especially Mycobacterium tuberculosis.     

Use of Glycated Hemoglobin to Assess Glycemic Control in Wistar Diabetic Fatty Rats and Zucker Fatty Rats

The Wistar Diabetic Fatty rat (WDF fafa) is a con-genic strain of the Wistar Kyoto rat. Studies using blood glucose reveal that only fatty male (not female) WDF rats spontaneously develop hyperglycemia when fed a stock diet Blood glucose values have not provided consistent results for evaluation of glycemic status in fatty male WDF rats. Zucker fatty (fafa) rats, while sharing the fa gene and the development of hyperinsulinemia and hyperlipemia, do not spontaneously become hyperglycemic. In order to examine strain differences and the effects of age on long-term average glycemic status in WDF and Zucker rats, glycated hemoglobin (GHb) was analyzed. Glycated hemoglobin was measured in male lean and obese WDF and Zucker rats at 2,3,6, and 12 months of age. Nonfasted plasma glucose was measured in male lean and obese WDF rats at 2, 3, 6, and 12 months of age and in lean and obese Zucker rats at 3, 6, and 12 months of age. Plasma insulin was measured in lean and obese WDF and Zucker rats at 3, 6, and 12 months of age. Obese WDF rats had significantly elevated GHb compared to lean controls at 3, 6, and 12 months of age. Glycated hemoglobin was substantially above the normal range (3.8-6.5%) at 3 months of age (14.1%). Glycated hemoglobin significantly declined in the obese WDF rats between 6 and 12 months of age. Nonfasted plasma glucose was significantly elevated in the obese WDF rats at 3 months (14.1 ± 2.1 mM/L) and 6 months of age (16.2 ± 2.3 mM/L) compared to lean controls. At 12 months of age there was no difference in plasma glucose between obese and lean WDF rats. Obese and lean Zucker rats had similar levels of GHb and plasma glucose at all ages. In conclusion, GHb provides more integrated data for classifying disease status of WDF rats and evaluation of potential long-term complications associated with hyperglycemia.     

Reinvestigation on the causes of genomic GC variation between the orthologous genes of Mycobacterium tuberculosis and Mycobacterium leprae

Genomic GC (overall G+C content of the coding sequences) variations were reinvestigated between the orthologous genes of Mycobacterium tuberculosis and Mycobacterium leprae species. It was observed that overall genomic GC variation between the species mainly originates from the combined effects GC(1) and GC(2) variations. But codons having identical amino acids with different codons (IA) (between the orthologous codon pairs) are responsible for the genomic GC(3) variation between the organisms, whereas orthologous codons having different amino acids (DA) between the two organisms are responsible for the variation of GC(1) levels. Further analyses indicate that duets and quartets are going in the same direction with same magnitude in changing the GC(3) levels for IA category, whereas GC(1) levels of duets of DA category decreases significantly from the overall GC(1) levels but GC(1) levels of quartets increases significantly from the overall GC(1) levels. GC(3) levels of informational genes for the IA category decrease more rapidly than the other functional categories of genes. The biological implications of these results have been discussed in this paper.