Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “unsupervised clustering”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11Linked to original sources

Unsupervised learning with independent component analysis can identify patterns of glaucomatous visual field defects.

PURPOSE: We previously reported the use of clustering by unsupervised learning with machine learning classifiers to segment clusters of patterns in standard automated perimetry (SAP) for glaucoma. In this study, the process of unsupervised learning by independent component analysis decomposed SAP field patterns into axes, and the information represented by these axes was evaluated. METHODS: SAP fields were obtained with the Humphrey Visual Field Analyzer on 189 normal eyes and 156 eyes with glaucomatous optic neuropathy (GON) determined by masked review with stereoscopic optic disc photos. The variational Bayesian independent component analysis mixture model (vB-ICA-mm) partitioned the SAP fields into the most informative number of clusters. Simultaneously, it learned an optimal number of maximally independent axes for each cluster. RESULTS: The most informative number of clusters was two. vB-ICA-mm placed 68.6% of the SAP fields from eyes with GON in a cluster labeled G and 98.4% of the fields from eyes with normal optic discs in a cluster labeled N. Cluster G optimally contained six axes. Post hoc analysis of patterns generated at -1 SD and +2 SD from the cluster G mean on the six axes revealed defects similar to those identified by experts as indicative of glaucoma. SAP fields associated with an axis showed increasing severity as they were located farther in the positive direction from the cluster G mean. CONCLUSIONS: vB-ICA-mm represented the SAP fields with patterns that were meaningful for glaucoma experts. This process also captured severity in the patterns uncovered. These findings should validate vB-ICA-mm as a data mining technique for new and unfamiliar complex tests.

Artificial Intelligence↗

Altered expression of TFF-1 and CES-2 in Barrett's Esophagus and associated adenocarcinomas.

Identification of biomarkers to recognize individuals with Barrett's esophagus (BE) predisposed to develop malignancy is currently a pressing issue. We utilized gene expression profiling to compare molecular signatures of normal esophagus and stomach, BE, and adenocarcinoma (AC) to identify such potential biomarkers. Over 22,000 genes were analyzed by oligonucleotide microarrays on 38 unique RNA Unsupervised and supervised clusterings were performed on a subset of 2849 genes that varied most significantly across the specimens. Immunohistochemistry (IHC) for two of the significantly differentially expressed gene products was performed on tissue microarrays. Unsupervised clustering identified two discernable molecular BE profiles, one of which was similar to normal gastric tissue ("BE1"), and another that was shared by several of the AC specimens ("BE2"). The BE1 profile included expression of several genes that have been described as tumor-suppressor genes, most notably trefoil factor 1 (TFF-1). The BE2 profile included expression of genes previously found overexpressed in cancers, such as carboxylesterase-2 (CES-2). IHC demonstrated the loss of TFF-1 late in the progression of BE to AC. It also revealed CES-2 as being upregulated in AC documented to have arisen in the presence of BE. These potential biomarkers, as well as the relative expression of genes from BE1 versus those from BE2, may be validated in the future to aid in risk stratification and guide treatment protocols in patients with BE and associated AC.

Adenocarcinoma↗

Exploring supervised and unsupervised methods to detect topics in biomedical text.

BACKGROUND: Topic detection is a task that automatically identifies topics (e.g., "biochemistry" and "protein structure") in scientific articles based on information content. Topic detection will benefit many other natural language processing tasks including information retrieval, text summarization and question answering; and is a necessary step towards the building of an information system that provides an efficient way for biologists to seek information from an ocean of literature. RESULTS: We have explored the methods of Topic Spotting, a task of text categorization that applies the supervised machine-learning technique naïve Bayes to assign automatically a document into one or more predefined topics; and Topic Clustering, which apply unsupervised hierarchical clustering algorithms to aggregate documents into clusters such that each cluster represents a topic. We have applied our methods to detect topics of more than fifteen thousand of articles that represent over sixteen thousand entries in the Online Mendelian Inheritance in Man (OMIM) database. We have explored bag of words as the features. Additionally, we have explored semantic features; namely, the Medical Subject Headings (MeSH) that are assigned to the MEDLINE records, and the Unified Medical Language System (UMLS) semantic types that correspond to the MeSH terms, in addition to bag of words, to facilitate the tasks of topic detection. Our results indicate that incorporating the MeSH terms and the UMLS semantic types as additional features enhances the performance of topic detection and the naïve Bayes has the highest accuracy, 66.4%, for predicting the topic of an OMIM article as one of the total twenty-five topics. CONCLUSION: Our results indicate that the supervised topic spotting methods outperformed the unsupervised topic clustering; on the other hand, the unsupervised topic clustering methods have the advantages of being robust and applicable in real world settings.

