Application of Volterra series in nonlinear mechanical system identification and in structural health monitoring problems

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

DNA damage defines sites of recurrent chromosomal translocations in B lymphocytes

Recurrent chromosomal translocations underlie both haematopoietic and solid tumours. Their origin has been ascribed to selection of random rearrangements, targeted DNA damage, or frequent nuclear interactions between translocation partners; however, the relative contribution of each of these elements has not been measured directly or on a large scale. Here we examine the role of nuclear architecture and frequency of DNA damage in the genesis of chromosomal translocations by measuring these parameters simultaneously in cultured mouse B lymphocytes. In the absence of recurrent DNA damage, translocations between Igh or Myc and all other genes are directly related to their contact frequency. Conversely, translocations associated with recurrent site-directed DNA damage are proportional to the rate of DNA break formation, as measured by replication protein A accumulation at the site of damage. Thus, non-targeted rearrangements reflect nuclear organization whereas DNA break formation governs the location and frequency of recurrent translocations, including those driving B-cell malignancies.

Optimization of anti-Trypanosoma cruzi oxadiazoles leads to identification of compounds with efficacy in infected mice

We recently showed that oxadiazoles have anti-Trypanosoma cruzi activity at micromolar concentrations. These compounds are easy to synthesize and show a number of clear and interpretable structure-activity relationships (SAR), features that make them attractive to pursue potency enhancement. We present here the structural design, synthesis, and anti-T. cruzi evaluation of new oxadiazoles denoted 5a-h and 6a-h. The design of these compounds was based on a previous model of computational docking of oxadiazoles on the T. cruzi protease cruzain. We tested the ability of these compounds to inhibit catalytic activity of cruzain, but we found no correlation between the enzyme inhibition and the antiparasitic activity of the compounds. However, we found reliable SAR data when we tested these compounds against the whole parasite. While none of these oxadiazoles showed toxicity for mammalian cells, oxadiazoles 6c (fluorine), 6d (chlorine), and 6e (bromine) reduced epimastigote proliferation and were cidal for trypomastigotes of T. cruzi Y strain. Oxadiazoles 6c and 6d have IC50 of 9.5 +/- 2.8 and 3.5 +/- 1.8 mu M for trypomastigotes, while Benznidazole, which is the currently used drug for Chagas disease treatment, showed an IC50 of 11.3 +/- 2.8 mu M. Compounds 6c and 6d impair trypomastigote development and invasion in macrophages, and also induce ultrastructural alterations in trypomastigotes. Finally, compound 6d given orally at 50 mg/kg substantially reduces the parasitemia in T. cruzi-infected BALB/c mice. Our drug design resulted in potency enhancement of oxadiazoles as anti-Chagas disease agents, and culminated with the identification of oxadiazole 6d, a trypanosomicidal compound in an animal model of infection. (C) 2012 Elsevier Ltd. All rights reserved.

Substructuring approache in state space models for dynamic system parameters identification

In the recent decade, the request for structural health monitoring expertise increased exponentially in the United States. The aging issues that most of the transportation structures are experiencing can put in serious jeopardy the economic system of a region as well as of a country. At the same time, the monitoring of structures is a central topic of discussion in Europe, where the preservation of historical buildings has been addressed over the last four centuries. More recently, various concerns arose about security performance of civil structures after tragic events such the 9/11 or the 2011 Japan earthquake: engineers looks for a design able to resist exceptional loadings due to earthquakes, hurricanes and terrorist attacks. After events of such a kind, the assessment of the remaining life of the structure is at least as important as the initial performance design. Consequently, it appears very clear that the introduction of reliable and accessible damage assessment techniques is crucial for the localization of issues and for a correct and immediate rehabilitation. The System Identification is a branch of the more general Control Theory. In Civil Engineering, this field addresses the techniques needed to find mechanical characteristics as the stiffness or the mass starting from the signals captured by sensors. The objective of the Dynamic Structural Identification (DSI) is to define, starting from experimental measurements, the modal fundamental parameters of a generic structure in order to characterize, via a mathematical model, the dynamic behavior. The knowledge of these parameters is helpful in the Model Updating procedure, that permits to define corrected theoretical models through experimental validation. The main aim of this technique is to minimize the differences between the theoretical model results and in situ measurements of dynamic data. Therefore, the new model becomes a very effective control practice when it comes to rehabilitation of structures or damage assessment. The instrumentation of a whole structure is an unfeasible procedure sometimes because of the high cost involved or, sometimes, because it’s not possible to physically reach each point of the structure. Therefore, numerous scholars have been trying to address this problem. In general two are the main involved methods. Since the limited number of sensors, in a first case, it’s possible to gather time histories only for some locations, then to move the instruments to another location and replay the procedure. Otherwise, if the number of sensors is enough and the structure does not present a complicate geometry, it’s usually sufficient to detect only the principal first modes. This two problems are well presented in the works of Balsamo [1] for the application to a simple system and Jun [2] for the analysis of system with a limited number of sensors. Once the system identification has been carried, it is possible to access the actual system characteristics. A frequent practice is to create an updated FEM model and assess whether the structure fulfills or not the requested functions. Once again the objective of this work is to present a general methodology to analyze big structure using a limited number of instrumentation and at the same time, obtaining the most information about an identified structure without recalling methodologies of difficult interpretation. A general framework of the state space identification procedure via OKID/ERA algorithm is developed and implemented in Matlab. Then, some simple examples are proposed to highlight the principal characteristics and advantage of this methodology. A new algebraic manipulation for a prolific use of substructuring results is developed and implemented.

