Nonlinear dynamic model and simulation of morphing wing profile actuated by shape memory alloys

Morphing aircraft have the ability to actively adapt and change their shape to achieve different missions efficiently. The development of morphing structures is deeply related with the ability to model precisely different designs in order to evaluate its characteristics. This paper addresses the dynamic modeling of a sectioned wing profile (morphing airfoil) connected by rotational joints (hinges). In this proposal, a pair of shape memory alloy (SMA) wires are connected to subsequent sections providing torque by reducing its length (changing airfoil camber). The dynamic model of the structure is presented for one pair of sections considering the system with one degree of freedom. The motion equations are solved using numerical techniques due the nonlinearities of the model. The numerical results are compared with experimental data and a discussion of how good this approach captures the physical phenomena associated with this problem. © The Society for Experimental Mechanics, Inc. 2012.

Numerical Error Prediction and its applications in CFD using tau-estimation

Nowadays, Computational Fluid Dynamics (CFD) solvers are widely used within the industry to model fluid flow phenomenons. Several fluid flow model equations have been employed in the last decades to simulate and predict forces acting, for example, on different aircraft configurations. Computational time and accuracy are strongly dependent on the fluid flow model equation and the spatial dimension of the problem considered. While simple models based on perfect flows, like panel methods or potential flow models can be very fast to solve, they usually suffer from a poor accuracy in order to simulate real flows (transonic, viscous). On the other hand, more complex models such as the full Navier- Stokes equations provide high fidelity predictions but at a much higher computational cost. Thus, a good compromise between accuracy and computational time has to be fixed for engineering applications. A discretisation technique widely used within the industry is the so-called Finite Volume approach on unstructured meshes. This technique spatially discretises the flow motion equations onto a set of elements which form a mesh, a discrete representation of the continuous domain. Using this approach, for a given flow model equation, the accuracy and computational time mainly depend on the distribution of nodes forming the mesh. Therefore, a good compromise between accuracy and computational time might be obtained by carefully defining the mesh. However, defining an optimal mesh for complex flows and geometries requires a very high level expertize in fluid mechanics and numerical analysis, and in most cases a simple guess of regions of the computational domain which might affect the most the accuracy is impossible. Thus, it is desirable to have an automatized remeshing tool, which is more flexible with unstructured meshes than its structured counterpart. However, adaptive methods currently in use still have an opened question: how to efficiently drive the adaptation ? Pioneering sensors based on flow features generally suffer from a lack of reliability, so in the last decade more effort has been made in developing numerical error-based sensors, like for instance the adjoint-based adaptation sensors. While very efficient at adapting meshes for a given functional output, the latter method is very expensive as it requires to solve a dual set of equations and computes the sensor on an embedded mesh. Therefore, it would be desirable to develop a more affordable numerical error estimation method. The current work aims at estimating the truncation error, which arises when discretising a partial differential equation. These are the higher order terms neglected in the construction of the numerical scheme. The truncation error provides very useful information as it is strongly related to the flow model equation and its discretisation. On one hand, it is a very reliable measure of the quality of the mesh, therefore very useful in order to drive a mesh adaptation procedure. On the other hand, it is strongly linked to the flow model equation, so that a careful estimation actually gives information on how well a given equation is solved, which may be useful in the context of _ -extrapolation or zonal modelling. The following work is organized as follows: Chap. 1 contains a short review of mesh adaptation techniques as well as numerical error prediction. In the first section, Sec. 1.1, the basic refinement strategies are reviewed and the main contribution to structured and unstructured mesh adaptation are presented. Sec. 1.2 introduces the definitions of errors encountered when solving Computational Fluid Dynamics problems and reviews the most common approaches to predict them. Chap. 2 is devoted to the mathematical formulation of truncation error estimation in the context of finite volume methodology, as well as a complete verification procedure. Several features are studied, such as the influence of grid non-uniformities, non-linearity, boundary conditions and non-converged numerical solutions. This verification part has been submitted and accepted for publication in the Journal of Computational Physics. Chap. 3 presents a mesh adaptation algorithm based on truncation error estimates and compares the results to a feature-based and an adjoint-based sensor (in collaboration with Jorge Ponsín, INTA). Two- and three-dimensional cases relevant for validation in the aeronautical industry are considered. This part has been submitted and accepted in the AIAA Journal. An extension to Reynolds Averaged Navier- Stokes equations is also included, where _ -estimation-based mesh adaptation and _ -extrapolation are applied to viscous wing profiles. The latter has been submitted in the Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering. Keywords: mesh adaptation, numerical error prediction, finite volume Hoy en día, la Dinámica de Fluidos Computacional (CFD) es ampliamente utilizada dentro de la industria para obtener información sobre fenómenos fluidos. La Dinámica de Fluidos Computacional considera distintas modelizaciones de las ecuaciones fluidas (Potencial, Euler, Navier-Stokes, etc) para simular y predecir las fuerzas que actúan, por ejemplo, sobre una configuración de aeronave. El tiempo de cálculo y la precisión en la solución depende en gran medida de los modelos utilizados, así como de la dimensión espacial del problema considerado. Mientras que modelos simples basados en flujos perfectos, como modelos de flujos potenciales, se pueden resolver rápidamente, por lo general aducen de una baja precisión a la hora de simular flujos reales (viscosos, transónicos, etc). Por otro lado, modelos más complejos tales como el conjunto de ecuaciones de Navier-Stokes proporcionan predicciones de alta fidelidad, a expensas de un coste computacional mucho más elevado. Por lo tanto, en términos de aplicaciones de ingeniería se debe fijar un buen compromiso entre precisión y tiempo de cálculo. Una técnica de discretización ampliamente utilizada en la industria es el método de los Volúmenes Finitos en mallas no estructuradas. Esta técnica discretiza espacialmente las ecuaciones del movimiento del flujo sobre un conjunto de elementos que forman una malla, una representación discreta del dominio continuo. Utilizando este enfoque, para una ecuación de flujo dado, la precisión y el tiempo computacional dependen principalmente de la distribución de los nodos que forman la malla. Por consiguiente, un buen compromiso entre precisión y tiempo de cálculo se podría obtener definiendo cuidadosamente la malla, concentrando sus elementos en aquellas zonas donde sea estrictamente necesario. Sin embargo, la definición de una malla óptima para corrientes y geometrías complejas requiere un nivel muy alto de experiencia en la mecánica de fluidos y el análisis numérico, así como un conocimiento previo de la solución. Aspecto que en la mayoría de los casos no está disponible. Por tanto, es deseable tener una herramienta que permita adaptar los elementos de malla de forma automática, acorde a la solución fluida (remallado). Esta herramienta es generalmente más flexible en mallas no estructuradas que con su homóloga estructurada. No obstante, los métodos de adaptación actualmente en uso todavía dejan una pregunta abierta: cómo conducir de manera eficiente la adaptación. Sensores pioneros basados en las características del flujo en general, adolecen de una falta de fiabilidad, por lo que en la última década se han realizado grandes esfuerzos en el desarrollo numérico de sensores basados en el error, como por ejemplo los sensores basados en el adjunto. A pesar de ser muy eficientes en la adaptación de mallas para un determinado funcional, este último método resulta muy costoso, pues requiere resolver un doble conjunto de ecuaciones: la solución y su adjunta. Por tanto, es deseable desarrollar un método numérico de estimación de error más asequible. El presente trabajo tiene como objetivo estimar el error local de truncación, que aparece cuando se discretiza una ecuación en derivadas parciales. Estos son los términos de orden superior olvidados en la construcción del esquema numérico. El error de truncación proporciona una información muy útil sobre la solución: es una medida muy fiable de la calidad de la malla, obteniendo información que permite llevar a cabo un procedimiento de adaptación de malla. Está fuertemente relacionado al modelo matemático fluido, de modo que una estimación precisa garantiza la idoneidad de dicho modelo en un campo fluido, lo que puede ser útil en el contexto de modelado zonal. Por último, permite mejorar la precisión de la solución resolviendo un nuevo sistema donde el error local actúa como término fuente (_ -extrapolación). El presenta trabajo se organiza de la siguiente manera: Cap. 1 contiene una breve reseña de las técnicas de adaptación de malla, así como de los métodos de predicción de los errores numéricos. En la primera sección, Sec. 1.1, se examinan las estrategias básicas de refinamiento y se presenta la principal contribución a la adaptación de malla estructurada y no estructurada. Sec 1.2 introduce las definiciones de los errores encontrados en la resolución de problemas de Dinámica Computacional de Fluidos y se examinan los enfoques más comunes para predecirlos. Cap. 2 está dedicado a la formulación matemática de la estimación del error de truncación en el contexto de la metodología de Volúmenes Finitos, así como a un procedimiento de verificación completo. Se estudian varias características que influyen en su estimación: la influencia de la falta de uniformidad de la malla, el efecto de las no linealidades del modelo matemático, diferentes condiciones de contorno y soluciones numéricas no convergidas. Esta parte de verificación ha sido presentada y aceptada para su publicación en el Journal of Computational Physics. Cap. 3 presenta un algoritmo de adaptación de malla basado en la estimación del error de truncación y compara los resultados con sensores de featured-based y adjointbased (en colaboración con Jorge Ponsín del INTA). Se consideran casos en dos y tres dimensiones, relevantes para la validación en la industria aeronáutica. Este trabajo ha sido presentado y aceptado en el AIAA Journal. También se incluye una extensión de estos métodos a las ecuaciones RANS (Reynolds Average Navier- Stokes), en donde adaptación de malla basada en _ y _ -extrapolación son aplicados a perfiles con viscosidad de alas. Este último trabajo se ha presentado en los Actas de la Institución de Ingenieros Mecánicos, Parte G: Journal of Aerospace Engineering. Palabras clave: adaptación de malla, predicción del error numérico, volúmenes finitos