Abstracting and Indexing↗

A hierarchical unsupervised growing neural network for clustering gene expression patterns.

MOTIVATION: We describe a new approach to the analysis of gene expression data coming from DNA array experiments, using an unsupervised neural network. DNA array technologies allow monitoring thousands of genes rapidly and efficiently. One of the interests of these studies is the search for correlated gene expression patterns, and this is usually achieved by clustering them. The Self-Organising Tree Algorithm, (SOTA) (Dopazo,J. and Carazo,J.M. (1997) J. Mol. Evol., 44, 226-233), is a neural network that grows adopting the topology of a binary tree. The result of the algorithm is a hierarchical cluster obtained with the accuracy and robustness of a neural network. RESULTS: SOTA clustering confers several advantages over classical hierarchical clustering methods. SOTA is a divisive method: the clustering process is performed from top to bottom, i.e. the highest hierarchical levels are resolved before going to the details of the lowest levels. The growing can be stopped at the desired hierarchical level. Moreover, a criterion to stop the growing of the tree, based on the approximate distribution of probability obtained by randomisation of the original data set, is provided. By means of this criterion, a statistical support for the definition of clusters is proposed. In addition, obtaining average gene expression patterns is a built-in feature of the algorithm. Different neurons defining the different hierarchical levels represent the averages of the gene expression patterns contained in the clusters. Since SOTA runtimes are approximately linear with the number of items to be classified, it is especially suitable for dealing with huge amounts of data. The method proposed is very general and applies to any data providing that they can be coded as a series of numbers and that a computable measure of similarity between data items can be used. AVAILABILITY: A server running the program can be found at: http://bioinfo.cnio.es/sotarray.

Algorithms↗

Pitfalls in the clustering of neuroimage data and improvements by global optimization strategies.

In this paper, we examined three vector quantization (VQ) methods used for the unsupervised classification (clustering) of functional magnetic resonance imaging (fMRI) data. Classification means that each brain volume element (voxel), according to a given scanning raster, was assigned to one group of voxels based on similarity of the fMRI signal patterns. It was investigated how the VQ methods can isolate a cluster that describes the region involved in a particular brain function. As an example, word processing was stimulated by a word comparison task. VQ analysis methodology was verified using simulated fMRI response patterns. It was demonstrated in detail that VQ based on global rather than local optimization of the objective function yielded a higher performance. Performance was measured in statistically relevant series of VQ attempts using several indices for goodness, reliability and efficiency of VQ solutions. Furthermore, it was shown that a poor local optimization caused either an underestimation or an overestimation of the stimulus-induced brain activation. However, this was not observed if the cluster analysis was based upon a global optimization strategy.

Adult↗

TreeSOM: Cluster analysis in the self-organizing map.

Clustering problems arise in various domains of science and engineering. A large number of methods have been developed to date. The Kohonen self-organizing map (SOM) is a popular tool that maps a high-dimensional space onto a small number of dimensions by placing similar elements close together, forming clusters. Cluster analysis is often left to the user. In this paper we present the method TreeSOM and a set of tools to perform unsupervised SOM cluster analysis, determine cluster confidence and visualize the result as a tree facilitating comparison with existing hierarchical classifiers. We also introduce a distance measure for cluster trees that allows one to select a SOM with the most confident clusters.

Age Factors↗

A DNA microarray survey of gene expression in normal human tissues.

