The organisation and expression of the genes encoding the mitochondrial glycine decarboxylase complex and serine hydroxymethyltransferase in pea (Pisum sativum)

Restriction fragment length polymorphisms have been used to determine the chromosomal location of the genes encoding the glycine decarboxylase complex (GDC) and serine hydroxymethyltransferase (SHMT) of pea leaf mitochondria. The genes encoding the H subunit of GDC and the genes encoding SHMT both show linkage to the classical group I marker i. In addition, the genes for the P protein of GDC show linkage to the classic group I marker a. The genes for the L and T proteins of GDC are linked to one another and are probably situated on the satellite of chromosome 7. The mRNAs encoding the five polypeptides that make up GDC and SHMT are strongly induced when dark-grown etiolated pea seedlings are placed in the light. Similarly, when mature plants are placed in the dark for 48 h, the levels of both GDC protein and SHMT mRNAs decline dramatically and then are induced strongly when these plants are returned to the light. During both treatments a similar pattern of mRNA induction is observed, with the mRNA encoding the P protein of GDC being the most rapidly induced and the mRNA for the H protein the slowest. Whereas during the greening of etiolated seedlings the polypeptides of GDC and SHMT show patterns of accumulation similar to those of the corresponding mRNAs, very little change in the level of the polypeptides is seen when mature plants are placed in the dark and then re-exposed to the light.

Single nucleotide polymorphisms and mRNA expression of VEGF-A in papillary thyroid carcinoma: Potential markers for aggressive phenotypes

BACKGROUND AND OBJECTIVES Polymorphisms of the VEGF gene are known to affect the biological behaviour of cancers but have seldom been studied in thyroid cancer. The aim of the current study is to evaluate the prevalence and relevance of VEGF-A polymorphisms and mRNA expression in papillary thyroid carcinoma (PTC). MATERIALS AND METHODS Genomic DNA and total RNA were isolated from paraffin-embedded tissue from 91 PTC (51 conventional PTC and 40 follicular variant) and 78 control thyroid tissues. Three DNA polymorphisms (+936C > T, +405C > G and -141A > C) in the 3' and 5' untranslated region (3'-UTR, 5'-UTR) of VEGF-A were studied using PCR and RFLP. Also, the mRNA expression of VEGF-A in these tissues was studied by real-time PCR. RESULTS Distribution of polymorphisms in the 5'-UTR (VEGF-VEGF -141A > C and +405C > G) and 3'-UTR (VEGF +936C > T) were all significantly different in PTC and benign thyroid tissue (p = 0.0001, 0.001 and 0.028 respectively). The VEGF -141 C allele was more common in PTC with lymph node metastases (p = 0.026). VEGF + 405 Galleles andVEGF +936 CC genotype were more common in PTC of advanced pathological staging (p = 0.018 and 0.017 respectively). Also, increased expression of VEGF-A mRNA was noted in PTC compared to control (p = 0.009). Within the group of patients with conventional PTC, those with lymph nodal metastases had a higher level of VEGF-A mRNA expression than other patients (p = 0.0003). CONCLUSION These findings suggest that VEGF polymorphisms and mRNA expression may predict the aggressiveness behaviour of thyroid cancer.

Association of the SNP rs2623047 in the HSPG modification enzyme SULF1 with an Australian Caucasian Breast Cancer Cohort

Breast cancer is the second most common cancer worldwide and the most common cancer reported in women. This malignant tumour is characterised by a number of specific features including uncontrolled cell proliferation. It ranks fifth in the world as a cause of cancer death in women. Early diagnosis increases 5 year survival rates up to 95%. Heparan sulfate proteoglycans (HSPGs) are complex proteins composed of a core protein to which a number of highly sulfated side chains are synthesised by a highly co-ordinated process resulting in distinct sulfation patterns, which determine specific interations with cell-signaling partners including growth factors, their receptors, ligands and morphogens. The enzymes responsible for chain initiation, elongation and sulfation are critical for creating HS chain variability conferring biological functionality. This study investigated single nucleotide polymorphism in SULF1, the enzyme responsible for the 6-0 desulfation of heparan sulfate side chains. We investigated this SNP in an Australian Caucasian case-control breast cancer population and found a significant association between SULF1 and breast cancer at both the allelic and genotypic level (allele, p=0.016; genotype, p=0.032). Our results suggest the res2623047 SNP in SULF1 may impact breast cancer susceptibility. Specifically, the T allele of rs2623047 in SULF1 is associated with a increased risk of developing breast cancer in our cohort. The identification of markers including SULF1 may improve detection of this disease at its earliest stages improving patient treatment and prognosis.