Comparison of fatty acid profiles of edible meat, adipose tissues and muscles between cocks and capons

The effect of caponisation on fat composition by parts (wing, breast, thigh, and drumstick) and tissues (skin, subcutaneous adipose tissue, intermuscular adipose tissue and muscle) was examined in the present study and fatty acid profiles of abdominal fat and edible meat by parts and tissue components were determined. The sample was made up of twenty-eight castrated and twenty male Penedesenca Negra chicks reared under free-range conditions and slaughtered at 28 wk of age; the birds were castrated at four or eight weeks. Caponisation significantly increased (P < 0.01) the chemical fat content in all parts (16.31% to 37.98% in breast; 21.98% to 34.13% in wing; 21.09% to 49.57% in thigh; 14.33% to 24.82% in drumstick) and led to minor modifications in fat haracteristics, particularly in the thigh and the drumstick, where the unsaturated vs. saturated fatty acid ratio increased from 1.31 to 1.76 ( P < 0.01) and from 1.48 to 2.07 (P < 0.01), respectively. Delaying the age of castration from 4 to 8 weeks increased this ratio by 0.35 in the edible meat. Even though the profile of the abdominal fat is less saturated in capons, all changes occurring on fat quality after caponisation indicate that increased fatness after castration does not imply worse fat nutritional properties.

