Development of a high performance parallel computing platform and its use in the study of nanostructures : clusters, sheets and tubes .

Small clusters of gallium oxide, technologically important high temperature ceramic, together with interaction of nucleic acid bases with graphene and small-diameter carbon nanotube are focus of first principles calculations in this work. A high performance parallel computing platform is also developed to perform these calculations at Michigan Tech. First principles calculations are based on density functional theory employing either local density or gradient-corrected approximation together with plane wave and gaussian basis sets. The bulk Ga2O3 is known to be a very good candidate for fabricating electronic devices that operate at high temperatures. To explore the properties of Ga2O3 at nonoscale, we have performed a systematic theoretical study on the small polyatomic gallium oxide clusters. The calculated results find that all lowest energy isomers of GamOn clusters are dominated by the Ga-O bonds over the metal-metal or the oxygen-oxygen bonds. Analysis of atomic charges suggest the clusters to be highly ionic similar to the case of bulk Ga2O3. In the study of sequential oxidation of these slusters starting from Ga2O, it is found that the most stable isomers display up to four different backbones of constituent atoms. Furthermore, the predicted configuration of the ground state of Ga2O is recently confirmed by the experimental result of Neumark's group. Guided by the results of calculations the study of gallium oxide clusters, performance related challenge of computational simulations, of producing high performance computers/platforms, has been addressed. Several engineering aspects were thoroughly studied during the design, development and implementation of the high performance parallel computing platform, rama, at Michigan Tech. In an attempt to stay true to the principles of Beowulf revolutioni, the rama cluster was extensively customized to make it easy to understand, and use - for administrators as well as end-users. Following the results of benchmark calculations and to keep up with the complexity of systems under study, rama has been expanded to a total of sixty four processors. Interest in the non-covalent intereaction of DNA with carbon nanotubes has steadily increased during past several years. This hybrid system, at the junction of the biological regime and the nanomaterials world, possesses features which make it very attractive for a wide range of applicatioins. Using the in-house computational power available, we have studied details of the interaction between nucleic acid bases with graphene sheet as well as high-curvature small-diameter carbon nanotube. The calculated trend in the binding energies strongly suggests that the polarizability of the base molecules determines the interaction strength of the nucleic acid bases with graphene. When comparing the results obtained here for physisorption on the small diameter nanotube considered with those from the study on graphene, it is observed that the interaction strength of nucleic acid bases is smaller for the tube. Thus, these results show that the effect of introducing curvature is to reduce the binding energy. The binding energies for the two extreme cases of negligible curvature (i.e. flat graphene sheet) and of very high curvature (i.e. small diameter nanotube) may be considered as upper and lower bounds. This finding represents an important step towards a better understanding of experimentally observed sequence-dependent interaction of DNA with Carbon nanotubes.

ANALYSIS AND VISUALIZATION OF FLOW FIELDS USING INFORMATION-THEORETIC TECHNIQUES AND GRAPH-BASED REPRESENTATIONS

Three-dimensional flow visualization plays an essential role in many areas of science and engineering, such as aero- and hydro-dynamical systems which dominate various physical and natural phenomena. For popular methods such as the streamline visualization to be effective, they should capture the underlying flow features while facilitating user observation and understanding of the flow field in a clear manner. My research mainly focuses on the analysis and visualization of flow fields using various techniques, e.g. information-theoretic techniques and graph-based representations. Since the streamline visualization is a popular technique in flow field visualization, how to select good streamlines to capture flow patterns and how to pick good viewpoints to observe flow fields become critical. We treat streamline selection and viewpoint selection as symmetric problems and solve them simultaneously using the dual information channel [81]. To the best of my knowledge, this is the first attempt in flow visualization to combine these two selection problems in a unified approach. This work selects streamline in a view-independent manner and the selected streamlines will not change for all viewpoints. My another work [56] uses an information-theoretic approach to evaluate the importance of each streamline under various sample viewpoints and presents a solution for view-dependent streamline selection that guarantees coherent streamline update when the view changes gradually. When projecting 3D streamlines to 2D images for viewing, occlusion and clutter become inevitable. To address this challenge, we design FlowGraph [57, 58], a novel compound graph representation that organizes field line clusters and spatiotemporal regions hierarchically for occlusion-free and controllable visual exploration. We enable observation and exploration of the relationships among field line clusters, spatiotemporal regions and their interconnection in the transformed space. Most viewpoint selection methods only consider the external viewpoints outside of the flow field. This will not convey a clear observation when the flow field is clutter on the boundary side. Therefore, we propose a new way to explore flow fields by selecting several internal viewpoints around the flow features inside of the flow field and then generating a B-Spline curve path traversing these viewpoints to provide users with closeup views of the flow field for detailed observation of hidden or occluded internal flow features [54]. This work is also extended to deal with unsteady flow fields. Besides flow field visualization, some other topics relevant to visualization also attract my attention. In iGraph [31], we leverage a distributed system along with a tiled display wall to provide users with high-resolution visual analytics of big image and text collections in real time. Developing pedagogical visualization tools forms my other research focus. Since most cryptography algorithms use sophisticated mathematics, it is difficult for beginners to understand both what the algorithm does and how the algorithm does that. Therefore, we develop a set of visualization tools to provide users with an intuitive way to learn and understand these algorithms.

