Hierarchical clustering using the arithmetic-harmonic cut: complexity and experiments

Clustering, particularly hierarchical clustering, is an important method for understanding and analysing data across a wide variety of knowledge domains with notable utility in systems where the data can be classified in an evolutionary context. This paper introduces a new hierarchical clustering problem defined by a novel objective function we call the arithmetic-harmonic cut. We show that the problem of finding such a cut is NP-hard and APX-hard but is fixed-parameter tractable, which indicates that although the problem is unlikely to have a polynomial time algorithm (even for approximation), exact parameterized and local search based techniques may produce workable algorithms. To this end, we implement a memetic algorithm for the problem and demonstrate the effectiveness of the arithmetic-harmonic cut on a number of datasets including a cancer type dataset and a corona virus dataset. We show favorable performance compared to currently used hierarchical clustering techniques such as k-Means, Graclus and Normalized-Cut. The arithmetic-harmonic cut metric overcoming difficulties other hierarchical methods have in representing both intercluster differences and intracluster similarities.     

Aptamer-based multiplexed proteomic technology for biomarker discovery

Background

The interrogation of proteomes (“proteomics”) in a highly multiplexed and efficient manner remains a coveted and challenging goal in biology and medicine.

Methodology/Principal Findings

We present a new aptamer-based proteomic technology for biomarker discovery capable of simultaneously measuring thousands of proteins from small sample volumes (15 µL of serum or plasma). Our current assay measures 813 proteins with low limits of detection (1 pM median), 7 logs of overall dynamic range (∼100 fM–1 µM), and 5% median coefficient of variation. This technology is enabled by a new generation of aptamers that contain chemically modified nucleotides, which greatly expand the physicochemical diversity of the large randomized nucleic acid libraries from which the aptamers are selected. Proteins in complex matrices such as plasma are measured with a process that transforms a signature of protein concentrations into a corresponding signature of DNA aptamer concentrations, which is quantified on a DNA microarray. Our assay takes advantage of the dual nature of aptamers as both folded protein-binding entities with defined shapes and unique nucleotide sequences recognizable by specific hybridization probes. To demonstrate the utility of our proteomics biomarker discovery technology, we applied it to a clinical study of chronic kidney disease (CKD). We identified two well known CKD biomarkers as well as an additional 58 potential CKD biomarkers. These results demonstrate the potential utility of our technology to rapidly discover unique protein signatures characteristic of various disease states.

Conclusions/Significance

We describe a versatile and powerful tool that allows large-scale comparison of proteome profiles among discrete populations. This unbiased and highly multiplexed search engine will enable the discovery of novel biomarkers in a manner that is unencumbered by our incomplete knowledge of biology, thereby helping to advance the next generation of evidence-based medicine.     

Tempo and mode of performance evolution across multiple independent origins of adhesive toe pads in lizards

Understanding macroevolutionary dynamics of trait evolution is an important endeavor in evolutionary biology. Ecological opportunity can liberate a trait as it diversifies through trait space, while genetic and selective constraints can limit diversification. While many studies have examined the dynamics of morphological traits, diverse morphological traits may yield the same or similar performance and as performance is often more proximately the target of selection, examining only morphology may give an incomplete understanding of evolutionary dynamics. Here, we ask whether convergent evolution of pad‐bearing lizards has followed similar evolutionary dynamics, or whether independent origins are accompanied by unique constraints and selective pressures over macroevolutionary time. We hypothesized that geckos and anoles each have unique evolutionary tempos and modes. Using performance data from 59 species, we modified Brownian motion (BM) and Ornstein–Uhlenbeck (OU) models to account for repeated origins estimated using Bayesian ancestral state reconstructions. We discovered that adhesive performance in geckos evolved in a fashion consistent with Brownian motion with a trend, whereas anoles evolved in bounded performance space consistent with more constrained evolution (an Ornstein–Uhlenbeck model). Our results suggest that convergent phenotypes can have quite distinctive evolutionary patterns, likely as a result of idiosyncratic constraints or ecological opportunities.     

Using consecutive prescribed fires to reduce shrub encroachment in grassland by increasing shrub mortality

Woody plant encroachment into open grasslands occurs worldwide and causes multiple ecological and management impacts. Prescribed fire could be used to conserve grassland habitat but often has limited efficacy because many woody plants resprout after fire and rapidly reestablish abundance. If fire‐induced mortality could be increased, prescribed fire would be a more effective management tool. In California's central coast, shrub encroachment, especially of Baccharis pilularis (coyote brush), is converting coastal prairie into shrub‐dominated communities, with a consequent loss of native herbaceous species and open grassland habitat. B. pilularis has not been successfully controlled with single prescribed fire events because the shrub resprouts and reestablishes cover within a few years. We investigated whether two consecutive annual burns would control B. pilularis by killing resprouting shrubs, without reducing native herbaceous species or encouraging invasive plants. As expected, resprouting did occur; however, 2 years after the second burn, B. pilularis cover on burned plots was only 41% of the cover on unburned plots. Mortality of B. pilularis more than doubled following the second burn, likely maintaining a reduction in B. pilularis cover for longer than a single burn would have. Three native coastal prairie perennial grasses did not appear to be adversely affected by the two burns, nor did the burns result in increased cover of invasive species. Managers wanting to restore coastal prairie following B. pilularis encroachment should consider two consecutive annual burns, especially if moderate fire intensity is achievable.     