Identification of noninvasive biomarkers for alcohol-induced liver disease using urinary metabolomics and the Ppara-null mouse

Alcohol-induced liver disease (ALD) is a leading cause of nonaccident-related deaths in the United States. Although liver damage caused by ALD is reversible when discovered at the earlier stages, current risk assessment tools are relatively nonspecific. Identification of an early specific signature of ALD would aid in therapeutic intervention and recovery. In this study, the metabolic changes associated with ALD were examined using alcohol-fed male Ppara-null mouse as a model of ALD. Principal components analysis of the mass spectrometry-based urinary metabolic profile showed that alcohol-treated wild-type and Ppara-null mice could be distinguished from control animals without information on history of alcohol consumption. The urinary excretion of ethyl-sulfate, ethyl-beta-d-glucuronide, 4-hydroxyphenylacetic acid, and 4-hydroxyphenylacetic acid sulfate was elevated and that of the 2-hydroxyphenylacetic acid, adipic acid, and pimelic acid was depleted during alcohol treatment in both wild-type and the Ppara-null mice albeit to different extents. However, indole-3-lactic acid was exclusively elevated by alcohol exposure in Ppara-null mice. The elevation of indole-3-lactic acid is mechanistically related to the molecular events associated with development of ALD in alcohol-treated Ppara-null mice. This study demonstrated the ability of a metabolomics approach to identify early, noninvasive biomarkers of ALD pathogenesis in Ppara-null mouse model.

Identification of the p53 family-responsive element in the promoter region of the tumor suppressor gene hypermethylated in cancer 1

The tumor suppressor gene hypermethylated in cancer 1 (HIC1), located on human chromosome 17p13.3, is frequently silenced in cancer by epigenetic mechanisms. Hypermethylated in cancer 1 belongs to the bric à brac/poxviruses and zinc-finger family of transcription factors and acts by repressing target gene expression. It has been shown that enforced p53 expression leads to increased HIC1 mRNA, and recent data suggest that p53 and Hic1 cooperate in tumorigenesis. In order to elucidate the regulation of HIC1 expression, we have analysed the HIC1 promoter region for p53-dependent induction of gene expression. Using progressively truncated luciferase reporter gene constructs, we have identified a p53-responsive element (PRE) 500 bp upstream of the TATA-box containing promoter P0 of HIC1, which is sequence specifically bound by p53 in vitro as assessed by electrophoretic mobility shift assays. We demonstrate that this HIC1 p53-responsive element (HIC1.PRE) is necessary and sufficient to mediate induction of transcription by p53. This result is supported by the observation that abolishing endogenous wild-type p53 function prevents HIC1 mRNA induction in response to UV-induced DNA damage. Other members of the p53 family, notably TAp73beta and DeltaNp63alpha, can also act through this HIC1.PRE to induce transcription of HIC1, and finally, hypermethylation of the HIC1 promoter attenuates inducibility by p53.

Characterization of retinal damage in the episcleral vein cauterization rat glaucoma model

