SNP genotyping reveals genetic diversity between cultivated landraces and contemporary varieties of tomato

Background

The tomato (Solanum lycopersium L.) is the most widely grown vegetable in the world. It was domesticated in Latin America and Italy and Spain are considered secondary centers of diversification. This food crop has experienced severe genetic bottlenecks and modern breeding activities have been characterized by trait introgression from wild species and divergence in different market classes.

Results

With the aim to examine patterns of polymorphism, characterize population structure and identify putative loci under positive selection, we genotyped 214 tomato accessions (which include cultivated landraces, commercial varieties and wild relatives) using a custom-made Illumina SNP-panel. Most of the 175 successfully scored SNP loci were found to be polymorphic. Population structure analysis and estimates of genetic differentiation indicated that landraces constitute distinct sub-populations. Furthermore, contemporary varieties could be separated in groups (processing, fresh and cherry) that are consistent with the recent breeding aimed at market-class specialization. In addition, at the 95% confidence level, we identified 30, 34 and 37 loci under positive selection between landraces and each of the groups of commercial variety (cherry, processing and fresh market, respectively). Their number and genomic locations imply the presence of some extended regions with high genetic variation between landraces and contemporary varieties.

Conclusions

Our work provides knowledge concerning the level and distribution of genetic variation within cultivated tomato landraces and increases our understanding of the genetic subdivision of contemporary varieties. The data indicate that adaptation and selection have led to a genomic signature in cultivated landraces and that the subpopulation structure of contemporary varieties is shaped by directed breeding and largely of recent origin. The genomic characterization presented here is an essential step towards a future exploitation of the available tomato genetic resources in research and breeding programs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-835) contains supplementary material, which is available to authorized users.     

Comprehensive genotyping of the USA national maize inbred seed bank

Background

Genotyping by sequencing, a new low-cost, high-throughput sequencing technology was used to genotype 2,815 maize inbred accessions, preserved mostly at the National Plant Germplasm System in the USA. The collection includes inbred lines from breeding programs all over the world.

Results

The method produced 681,257 single-nucleotide polymorphism (SNP) markers distributed across the entire genome, with the ability to detect rare alleles at high confidence levels. More than half of the SNPs in the collection are rare. Although most rare alleles have been incorporated into public temperate breeding programs, only a modest amount of the available diversity is present in the commercial germplasm. Analysis of genetic distances shows population stratification, including a small number of large clusters centered on key lines. Nevertheless, an average fixation index of 0.06 indicates moderate differentiation between the three major maize subpopulations. Linkage disequilibrium (LD) decays very rapidly, but the extent of LD is highly dependent on the particular group of germplasm and region of the genome. The utility of these data for performing genome-wide association studies was tested with two simply inherited traits and one complex trait. We identified trait associations at SNPs very close to known candidate genes for kernel color, sweet corn, and flowering time; however, results suggest that more SNPs are needed to better explore the genetic architecture of complex traits.

Conclusions

The genotypic information described here allows this publicly available panel to be exploited by researchers facing the challenges of sustainable agriculture through better knowledge of the nature of genetic diversity.     

Whole genome resequencing in tomato reveals variation associated with introgression and breeding events

Background

One of the goals of genomics is to identify the genetic loci responsible for variation in phenotypic traits. The completion of the tomato genome sequence and recent advances in DNA sequencing technology allow for in-depth characterization of genetic variation present in the tomato genome. Like many self-pollinated crops, cultivated tomato accessions show a low molecular but high phenotypic diversity. Here we describe the whole-genome resequencing of eight accessions (four cherry-type and four large fruited lines) chosen to represent a large range of intra-specific variability and the identification and annotation of novel polymorphisms.

Results

The eight genomes were sequenced using the GAII Illumina platform. Comparison of the sequences with the reference genome yielded more than 4 million single nucleotide polymorphisms (SNPs). This number varied from 80,000 to 1.5 million according to the accessions. Almost 128,000 InDels were detected. The distribution of SNPs and InDels across and within chromosomes was highly heterogeneous revealing introgressions from wild species and the mosaic structure of the genomes of the cherry tomato accessions. In-depth annotation of the polymorphisms identified more than 16,000 unique non-synonymous SNPs. In addition 1,686 putative copy-number variations (CNVs) were identified.

Conclusions