BACKGROUND: Numerous studies have used DNA microarrays to survey gene expression in cancer and other disease states. Comparatively little is known about the genes expressed across the gamut of normal human tissues. Systematic studies of global gene-expression patterns, by linking variation in the expression of specific genes to phenotypic variation in the cells or tissues in which they are expressed, provide clues to the molecular organization of diverse cells and to the potential roles of the genes. RESULTS: Here we describe a systematic survey of gene expression in 115 human tissue samples representing 35 different tissue types, using cDNA microarrays representing approximately 26,000 different human genes. Unsupervised hierarchical cluster analysis of the gene-expression patterns in these tissues identified clusters of genes with related biological functions and grouped the tissue specimens in a pattern that reflected their anatomic locations, cellular compositions or physiologic functions. In unsupervised and supervised analyses, tissue-specific patterns of gene expression were readily discernable. By comparative hybridization to normal genomic DNA, we were also able to estimate transcript abundances for expressed genes. CONCLUSIONS: Our dataset provides a baseline for comparison to diseased tissues, and will aid in the identification of tissue-specific functions. In addition, our analysis identifies potential molecular markers for detection of injury to specific organs and tissues, and provides a foundation for selection of potential targets for selective anticancer therapy.

Cluster Analysis↗

SOPHIA: an interactive cluster-based retrieval system for the OHSUMED collection.

The ability to perform an exploratory search and retrieval of relevant documents from a large collection of domain-specific documents is an important requirement both in the field of medicine and other areas. In this paper, we present a unsupervised distributional clustering technique called SOPHIA. SOPHIA provides a semantically meaningful visual clustering of the document corpus in conjunction with an intuitive interactive search facility. We assess the effectiveness of SOPHIA's cluster-based information retrieval for the MEDLINE testset collection known as OHSUMED.

Artificial Intelligence↗

Gene expression profiles of normal proliferating and differentiating human intestinal epithelial cells: a comparison with the Caco-2 cell model.

cDNA microarray technology enables detailed analysis of gene expression throughout complex processes such as differentiation. The aim of this study was to analyze the gene expression profile of normal human intestinal epithelial cells using cell models that recapitulate the crypt-villus axis of intestinal differentiation in comparison with the widely used Caco-2 cell model. cDNA microarrays (19,200 human genes) and a clustering algorithm were used to identify patterns of gene expression in the crypt-like proliferative HIEC and tsFHI cells, and villus epithelial cells as well as Caco-2/15 cells at two distinct stages of differentiation. Unsupervised hierarchical clustering analysis of global gene expression among the cell lines identified two branches: one for the HIEC cells versus a second comprised of two sub-groups: (a) the proliferative Caco-2 cells and (b) the differentiated Caco-2 cells and closely related villus epithelial cells. At the gene level, supervised hierarchical clustering with 272 differentially expressed genes revealed distinct expression patterns specific to each cell phenotype. We identified several upregulated genes that could lead to the identification of new regulatory pathways involved in cell differentiation and carcinogenesis. The combined use of microarray analysis and human intestinal cell models thus provides a powerful tool for establishing detailed gene expression profiles of proliferative to terminally differentiated intestinal cells. Furthermore, the molecular differences between the normal human intestinal cell models and Caco-2 cells clearly point out the strengths and limitations of this widely used experimental model for studying intestinal cell proliferation and differentiation.

Caco-2 Cells↗

Cytokeratin profiles identify diagnostic signatures in colorectal cancer using multiplex analysis of tissue microarrays.

BACKGROUND AND AIMS: Recent cDNA expression profiling analyses indicate that within specific organ cancers Cytokeratins (CKs) dysregulation may identify subgroups with distinct biological phenotypes. Our objectives in this study were (1) to test whether cytokeratins were also distinct on the protein level, (2) to evaluate these biomarkers in a series of well-characterised CRCs, (3) to apply hierarchical cluster analysis to immunohistochemical data. METHODS: Tissue microarrays (TMA) comprising 468 CRC specimens from 203 patients were constructed to evaluate CK5, CK7, CK8, CK13, CK14, CK16, CK17, CK18, CK19 and CK20. In total, 2919 samples were analyzed. RESULTS: Unsupervised hierarchical clustering discovered subgroups represented by reduced CK8 and CK20 expression, that differed by a shorter patients survival. The evaluation of the specific biomarkers by Kaplan-Meier analysis showed that reduced CK8 expression (p<0.01) was significantly associated with shorter patients' survival, but was not an independent factor correlated with tumour stage (pT), grading (G) and nodal stage (pN). CONCLUSIONS: Reduced coexpression of CK8 and CK20 may indicate an epithelial-mesenchymal transition (EMT) representing an important step in the development of more aggressive CRCs. In addition, multiplex analysis of TMAs together with immunohistochemistry (IHC) supplemented by hierarchical clustering are a useful, promising and very powerful tool for the identification of tumour subgroups with diagnostic and prognostic signatures.