Episcleral vein cauterization (EVC) is used in rats to generate a glaucoma model with high intraocular pressure (IOP). The long-term retinal damage in this glaucoma model, however, has not been accurately quantified. We report the location and amount of retinal ganglion cell (RGC) damage caused by (EVC) induced IOP elevation in two rat strains. IOP was raised in one eye of Wistar (N = 5) and Brown-Norway(B-N)(N = 7) rats by EVC and monitored monthly until IOP in contralateral eyes equalized at 5 months post-surgery. Animals were maintained for 3.5-4.5 additional months. B-N rats (N = 7) that had no EVC served as controls for this strain. Scotopic flash ERGs were recorded at baseline and just prior to euthanasia. Automated counts of all retrogradely labeled RGCs in retinal flat-mounts were determined and compared between contralateral eyes. RGC density maps were constructed and RGC size distribution was determined. Oscillatory potentials in the group of eyes which had elevated IOP were decreased at the time of euthanasia, when IOP had returned to normal. The group of normal B-N rats had similar RGC counts between contralateral eyes. In the experimental group the mean number of RGCs was not significantly different between control and experimental eyes, but 1 of 5 Wistar and 2 of 7 B-N experimental eyes had at least 30% fewer RGCs than contralateral control eyes. Total retinal area in B-N experimental eyes was higher compared to contralateral eyes. Cumulative IOP exposure of the experimental eyes was modestly correlated with RGC loss while oscillatory potentials appeared to be inversely related to RGC loss. In retinas with extensive (> 30% RGC loss) but not complete damage, smaller cells were preserved better than larger ones. The above results indicate that RGC loss in both Wistar and B-N strains is variable after a prolonged elevation of IOP via EVC. Such variability despite equivalent IOP levels and ERG abnormalities, suggests unknown factors that can protect IOP-stressed RGCs. Identification and enhancement of such factors could prove useful for glaucoma therapy.

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain D(κ) , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.

The influence of yield surface shape and damage in the depth-dependent response of bone tissue to nanoindentation using spherical and Berkovich indenters

Prevention and treatment of osteoporosis rely on understanding of the micromechanical behaviour of bone and its influence on fracture toughness and cell-mediated adaptation processes. Postyield properties may be assessed by nonlinear finite element simulations of nanoindentation using elastoplastic and damage models. This computational study aims at determining the influence of yield surface shape and damage on the depth-dependent response of bone to nanoindentation using spherical and conical tips. Yield surface shape and damage were shown to have a major impact on the indentation curves. Their influence on indentation modulus, hardness, their ratio as well as the elastic-to-total work ratio is well described by multilinear regressions for both tip shapes. For conical tips, indentation depth was not statistically significant (p<0.0001). For spherical tips, damage was not a significant parameter (p<0.0001). The gained knowledge can be used for developing an inverse method for identification of postelastic properties of bone from nanoindentation.

Identification of cadmium-induced mutations in a mammalian cell line

Epidemiological studies have shown cadmium to induce cancer in humans, while experimental studies have proven this metal to be a potent tumor inducer in animals. However, cadmium appears nonmutagenic in most prokaryotic and eukaryotic mutagenesis assays. In this study, we present the identification of mutations in normal rat kidney cells infected with the mutant MuSVts110 retrovirus (6m2 cells) as a result of treatment with cadmium chloride. The detection of these mutations was facilitated by the use of a novel mutagenesis assay established in this laboratory. The 6m2 reversion assay is a positive selection system based on the conditional expression of the MuSVts110 v-mos gene. In MuSVts110 the gag and mos genes are fused out of frame, thus the translation of the v-mos sequence requires a frameshift in the genomic RNA. In 6m2 cells this frameshift is accomplished by the temperature-dependent splicing of the primary MuSVts110 transcript. Splicing of MuSVts110, which is mediated by cis-acting sequences, occurs when 6m2 cells are grown at 33$\sp\circ$C and below, but not at 39$\sp\circ$C. Therefore, 6m2 cells appear transformed at low growth temperatures, but take on a morphologically normal appearance when grown at high temperatures. The treatment of 6m2 cells with cadmium chloride resulted in the outgrowth of a number of cells that reverted to the transformed state at high growth temperatures. Analysis of the viral proteins expressed in these cadmium-induced 6m2 revertants suggested that they contained mutations in their MuSVts110 DNA. Sequencing of the viral DNA from three revertants that constitutively expressed the P85$\sp{gag{-}mos}$ transforming protein revealed five different mutations. The Cd-B2 revertant contained three of those mutations: an A-to-G transition 48 bases downstream of the MuSVts110 3$\sp\prime$ splice site, plus a G-to-T and an A-to-T transversion 84 and 100 bases downstream of the 5$\sp\prime$ splice site, respectively. The Cd-15-5 revertant also contained a point mutation, a T-to-C transition 46 bases downstream of the 5$\sp\prime$ splice site, while Cd-10-5 contained a three base deletion of MuSVts110 11 bases upstream of the 3$\sp\prime$ splice site. A fourth revertant, Cd-10, expressed a P100$\sp{gag{-}mos}$ transforming protein, and was found to have a two base deletion. This deletion accomplished the frameshift necessary for v-mos expression, but did not alter MuSVts110 RNA splicing and the expression of p85$\sp{gag{-}mos}.$ Lastly, sequencing of the MuSVts110 DNA from three spontaneous revertants revealed the same G to T transversion in each one. This was the same mutation that was found in the Cd-B2 revertant. These findings provide the first example of mutations resulting from exposure to cadmium and suggest, by the difference in each mutation, the complexity of the mechanism utilized by cadmium to induce DNA damage. ^