Genetic testing for haemochromatosis in patients with chondrocalcinosis

Hereditary haemochromatosis (HH) is the most common lethal monogenic human disease, affecting roughly 1 in 300 white northern Europeans. Homozygosity for the C282Y polymorphism within the HFE gene causes more than 80% of cases, with compound heterozygosity of the C282Y and H63D polymorphism also increasing susceptibility to disease. The aim of this study was to determine the frequency of the C282Y and H63D polymorphisms in the disease, and to assess the risk of HH in heterozygotes for the C282Y polymorphism. 128 patients were recruited because of either radiographic chondrocalcinosis (at least bicompartmental knee disease or joints other than the knee involved) or CPPD pseudogout. Genotyping of the HFE C282Y and H63D mutations was performed using PCR/SSP and genotypes for the C282Y polymorphism confirmed by PCR/RFLP. Historical white European control data were used for comparison. Two previously undiagnosed C282Y homozygotes (1.6%), and 16 C282Y heterozygotes (12.5%), including four (3.1%) C282Y/ H63D compound heterozygotes were identified. This represents a significant overrepresentation of C282Y homozygotes (relative risk 3.4, p-0.037), but the number of heterozygotes was not significantly increased. At a cost per test of £1 for each subject, screening all patients with chondrocalcinosis using the above ascertainment criteria costs only £64 for each case of haemochromatosis identified, clearly a highly cost effective test given the early mortality associated with untreated haemochromatosis. Routine screening for haemochromatosis in patients with appreciable chondrocatcinosis is recommended.

Mapping genes for osteoporosis-Old dogs and new tricks

In stark contrast to its horticultural origins, modern genetics is an extremely technology-driven field. Almost all the major advances in the field over the past 20 years have followed technological developments that have permitted change in study designs. The development of PCR in the 1980s led to RFLP mapping of monogenic diseases. The development of fluorescent-tagged genotyping methods led to linkage mapping approaches for common diseases that dominated the 1990s. The development of microarray SNP genotyping has led to the genome-wide association study era of the new millennium. And now the development of next-generation sequencing technologies is about to open up a new era of gene-mapping, enabling many potential new study designs. This review aims to present the strengths and weaknesses of the current approaches, and present some new ideas about gene-mapping approaches that are likely to advance our knowledge of the genes involved in heritable bone traits such as bone mineral density (BMD) and fracture.

Assessing the relationship between patch type and soil mites: A molecular approach

An urgent need exists for indicators of soil health and patch functionality in extensive rangelands that can be measured efficiently and at low cost. Soil mites are candidate indicators, but their identification and handling is so specialised and time-consuming that their inclusion in routine monitoring is unlikely. The aim of this study was to measure the relationship between patch type and mite assemblages using a conventional approach. An additional aim was to determine if a molecular approach traditionally used for soil microbes could be adapted for soil mites to overcome some of the bottlenecks associated with soil fauna diversity assessment. Soil mite species abundance and diversity were measured using conventional ecological methods in soil from patches with perennial grass and litter cover (PGL), and compared to soil from bare patches with annual grasses and/or litter cover (BAL). Soil mite assemblages were also assessed using a molecular method called terminal-restriction fragment length polymorphism (T-RFLP) analysis. The conventional data showed a relationship between patch type and mite assemblage. The Prostigmata and Oribatida were well represented in the PGL sites, particularly the Aphelacaridae (Oribatida). For T-RFLP analysis, the mite community was represented by a series of DNA fragment lengths that reflected mite sequence diversity. The T-RFLP data showed a distinct difference in the mite assemblage between the patch types. Where possible, T-RFLP peaks were matched to mite families using a reference 18S rDNA database, and the Aphelacaridae prevalent in the conventional samples at PGL sites were identified, as were prostigmatids and oribatids. We identified limits to the T-RFLP approach and this included an inability to distinguish some species whose DNA sequences were similar. Despite these limitations, the data still showed a clear difference between sites, and the molecular taxonomic inferences also compared well with the conventional ecological data. The results from this study indicated that the T-RFLP approach was effective in measuring mite assemblages in this system. The power of this technique lies in the fact that species diversity and abundance data can be obtained quickly because of the time taken to process hundreds of samples, from soil DNA extraction to data output on the gene analyser, can be as little as 4 days.