This study represents the first whole genome resequencing experiment in cultivated tomato. Substantial genetic differences exist between the sequenced tomato accessions and the reference sequence. The heterogeneous distribution of the polymorphisms may be related to introgressions that occurred during domestication or breeding. The annotated SNPs, InDels and CNVs identified in this resequencing study will serve as useful genetic tools, and as candidate polymorphisms in the search for phenotype-altering DNA variations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-791) contains supplementary material, which is available to authorized users.     

Genetic Variation and Recent Positive Selection in Worldwide Human Populations: Evidence from Nearly 1 Million SNPs

Background

Genome-wide scans of hundreds of thousands of single-nucleotide polymorphisms (SNPs) have resulted in the identification of new susceptibility variants to common diseases and are providing new insights into the genetic structure and relationships of human populations. Moreover, genome-wide data can be used to search for signals of recent positive selection, thereby providing new insights into the genetic adaptations that occurred as modern humans spread out of Africa and around the world.

Methodology

We genotyped approximately 500,000 SNPs in 255 individuals (5 individuals from each of 51 worldwide populations) from the Human Genome Diversity Panel (HGDP-CEPH). When merged with non-overlapping SNPs typed previously in 250 of these same individuals, the resulting data consist of over 950,000 SNPs. We then analyzed the genetic relationships and ancestry of individuals without assigning them to populations, and we also identified candidate regions of recent positive selection at both the population and regional (continental) level.

Conclusions

Our analyses both confirm and extend previous studies; in particular, we highlight the impact of various dispersals, and the role of substructure in Africa, on human genetic diversity. We also identified several novel candidate regions for recent positive selection, and a gene ontology (GO) analysis identified several GO groups that were significantly enriched for such candidate genes, including immunity and defense related genes, sensory perception genes, membrane proteins, signal receptors, lipid binding/metabolism genes, and genes involved in the nervous system. Among the novel candidate genes identified are two genes involved in the thyroid hormone pathway that show signals of selection in African Pygmies that may be related to their short stature.     

Spatial Genetic Heterogeneity in Populations of a Newly Invasive Whitefly in China Revealed by a Nation-Wide Field Survey

Background

Even though introductions of exotic species provide ready-made experiments of rapid evolution, few studies have examined the genetic structure of an exotic species shortly after its initial introduction and subsequent spread. To determine the genetic structure of its populations during the initial introduction, we investigated the invasive sweet potato whitefly (Bemisia tabaci Q, commonly known as B. tabaci biotype Q) in China, which was introduced in approximately 2003. A total of 619 B. tabaci Q individuals in 20 provinces throughout China were collected and analyzed using five microsatellite loci.

Results

The introduced populations of B. tabaci Q in China represent eight genetic clusters with different geographic distributions. The populations in Yunnan Province, where B. tabaci Q was first detected, are genetically different from the other populations in China.

Conclusion

The introduced populations of B. tabaci Q in China have high spatial genetic heterogeneity. Additional research is required to determine whether the heterogeneity results from multiple introductions, rapid evolution following one or few introductions, or some combination of multiple introductions and rapid evolution. The heterogeneity, however, is inconsistent with a single introduction at Yunnan Province, where B. tabaci Q was first detected, followed by spread.     

A powerful tool for genome analysis in maize: development and evaluation of the high density 600 k SNP genotyping array

Background

High density genotyping data are indispensable for genomic analyses of complex traits in animal and crop species. Maize is one of the most important crop plants worldwide, however a high density SNP genotyping array for analysis of its large and highly dynamic genome was not available so far.

Results

We developed a high density maize SNP array composed of 616,201 variants (SNPs and small indels). Initially, 57 M variants were discovered by sequencing 30 representative temperate maize lines and then stringently filtered for sequence quality scores and predicted conversion performance on the array resulting in the selection of 1.2 M polymorphic variants assayed on two screening arrays. To identify high-confidence variants, 285 DNA samples from a broad genetic diversity panel of worldwide maize lines including the samples used for sequencing, important founder lines for European maize breeding, hybrids, and proprietary samples with European, US, semi-tropical, and tropical origin were used for experimental validation. We selected 616 k variants according to their performance during validation, support of genotype calls through sequencing data, and physical distribution for further analysis and for the design of the commercially available Affymetrix® Axiom® Maize Genotyping Array. This array is composed of 609,442 SNPs and 6,759 indels. Among these are 116,224 variants in coding regions and 45,655 SNPs of the Illumina® MaizeSNP50 BeadChip for study comparison. In a subset of 45,974 variants, apart from the target SNP additional off-target variants are detected, which show only a minor bias towards intermediate allele frequencies. We performed principal coordinate and admixture analyses to determine the ability of the array to detect and resolve population structure and investigated the extent of LD within a worldwide validation panel.