Identification of the signal pathways modulated by expression of p210$\rm\sp{bcr-abl}$ and the inhibition of p210$\rm\sp{bcr-abl}$ transforming activity by polypeptides interacting with p210$\rm\sp{bcr-abl}$ in murine myeloid 32D cells

The Bcr-Abl fusion oncogene which resulted from a balanced reciprocal translocation between chromosome 9 and 22, t(9;22)(q11, q34), encodes a 210 KD elevated tyrosine specific protein kinase that is found in more than 95 percent of chronic myelogenous leukemia patients (CML). Increase of level of phosphorylation of tyrosine is observed on cell cycle regulatory proteins in cells overexpressing the Bcr-Abl oncogene, which activates multiple signaling pathways. In addition, distinct signals are required for transforming susceptible fibroblast and hematopoietic cells, and the minimal signals essential for transforming hematopoietic cells are yet to be defined. In the present study, we first established a tetracycline repressible p210$\rm\sp{bcr-abl}$ expression system in a murine myeloid cell line 32D c13, which depends on IL3 to grow in the presence of tetracycline and proliferate independent of IL3 in the absence of tetracycline. Interestingly, one of these sublines does not form tumors in athymic nude mice suggesting that these cells may not be completely transformed. These cells also exhibit a dose-dependent growth and expression of p210$\rm\sp{bcr-abl}$ at varying concentrations of tetracycline in the culture. However, p210$\rm\sp{bcr-abl}$ rescues IL3 deprivation induced apoptosis in a non-dose dependent fashion. DNA genotoxic damage induced by gamma-irradiation activates c-Abl tyrosine kinase, the cellular homologue of p210$\rm\sp{bcr-abl},$ and leads to activation of p38 MAP kinase in the cells. However, in the presence of p210$\rm\sp{bcr-abl}$ the irradiation failed to activate the p38 MAP kinase as examined by an antibody against phosphorylated p38 MAP kinase. Similarly, an altered tyrosine phosphorylation of the JAK1-STAT1 pathways was identified in cells constitutively overexpressing p210$\rm\sp{bcr-abl}.$ This may provided a molecular mechanism for altered therapeutic response of CML patients to IFN-$\alpha.$^ Bcr-Abl oncoprotein has multiple functional domains which have been identified by the work of others. The Bcr tetramerization domain, which may function to stabilize the association of the Bcr-Abl with actin filaments in p210$\rm\sp{bcr-abl}$ susceptible cells, are essential for transforming both fibroblast and hematopoietic cells. We designed a transcription unit encoding first 160 amino acids polypeptide of Bcr protein to test if this polypeptide can inhibit the transforming activity of the p210$\rm\sp{bcr-abl}$ oncoprotein in the 32D c13 cells. When this vector was transfected transiently along with the p210$\rm\sp{bcr-abl}$ expression vector, it can block the transforming activity of p210$\rm\sp{bcr-abl}.$ On the other hand, the retinoblastoma tumor suppressor protein (Rb), a naturally occurring negative regulator of the c-Abl kinase, the cellular homologue of Bcr-Abl oncoprotein, binds to and inhibits the c-Abl kinase in a cell cycle dependent manner. A polypeptide obtained from the carboxyl terminal end of the retinoblastoma tumor suppressor protein, in which the nuclear localization signal was mutated, was used to inhibit the kinase activity of the p210$\rm\sp{bcr-abl}$ in the cytoplasm. This polypeptide, called Rb MC-box, and its wild type form, Rb C-box, when overexpressed in the 32D cells are mainly localized in the cytoplasm. Cotransfection of a plasmid transcription unit coding for this polypeptide and the gene for the p210$\rm\sp{bcr-abl}$ resulted in reduced plating efficiency of p210$\rm\sp{bcr-abl}$ transfected IL3 independent 32D cells. Together, these results may lead to a molecular approach to therapy of CML and an in vitro assay system to identify new targets to which an inhibitory polypeptide transcription unit may be directed. ^