Biomarkers, Tumor↗

A fuzzy clustering based segmentation system as support to diagnosis in medical imaging.

In medical imaging uncertainty is widely present in data, because of the noise in acquisition and of the partial volume effects originating from the low resolution of sensors. In particular, borders between tissues are not exactly defined and memberships in the boundary regions are intrinsically fuzzy. Therefore, computer assisted unsupervised fuzzy clustering methods turn out to be particularly suitable for handling a decision making process concerning segmentation of multimodal medical images. By using the possibilistic c-means algorithm as a refinement of a neural network based clustering algorithm named capture effect neural network, we developed the possibilistic neuro fuzzy c-means algorithm (PNFCM). In this paper the PNFCM has been applied to two different multimodal data sets and the results have been compared to those obtained by using the classical fuzzy c-means algorithm. Furthermore, a discussion is presented about the role of fuzzy clustering as a support to diagnosis in medical imaging.

Algorithms↗

Functional genomics and proteomics in the clinical neurosciences: data mining and bioinformatics.

The goal of this chapter is to introduce some of the available computational methods for expression analysis. Genomic and proteomic experimental techniques are briefly discussed to help the reader understand these methods and results better in context with the biological significance. Furthermore, a case study is presented that will illustrate the use of these analytical methods to extract significant biomarkers from high-throughput microarray data. Genomic and proteomic data analysis is essential for understanding the underlying factors that are involved in human disease. Currently, such experimental data are generally obtained by high-throughput microarray or mass spectrometry technologies among others. The sheer amount of raw data obtained using these methods warrants specialized computational methods for data analysis. Biomarker discovery for neurological diagnosis and prognosis is one such example. By extracting significant genomic and proteomic biomarkers in controlled experiments, we come closer to understanding how biological mechanisms contribute to neural degenerative diseases such as Alzheimers' and how drug treatments interact with the nervous system. In the biomarker discovery process, there are several computational methods that must be carefully considered to accurately analyze genomic or proteomic data. These methods include quality control, clustering, classification, feature ranking, and validation. Data quality control and normalization methods reduce technical variability and ensure that discovered biomarkers are statistically significant. Preprocessing steps must be carefully selected since they may adversely affect the results of the following expression analysis steps, which generally fall into two categories: unsupervised and supervised. Unsupervised or clustering methods can be used to group similar genomic or proteomic profiles and therefore can elucidate relationships within sample groups. These methods can also assign biomarkers to sub-groups based on their expression profiles across patient samples. Although clustering is useful for exploratory analysis, it is limited due to its inability to incorporate expert knowledge. On the other hand, classification and feature ranking are supervised, knowledge-based machine learning methods that estimate the distribution of biological expression data and, in doing so, can extract important information about these experiments. Classification is closely coupled with feature ranking, which is essentially a data reduction method that uses classification error estimation or other statistical tests to score features. Biomarkers can subsequently be extracted by eliminating insignificantly ranked features. These analytical methods may be equally applied to genetic and proteomic data. However, because of both biological differences between the data sources and technical differences between the experimental methods used to obtain these data, it is important to have a firm understanding of the data sources and experimental methods. At the same time, regardless of the data quality, it is inevitable that some discovered biomarkers are false positives. Thus, it is important to validate discovered biomarkers. The validation process may be slow; yet, the overall biomarker discovery process is significantly accelerated due to initial feature ranking and data reduction steps. Information obtained from the validation process may also be used to refine data analysis procedures for future iteration. Biomarker validation may be performed in a number of ways - bench-side in traditional labs, web-based electronic resources such as gene ontology and literature databases, and clinical trials.

Animals↗

Use of gene-expression profiling to identify prognostic subclasses in adult acute myeloid leukemia.

