A novel harmony search-K means hybrid algorithm for clustering gene expression data

Recent progress in bioinformatics research has led to the accumulation of huge quantities of biological data at various data sources. The DNA microarray technology makes it possible to simultaneously analyze large number of genes across different samples. Clustering of microarray data can reveal the hidden gene expression patterns from large quantities of expression data that in turn offers tremendous possibilities in functional genomics, comparative genomics, disease diagnosis and drug development. The k- ¬means clustering algorithm is widely used for many practical applications. But the original k-¬means algorithm has several drawbacks. It is computationally expensive and generates locally optimal solutions based on the random choice of the initial centroids. Several methods have been proposed in the literature for improving the performance of the k-¬means algorithm. A meta-heuristic optimization algorithm named harmony search helps find out near-global optimal solutions by searching the entire solution space. Low clustering accuracy of the existing algorithms limits their use in many crucial applications of life sciences. In this paper we propose a novel Harmony Search-K means Hybrid (HSKH) algorithm for clustering the gene expression data. Experimental results show that the proposed algorithm produces clusters with better accuracy in comparison with the existing algorithms.     

Genome-wide linkage and peak-wide association study of obesity-related quantitative traits in Caribbean Hispanics

Although obesity is more prevalent in Hispanics than non-Hispanic whites in the United States, little is known about the genetic etiology of the related traits in this population. To identify genetic loci influencing obesity in non-Mexican Hispanics, we performed a genome-wide linkage scan in 1,390 subjects from 100 Caribbean Hispanic families on six obesity-related quantitative traits: body mass index (BMI), body weight, waist circumference, waist-to-hip ratio, abdominal and average triceps skinfold thickness after adjusting for significant demographic and lifestyle factors. We then carried out an association analysis of the linkage peaks and the FTO gene in an independent community-based Hispanic subcohort (N = 652, 64% Caribbean Hispanics) from the Northern Manhattan Study. Evidence of linkage was strongest on 1q43 with multipoint LOD score of 2.45 (p = 0.0004) for body weight. Suggestive linkage evidence of LOD > 2.0 was also identified on 1q43 for BMI (LOD = 2.03), 14q32 for abdominal skinfold thickness (LOD = 2.17), 16p12 for BMI (LOD = 2.27) and weight (LOD = 2.26), and 16q23–24 for average triceps skinfold thickness (LOD = 2.32). In the association analysis of 6,440 single nucleotide polymorphisms (SNPs) under 1-LOD unit down regions of our linkage peaks on chromosome 1q43 and 16p12 as well as in the FTO gene, we found that two SNPs (rs6665519 and rs669231) on 1q43 and one FTO SNP (rs12447427) were significantly associated with BMI or body weight after adjustment for multiple testing. Our results suggest that in addition to FTO, multiple genetic loci, particularly those on 1q43 region, may contribute to the variations in obesity-related quantitative traits in Caribbean Hispanics.     

Characteristics of the Human Host Have Little Influence on Which Local Schistosoma mansoni Populations Are Acquired

Background

Brazil remains the country in the Americas with the highest prevalence of schistosomiasis. A combination of control efforts and development, however, has sharply reduced its intensity and distribution. The acquisition of specific schistosome populations may be dependent on host characteristics such as sex, age, geography, work, habits and culture. How these and other host characteristics align with parasite subpopulations may guide approaches to improve control.

Methodology

A cohort of more than 90% of the residents in two rural communities in Brazil participated in an epidemiologic survey of demographic, socio-economic and behavioral characteristics. The variables sex, age, intensity of infection, socio-economic index, % lifetime spent on site, previous infection, and trips outside the district were used to group parasites infecting individuals. Schistosoma mansoni infection status was determined by examination of stools submitted on 3 different days. The aggregate of eggs collected from the whole stool was used to determine degree of population differentiation from allele frequencies for 15 microsatellites.

Conclusions/Significance

Infection prevalence was 41% for these communities, and the epidemiologic characteristics were similar to many of the endemic areas of Brazil and the world. Parasite population structuring was observed between the two communities (Jost''s D 0.046, CI95% 0.042–0.051), although separated by only 8 km and connected by a highway. No structuring was observed when infected individuals were stratified by host''s biologic, demographic or epidemiologic characteristics. Those most heavily infected best reflected the communities'' overall parasite diversity. The lack of differentiation within villages suggests that individuals are likely to get infected at the same sites or that the same parasite multilocus genotypes can be found at most sites. The geographic structuring between villages and the lack of structuring by age of the host further supports the impression of a population little affected by migration or drift.     