Conclusions

The high density Affymetrix® Axiom® Maize Genotyping Array is optimized for European and American temperate maize and was developed based on a diverse sample panel by applying stringent quality filter criteria to ensure its suitability for a broad range of applications. With 600 k variants it is the largest currently publically available genotyping array in crop species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-823) contains supplementary material, which is available to authorized users.     

Linkage disequilibrium and genome-wide association analysis for anthocyanin pigmentation and fruit color in eggplant

Background

The genome-wide association (GWA) approach represents an alternative to biparental linkage mapping for determining the genetic basis of trait variation. Both approaches rely on recombination to re-arrange the genome, and seek to establish correlations between phenotype and genotype. The major advantages of GWA lie in being able to sample a much wider range of the phenotypic and genotypic variation present, in being able to exploit multiple rounds of historical recombination in many different lineages and to include multiple accessions of direct relevance to crop improvement.

Results

A 191 accessions eggplant (Solanum melongena L.) association panel, comprising a mixture of breeding lines, old varieties and landrace selections originating from Asia and the Mediterranean Basin, was SNP genotyped and scored for anthocyanin pigmentation and fruit color at two locations over two years. The panel formed two major clusters, reflecting geographical provenance and fruit type. The global level of linkage disequilibrium was 3.4 cM. A mixed linear model appeared to be the most appropriate for GWA. A set of 56 SNP locus/phenotype associations was identified and the genomic regions harboring these loci were distributed over nine of the 12 eggplant chromosomes. The associations were compared with the location of known QTL for the same traits.

Conclusion

The GWA mapping approach was effective in validating a number of established QTL and, thanks to the wide diversity captured by the panel, was able to detect a series of novel marker/trait associations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-896) contains supplementary material, which is available to authorized users.     

Multiple introductions boosted genetic diversity in the invasive range of black cherry (Prunus serotina; Rosaceae)

Background and Aims

Black cherry (Prunus serotina) is a North American tree that is rapidly invading European forests. This species was introduced first as an ornamental plant then it was massively planted by foresters in many countries but its origins and the process of invasion remain poorly documented. Based on a genetic survey of both native and invasive ranges, the invasion history of black cherry was investigated by identifying putative source populations and then assessing the importance of multiple introductions on the maintenance of gene diversity.

Methods

Genetic variability and structure of 23 populations from the invasive range and 22 populations from the native range were analysed using eight nuclear microsatellite loci and five chloroplast DNA regions.

Key Results

Chloroplast DNA diversity suggests there were multiple introductions from a single geographic region (the north-eastern United States). A low reduction of genetic diversity was observed in the invasive range for both nuclear and plastid genomes. High propagule pressure including both the size and number of introductions shaped the genetic structure in Europe and boosted genetic diversity. Populations from Denmark, The Netherlands, Belgium and Germany showed high genetic diversity and low differentiation among populations, supporting the hypothesis that numerous introduction events, including multiple individuals and exchanges between sites, have taken place during two centuries of plantation.

Conclusions

This study postulates that the invasive black cherry has originated from east of the Appalachian Mountains (mainly the Allegheny plateau) and its invasiveness in north-western Europe is mainly due to multiple introductions containing high numbers of individuals.     

Genetic characterization of human T-cell lymphotropic virus type 1 in Mozambique: transcontinental lineages drive the HTLV-1 endemic

Background

Human T-Cell Lymphotropic Virus Type 1 (HTLV-1) is the etiological agent of adult T-cell leukemia (ATL) and HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP). It has been estimated that 10–20 million people are infected worldwide, but no successful treatment is available. Recently, the epidemiology of this virus was addressed in blood donors from Maputo, showing rates from 0.9 to 1.2%. However, the origin and impact of HTLV endemic in this population is unknown.

Objective

To assess the HTLV-1 molecular epidemiology in Mozambique and to investigate their relationship with HTLV-1 lineages circulating worldwide.

Methods