BACKGROUND: In patients with acute myeloid leukemia (AML), the presence or absence of recurrent cytogenetic aberrations is used to identify the appropriate therapy. However, the current classification system does not fully reflect the molecular heterogeneity of the disease, and treatment stratification is difficult, especially for patients with intermediate-risk AML with a normal karyotype. METHODS: We used complementary-DNA microarrays to determine the levels of gene expression in peripheral-blood samples or bone marrow samples from 116 adults with AML (including 45 with a normal karyotype). We used unsupervised hierarchical clustering analysis to identify molecular subgroups with distinct gene-expression signatures. Using a training set of samples from 59 patients, we applied a novel supervised learning algorithm to devise a gene-expression-based clinical-outcome predictor, which we then tested using an independent validation group comprising the 57 remaining patients. RESULTS: Unsupervised analysis identified new molecular subtypes of AML, including two prognostically relevant subgroups in AML with a normal karyotype. Using the supervised learning algorithm, we constructed an optimal 133-gene clinical-outcome predictor, which accurately predicted overall survival among patients in the independent validation group (P=0.006), including the subgroup of patients with AML with a normal karyotype (P=0.046). In multivariate analysis, the gene-expression predictor was a strong independent prognostic factor (odds ratio, 8.8; 95 percent confidence interval, 2.6 to 29.3; P<0.001). CONCLUSIONS: The use of gene-expression profiling improves the molecular classification of adult AML.

Acute Disease↗

A fuzzy relational clustering algorithm based on a dissimilarity measure extracted from data.

One of the critical aspects of clustering algorithms is the correct identification of the dissimilarity measure used to drive the partitioning of the data set. The dissimilarity measure induces the cluster shape and therefore determines the success of clustering algorithms. As cluster shapes change from a data set to another, dissimilarity measures should be extracted from data. To this aim, we exploit some pairs of points with known dissimilarity value to teach a dissimilarity relation to a feed-forward neural network. Then, we use the neural dissimilarity measure to guide an unsupervised relational clustering algorithm. Experiments on synthetic data sets and on the Iris data set show that the relational clustering algorithm based on the neural dissimilarity outperforms some popular clustering algorithms (with possible partial supervision) based on spatial dissimilarity.

Journal Article↗

Immunoprofiles of 11 biomarkers using tissue microarrays identify prognostic subgroups in colorectal cancer.

BACKGROUND AND AIMS: Genomewide expression profiling has identified a number of genes expressed at higher levels in colorectal cancer (CRC) than in normal tissues. Our objectives in this study were: 1) to test whether genes were also distinct on the protein level; 2) to evaluate these biomarkers in a series of well-characterized CRCs; and 3) to apply hierarchical cluster analysis to the immunohistochemical data. METHODS: Tissue microarrays (TMAs) comprising 351 CRC specimens from 270 patients were constructed to evaluate the genes Adam10, Cyclin D1, Annexin II, NFKB, Casein kinase 2 beta (CK2B), YB-1, P32, Rad51, c-fos, IGFBP4, and Connexin26 (Cx26). In total, 3,797 samples were analyzed. RESULTS: Unsupervised hierarchical clustering discovered subgroups of CRC that differed by tumor stage and survival. Kaplan-Meier analysis showed that reduced Cx26 expression was significantly associated with shorter patient survival and higher tumor grade (G1/G2 vs G3, P = .02), and Adam10 expression with a higher tumor stage (pT1/2 vs pT3/4, P = .04). CONCLUSIONS: Our study highlights the potential of TMAs for a higher-dimensional analysis by evaluating serial sections of the same tissue core (three-dimensional TMA analysis). In addition, it endorses the use of immunohistochemistry supplemented by hierarchical clustering for the identification of tumor subgroups with diagnostic and prognostic signatures.

Biomarkers, Tumor↗

Novel potential ALL low-risk markers revealed by gene expression profiling with new high-throughput SSH-CCS-PCR.