Whole-Exome Sequencing Efficiently Detects Rare Mutations in Autosomal Recessive Nonsyndromic Hearing Loss

Identification of the pathogenic mutations underlying autosomal recessive nonsyndromic hearing loss (ARNSHL) is difficult, since causative mutations in 39 different genes have so far been reported. After excluding mutations in the most common ARNSHL gene, GJB2, via Sanger sequencing, we performed whole-exome sequencing (WES) in 30 individuals from 20 unrelated multiplex consanguineous families with ARNSHL. Agilent SureSelect Human All Exon 50 Mb kits and an Illumina Hiseq2000 instrument were used. An average of 93%, 84% and 73% of bases were covered to 1X, 10X and 20X within the ARNSHL-related coding RefSeq exons, respectively. Uncovered regions with WES included those that are not targeted by the exome capture kit and regions with high GC content. Twelve homozygous mutations in known deafness genes, of which eight are novel, were identified in 12 families: MYO15A-p.Q1425X, -p.S1481P, -p.A1551D; LOXHD1-p.R1494X, -p.E955X; GIPC3-p.H170N; ILDR1-p.Q274X; MYO7A-p.G2163S; TECTA-p.Y1737C; TMC1-p.S530X; TMPRSS3-p.F13Lfs*10; TRIOBP-p.R785Sfs*50. Each mutation was within a homozygous run documented via WES. Sanger sequencing confirmed co-segregation of the mutation with deafness in each family. Four rare heterozygous variants, predicted to be pathogenic, in known deafness genes were detected in 12 families where homozygous causative variants were already identified. Six heterozygous variants that had similar characteristics to those abovementioned variants were present in 15 ethnically-matched individuals with normal hearing. Our results show that rare causative mutations in known ARNSHL genes can be reliably identified via WES. The excess of heterozygous variants should be considered during search for causative mutations in ARNSHL genes, especially in small-sized families.     

Southeastern U.S.A. continental shelf respiratory rates revisited

Respiratory rates on the U. S. southeastern continental shelf have been estimated several times by different investigators, most recently by Jiang et al. (Biogeochemistry 98:101–113, 2010) who report lower mean rates than were found in earlier work and attribute the differences to analytical error in all methods used in earlier studies. The differences are, instead, attributable to the differences in the geographical scope of the studies. The lower estimates of regional organic carbon flux of Jiang et al. (Biogeochemistry 98:101–113, 2010) are a consequence of their extrapolation of data from a small portion of the shelf to the entire South Atlantic Bight. This comment examines the methodologies used as well as the variability of respiratory rates in this region over space and time.     

A computational approach to characterization of bovine sperm chromatin alterations

We describe here a computational morphology-based approach to the investigation of possible causes of chromatin alterations in sperm. A comprehensive set of state-of-the-art and geometric measures are computationally extracted from toluidine blue stained images and analyzed to infer the possible processes leading to normal and abnormal chromatin formation while seeking a possible taxonomy of chromatin alterations and their influence on sperm head morphology. Using this methodology, we have identified higher chromatin fragility at some specific points of the sperm head. Despite the lack of correlation between morphologies of sperm head and chromatin structure, four main morphological types of chromatin alterations in bull spermatozoa have been identified and their possible causes discussed.     

A comparative analysis of the genetic epidemiology of deafness in the United States in two sets of pedigrees collected more than a century apart

In 1898, E.A. Fay published an analysis of nearly 5000 marriages among deaf individuals in America collected during the 19^th century. Each pedigree included three-generation data on marriage partners that included at least one deaf proband, who were ascertained by complete selection. We recently proposed that the intense phenotypic assortative mating among the deaf might have greatly accelerated the normally slow response to relaxed genetic selection against deafness that began in many Western countries with the introduction of sign language and the establishment of residential schools. Simulation studies suggest that this mechanism might have doubled the frequency of the commonest forms of recessive deafness (DFNB1) in this country during the past 200 years. To test this prediction, we collected pedigree data on 311 contemporary marriages among deaf individuals that were comparable to those collected by Fay. Segregation analysis of the resulting data revealed that the estimated proportion of noncomplementary matings that can produce only deaf children has increased by a factor of more than five in the past 100 years. Additional analysis within our sample of contemporary pedigrees showed that there was a statistically significant linear increase in the prevalence of pathologic GJB2 mutations when the data on 441 probands were partitioned into three 20-year birth cohorts (1920 through 1980). These data are consistent with the increase in the frequency of DFNB1 predicted by our previous simulation studies and provide convincing evidence for the important influence that assortative mating can have on the frequency of common genes for deafness.