Phylogenetic relationship of equine Actinobacillus species and distribution of RTX toxin genes among clusters

Equine Actinobacillus species were analysed phylogenetically by 16S rRNA gene (rrs) sequencing focusing on the species Actinobacillus equuli, which has recently been subdivided into the non-haemolytic A. equuli subsp. equuli and the haemolytic A. equuli subsp. haemolyticus. In parallel we determined the profile for RTX toxin genes of the sample of strains by PCR testing for the presence of the A. equuli haemolysin gene aqx, and the toxin genes apxI, apxII, apxIII and apxIV, which are known in porcine pathogens such as Actinobacillus pleuropneumoniae and Actinobacillus suis. The rrs-based phylogenetic analysis revealed two distinct subclusters containing both A. equuli subsp. equuli and A. equuli subsp. haemolyticus distributed through both subclusters with no correlation to taxonomic classification. Within one of the rrs-based subclusters containing the A. equuli subsp. equuli type strain, clustered as well the porcine Actinobacillus suis strains. This latter is known to be also phenotypically closely related to A. equuli. The toxin gene analysis revealed that all A. equuli subsp. haemolyticus strains from both rrs subclusters specifically contained the aqx gene while the A. suis strains harboured the genes apxI and apxII. The aqx gene was found to be specific for A. equuli subsp. haemolyticus, since A. equuli subsp. equuli contained no aqx nor any of the other RTX genes tested. The specificity of aqx for the haemolytic equine A. equuli and ApxI and ApxII for the porcine A. suis indicates a role of these RTX toxins in host species predilection of the two closely related species of bacterial pathogens and allows PCR based diagnostic differentiation of the two.

Cluster randomized clinical trials in orthodontics: design, analysis and reporting issues

Cluster randomized trials (CRTs) use as the unit of randomization clusters, which are usually defined as a collection of individuals sharing some common characteristics. Common examples of clusters include entire dental practices, hospitals, schools, school classes, villages, and towns. Additionally, several measurements (repeated measurements) taken on the same individual at different time points are also considered to be clusters. In dentistry, CRTs are applicable as patients may be treated as clusters containing several individual teeth. CRTs require certain methodological procedures during sample calculation, randomization, data analysis, and reporting, which are often ignored in dental research publications. In general, due to similarity of the observations within clusters, each individual within a cluster provides less information compared with an individual in a non-clustered trial. Therefore, clustered designs require larger sample sizes compared with non-clustered randomized designs, and special statistical analyses that account for the fact that observations within clusters are correlated. It is the purpose of this article to highlight with relevant examples the important methodological characteristics of cluster randomized designs as they may be applied in orthodontics and to explain the problems that may arise if clustered observations are erroneously treated and analysed as independent (non-clustered).

Analysis of close associations of uropod-associated proteins in human T-cells using the proximity ligation assay