Molting dynamics and juvenile hormone titer profiles in the nymphal stages of a lower termite, Cryptotermes secundus (Kalotermitidae) - signatures of developmental plasticity

Termites are social cockroaches and this sociality is founded on a high plasticity during development. Three molting types (progressive, stationary and regressive molts) are fundamental to achieve plasticity during alate/sexual development, and they make termites a major challenge to any model on endocrine regulation in insect development. As the endocrine signatures underpinning this plasticity are barely understood, we studied the developmental dynamics and their underlying juvenile hormone OH) titers in a wood-dwelling termite. Cryptotermes secundus, which is characterized by an ancestral life style of living in dead wood and individuals being totipotent in development. The following general pattern elements could be identified during winged sexual development (i) regressive molts were accompanied by longer intermolt periods than other molting types, (ii) JH titers decreased gradually during the developmental transition from larva (immatures without wing buds), to nymph (immatures with wing buds), to winged adult, (iii) in all nymphal stages, the JH titer rose before the next molt and dropped thereafter within the first week, (iv) considerable variation in JH titers occurred in the midphase of the molting cycle of the 2nd and 3rd nymphal instar, inferring that this variation may reflect the underlying endocrine signature of each of the three molting types, (v) the 4th nymphal instar, the shortest of all, seems to be a switch point in development, as nymphs in this stage mainly developed progressively. When comparing these patterns with endocrine signatures seen in cockroaches, the developmental program of Cryprotermes can be interpreted as a co-option and repetitive use of hormonal dynamics of the post dorsal-closure phase of cockroach embryonic development. (C) 2012 Elsevier Ltd. All tights reserved.