Genetic population structure of red snappers (Lutjanus malabaricus Bloch & Schneider, 1801 and Lutjanus erythropterus Bloch, 1790) in central and eastern Indonesia and northern Australia

The genetic population structure of red snapper Lutjanus malabaricus and Lutjanus erythropterus in eastern Indonesia and northern Australia was investigated by allozyme electrophoresis and sequence variation in the control region of mtDNA. Samples were collected from eight sites in Indonesia and four sites in northern Australia for both species. A total of 13 allozyme loci were scored. More variable loci were observed in L. malabaricus than in L. erythropterus. Sequence variation in the control region (left domain) of the mitochondrial genome was assessed by RFLP and direct sequencing. MtDNA haplotype diversity was high (L. erythropterus, 0.95 and L. malabaricus, 0.97), as was intraspecific sequence divergence, (L. erythropterus, 0.0-12.5% and L. malabaricus, 0.0-9.5%). The pattern of mtDNA haplotype frequencies grouped both species into two broad fisheries stocks with a genetic boundary either between Kupang and Sape (L. malabaricus) or between Kupang and Australian Timor Sea (L. erythropertus). The allozyme analyses revealed similar boundaries for L. erythropterus. Seven allozymes stocks compared to two mtDNA stocks of L. malabaricus including Ambon, which was not sampled with mtDNA, however, were reported. Possible reasons for differences in discrimination between the methods include: i) increased power of multiple allozyme loci over the single mtDNA locus, ii) insufficient gene sampling in the mtDNA control region and iii) relative evolutionary dynamics of nuclear (allozyme loci) and mitochondrial DNA in these taxa. Allozyme and haplotype data did not distinguish separate stocks among the four Australian locations nor the central Indonesian (Bali and Sape locations) for both L. malabaricus and L. erythropterus.

Characterisation of chicken Campylobacter jejuni isolates using resolution optimised single nucleotide polymorphisms and binary gene markers

The principal objective of this study was to determine if Campylobacter jejuni genotyping methods based upon resolution optimised sets of single nucleotide polymorphisms (SNPs) and binary genetic markers were capable of identifying epidemiologically linked clusters of chicken-derived isolates. Eighty-eight C. jejuni isolates of known flaA RFLP type were included in the study. They encompassed three groups of ten isolates that were obtained at the same time and place and possessed the same flaA type. These were regarded as being epidemiologically linked. Twenty-six unlinked C. jejuni flaA type I isolates were included to test the ability of SNP and binary typing to resolve isolates that were not resolved by flaA RFLP. The remaining isolates were of different flaA types. All isolates were typed by real-time PCR interrogation of the resolution optimised sets of SNPs and binary markers. According to each typing method, the three epidemiologically linked clusters were three different clones that were well resolved from the other isolates. The 26 unlinked C. jejuni flaA type I isolates were resolved into 14 SNP-binary types, indicating that flaA typing can be unreliable for revealing epidemiological linkage. Comparison of the data with data from a fully typed set of isolates associated with human infection revealed that abundant lineages in the chicken isolates that were also found in the human isolates belonged to clonal complex (CC) -21 and CC-353, with the usually rare C-353 member ST-524 being especially abundant in the chicken collection. The chicken isolates selected to be diverse according to flaA were also diverse according to SNP and binary typing. It was observed that CC-48 was absent in the chicken isolates, despite being very common in Australian human infection isolates, indicating that this may be a major cause of human disease that is not chicken associated.

Mapping quantitative trait loci for resistance to Pratylenchus thornei from synthetic hexaploid wheat in the International Triticeae Mapping Initiative (ITMI) population