We have shown previously that the raft-associated proteins flotillin-1 and -2 are rapidly recruited to the uropods of chemoattractant-stimulated human neutrophils and T-cells and are involved in cell polarization. Other proteins such as the adhesion receptor PSGL-1, the actin-membrane linker proteins ezrin/radixin/moesin (ERM) and the signaling enzyme phosphatidylinositol-4-phosphate 5-kinase type Iγ90 (PIPKIγ90) also accumulate in the T-cell uropod. Using the in situ proximity ligation assay (PLA) we now have investigated putative close associations of these proteins in human freshly isolated T-cells before and after chemokine addition. The PLA allows in situ subcellular localization of close proximity of endogenous proteins at single-molecule resolution in fixed cells. It allows detection also of weaker and transient complexes that would not be revealed with co-immunoprecipitation approaches. We previously provided evidence for heterodimer formation of tagged flotillin-1 and -2 in T-cells before and after chemokine addition using fluorescence resonance energy transfer (FRET). We now confirm these findings using PLA for the endogenous flotillins in fixed human T-cells. Moreover, in agreement with the literature, our PLA findings confirm a close association of endogenous PSGL-1 and ERM proteins both in resting and chemokine-activated human T-cells. In addition, we provide novel evidence using the PLA for close associations of endogenous activated ERM proteins with PIPKIγ90 and of endogenous flotillins with PSGL-1 in human T-cells, before and after chemokine addition. Our findings suggest that preformed clusters of these proteins coalesce in the uropod upon cell stimulation.

Evaluating the disparity of female breast cancer mortality among racial groups - a spatiotemporal analysis

Background The literature suggests that the distribution of female breast cancer mortality demonstrates spatial concentration. There remains a lack of studies on how the mortality burden may impact racial groups across space and over time. The present study evaluated the geographic variations in breast cancer mortality in Texas females according to three predominant racial groups (non-Hispanic White, Black, and Hispanic females) over a twelve-year period. It sought to clarify whether the spatiotemporal trend might place an uneven burden on particular racial groups, and whether the excess trend has persisted into the current decade. Methods The Spatial Scan Statistic was employed to examine the geographic excess of breast cancer mortality by race in Texas counties between 1990 and 2001. The statistic was conducted with a scan window of a maximum of 90% of the study period and a spatial cluster size of 50% of the population at risk. The next scan was conducted with a purely spatial option to verify whether the excess mortality persisted further. Spatial queries were performed to locate the regions of excess mortality affecting multiple racial groups. Results The first scan identified 4 regions with breast cancer mortality excess in both non-Hispanic White and Hispanic female populations. The most likely excess mortality with a relative risk of 1.12 (p = 0.001) occurred between 1990 and 1996 for non-Hispanic Whites, including 42 Texas counties along Gulf Coast and Central Texas. For Hispanics, West Texas with a relative risk of 1.18 was the most probable region of excess mortality (p = 0.001). Results of the second scan were identical to the first. This suggested that the excess mortality might not persist to the present decade. Spatial queries found that 3 counties in Southeast and 9 counties in Central Texas had excess mortality involving multiple racial groups. Conclusion Spatiotemporal variations in breast cancer mortality affected racial groups at varying levels. There was neither evidence of hot-spot clusters nor persistent spatiotemporal trends of excess mortality into the present decade. Non-Hispanic Whites in the Gulf Coast and Hispanics in West Texas carried the highest burden of mortality, as evidenced by spatial concentration and temporal persistence.

Identification and phenotypic characterization of a second collagen adhesin, Scm, and genome-based identification and analysis of 13 other predicted MSCRAMMs, including four distinct pilus loci, in Enterococcus faecium.

Attention has recently been drawn to Enterococcus faecium because of an increasing number of nosocomial infections caused by this species and its resistance to multiple antibacterial agents. However, relatively little is known about the pathogenic determinants of this organism. We have previously identified a cell-wall-anchored collagen adhesin, Acm, produced by some isolates of E. faecium, and a secreted antigen, SagA, exhibiting broad-spectrum binding to extracellular matrix proteins. Here, we analysed the draft genome of strain TX0016 for potential microbial surface components recognizing adhesive matrix molecules (MSCRAMMs). Genome-based bioinformatics identified 22 predicted cell-wall-anchored E. faecium surface proteins (Fms), of which 15 (including Acm) had characteristics typical of MSCRAMMs, including predicted folding into a modular architecture with multiple immunoglobulin-like domains. Functional characterization of one [Fms10; redesignated second collagen adhesin of E. faecium (Scm)] revealed that recombinant Scm(65) (A- and B-domains) and Scm(36) (A-domain) bound to collagen type V efficiently in a concentration-dependent manner, bound considerably less to collagen type I and fibrinogen, and differed from Acm in their binding specificities to collagen types IV and V. Results from far-UV circular dichroism measurements of recombinant Scm(36) and of Acm(37) indicated that these proteins were rich in beta-sheets, supporting our folding predictions. Whole-cell ELISA and FACS analyses unambiguously demonstrated surface expression of Scm in most E. faecium isolates. Strikingly, 11 of the 15 predicted MSCRAMMs clustered in four loci, each with a class C sortase gene; nine of these showed similarity to Enterococcus faecalis Ebp pilus subunits and also contained motifs essential for pilus assembly. Antibodies against one of the predicted major pilus proteins, Fms9 (redesignated EbpC(fm)), detected a 'ladder' pattern of high-molecular-mass protein bands in a Western blot analysis of cell surface extracts from E. faecium, suggesting that EbpC(fm) is polymerized into a pilus structure. Further analysis of the transcripts of the corresponding gene cluster indicated that fms1 (ebpA(fm)), fms5 (ebpB(fm)) and ebpC(fm) are co-transcribed, a result consistent with those for pilus-encoding gene clusters of other Gram-positive bacteria. All 15 genes occurred frequently in 30 clinically derived diverse E. faecium isolates tested. The common occurrence of MSCRAMM- and pilus-encoding genes and the presence of a second collagen-binding protein may have important implications for our understanding of this emerging pathogen.