The current systems of risk grouping in pediatric acute lymphoblastic leukemia (ALL) fail to predict therapeutic success in 10-35% of patients. To identify better predictive markers of clinical behavior in ALL, we have developed an integrated approach for gene expression profiling that couples suppression subtractive hybridization, concatenated cDNA sequencing, and reverse transcriptase real-time quantitative PCR. Using this approach, a total of 600 differentially expressed genes were identified between t(4;11) ALL and pre-B ALL with no determinant chromosomal translocation. The expression of 67 genes was analyzed in different cytogenetic ALL subgroups and B lymphocytes isolated from healthy donors. Three genes, BACH1, TP53BPL, and H2B/S, were consistently expressed as a significant cluster associated with the low-risk ALL subgroups. A total of 42 genes were differentially expressed in ALL vs normal B lymphocytes, with no specific association with any particular ALL subgroups. The remaining 22 genes were part of a specific expression profile associated with the hyperdiploid, t(12;21), or t(4;11) subgroups. Using an unsupervised hierarchical cluster analysis, the discriminating power of these specific expression profiles allowed the clustering of patients according to their subgroups. These genes could help to understand the difference in treatment response and become therapeutical targets to improve ALL clinical outcomes.

B-Lymphocytes↗

Cluster analysis of p-glycoprotein, c-erb-B2 and P53 in relation to tumor histology strongly indicates prognosis in patients with operable non-small cell lung cancer.

BACKGROUND: Various biomarkers have prognostic value in non-small cell lung cancer (NSCLC). We aimed to identify the roles of P53, c-erb- B2 and p-glycoprotein (pgp) as prognostic factors, independently or in conjunction with each other, in operable NSCLC. MATERIAL/METHODS: Seventy operable NSCLC cases were retrospectively evaluated for P53, c-erb-B2 and pgp expression patterns by immunohistochemistry. An unsupervised hierarchical cluster analysis of the 3 biomarkers was conducted. Univariate and multivariate survival analyses were made in relation to cluster affiliation. RESULTS: Cluster analysis yielded two distinct subgroups; group A of high biomarker expressors (n=26, 37%), and group B (n=44, 63%) of low expressors. Cluster affiliation with regard to tumor histology (interaction term) was independently associated with Recurrence- free survival (RFS) and Overall survival (OAS) with a Hazard Ratio (HR) of 5.88, P=0.003, and HR=4.68, P=0.012, respectively. The median OAS times for cluster A and B in the squamous cell carcinoma subgroup were 328 and 596 days, whereas the corresponding figures in the non-squamous cell carcinoma subgroup were non-measurable and 298 days. CONCLUSIONS: In operable NSCLC there may be different relationships of P53, c-erb-B2 and pgp with patient outcome for different tumor histologies. The prognostic utility of cluster affiliation with regard to these biomarkers, and in relation to tumor histology, deserves further testing.

ATP Binding Cassette Transporter, Subfamily B, Mem↗

Interleukin-1beta is prominent in the early pulmonary inflammatory response after hepatic injury.

BACKGROUND: Acute lung injury and inflammation can occur after hepatic ischemia/reperfusion or cryoablation. The etiology of this response is uncertain although it involves NF-kappaB-mediated cytokine release from the liver. METHODS: Inflammation-specific complementary DNA microarrays were utilized to evaluate cytokine upregulation in mouse lung at 4 hours after partial-volume hepatic cryoablation with a recirculating liquid N(2) probe. Hierarchical cluster analysis was performed to identify candidate genes. On the basis of these results, an enzyme-linked immunosorbent assay for interleukin-1beta (IL-1beta) was conducted on serum and pulmonary parenchymal specimens. The time course of IL-1beta transcriptional upregulation in the liver and lungs was evaluated by quantitative reverse transcription/real-time polymerase chain reaction. RESULTS: Starting with a pool of 35 genes generated from normalization and variation filtration, unsupervised hierarchical clustering clearly distinguished lungs of hepatic cryo-injured mice from controls. Genes from the IL-1-family were prominent in the signature. IL-1beta was demonstrable in serum within 2 hours postinjury (218 +/- 89 pg/mL vs 0 at baseline, P = .01). In the lung, IL-1beta was more than 4-fold greater at 4 hours than at baseline. Real-time polymerase chain reaction showed a transcription peak of IL-1beta at 30 minutes in the liver, whereas expression in the lungs remained low until 60 minutes, then continued to increase through 4 hours. CONCLUSIONS: Microarray analysis identified cytokines of the IL-1 family as central components of acute lung injury after hepatic cryoablation. IL-1beta levels increased in both serum and lung tissue over 4 hours after liver injury. Expression of IL-1beta peaked early in the injured liver remnant, followed by subsequent increases in the lungs. Targeted intervention against IL-1beta may ameliorate liver-mediated lung injury.

Animals↗