Characterization of the potential role for RNA-binding protein FUS/TLS in DNA damage response: A quantitative proteomic approach

FUS/TLS (fused in sarcoma/translocated in liposarcoma) is a ubiquitously expressed RNA-binding protein of the hnRNP family, that has been discovered as fused to transcription factors, through chromosomal translocations, in several human sarcomas and found in protein aggregates in neurons of patients with an inherited form of Amyotrophic Lateral Sclerosis (ALS) [1]. To date, FUS/TLS has been implicated in a variety of cellular processes such as gene expression control, transcriptional regulation, pre-mRNA splicing and miRNA processing [2]. In addition, some evidences link FUS/TLS to genome stability control and DNA damage response. In fact, mice lacking FUS/TLS are hypersensitive to ionizing radiation (IR) and show high levels of chromosome instability and in response to double-strand breaks, FUS/TLS gets phosphorylated by the protein kinase ATM [3,4,5]. Furthermore, the inducible depletion of FUS/TLS in a neuroblastoma cell line (SH-SY5Y FUS/TLS TET-off iKD) subjected to genotoxic stress (IR) resulted in an increased phosphorylation of γH2AX respect to control cells, suggesting an higher activation of the DNA damage response. The study aims to investigate the specific role of FUS/TLS in DNA damage response through the characterization of the proteomic profile of SH-SY5Y FUS/TLS iKD cells subjected to DNA damage stress, by mass spectrometry-based quantitative proteomics (e.g. SILAC). Preliminary results of mass spectrometric identification of FUS/TLS interacting proteins in HEK293 cells, expressing a recombinant flag-tagged FUS/TLS protein, highlighted the interactions with several proteins involved in DNA damage response, such as DNA-PK, XRCC-5/-6, and ERCC-6, raising the possibilities that FUS/TLS is involved in this pathway, even thou its exact role still need to be addressed.

Characterization of the potential role for RNA-binding protein FUS in DNA damage response: A quantitative proteomic approach

FUS/TLS (fused in sarcoma/translocated in liposarcoma) is a ubiquitously expressed protein of the hnRNP family, that has been discovered as fused to transcription factors in several human sarcomas and found in protein aggregates in neurons of patients with an inherited form of Amyotrophic Lateral Sclerosis [Vance C. et al., 2009]. FUS is a 53 kDa nuclear protein that contains structural domains, such as a RNA Recognition Motif (RRM) and a zinc finger motif, that give to FUS the ability to bind to both RNA and DNA sequences. It has been implicated in a variety of cellular processes, such as pre-mRNA splicing, miRNA processing, gene expression control and transcriptional regulation [Fiesel FC. and Kahle PJ., 2011]. Moreover, some evidences link FUS to genome stability control and DNA damage response: mice lacking FUS are hypersensitive to ionizing radiation (IR) and show high levels of chromosome instability and, in response to double-strand breaks, FUS is phosphorylated by the protein kinase ATM [Kuroda M. et al., 2000; Hicks GG. et al., 2000; Gardiner M. et al., 2008]. Furthermore, preliminary results of mass spectrometric identification of FUS interacting proteins in HEK293 cells, expressing a recombinant flag-tagged FUS protein, highlighted the interactions with proteins involved in DNA damage response, such as DNA-PK, XRCC-5/-6, and ERCC-6, raising the possibilities that FUS is involved in this pathway, even though its role still needs to be clarified. This study aims to investigate the biological roles of FUS in human cells and in particular the putative role in DNA damage response through the characterization of the proteomic profile of the neuroblastoma cell line SH-SY5Y upon FUS inducible depletion, by a quantitative proteomic approach. The SH-SY5Y cell line that will be used in this study expresses, in presence of tetracycline, a shRNA that targets FUS mRNA, leading to FUS protein depletion (SH-SY5Y FUS iKD cells). To quantify changes in proteins expression levels a SILAC strategy (Stable Isotope Labeling by Amino acids in Cell culture) will be conducted on SH-SY5Y FUS iKD cells and a control SH-SY5Y cell line (that expresses a mock shRNA) and the relative changes in proteins levels will be evaluated after five and seven days upon FUS depletion, by nanoliquid chromatography coupled to tandem mass spectrometry (nLC-MS/MS) and bioinformatics analysis. Preliminary experiments demonstrated that the SH-SY5Y FUS iKD cells, when subjected to genotoxic stress (high dose of IR), upon inducible depletion of FUS, showed a increased phosphorylation of gH2AX with respect to control cells, suggesting an higher activation of the DNA damage response.