Preliminary ichthyolith biostratigraphy and age profiles at DSDP Sites 91-595 and 91-596

Seventy meters of Cenozoic and Mesozoic pelagic clay cored at DSDP Sites 595 and 596 provide the basis for a preliminary analysis of ichthyolith biostratigraphy in the southwest Pacific. A most likely order of the more reliable ichthyolith events is compared with a synthesis of ichthyolith biostratigraphy in the North Pacific and with dated composite ranges. The resultant preliminary ichthyolith stratigraphy suggests that the Cenozoic is represented by the upper 20 m at Site 596 and 16 to 22 m at Site 595. Mixing of taxa precludes a clear recognition of the Cretaceous/Tertiary boundary at Site 595. The occurrence of 13 newly described subtypes is recorded in Mesozoic sediments at Sites 595 and 596. These new subtypes and previously described Mesozoic forms may be useful for recognizing Mesozoic subdivisions when their occurrences in sequences dated by other microfossils are investigated.

Two pollen profiles of the middle Pleistocene at Samerberg/Upper Bavaria, Germany

Previous pollen analytical studies on sediments from the pleistocene lake basin at Samerberg, situated on the northern edge of the Bavarian Alps (47°45' N, 12°12' E, 607 m a.s.l.) had been performed on samples taken from cores and exposures close to the southern shore of the former lake. After geoelectric and refraction-seismic measurements had shown that the lake basin had been much deeper in its northern part, another core was taken where maximum depth could be expected. The corer penetrated three moraines, two of them lying above pollen-bearing sediments, and one below them, and reached the hard rock (Kössener Kalk) at a depth of 93 m. Two forest phases could be identified by pollen analysis. The pollen record begins abruptly in a forest phase at the end of a spruce-dominated period when fir started to spread (DA 1, DA = pollen zone). Following this, Abies (fir) was the main tree species at Samerberg, Picea being second, and deciduous trees were almost non-existent. First box (Buxus) was of major importance in the fir forests (DA 2), but later on beech (Fagus) and wing-nut (Pterocarya) spread (DA 3). Finally this forest gave way to a spruce forest with pine (DA 4). The beginning and the end of this interglacial cycle are not recorded. Its vegetational development is different from the eemian one known from earlier studies at Samerberg. It is characterized by the occurrence of Abies together with Buxus, Pterocarya and Fagus. A similar association of woody species is known only from the Holsteinian age deposits in an area ranging from England to Poland, though at no other place these species were such important constituents of the vegetation as at Samerberg. Therefore zone 1 to 4 are attributed to the Holsteinian interglacial period. The younger forest phase, separated from the interglacial by a stadial with open vegetation (DA 5), seems to be completely represented, though its sediments are disturbed, apparently by sliding which caused repetition of same-age-sediments in the core (DA 7a, b, c) The vegetational development is simple. A juniper phase (DA 6) was followed by reforestation with spruce, accompanied by some fir (DA 7, 9). Finally pine became the dominant species (DA 9). The simple vegetational development of this younger forest phase does not allow a safe correlation with one of the known pre-eemian interstadials, but for stratigraphical reasons it can be related best to the Dömnitz-interglacial, which among others is also known as Wacken- or Holstein-II-interglacial. Possibly another phase of reforestation is indicated at the end of the following stadial (DA 10). But due to an erosional unconformity nothing than the rise of the juniper curve can be stated. It was only after this sequence of forest phases and periods with open vegetation that glaciers reached the Samerberg area again.