Root-lesion nematode (Pratylenchus thornei) is a serious pathogen of wheat in many countries. The International Triticeae Mapping Initiative (ITMI) population of recombinant inbred lines (RILs) was assessed for resistance to P. thornei to determine the chromosome locations of the resistance genes. The ITMI population is derived from a cross between the resistant synthetic hexaploid wheat W-7984 and a susceptible bread wheat cultivar Opata 85. Two years of phenotypic data for resistance to P. thornei were obtained in replicated glasshouse trials. Quantitative trait locus (QTL) analysis was performed using available segregation and map data for 114 RILs. A QTL on chromosome 6DS showed consistent effects for reduced nematode numbers (partial resistance) across years and accounted for 11% and 23% of the phenotypic variation. A second QTL for P. thornei resistance on chromosome 2BS accounted for an additional 19% and 5%. Restriction fragment length polymorphism (RFLP) and simple sequence repeat (SSR) markers associated with the QTLs are physically located in regions rich in major genes at the distal ends of the short chromosome arms of 6D and 2B. SSR markers with potential for marker-assisted selection of P. thornei resistance effective in different genetic backgrounds have been identified.

The assessment of TGGE for the detection of interspecific and intergeneric DNA-marker polymorphism within Solanaceae

RFLP markers are currently the most appropriate marker system for the identification of uncharacterised polymorphism at the interspecific and intergeneric level. Given the benefits of a PCR-based marker system and the availability of sequence information for many Solanaceous cDNA clones, it is now possible to target conserved fragments, for primer development, that flank sequences possessing interspecific polymorphism. The potential outcome is the development of a suite of markers that amplify widely in Solanaceae. Temperature gradient gel electrophoresis (TGGE) is a relatively inexpensive gel-based system that is suitable for the detection of most single-base changes. TGGE can be used to screen for both known and unknown polymorphisms, and has been assessed here, for the development of PCR-based markers that are useful for the detection of interspecific variation within Solanaceae. Fifteen markers are presented where differences between Lycopersicon esculentum and L. pennellii have been detected by TGGE. The markers were assessed on a wider selection of plant species and found to be potentially useful for the identification of interspecific and intergeneric polymorphism in Solanaceous plants.

Molecular characterisation, pathogenesis and fungicide sensitivity of Pythium spp. from table beet (Beta vulgaris var. vulgaris) grown in the Lockyer Valley, Queensland

Table beet production in the Lockyer Valley of south-eastern Queensland is known to be adversely affected by soilborne root disease from infection by Pythium spp. However, little is known regarding the species or genotypes that are the causal agents of both pre- and post-emergence damping off. Based on RFLP analysis with HhaI, HinfI and MboI of the PCR amplified ITS region DNA from soil and diseased plant samples, the majority of 130 Pythium isolates could be grouped into three genotypes, designated LVP A, LVP B and LVP C. These groups comprised 43, 41 and 7% of all isolates, respectively. Deoxyribonucleic acid sequence analysis of the ITS region indicated that LVP A was a strain of Pythium aphanidermatum, with greater than 99% similarity to the corresponding P. aphanidermatum sequences from the publicly accessible databases. The DNA sequences from LVP B and LVP C were most closely related to P. ultimum and P. dissotocum, respectively. Lower frequencies of other distinct isolates with unique RFLP patterns were also obtained with high levels of similarity (>97%) to P. heterothallicum, P. periplocum and genotypes of P. ultimum other than LVP B. Inoculation trials of 1- and 4-week-old beet seedlings indicated that compared with isolates of the LVP B genotype, a higher frequency of LVP A isolates caused disease. Isolates with the LVP A, LVP B and LVP C genotypes were highly sensitive to the fungicide Ridomil MZ, which suppressed radial growth on V8 agar between approximately four and thirty fold at 5 μg/mL metalaxyl and 40 μg/mL mancozeb, a concentration far lower than the recommended field application rate.

A consensus genetic map of sorghum that integrates multiple component maps and high-throughput Diversity Array Technology (DArT) markers