"When we like what we know – A parametric fMRI analysis of beauty and familiarity"

This paper presents a neuroscientific study of aesthetic judgments on written texts. In an fMRI experiment participants read a number of proverbs without explicitly evaluating them. In a post-scan rating they rated each item for familiarity and beauty. These individual ratings were correlated with the functional data to investigate the neural correlates of implicit aesthetic judgments. We identified clusters in which BOLD activity was correlated with individual post-scan beauty ratings. This indicates that some spontaneous aesthetic evaluation takes place during reading, even if not required by the task. Positive correlations were found in the ventral striatum and in medial prefrontal cortex, likely reflecting the rewarding nature of sentences that are aesthetically pleasing. On the contrary, negative correlations were observed in the classic left frontotemporal reading network. Midline structures and bilateral temporo-parietal regions correlated positively with familiarity, suggesting a shift from the task-network towards the default network with increasing familiarity.

A cladistic analysis of apolipoprotein B gene variation and plasma lipid and apolipoprotein B levels, using familial data

Any functionally important mutation is embedded in an evolutionary matrix of other mutations. Cladistic analysis, based on this, is a method of investigating gene effects using a haplotype phylogeny to define a set of tests which localize causal mutations to branches of the phylogeny. Previous implementations of cladistic analysis have not addressed the issue of analyzing data from related individuals, though in human studies, family data are usually needed to obtain unambiguous haplotypes. In this study, a method of cladistic analysis is described in which haplotype effects are parameterized in a linear model which accounts for familial correlations. The method was used to study the effect of apolipoprotein (Apo) B gene variation on total-, LDL-, and HDL-cholesterol, triglyceride, and Apo B levels in 121 French families. Five polymorphisms defined Apo B haplotypes: the signal peptide Insertion/deletion, Bsp 1286I, XbaI, MspI, and EcoRI. Eleven haplotypes were found, and a haplotype phylogeny was constructed and used to define a set of tests of haplotype effects on lipid and apo B levels.^ This new method of cladistic analysis, the parametric method, found significant effects for single haplotypes for all variables. For HDL-cholesterol, 3 clusters of evolutionarily-related haplotypes affecting levels were found. Haplotype effects accounted for about 10% of the genetic variance of triglyceride and HDL-cholesterol levels. The results of the parametric method were compared to those of a method of cladistic analysis based on permutational testing. The permutational method detected fewer haplotype effects, even when modified to account for correlations within families. Simulation studies exploring these differences found evidence of systematic errors in the permutational method due to the process by which haplotype groups were selected for testing.^ The applicability of cladistic analysis to human data was shown. The parametric method is suggested as an improvement over the permutational method. This study has identified candidate haplotypes for sequence comparisons in order to locate the functional mutations in the Apo B gene which may influence plasma lipid levels. ^

{\it Hox\/} genes and specification of regional identity in the axial skeleton: Analysis of the individual and combined mutant phenotypes of the paralogous group 4 murine {\it Hox\/} genes; {\it hoxa\/}-4, {\it hoxb\/}-4 and {\it hoxd\/}-4