Cluster Analysis To Identify Elderly People's Profiles: A Healthcare Strategy Based On Frailty Characteristics.

The new social panorama resulting from aging of the Brazilian population is leading to significant transformations within healthcare. Through the cluster analysis strategy, it was sought to describe the specific care demands of the elderly population, using frailty components. Cross-sectional study based on reviewing medical records, conducted in the geriatric outpatient clinic, Hospital de Clínicas, Universidade Estadual de Campinas (Unicamp). Ninety-eight elderly users of this clinic were evaluated using cluster analysis and instruments for assessing their overall geriatric status and frailty characteristics. The variables that most strongly influenced the formation of clusters were age, functional capacities, cognitive capacity, presence of comorbidities and number of medications used. Three main groups of elderly people could be identified: one with good cognitive and functional performance but with high prevalence of comorbidities (mean age 77.9 years, cognitive impairment in 28.6% and mean of 7.4 comorbidities); a second with more advanced age, greater cognitive impairment and greater dependence (mean age 88.5 years old, cognitive impairment in 84.6% and mean of 7.1 comorbidities); and a third younger group with poor cognitive performance and greater number of comorbidities but functionally independent (mean age 78.5 years old, cognitive impairment in 89.6% and mean of 7.4 comorbidities). These data characterize the profile of this population and can be used as the basis for developing efficient strategies aimed at diminishing functional dependence, poor self-rated health and impaired quality of life.

Hydrometeor and latent heat profiles of tropical cyclones Conson, Ivan and Catarina using PR/TRMM data

ABSTRACT Microphysical and thermodynamical features of two tropical systems, namely Hurricane Ivan and Typhoon Conson, and one sub-tropical, Catarina, have been analyzed based on space-born radar PR measurements available on the TRMM satellite. The procedure to classify the reflectivity profiles followed the Heymsfield et al (2000) and Steiner et al (1995) methodologies. The water and ice content have been calculated using a relationship obtained with data of the surface SPOL radar and PR in Rondonia State in Brazil. The diabatic heating rate due to latent heat release has been estimated using the methodology developed by Tao et al (1990). A more detailed analysis has been performed for Hurricane Catarina, the first of its kind in South Atlantic. High water content mean value has been found in Conson and Ivan at low levels and close to their centers. Results indicate that hurricane Catarina was shallower than the other two systems, with less water and the water was concentrated closer to its center. The mean ice content in Catarina was about 0.05 g kg-1 while in Conson it was 0.06 g kg-1 and in Ivan 0.08 g kg-1. Conson and Ivan had water content up to 0.3 g kg-1 above the 0ºC layer, while Catarina had less than 0.15 g kg-1. The latent heat released by Catarina showed to be very similar to the other two systems, except in the regions closer to the center.

Wing diagnostic characters for Culex quinquefasciatus and Culex nigripalpus (Diptera, Culicidae)

Wing diagnostic characters for Culex quinquefasciatus and Culex nigripalpus (Diptera, Culicidae). Culex quinquefasciatus and Culex nigripalpus are mosquitoes of public health interest, which can occur sympatrically in urban and semi-urban localities. Morphological identification of these species may be difficult when specimens are not perfectly preserved. In order to suggest an alternative taxonomical diagnosis, wings of these species were comparatively characterized using geometric morphometrics. Both species could be distinguished by wing shape with accuracy rates ranging from 85-100%. Present results indicate that one can identify these species relying only on wing characters when traditional taxonomical characters are not visible.