Blood donors and HIV patients were screened for HTLV antibodies by using enzyme immunoassay, followed by Western Blot. PCR and sequencing of HTLV-1 LTR region were applied and genetic HTLV-1 subtypes were assigned by the neighbor-joining method. The mean genetic distance of Mozambican HTLV-1 lineages among the genetic clusters were determined. Human mitochondrial (mt) DNA analysis was performed and individuals classified in mtDNA haplogroups.

Results

LTR HTLV-1 analysis demonstrated that all isolates belong to the Transcontinental subgroup of the Cosmopolitan subtype. Mozambican HTLV-1 sequences had a high inter-strain genetic distance, reflecting in three major clusters. One cluster is associated with the South Africa sequences, one is related with Middle East and India strains and the third is a specific Mozambican cluster. Interestingly, 83.3% of HIV/HTLV-1 co-infection was observed in the Mozambican cluster. The human mtDNA haplotypes revealed that all belong to the African macrohaplogroup L with frequencies representatives of the country.

Conclusions

The Mozambican HTLV-1 genetic diversity detected in this study reveals that although the strains belong to the most prevalent and worldwide distributed Transcontinental subgroup of the Cosmopolitan subtype, there is a high HTLV diversity that could be correlated with at least 3 different HTLV-1 introductions in the country. The significant rate of HTLV-1a/HIV-1C co-infection, particularly in the Mozambican cluster, has important implications for the controls programs of both viruses.     

The origin and evolutionary history of HIV-1 subtype C in Senegal

Background

The classification of HIV-1 strains in subtypes and Circulating Recombinant Forms (CRFs) has helped in tracking the course of the HIV pandemic. In Senegal, which is located at the tip of West Africa, CRF02_AG predominates in the general population and Female Sex Workers (FSWs). In contrast, 40% of Men having Sex with Men (MSM) in Senegal are infected with subtype C. In this study we analyzed the geographical origins and introduction dates of HIV-1 C in Senegal in order to better understand the evolutionary history of this subtype, which predominates today in the MSM population

Methodology/Principal Findings

We used a combination of phylogenetic analyses and a Bayesian coalescent-based approach, to study the phylogenetic relationships in pol of 56 subtype C isolates from Senegal with 3,025 subtype C strains that were sampled worldwide. Our analysis shows a significantly well supported cluster which contains all subtype C strains that circulate among MSM in Senegal. The MSM cluster and other strains from Senegal are widely dispersed among the different subclusters of African HIV-1 C strains, suggesting multiple introductions of subtype C in Senegal from many different southern and east African countries. More detailed analyses show that HIV-1 C strains from MSM are more closely related to those from southern Africa. The estimated date of the MRCA of subtype C in the MSM population in Senegal is estimated to be in the early 80''s.

Conclusions/Significance

Our evolutionary reconstructions suggest that multiple subtype C viruses with a common ancestor originating in the early 1970s entered Senegal. There was only one efficient spread in the MSM population, which most likely resulted from a single introduction, underlining the importance of high-risk behavior in spread of viruses.     

Single nucleotide polymorphisms generated by genotyping by sequencing to characterize genome-wide diversity,linkage disequilibrium,and selective sweeps in cultivated watermelon

Background

A large single nucleotide polymorphism (SNP) dataset was used to analyze genome-wide diversity in a diverse collection of watermelon cultivars representing globally cultivated, watermelon genetic diversity. The marker density required for conducting successful association mapping depends on the extent of linkage disequilibrium (LD) within a population. Use of genotyping by sequencing reveals large numbers of SNPs that in turn generate opportunities in genome-wide association mapping and marker-assisted selection, even in crops such as watermelon for which few genomic resources are available. In this paper, we used genome-wide genetic diversity to study LD, selective sweeps, and pairwise F_ST distributions among worldwide cultivated watermelons to track signals of domestication.

Results

We examined 183 Citrullus lanatus var. lanatus accessions representing domesticated watermelon and generated a set of 11,485 SNP markers using genotyping by sequencing. With a diverse panel of worldwide cultivated watermelons, we identified a set of 5,254 SNPs with a minor allele frequency of ≥ 0.05, distributed across the genome. All ancestries were traced to Africa and an admixture of various ancestries constituted secondary gene pools across various continents. A sliding window analysis using pairwise F_ST values was used to resolve selective sweeps. We identified strong selection on chromosomes 3 and 9 that might have contributed to the domestication process. Pairwise analysis of adjacent SNPs within a chromosome as well as within a haplotype allowed us to estimate genome-wide LD decay. LD was also detected within individual genes on various chromosomes. Principal component and ancestry analyses were used to account for population structure in a genome-wide association study. We further mapped important genes for soluble solid content using a mixed linear model.