The Hox gene products are transcription factors involved in specifying regional identity along the anteroposterior body axis. In Drosophila, where these genes are known as HOM-C (Homeotic-complex) genes and where they have been most extensively studied, they are expressed in restricted domains along the anteroposterior axis with different anterior limits. Genetic analysis of a large number of gain- and loss-of-function alleles of these genes has revealed that these genes are important in specifying segmental identity at their anterior limits of expression. Furthermore, there is a functional dominance of posterior genes over anterior genes, such that posterior genes can dominantly specify their developmental programs in spite of the expression of more anterior genes in the same segment. In the mouse, there are four clusters of HOM-C genes, called Hox genes. Thus, there may be up to four genes, called paralogs, that are more highly homologous to each other and to their Drosophila homolog than they are to the other mouse Hox genes. The single mutants for two paralogous genes, hoxa-4 and hoxd-4, presented in this dissertation, are similar to several other mouse Hox mutants in that they show partial, incompletely penetrant homeotic transformations of vertebrae at their anterior limit of expression. These mutants were then bred with hoxb-4 mutants (Ramirez-Solis, et al. 1993) to generate the three possible double mutant combinations as well as the triple mutant. The skeletal phenotypes of these group 4 Hox compound mutants displayed clear alterations in regional identity, such that a nearly complete transformation towards the morphology of the first cervical vertebra occurs. These results suggest a certain degree of functional redundancy among paralogous genes in specifying regional identity. Furthermore, there was a remarkable dose-dependent increase in the number of vertebrae transformed to a first cervical vertebra identity, including the second through the fifth cervical vertebrae in the triple mutant. Thus, these genes are required in a larger anteroposterior domain than is revealed by the single mutant phenotypes alone, such that multiple mutations in these genes result in transformations of vertebrae that are not at their anterior limit of expression. ^

Design and analysis for three level hierarchical structures with count response

In numerous intervention studies and education field trials, random assignment to treatment occurs in clusters rather than at the level of observation. This departure of random assignment of units may be due to logistics, political feasibility, or ecological validity. Data within the same cluster or grouping are often correlated. Application of traditional regression techniques, which assume independence between observations, to clustered data produce consistent parameter estimates. However such estimators are often inefficient as compared to methods which incorporate the clustered nature of the data into the estimation procedure (Neuhaus 1993).1 Multilevel models, also known as random effects or random components models, can be used to account for the clustering of data by estimating higher level, or group, as well as lower level, or individual variation. Designing a study, in which the unit of observation is nested within higher level groupings, requires the determination of sample sizes at each level. This study investigates the design and analysis of various sampling strategies for a 3-level repeated measures design on the parameter estimates when the outcome variable of interest follows a Poisson distribution. ^ Results study suggest that second order PQL estimation produces the least biased estimates in the 3-level multilevel Poisson model followed by first order PQL and then second and first order MQL. The MQL estimates of both fixed and random parameters are generally satisfactory when the level 2 and level 3 variation is less than 0.10. However, as the higher level error variance increases, the MQL estimates become increasingly biased. If convergence of the estimation algorithm is not obtained by PQL procedure and higher level error variance is large, the estimates may be significantly biased. In this case bias correction techniques such as bootstrapping should be considered as an alternative procedure. For larger sample sizes, those structures with 20 or more units sampled at levels with normally distributed random errors produced more stable estimates with less sampling variance than structures with an increased number of level 1 units. For small sample sizes, sampling fewer units at the level with Poisson variation produces less sampling variation, however this criterion is no longer important when sample sizes are large. ^ 1Neuhaus J (1993). “Estimation efficiency and Tests of Covariate Effects with Clustered Binary Data”. Biometrics , 49, 989–996^

Transcriptional analysis of left-sided colitis, pancolitis, and ulcerative colitis-associated dysplasia.