Geometric morphometrics of the wing as a tool for assigning genetic lineages and geographic origin to Melipona beecheii (Hymenoptera: Meliponini)

The stingless bee Melipona beecheii presents great variability and is considered a complex of species. In order to better understand this species complex, we need to evaluate its diversity and develop methods that allow geographic traceability of the populations. Here we present a fast, efficient, and inexpensive means to accomplish this using geometric morphometrics of wings. We collected samples from Mexico, Guatemala, El Salvador, Nicaragua, and Costa Rica and we were able to correctly assign 87.1% of the colonies to their sampling sites and 92.4% to their haplotype. We propose that geometric morphometrics of the wing could be used as a first step analysis leaving the more expensive molecular analysis only to doubtful cases.

Morphometric and genetic changes in a population of Apis mellifera after 34 years of Africanization

Though the replacement of European bees by Africanized honey bees in tropical America has attracted considerable attention, little is known about the temporal changes in morphological and genetic characteristics in these bee populations. We examined the changes in the morphometric and genetic profiles of an Africanized honey bee population collected near where the original African swarms escaped, after 34 years of Africanization. Workers from colonies sampled in 1968 and in 2002 were morphometrically analyzed using relative warps analysis and an Automatic Bee Identification System (ABIS). All the colonies had their mitochondrial DNA identified. The subspecies that mixed to form the Africanized honey bees were used as a comparison for the morphometric analysis. The two morphometric approaches showed great similarity of Africanized bees with the African subspecies, Apis mellifera scutellata, corroborating with other markers. We also found the population of 1968 to have the pattern of wing venation to be more similar to A. m. scutellata than the current population. The mitochondrial DNA of European origin, which was very common in the 1968 population, was not found in the current population, indicating selective pressure replacing the European with the African genome in this tropical region. Both morphometric methodologies were very effective in discriminating the A. mellifera groups; the non-linear analysis of ABIS was the most successful in identifying the bees, with more than 94% correct classifications.

Gender identification of five genera of stingless bees (Apidae, Meliponini) based on wing morphology

Currently, the identification of pollinators is a critical necessity of conservation programs. After it was found that features extracted from patterns of wing venation are sufficient to discriminate among insect species, various studies have focused on this structure. We examined wing venation patterns of males and workers of five stingless bee species in order to determine if there are differences between sexes and if these differences are greater within than between species. Geometric morphometric analyses were made of the forewings of males and workers of Nannotrigona testaceicornis, Melipona quadrifasciata, Frieseomelitta varia, and Scaptotrigona aff. depilis and Plebeia remota. The patterns of males and workers from the same species were more similar than the patterns of individuals of the same sex from different species, and the patterns of both males and workers, when analyzed alone, were sufficiently different to distinguish among these five species. This demonstrates that we can use this kind of analysis for the identification of stingless bee species and that the sex of the individual does not impede identification. Computer-assisted morphometric analysis of bee wing images can be a useful tool for biodiversity studies and conservation programs.

MicroRNA transcriptome profiles during swine skeletal muscle development

Background: MicroRNA (miR) are a class of small RNAs that regulate gene expression by inhibiting translation of protein encoding transcripts. To evaluate the role of miR in skeletal muscle of swine, global microRNA abundance was measured at specific developmental stages including proliferating satellite cells, three stages of fetal growth, day-old neonate, and the adult. Results: Twelve potential novel miR were detected that did not match previously reported sequences. In addition, a number of miR previously reported to be expressed in mammalian muscle were detected, having a variety of abundance patterns through muscle development. Muscle-specific miR-206 was nearly absent in proliferating satellite cells in culture, but was the highest abundant miR at other time points evaluated. In addition, miR-1 was moderately abundant throughout developmental stages with highest abundance in the adult. In contrast, miR-133 was moderately abundant in adult muscle and either not detectable or lowly abundant throughout fetal and neonate development. Changes in abundance of ubiquitously expressed miR were also observed. MiR-432 abundance was highest at the earliest stage of fetal development tested (60 day-old fetus) and decreased throughout development to the adult. Conversely, miR-24 and miR-27 exhibited greatest abundance in proliferating satellite cells and the adult, while abundance of miR-368, miR-376, and miR-423-5p was greatest in the neonate. Conclusion: These data present a complete set of transcriptome profiles to evaluate miR abundance at specific stages of skeletal muscle growth in swine. Identification of these miR provides an initial group of miR that may play a vital role in muscle development and growth.