Conclusions

Information concerning the SNP resources, population structure, and LD developed in this study will help in identifying agronomically important candidate genes from the genomic regions underlying selection and for mapping quantitative trait loci using a genome-wide association study in sweet watermelon.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-767) contains supplementary material, which is available to authorized users.     

Genetic structure of the tiger mosquito, Aedes albopictus, in Cameroon (Central Africa)

Background

Aedes albopictus (Skuse, 1884) (Diptera: Culicidae), a mosquito native to Asia, has recently invaded all five continents. In Central Africa it was first reported in the early 2000s, and has since been implicated in the emergence of arboviruses such as dengue and chikungunya in this region. Recent genetic studies of invasive species have shown that multiple introductions are a key factor for successful expansion in new areas. As a result, phenotypic characters such as vector competence and insecticide susceptibility may vary within invasive pest species, potentially affecting vector efficiency and pest management. Here we assessed the genetic variability and population genetics of Ae. albopictus isolates in Cameroon (Central Africa), thereby deducing their likely geographic origin.

Methods and Results

Mosquitoes were sampled in 2007 in 12 localities in southern Cameroon and analyzed for polymorphism at six microsatellite loci and in two mitochondrial DNA regions (ND5 and COI). All the microsatellite markers were successfully amplified and were polymorphic, showing moderate genetic structureamong geographic populations (F_ST = 0.068, P<0.0001). Analysis of mtDNA sequences revealed four haplotypes each for the COI and ND5 genes, with a dominant haplotype shared by all Cameroonian samples. The weak genetic variation estimated from the mtDNA genes is consistent with the recent arrival of Ae. albopictus in Cameroon. Phylogeographic analysis based on COI polymorphism indicated that Ae. albopictus populations from Cameroon are related to tropical rather than temperate or subtropical outgroups.

Conclusion

The moderate genetic diversity observed among Cameroonian Ae. albopictus isolates is in keeping with recent introduction and spread in this country. The genetic structure of natural populations points to multiple introductions from tropical regions.     

Unravelling the complex trait of harvest index in rapeseed (Brassica napus L.) with association mapping

Background

Harvest index (HI), the ratio of grain yield to total biomass, is considered as a measure of biological success in partitioning assimilated photosynthate to the harvestable product. While crop production can be dramatically improved by increasing HI, the underlying molecular genetic mechanism of HI in rapeseed remains to be shown.

Results

In this study, we examined the genetic architecture of HI using 35,791 high-throughput single nucleotide polymorphisms (SNPs) genotyped by the Illumina BrassicaSNP60 Bead Chip in an association panel with 155 accessions. Five traits including plant height (PH), branch number (BN), biomass yield per plant (BY), harvest index (HI) and seed yield per plant (SY), were phenotyped in four environments. HI was found to be strongly positively correlated with SY, but negatively or not strongly correlated with PH. Model comparisons revealed that the A–D test (ADGWAS model) could perfectly balance false positives and statistical power for HI and associated traits. A total of nine SNPs on the C genome were identified to be significantly associated with HI, and five of them were identified to be simultaneously associated with HI and SY. These nine SNPs explained 3.42 % of the phenotypic variance in HI.

Conclusions

Our results showed that HI is a complex polygenic phenomenon that is strongly influenced by both environmental and genotype factors. The implications of these results are that HI can be increased by decreasing PH or reducing inefficient transport from pods to seeds in rapeseed. The results from this association mapping study can contribute to a better understanding of natural variations of HI, and facilitate marker-based breeding for HI.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1607-0) contains supplementary material, which is available to authorized users.     

Why Africa matters: evolution of Old World Salvia (Lamiaceae) in Africa

Background and Aims

Salvia is the largest genus in Lamiaceae and it has recently been found to be non-monophyletic. Molecular data on Old World Salvia are largely lacking. In this study, we present data concerning Salvia in Africa. The focus is on the colonization of the continent, character evolution and the switch of pollination systems in the genus.