BACKGROUND It is unknown why patients with extensive ulcerative colitis (UC) have a higher risk of colorectal cancer compared with patients with left-sided UC. This study characterizes the inflammatory processes in left-sided UC, pancolitis, and UC-associated dysplasia at the transcriptional level to identify potential biomarkers and transcripts of importance for the carcinogenic behavior of chronic inflammation. METHODS The Affymetrix GeneChip Human Genome U133 Plus 2.0 was applied on colonic biopsies from UC patients with left-sided UC, pancolitis, dysplasia, and controls. Reverse transcription polymerase chain reaction and immunohistochemistry were performed for validating selected transcripts in the initial cohort and in 2 independent cohorts of patients with UC. Microarray data were analyzed by principal component analysis, and reverse transcription polymerase chain reaction and immunohistochemistry data by the Wilcoxon's rank-sum test. RESULTS The principal component analysis results revealed separate clusters for left-sided UC, pancolitis, dysplasia, and controls. Close clustering of dysplastic and pancolitic samples indicated similarities in gene expression. Indeed, 101 and 656 parallel upregulated and downregulated transcripts, respectively, were identified in specimens from dysplasia and pancolitis. Validation of selected transcripts hereof identified insulin receptor alpha (INSRA) and MAP kinase interacting serine/threonine kinase 2 (MKNK2) with an enhanced expression in dysplasia compared with left-sided UC and controls, whereas laminin γ2 (LAMC2) was found with a lower expression in dysplasia compared with the remaining 3 groups. CONCLUSIONS This study demonstrates pancolitis and left-sided UC as distinct inflammatory processes at the transcriptional level, and identifies INSRA, MKNK2, and LAMC2 as potential critical transcripts in the inflammation-driven preneoplastic process of UC.

Bayesian models for the analysis of genetic structure when populations are correlated

Motivation: Population allele frequencies are correlated when populations have a shared history or when they exchange genes. Unfortunately, most models for allele frequency and inference about population structure ignore this correlation. Recent analytical results show that among populations, correlations can be very high, which could affect estimates of population genetic structure. In this study, we propose a mixture beta model to characterize the allele frequency distribution among populations. This formulation incorporates the correlation among populations as well as extending the model to data with different clusters of populations. Results: Using simulated data, we show that in general, the mixture model provides a good approximation of the among-population allele frequency distribution and a good estimate of correlation among populations. Results from fitting the mixture model to a dataset of genotypes at 377 autosomal microsatellite loci from human populations indicate high correlation among populations, which may not be appropriate to neglect. Traditional measures of population structure tend to over-estimate the amount of genetic differentiation when correlation is neglected. Inference is performed in a Bayesian framework.

AN ANALYSIS OF THE GENERAL WELL-BEING OF BLACKS AND WHITES

The importance of race as a factor in mental health status has been a topic of controversy. This study reviews the history of research in this area and examines racial variances in the relationship between selected socio-demographic variables and general well-being. The study also examines the appropriateness of an additive versus an interactive statistical model for this investigation.^ The sample consists of 6,913 persons who completed the General Well-Being Schedules as administered in the detailed component of the first National Health and Nutrition Examination Survey (NHANES I) conducted by the National Center for Health Statistics between April, 1971 and October, 1975. The sampling design is a multistage, probability sample of clusters of persons in area based segments. Of the 6,913 persons, 873 are Black.^ Unlike other recent community based mental health studies, this study revealed significant differences between the general well-being of Blacks and Whites. Blacks continued to exhibit significantly lower levels of well-being even after adjustments were made for income, education, marital status, sex, age and place of residence. Statistical interaction was found between race and sex with Black females reporting lower levels of well-being than either Black or White males or their White female counterparts.^ The study includes a detailed review of the NHANES I sample design. It is shown that selected aspects of the design make it difficult to render appropriate national comparisons of Black-White differences. As a result conclusions pertaining to these differences based on NHANES I may be of questionable validity. ^

Spatial analysis of pertussis cases in Harris County from 2005-2010

Objective: This study examined the recent trends and characteristics of reported pertussis in Harris County from 2005-2010. ^ Methods: The study population included surveillance data from all reported pertussis cases from January 1, 2005 to December 31, 2010 to Harris County Public Health and Environmental Services (HCPHES). We calculated incidence and attack rates for varying age groups, race/ethnicity, and gender. Spatial analyses were conducted of hot spot and cluster of incident cases in Harris County census tracts. Maps were constructed using geographic information system. ^ Results: Age-specific incidence rates of reported cases of pertussis were highest among infants under a year of age and lowest among adults age 20 and older. Hispanics represented the most cases reported compared to any other race or ethnic group (42% of 483 cases). Age-adjusted rates were highest in 2009 at 9.81 cases per 100,000 population. Only 31.2% of people received at least four of the recommended five doses of vaccine. Spatial analyses revealed statistically significant clusters within the northeast region of Harris County. ^ Conclusions: Hispanic infants are the most at risk group for pertussis. Although 70% of cases had a history of immunization, 41.8% of infants were appropriately vaccinated for their age. Increased vaccination coverage may decrease the incidence of pertussis.^