Gene Expression Profiles in Radiation Workers Occupationally Exposed to Ionizing Radiation

Ionizing radiation OR) imposes risks to human health and the environment. IR at low doses and low (lose rates has the potency to initiate carcinogenesis. Genotoxic environmental agents such as IR trigger a cascade of signal transduction pathways for cellular protection. In this study, using cDNA microarray technique, we monitored the gene expression profiles in lymphocytes derived from radiation-ex posed individuals (radiation workers). Physical dosimetry records on these patients indicated that the absorbed dose ranged from 0.696 to 39.088 mSv. Gene expression analysis revealed statistically significant transcriptional changes in a total of 78 genes (21 up-regulated and 57 clown-regulated) involved in several biological processes such as ubiquitin cycle (UHRF2 and PIAS1), DNA repair (LIG3, XPA, ERCC5, RAD52, DCLRE1C), cell cycle regulation/proliferation (RHOA, CABLES2, TGFB2, IL16), and stress response (GSTP1, PPP2R5A, DUSP22). Some of the genes that showed altered expression profiles in this study call be used as biomarkers for monitoring the chronic low level exposure in humans. Additionally, alterations in gene expression patterns observed in chronically exposed radiation workers reinforces the need for defining the effective radiation dose that causes immediate genetic damage as well as the long-term effects on genomic instability, including cancer.

Suppression subtractive hybridization profiles of radial growth phase and metastatic melanoma cell lines reveal novel potential targets

Background: Melanoma progression occurs through three major stages: radial growth phase (RGP), confined to the epidermis; vertical growth phase (VGP), when the tumor has invaded into the dermis; and metastasis. In this work, we used suppression subtractive hybridization (SSH) to investigate the molecular signature of melanoma progression, by comparing a group of metastatic cell lines with an RGP-like cell line showing characteristics of early neoplastic lesions including expression of the metastasis suppressor KISS1, lack of alpha v beta 3-integrin and low levels of RHOC. Methods: Two subtracted cDNA collections were obtained, one (RGP library) by subtracting the RGP cell line (WM1552C) cDNA from a cDNA pool from four metastatic cell lines (WM9, WM852, 1205Lu and WM1617), and the other (Met library) by the reverse subtraction. Clones were sequenced and annotated, and expression validation was done by Northern blot and RT-PCR. Gene Ontology annotation and searches in large-scale melanoma expression studies were done for the genes identified. Results: We identified 367 clones from the RGP library and 386 from the Met library, of which 351 and 368, respectively, match human mRNA sequences, representing 288 and 217 annotated genes. We confirmed the differential expression of all genes selected for validation. In the Met library, we found an enrichment of genes in the growth factors/receptor, adhesion and motility categories whereas in the RGP library, enriched categories were nucleotide biosynthesis, DNA packing/repair, and macromolecular/vesicular trafficking. Interestingly, 19% of the genes from the RGP library map to chromosome 1 against 4% of the ones from Met library. Conclusion: This study identifies two populations of genes differentially expressed between melanoma cell lines from two tumor stages and suggests that these sets of genes represent profiles of less aggressive versus metastatic melanomas. A search for expression profiles of melanoma in available expression study databases allowed us to point to a great potential of involvement in tumor progression for several of the genes identified here. A few sequences obtained here may also contribute to extend annotated mRNAs or to the identification of novel transcripts.