Methods

Maximum likelihood and Bayesian inference were used for phylogenetic reconstruction. Analyses were based on two nuclear markers [internal transcribed spacer (ITS) and external transcribed spacer (ETS)] and one plastid marker (rpl32-trnL). Sequence data were generated for 41 of the 62 African taxa (66 %). Mesquite was used to reconstruct ancestral character states for distribution, life form, calyx shape, stamen type and pollination syndrome.

Key Results

Salvia in Africa is non-monophyletic. Each of the five major regions in Africa, except Madagascar, was colonized at least twice, and floristic links between North African, south-west Asian and European species are strongly supported. The large radiation in Sub-Saharan Africa (23 species) can be traced back to dispersal from North Africa via East Africa to the Cape Region. Adaptation to bird pollination in southern Africa and Madagascar reflects parallel evolution.

Conclusions

The phenotypic diversity in African Salvia is associated with repeated introductions to the continent. Many important evolutionary processes, such as colonization, adaptation, parallelism and character transformation, are reflected in this comparatively small group. The data presented in this study can help to understand the evolution of Salvia sensu lato and other large genera.     

Linked genetic variants on chromosome 10 control ear morphology and body mass among dog breeds

Background

The domestic dog is a rich resource for mapping the genetic components of phenotypic variation due to its unique population history involving strong artificial selection. Genome-wide association studies have revealed a number of chromosomal regions where genetic variation associates with morphological characters that typify dog breeds. A region on chromosome 10 is among those with the highest levels of genetic differentiation between dog breeds and is associated with body mass and ear morphology, a common motif of animal domestication. We characterised variation in this region to uncover haplotype structure and identify candidate functional variants.

Results

We first identified SNPs that strongly associate with body mass and ear type by comparing sequence variation in a 3 Mb region between 19 breeds with a variety of phenotypes. We next genotyped a subset of 123 candidate SNPs in 288 samples from 46 breeds to identify the variants most highly associated with phenotype and infer haplotype structure. A cluster of SNPs that associate strongly with the drop ear phenotype is located within a narrow interval downstream of the gene MSRB3, which is involved in human hearing. These SNPs are in strong genetic linkage with another set of variants that correlate with body mass within the gene HMGA2, which affects human height. In addition we find evidence that this region has been under selection during dog domestication, and identify a cluster of SNPs within MSRB3 that are highly differentiated between dogs and wolves.

Conclusions

We characterise genetically linked variants that potentially influence ear type and body mass in dog breeds, both key traits that have been modified by selective breeding that may also be important for domestication. The finding that variants on long haplotypes have effects on more than one trait suggests that genetic linkage can be an important determinant of the phenotypic response to selection in domestic animals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1702-2) contains supplementary material, which is available to authorized users.     

Cluster analysis in severe emphysema subjects using phenotype and genotype data: an exploratory investigation

Background

Numerous studies have demonstrated associations between genetic markers and COPD, but results have been inconsistent. One reason may be heterogeneity in disease definition. Unsupervised learning approaches may assist in understanding disease heterogeneity.

Methods

We selected 31 phenotypic variables and 12 SNPs from five candidate genes in 308 subjects in the National Emphysema Treatment Trial (NETT) Genetics Ancillary Study cohort. We used factor analysis to select a subset of phenotypic variables, and then used cluster analysis to identify subtypes of severe emphysema. We examined the phenotypic and genotypic characteristics of each cluster.

Results

We identified six factors accounting for 75% of the shared variability among our initial phenotypic variables. We selected four phenotypic variables from these factors for cluster analysis: 1) post-bronchodilator FEV₁ percent predicted, 2) percent bronchodilator responsiveness, and quantitative CT measurements of 3) apical emphysema and 4) airway wall thickness. K-means cluster analysis revealed four clusters, though separation between clusters was modest: 1) emphysema predominant, 2) bronchodilator responsive, with higher FEV₁; 3) discordant, with a lower FEV₁ despite less severe emphysema and lower airway wall thickness, and 4) airway predominant. Of the genotypes examined, membership in cluster 1 (emphysema-predominant) was associated with TGFB1 SNP rs1800470.

Conclusions

Cluster analysis may identify meaningful disease subtypes and/or groups of related phenotypic variables even in a highly selected group of severe emphysema subjects, and may be useful for genetic association studies.