Cerebellar Pathology Does Not Impair Performance on Identification or Categorization Tasks

In comparison to the basal ganglia, prefrontal cortex, and medial temporal lobes, the cerebellum has been absent from recent research on the neural substrates of categorization and identification, two prominent tasks in the learning and memory literature. To investigate the contribution of the cerebellum to these tasks, we tested patients with cerebellar pathology (seven with bilateral degeneration, six with unilateral lesions, and two with midline damage) on rule-based and information-integration categorization tasks and an identification task. In rule-based tasks, it is assumed that participants learn the categories through an explicit reasoning process. In information-integration tasks, optimal performance requires the integration of information from multiple stimulus dimensions, and participants are typically unaware of the decision strategy. The identification task, in contrast, required participants to learn arbitrary, color-word associations. The cerebellar patients performed similar to matched controls on all three tasks and performance did not vary with the extent of cerebellar pathology. Although the interpretation of these null results requires caution, these data contribute to the current debate on cerebellar contributions to cognition by providing boundary conditions on understanding the neural substrates of categorization and identification, and help define the functional domain of the cerebellum in learning and memory.

Identification of genetic variants that influence atherosclerosis in young persons

Atherosclerosis is widely accepted as a complex genetic phenotype and is the usual cause of cardiovascular disease, the world’s leading killer. Genetic factors have been proven to be important risk contributors for atherosclerosis and much work has been done to identify promising candidates that might play a role in the development of atherosclerosis. It is well known that many independent replications are needed to unequivocally establish a valid genotype-phenotype association across different populations before the findings are extended to clinical settings and to the expensive follow-up studies designed to identify causal genetic variants. Aiming to replicate the association with atherosclerosis in the Pathobiological Determinants of Atherosclerosis in Youth (PDAY) study, we assessed the relationship of 32 atherosclerosis candidate SNPs to atherosclerosis in the PDAY cohort, consisting of AA and EA young people aged 15-34 years who died of non-medical causes. Two association studies, a whole sample study and a 1:1 matched case control study were performed by use of multiple linear regression and logistic regression analyses, respectively. For the whole sample association study, 32 SNPs among 2,650 individuals (1,369 AA and 1,281 EA) were tested for the association with six early atherosclerosis phenotypes: abdominal aorta fatty streaks, abdominal aorta raised lesions, right coronary artery fatty streaks, right coronary artery raised lesions, thoracic aorta fatty streaks, and thoracic aorta raised lesions. For the matched case-control association study, 337 case-control paired samples were included; cases were chosen with the highest total raised lesion scores from the studied population, while controls were randomly selected from individuals that had no raised lesions and matched to cases by age, gender and race. Sixteen SNPs in 13 genes were found to be significantly associated with atherosclerosis in at least one of the PDAY association studies. Among these 16 findings: eight SNPs (rs9579646, rs6053733, rs3849150, rs10499903, rs2148079, rs5073691, rs10116277, and rs17228212) successfully replicated previous results, six SNPs (rs17222814, rs10811661, rs7028570, rs7291467, rs16996148 and rs10401969) were reported as new findings exclusive to our study, the last two of the 16 SNPs, rs501120 and rs6922269, showed either intriguing or conflicting result. SNP rs17222814 in ALOX5AP and SNP rs3849150 in LRRC18 were consistently associated with atherosclerosis in both prior and the two PDAY association studies. SNP rs3849150 was also identified to be highly correlated with a non-synonymous coding SNP, rs17772611, which may damage the protein (polyphen score = 0.996), suggesting that SNP rs17772611 may be the causal functional variant.^ In conclusion, our study added more support for the association of these candidate genes with atherosclerosis. SNPs rs3849150 and rs17772611 of LRRC18, as well as SNP rs17222814 of ALOX5AP, were the most significant findings from our study, and may be ranked among the best for further study.^