Background: Sorghum genome mapping based on DNA markers began in the early 1990s and numerous genetic linkage maps of sorghum have been published in the last decade, based initially on RFLP markers with more recent maps including AFLPs and SSRs and very recently, Diversity Array Technology (DArT) markers. It is essential to integrate the rapidly growing body of genetic linkage data produced through DArT with the multiple genetic linkage maps for sorghum generated through other marker technologies. Here, we report on the colinearity of six independent sorghum component maps and on the integration of these component maps into a single reference resource that contains commonly utilized SSRs, AFLPs, and high-throughput DArT markers. Results: The six component maps were constructed using the MultiPoint software. The lengths of the resulting maps varied between 910 and 1528 cM. The order of the 498 markers that segregated in more than one population was highly consistent between the six individual mapping data sets. The framework consensus map was constructed using a "Neighbours" approach and contained 251 integrated bridge markers on the 10 sorghum chromosomes spanning 1355.4 cM with an average density of one marker every 5.4 cM, and were used for the projection of the remaining markers. In total, the sorghum consensus map consisted of a total of 1997 markers mapped to 2029 unique loci ( 1190 DArT loci and 839 other loci) spanning 1603.5 cM and with an average marker density of 1 marker/0.79 cM. In addition, 35 multicopy markers were identified. On average, each chromosome on the consensus map contained 203 markers of which 58.6% were DArT markers. Non-random patterns of DNA marker distribution were observed, with some clear marker-dense regions and some marker-rare regions. Conclusion: The final consensus map has allowed us to map a larger number of markers than possible in any individual map, to obtain a more complete coverage of the sorghum genome and to fill a number of gaps on individual maps. In addition to overall general consistency of marker order across individual component maps, good agreement in overall distances between common marker pairs across the component maps used in this study was determined, using a difference ratio calculation. The obtained consensus map can be used as a reference resource for genetic studies in different genetic backgrounds, in addition to providing a framework for transferring genetic information between different marker technologies and for integrating DArT markers with other genomic resources. DArT markers represent an affordable, high throughput marker system with great utility in molecular breeding programs, especially in crops such as sorghum where SNP arrays are not publicly available.

Campylobacter jejuni and Campylobacter coli Genotyping by High-Resolution Melting Analysis of a flaA Fragment.

The highly variable flagellin-encoding flaA gene has long been used for genotyping Campylobacter jejuni and Campylobacter coli. High-resolution melting (HRM) analysis is emerging as an efficient and robust method for discriminating DNA sequence variants. The objective of this study was to apply HRM analysis to flaA-based genotyping. The initial aim was to identify a suitable flaA fragment. It was found that the PCR primers commonly used to amplify the flaA short variable repeat (SVR) yielded a mixed PCR product unsuitable for HRM analysis. However, a PCR primer set composed of the upstream primer used to amplify the fragment used for flaA restriction fragment length polymorphism (RFLP) analysis and the downstream primer used for flaA SVR amplification generated a very pure PCR product, and this primer set was used for the remainder of the study. Eighty-seven C. jejuni and 15 C. coli isolates were analyzed by flaA HRM and also partial flaA sequencing. There were 47 flaA sequence variants, and all were resolved by HRM analysis. The isolates used had previously also been genotyped using single-nucleotide polymorphisms (SNPs), binary markers, CRISPR HRM, and flaA RFLP.flaA HRM analysis provided resolving power multiplicative to the SNPs, binary markers, and CRISPR HRM and largely concordant with the flaA RFLP. It was concluded that HRM analysis is a promising approach to genotyping based on highly variable genes.

Genotyping for genetic association studies: Methods and applications