Telomeric repeat [TTAGGG]n sequences of human chromosomes are conserved in chimpanzee (Pan troglodytes)

Summary Using a series of genetic parameters, attempts have been made for more than two decades to establish the close kinship of human (Homo sapiens) with chimpanzee (Pan troglodytes). Molecular and cytogenetic data presently suggest that the two species are closely related. The recent isolation of a human telomeric probe (P5097-B.5) has prompted us to cross hybridize it to chimpanzee chromosomes in order to explore convergence and/or divergence of the telomeric repeat sequences (TTAGGG)_n. On hybridization, the human probe bound to both ends (telomeres) of chimpanzee chromosomes, suggesting a concerted evolution of tandemly repeated short simple sequences (TTAGGG)_n. Even the terminal heterochromatin of chimpanzee chromosomes was found to be endowed with telomeric repeats, suggesting that evolution of heterochromatin and capping with tandemly repeated short sequences are highly complex phenomena.     

Receptors and signaling pathway underlying relaxations to isoprostanes in canine and porcine airway smooth muscle

Using muscle bath techniques, we examined the inhibitory activities of several E- and F-ring isoprostanes in canine and porcine airway smooth muscle. 8-Isoprostaglandin E1 and 8-isoprostaglandin E2 (8-iso PGE2) reversed cholinergic tone in a concentration-dependent manner, whereas the F-ring isoprostanes were ineffective. Desensitization with 8-iso-PGE2 and PGE2 implicated isoprostane activity at the PGE2 receptor (EP). We found that the inhibitory E-ring isoprostane responses were significantly augmented by rolipram (a type IV phosphodiesterase inhibitor), while 1H-[1,2,4]-oxadiazolo[4,3-a]quinoxalin-1-one (a guanylate cyclase inhibitor) had no effect, suggesting a role for cAMP in isoprostane-mediated relaxations. 8-Iso-PGE2 did not reverse KCl tone, suggesting that voltage-dependent Ca2+ influx and myosin light chain kinase are not suppressed by isoprostanes. Patch-clamp studies showed marked suppression of K+ currents by 8-iso-PGE2. We conclude that E-ring isoprostanes exert PGE2 receptor-directed, cAMP-dependent relaxations in canine and porcine airway smooth muscle. This activity is not dependent on K+ channel activation or the direct inhibition of voltage-operated Ca2+ influx or myosin light chain kinase.     

Design,synthesis, and evaluation of 3,7-substituted coumarin derivatives as multifunctional Alzheimer’s disease agents

Multitarget directed ligands (MTDLs) are emerging as promising treatment options for Alzheimer’s disease (AD). Coumarin derivatives serve as a good starting point for designing MTDLs due to their inherent inhibition of monoamine oxidase (MAO) and cholinesterase enzymes, which are complicit in AD’s complex pathophysiology. A preliminary series of 3,7-substituted coumarin derivatives were synthesised and evaluated for enzyme inhibitory activity, cytotoxicity as well as neuroprotective ability. The results indicated that the compounds are weak cholinesterase inhibitors with five compounds demonstrating relatively potent inhibition and selectivity towards MAO-B with IC₅₀ values between 0.014 and 0.498 hx00B5;µM. Significant neuroprotective effects towards MPP⁺-compromised SH-SY5Y neuroblastoma cells were also observed, with no inherent cytotoxicity at 10 µM for all compounds. The overall results demonstrated that substitution of the phenylethyloxy moiety at the 7-position imparted superior general activity to the derivatives, with the propargylamine substitution at the 3-position, in particular, displaying the best MAO-B selectivity and neuroprotection.     

Rapid Parallel Adaptation to Anthropogenic Heavy Metal Pollution

The impact of human-mediated environmental change on the evolutionary trajectories of wild organisms is poorly understood. In particular, capacity of species to adapt rapidly (in hundreds of generations or less), reproducibly and predictably to extreme environmental change is unclear. Silene uniflora is predominantly a coastal species, but it has also colonized isolated, disused mines with phytotoxic, zinc-contaminated soils. To test whether rapid, parallel adaptation to anthropogenic pollution has taken place, we used reduced representation sequencing (ddRAD) to reconstruct the evolutionary history of geographically proximate mine and coastal population pairs and found largely independent colonization of mines from different coastal sites. Furthermore, our results show that parallel evolution of zinc tolerance has occurred without gene flow spreading adaptive alleles between mine populations. In genomic regions where signatures of selection were detected across multiple mine-coast pairs, we identified genes with functions linked to physiological differences between the putative ecotypes, although genetic differentiation at specific loci is only partially shared between mine populations. Our results are consistent with a complex, polygenic genetic architecture underpinning rapid adaptation. This shows that even under a scenario of strong selection and rapid adaptation, evolutionary responses to human activities (and other environmental challenges) may be idiosyncratic at the genetic level and, therefore, difficult to predict from genomic data.