In this thesis, two separate single nucleotide polymorphism (SNP) genotyping techniques were set up at the Finnish Genome Center, pooled genotyping was evaluated as a screening method for large-scale association studies, and finally, the former approaches were used to identify genetic factors predisposing to two distinct complex diseases by utilizing large epidemiological cohorts and also taking environmental factors into account. The first genotyping platform was based on traditional but improved restriction-fragment-length-polymorphism (RFLP) utilizing 384-microtiter well plates, multiplexing, small reaction volumes (5 µl), and automated genotype calling. We participated in the development of the second genotyping method, based on single nucleotide primer extension (SNuPeTM by Amersham Biosciences), by carrying out the alpha- and beta tests for the chemistry and the allele-calling software. Both techniques proved to be accurate, reliable, and suitable for projects with thousands of samples and tens of markers. Pooled genotyping (genotyping of pooled instead of individual DNA samples) was evaluated with Sequenom s MassArray MALDI-TOF, in addition to SNuPeTM and PCR-RFLP techniques. We used MassArray mainly as a point of comparison, because it is known to be well suited for pooled genotyping. All three methods were shown to be accurate, the standard deviations between measurements being 0.017 for the MassArray, 0.022 for the PCR-RFLP, and 0.026 for the SNuPeTM. The largest source of error in the process of pooled genotyping was shown to be the volumetric error, i.e., the preparation of pools. We also demonstrated that it would have been possible to narrow down the genetic locus underlying congenital chloride diarrhea (CLD), an autosomal recessive disorder, by using the pooling technique instead of genotyping individual samples. Although the approach seems to be well suited for traditional case-control studies, it is difficult to apply if any kind of stratification based on environmental factors is needed. Therefore we chose to continue with individual genotyping in the following association studies. Samples in the two separate large epidemiological cohorts were genotyped with the PCR-RFLP and SNuPeTM techniques. The first of these association studies concerned various pregnancy complications among 100,000 consecutive pregnancies in Finland, of which we genotyped 2292 patients and controls, in addition to a population sample of 644 blood donors, with 7 polymorphisms in the potentially thrombotic genes. In this thesis, the analysis of a sub-study of pregnancy-related venous thromboses was included. We showed that the impact of factor V Leiden polymorphism on pregnancy-related venous thrombosis, but not the other tested polymorphisms, was fairly large (odds ratio 11.6; 95% CI 3.6-33.6), and increased multiplicatively when combined with other risk factors such as obesity or advanced age. Owing to our study design, we were also able to estimate the risks at the population level. The second epidemiological cohort was the Helsinki Birth Cohort of men and women who were born during 1924-1933 in Helsinki. The aim was to identify genetic factors that might modify the well known link between small birth size and adult metabolic diseases, such as type 2 diabetes and impaired glucose tolerance. Among ~500 individuals with detailed birth measurements and current metabolic profile, we found that an insertion/deletion polymorphism of the angiotensin converting enzyme (ACE) gene was associated with the duration of gestation, and weight and length at birth. Interestingly, the ACE insertion allele was also associated with higher indices of insulin secretion (p=0.0004) in adult life, but only among individuals who were born small (those among the lowest third of birth weight). Likewise, low birth weight was associated with higher indices of insulin secretion (p=0.003), but only among carriers of the ACE insertion allele. The association with birth measurements was also found with a common haplotype of the glucocorticoid receptor (GR) gene. Furthermore, the association between short length at birth and adult impaired glucose tolerance was confined to carriers of this haplotype (p=0.007). These associations exemplify the interaction between environmental factors and genotype, which, possibly due to altered gene expression, predisposes to complex metabolic diseases. Indeed, we showed that the common GR gene haplotype associated with reduced mRNA expression in thymus of three individuals (p=0.0002).

Bayesian statistical analysis of bacterial diversity

Bacteria play an important role in many ecological systems. The molecular characterization of bacteria using either cultivation-dependent or cultivation-independent methods reveals the large scale of bacterial diversity in natural communities, and the vastness of subpopulations within a species or genus. Understanding how bacterial diversity varies across different environments and also within populations should provide insights into many important questions of bacterial evolution and population dynamics. This thesis presents novel statistical methods for analyzing bacterial diversity using widely employed molecular fingerprinting techniques. The first objective of this thesis was to develop Bayesian clustering models to identify bacterial population structures. Bacterial isolates were identified using multilous sequence typing (MLST), and Bayesian clustering models were used to explore the evolutionary relationships among isolates. Our method involves the inference of genetic population structures via an unsupervised clustering framework where the dependence between loci is represented using graphical models. The population dynamics that generate such a population stratification were investigated using a stochastic model, in which homologous recombination between subpopulations can be quantified within a gene flow network. The second part of the thesis focuses on cluster analysis of community compositional data produced by two different cultivation-independent analyses: terminal restriction fragment length polymorphism (T-RFLP) analysis, and fatty acid methyl ester (FAME) analysis. The cluster analysis aims to group bacterial communities that are similar in composition, which is an important step for understanding the overall influences of environmental and ecological perturbations on bacterial diversity. A common feature of T-RFLP and FAME data is zero-inflation, which indicates that the observation of a zero value is much more frequent than would be expected, for example, from a Poisson distribution in the discrete case, or a Gaussian distribution in the continuous case. We provided two strategies for modeling zero-inflation in the clustering framework, which were validated by both synthetic and empirical complex data sets. We show in the thesis that our model that takes into account dependencies between loci in MLST data can produce better clustering results than those methods which assume independent loci. Furthermore, computer algorithms that are efficient in analyzing large scale data were adopted for meeting the increasing computational need. Our method that detects homologous recombination in subpopulations may provide a theoretical criterion for defining bacterial species. The clustering of bacterial community data include T-RFLP and FAME provides an initial effort for discovering the evolutionary dynamics that structure and maintain bacterial diversity in the natural environment.