BRAIN ACTIVATION AND CONNECTIVITY IN NON-DISABLED MULTIPLE SCLEROSIS PATIENTS

Multiple sclerosis (MS) is the most common demyelinating disease affecting the central nervous system. There is no cure for MS and current therapies have limited efficacy. While the majority of individuals with MS develop significant clinical disability, a subset experiences a disease course with minimal impairment even in the presence of significant apparent tissue damage on magnetic resonance imaging (MRI). The current studies combined functional MRI and diffusion tensor imaging (DTI) to elucidate brain mechanisms associated with lack of clinical disability in patients with MS. Recent evidence has implicated cortical reorganization as a mechanism to limit the clinical manifestation of the disease. Functional MRI was used to test the hypothesis that non-disabled MS patients (Expanded Disability Status Scale ≤ 1.5) show increased recruitment of cognitive control regions (dorsolateral prefrontal and anterior cingulate cortex) while performing sensory, motor and cognitive tasks. Compared to matched healthy controls, patients increased activation of cognitive control brain regions when performing non-dominant hand movements and the 2-back working memory task. Using dynamic causal modeling, we tested whether increased cognitive control recruitment is associated with alterations in connectivity in the working memory functional network. Patients exhibited similar network connectivity to that of control subjects when performing working memory tasks. We subsequently investigated the integrity of major white matter tracts to assess structural connectivity and its relation to activation and functional integration of the cognitive control system. Patients showed substantial alterations in callosal, inferior and posterior white matter tracts and less pronounced involvement of the corticospinal tracts and superior longitudinal fasciculi (SLF). Decreased structural integrity within the right SLF in patients was associated with decreased performance, and decreased activation and connectivity of the cognitive control system when performing working memory tasks. These studies suggest that patient with MS without clinical disability increase cognitive control system recruitment across functional domains and rely on preserved functional and structural connectivity of brain regions associated with this network. Moreover, the current studies show the usefulness of combining brain activation data from functional MRI and structural connectivity data from DTI to improve our understanding of brain adaptation mechanisms to neurological disease.

3D MPRAGE improves classification of cortical lesions in multiple sclerosis.

BACKGROUND: Gray matter lesions are known to be common in multiple sclerosis (MS) and are suspected to play an important role in disease progression and clinical disability. A combination of magnetic resonance imaging (MRI) techniques, double-inversion recovery (DIR), and phase-sensitive inversion recovery (PSIR), has been used for detection and classification of cortical lesions. This study shows that high-resolution three-dimensional (3D) magnetization-prepared rapid acquisition with gradient echo (MPRAGE) improves the classification of cortical lesions by allowing more accurate anatomic localization of lesion morphology. METHODS: 11 patients with MS with previously identified cortical lesions were scanned using DIR, PSIR, and 3D MPRAGE. Lesions were identified on DIR and PSIR and classified as purely intracortical or mixed. MPRAGE images were then examined, and lesions were re-classified based on the new information. RESULTS: The high signal-to-noise ratio, fine anatomic detail, and clear gray-white matter tissue contrast seen in the MPRAGE images provided superior delineation of lesion borders and surrounding gray-white matter junction, improving classification accuracy. 119 lesions were identified as either intracortical or mixed on DIR/PSIR. In 89 cases, MPRAGE confirmed the classification by DIR/PSIR. In 30 cases, MPRAGE overturned the original classification. CONCLUSION: Improved classification of cortical lesions was realized by inclusion of high-spatial resolution 3D MPRAGE. This sequence provides unique detail on lesion morphology that is necessary for accurate classification.

Deep gray matter atrophy in multiple sclerosis: a tensor based morphometry.

Tensor based morphometry (TBM) was applied to determine the atrophy of deep gray matter (DGM) structures in 88 relapsing multiple sclerosis (MS) patients. For group analysis of atrophy, an unbiased atlas was constructed from 20 normal brains. The MS brain images were co-registered with the unbiased atlas using a symmetric inverse consistent nonlinear registration. These studies demonstrate significant atrophy of thalamus, caudate nucleus, and putamen even at a modest clinical disability, as assessed by the expanded disability status score (EDSS). A significant correlation between atrophy and EDSS was observed for different DGM structures: (thalamus: r=-0.51, p=3.85 x 10(-7); caudate nucleus: r=-0.43, p=2.35 x 10(-5); putamen: r=-0.36, p=6.12 x 10(-6)). Atrophy of these structures also correlated with 1) T2 hyperintense lesion volumes (thalamus: r=-0.56, p=9.96 x 10(-9); caudate nucleus: r=-0.31, p=3.10 x 10(-3); putamen: r=-0.50, p=6.06 x 10(-7)), 2) T1 hypointense lesion volumes (thalamus: r=-0.61, p=2.29 x 10(-10); caudate nucleus: r=-0.35, p=9.51 x 10(-4); putamen: r=-0.43, p=3.51 x 10(-5)), and 3) normalized CSF volume (thalamus: r=-0.66, p=3.55 x 10(-12); caudate nucleus: r=-0.52, p=2.31 x 10(-7), and putamen: r=-0.66, r=2.13 x 10(-12)). More severe atrophy was observed mainly in thalamus at higher EDSS. These studies appear to suggest a link between the white matter damage and DGM atrophy in MS.

Proton magnetic resonance spectroscopy in multiple sclerosis.

Proton magnetic resonance spectroscopy ((1)H-MRS) provides tissue metabolic information in vivo. This article reviews the role of MRS-determined metabolic alterations in lesions, normal-appearing white matter, gray matter, and spinal cord in advancing our knowledge of pathologic changes in multiple sclerosis (MS). In addition, the role of MRS in objectively evaluating therapeutic efficacy is reviewed. This potential metabolic information makes MRS a unique tool to follow MS disease evolution, understand its pathogenesis, evaluate the disease severity, establish a prognosis, and objectively evaluate the efficacy of therapeutic interventions.

Automatic quantitation of multiple sclerosis lesions on MR images

A two-pronged approach for the automatic quantitation of multiple sclerosis (MS) lesions on magnetic resonance (MR) images has been developed. This method includes the design and use of a pulse sequence for improved lesion-to-tissue contrast (LTC) and seeks to identify and minimize the sources of false lesion classifications in segmented images. The new pulse sequence, referred to as AFFIRMATIVE (Attenuation of Fluid by Fast Inversion Recovery with MAgnetization Transfer Imaging with Variable Echoes), improves the LTC, relative to spin-echo images, by combining Fluid-Attenuated Inversion Recovery (FLAIR) and Magnetization Transfer Contrast (MTC). In addition to acquiring fast FLAIR/MTC images, the AFFIRMATIVE sequence simultaneously acquires fast spin-echo (FSE) images for spatial registration of images, which is necessary for accurate lesion quantitation. Flow has been found to be a primary source of false lesion classifications. Therefore, an imaging protocol and reconstruction methods are developed to generate "flow images" which depict both coherent (vascular) and incoherent (CSF) flow. An automatic technique is designed for the removal of extra-meningeal tissues, since these are known to be sources of false lesion classifications. A retrospective, three-dimensional (3D) registration algorithm is implemented to correct for patient movement which may have occurred between AFFIRMATIVE and flow imaging scans. Following application of these pre-processing steps, images are segmented into white matter, gray matter, cerebrospinal fluid, and MS lesions based on AFFIRMATIVE and flow images using an automatic algorithm. All algorithms are seamlessly integrated into a single MR image analysis software package. Lesion quantitation has been performed on images from 15 patient volunteers. The total processing time is less than two hours per patient on a SPARCstation 20. The automated nature of this approach should provide an objective means of monitoring the progression, stabilization, and/or regression of MS lesions in large-scale, multi-center clinical trials. ^

In-vivo diffusion tensor imaging of rat spinal cord with a phased array coil at 7T

Magnetic resonance imaging (MRI) is a non-invasive technique that offers excellent soft tissue contrast for characterizing soft tissue pathologies. Diffusion tensor imaging (DTI) is an MRI technique that has shown to have the sensitivity to detect subtle pathology that is not evident on conventional MRI. ^ Rats are commonly used as animal models in characterizing the spinal cord pathologies including spinal cord injury (SCI), cancer, multiple sclerosis, etc. These pathologies could affect both thoracic and cervical regions and complete characterization of these pathologies using MRI requires DTI characterization in both the thoracic and cervical regions. Prior to the application of DTI for investigating the pathologic changes in the spinal cord, it is essential to establish DTI metrics in normal animals. ^ To date, in-vivo DTI studies of rat spinal cord have used implantable coils for high signal-to-noise ratio (SNR) and spin-echo pulse sequences for reduced geometric distortions. Implantable coils have several disadvantages including: (1) the invasive nature of implantation, (2) loss of SNR due to frequency shift with time in the longitudinal studies, and (3) difficulty in imaging the cervical region. While echo planar imaging (EPI) offers much shorter acquisition times compared to spin-echo imaging, EPI is very sensitive to static magnetic field inhomogeneities and the existing shimming techniques implemented on the MRI scanner do not perform well on spinal cord because of its geometry. ^ In this work, an integrated approach has been implemented for in-vivo DTI characterization of rat spinal cord in the thoracic and cervical regions. A three element phased array coil was developed for improved SNR and extended spatial coverage. A field-map shimming technique was developed for minimizing the geometric distortions in EPI images. Using these techniques, EPI based DWI images were acquired with optimized diffusion encoding scheme from 6 normal rats and the DTI-derived metrics were quantified. ^ The phantom studies indicated higher SNR and smaller bias in the estimated DTI metrics than the previous studies in the cervical region. In-vivo results indicated no statistical difference in the DTI characteristics of either gray matter or white matter between the thoracic and cervical regions. ^

Re-emerging Vitamin D deficiency epidemic and its implications on the public health in the US

Vitamin D is essential in maintaining the bone health and Calcium homeostasis in the body. These actions are mediated through the Vitamin D receptors (VDR) present in cells through which the activated vitamin D acts [1]. In the past, it was known that these receptors existed in the intestine and bone cell. However, recent discovery of VDR in other tissues as well, has broadened the action of Vitamin D and increased its adequate intake [1].^ In the past, Vitamin D deficiency was most common among institutionalized, elderly patients and children and thought to be extinct in the healthy population. However, recent evidence has shown that, prevalence of vitamin D deficiency is increasing into an epidemic status in the overall population of the United States, including the healthy individuals [2-3]. The increased daily-recommended requirement and other multiple factors are responsible for the re-emergence of this epidemic [4-5]. Some of these factors could be used to control the epidemic. Studies have also shown the association between vitamin D deficiency and increased risk for developing chronic diseases such as diabetes, hypertension, multiple sclerosis, arthritis, and some fatal cancers like prostate, colon and breast cancers [1, 4, 6-14]. This issue results in increased disease burden, morbidity and mortality in the community [15-20].^ Methods: The literature search was conducted using the University of Texas Health Science Center at Houston (UTHSC) and University of Texas Southwestern Medical Center (UTSW) online library. The key search terms used are “vitamin D deficiency And prevalence Or epidemiology”, “vitamin D deficiency And implication And public health” using PubMed and Mesh database and “vitamin D deficiency” using systematic reviews. The search is limited to Humans and the English language. The articles considered for the review are limited to Healthy US population to avoid health conditions that predispose the population to vitamin D deficiency. Only US population is considered to narrow down the study.^ Results: There is an increased prevalence of low levels of Vitamin D levels below the normal range in the US population regardless of age and health status. Vitamin D deficiency is also associated with increased risk of chronic illnesses and fatal cancers.^ Conclusion: This increased prevalence and the association of the deficiency with increased all-cause mortality has increased the economic burden and compromised the quality of life among the population. This necessitates the health care providers to routinely screen their patients for the Vitamin D status and counsel them to avoid the harmful effects of the Vitamin D deficiency. ^

Characterization of anti -sulfatide autoantibodies in demyelinating diseases

Multiple sclerosis (MS) is the most common autoimmune disease of the central nerve system and Guillain Barré Syndrome (GBS) is an inflammatory neuropathy involving the peripheral nerves. Anti-myelin immunoglobins may play a role in the demyelination processes of the both diseases. Sulfatide is an abundant glycolipid on myelin and is a candidate target antigen for disease related autoantibodies. The objective of this study was to characterize anti-sulfatide antibodies and compare antibodies from GBS and MS patients with fetal antibodies. Our hypothesis is that some B cells producing disease-associated autoantibodies are derived from or related to B cells of the fetal repertoire. Here we report that reactivity of plasma IgM against sulfatide was elevated in twelve MS patients compared with twelve normal subjects. This result implies that anti-sulfatide antibodies are disease-related. A total of sixteen human B lymphocyte clones producing anti-sulfatide autoantibodies were isolated from MS patients, GBS patients and a human fetus. Seven of the clones were from three MS patients, four of the clones were from three GBS patients and five were from the spleen of a twenty-week human fetus. Sequences have been obtained for the heavy and light chain variable regions (VDJ and VJ regions) of all of the anti-sulfatide immunoglobulins. Seven of the sixteen antibodies used VH3 for the variable region gene of the heavy chain consistent with the rate of VH3 usage in randomly selected B cells. Somatic mutations were significantly more frequent in the patient antibodies than in the fetus and somatic mutations in CDR's (Complementarity Determining Region) were significantly more frequent than in framework regions. No significant difference was found between patients and fetus in length of VH CDRIII. However, it is reported that antibodies from randomly selected normal adult B cells have longer CDRIII lengths than those of the fetus (Sanz I, 1991 Journal of Immunology Sep 1;147(5):1720-9). Our results are consistent with derivation of the precursors of B cells producing these autoantibodies from B cells related to those of the fetal repertoire. These findings are consistent with a model in which quiescent B cells from clones produced early in development undergo proliferation in dysregulated disease states, accumulating somatic mutations and increasing in reactivity toward self-antigens. ^ Epitope mapping and molecular modeling were done to elucidate the relationships between antibody structure and binding characteristics. The autoantibodies were tested for binding activity to three different antigens: sulfatide, galactoceramide and ceramide. Molecular modeling suggests that antibodies with positive charge surrounded by or adjacent to hydrophobic groups in the binding pocket bind to the head of sulfatide via the sulfate group through electrostatic interactions. However, the antibodies with hydrophobic groups separated from positive charges appear to bind to the hydrophobic tail of sulfatide. This observation was supported by a study of the effect of NaCl concentration on antigen binding. The result suggested that electrostatic interactions played a major role in sulfate group binding and that hydrophobic interactions were of greater importance for binding to the ceramide group. Our three-dimensional structure data indicated that epitope specificity of these antibodies is more predictable at the level of tertiary than primary structure and suggested positive selection based on structure occurred in the. formation of those autoantibodies. ^

Thoracic aortic disease in tuberous sclerosis complex: molecular pathogenesis and potential therapies in Tsc2+/- mice.

Tuberous sclerosis complex (TSC) is a genetic disorder with pleiotropic manifestations caused by heterozygous mutations in either TSC1 or TSC2. One of the less investigated complications of TSC is the formation of aneurysms of the descending aorta, which are characterized on pathologic examination by smooth muscle cell (SMC) proliferation in the aortic media. SMCs were explanted from Tsc2(+/-) mice to investigate the pathogenesis of aortic aneurysms caused by TSC2 mutations. Tsc2(+/-) SMCs demonstrated increased phosphorylation of mammalian target of rapamycin (mTOR), S6 and p70S6K and increased proliferation rates compared with wild-type (WT) SMCs. Tsc2(+/-) SMCs also had reduced expression of SMC contractile proteins compared with WT SMCs. An inhibitor of mTOR signaling, rapamycin, decreased SMC proliferation and increased contractile protein expression in the Tsc2(+/-) SMCs to levels similar to WT SMCs. Exposure to alpha-elastin fragments also decreased proliferation of Tsc2(+/-) SMCs and increased levels of p27(kip1), but failed to increase expression of contractile proteins. In response to artery injury using a carotid artery ligation model, Tsc2(+/-) mice significantly increased neointima formation compared with the control mice, and the neointima formation was inhibited by treatment with rapamycin. These results demonstrate that Tsc2 haploinsufficiency in SMCs increases proliferation and decreases contractile protein expression and suggest that the increased proliferative potential of the mutant cells may be suppressed in vivo by interaction with elastin. These findings provide insights into the molecular pathogenesis of aortic disease in TSC patients and identify a potential therapeutic target for treatment of this complication of the disease.

An Innovative Phase I Trial Design Allowing for the Identification of Multiple Potential Maximum Tolerated Doses with Combination Therapy of Targeted Agents

Treatment for cancer often involves combination therapies used both in medical practice and clinical trials. Korn and Simon listed three reasons for the utility of combinations: 1) biochemical synergism, 2) differential susceptibility of tumor cells to different agents, and 3) higher achievable dose intensity by exploiting non-overlapping toxicities to the host. Even if the toxicity profile of each agent of a given combination is known, the toxicity profile of the agents used in combination must be established. Thus, caution is required when designing and evaluating trials with combination therapies. Traditional clinical design is based on the consideration of a single drug. However, a trial of drugs in combination requires a dose-selection procedure that is vastly different than that needed for a single-drug trial. When two drugs are combined in a phase I trial, an important trial objective is to determine the maximum tolerated dose (MTD). The MTD is defined as the dose level below the dose at which two of six patients experience drug-related dose-limiting toxicity (DLT). In phase I trials that combine two agents, more than one MTD generally exists, although all are rarely determined. For example, there may be an MTD that includes high doses of drug A with lower doses of drug B, another one for high doses of drug B with lower doses of drug A, and yet another for intermediate doses of both drugs administered together. With classic phase I trial designs, only one MTD is identified. Our new trial design allows identification of more than one MTD efficiently, within the context of a single protocol. The two drugs combined in our phase I trial are temsirolimus and bevacizumab. Bevacizumab is a monoclonal antibody targeting the vascular endothelial growth factor (VEGF) pathway which is fundamental for tumor growth and metastasis. One mechanism of tumor resistance to antiangiogenic therapy is upregulation of hypoxia inducible factor 1α (HIF-1α) which mediates responses to hypoxic conditions. Temsirolimus has resulted in reduced levels of HIF-1α making this an ideal combination therapy. Dr. Donald Berry developed a trial design schema for evaluating low, intermediate and high dose levels of two drugs given in combination as illustrated in a recently published paper in Biometrics entitled “A Parallel Phase I/II Clinical Trial Design for Combination Therapies.” His trial design utilized cytotoxic chemotherapy. We adapted this design schema by incorporating greater numbers of dose levels for each drug. Additional dose levels are being examined because it has been the experience of phase I trials that targeted agents, when given in combination, are often effective at dosing levels lower than the FDA-approved dose of said drugs. A total of thirteen dose levels including representative high, intermediate and low dose levels of temsirolimus with representative high, intermediate, and low dose levels of bevacizumab will be evaluated. We hypothesize that our new trial design will facilitate identification of more than one MTD, if they exist, efficiently and within the context of a single protocol. Doses gleaned from this approach could potentially allow for a more personalized approach in dose selection from among the MTDs obtained that can be based upon a patient’s specific co-morbid conditions or anticipated toxicities.

CHARACTERIZING AND TREATING THE NEUROPATHOLOGY OF TUBEROUS SCLEROSIS COMPLEX IN THE MOUSE

Tuberous sclerosis complex (TSC) is a multisystem, autosomal dominant disorder affecting approximately 1 in 6000 births. Developmental brain abnormalities cause substantial morbidity and mortality and often lead to neurological disease including epilepsy, cognitive disabilities, and autism. TSC is caused by inactivating mutations in either TSC1 or TSC2, whose protein products are known inhibitors of mTORC1, an important kinase regulating translation and cell growth. Nonetheless, neither the pathophysiology of the neurological manifestations of TSC nor the extent of mTORC1 involvement in the development of these lesions is known. Murine models would greatly advance the study of this debilitating disorder. This thesis will describe the generation and characterization of a novel brain-specific mouse model of TSC, Tsc2flox/ko;hGFAP-Cre. In this model, the Tsc2 gene has been removed from most neurons and glia of the cortex and hippocampus by targeted Cre-mediated deletion in radial glial neuroprogenitor cells. The Tsc2flox/ko;hGFAP-Cre mice fail to thrive beginning postnatal day 8 and die from seizures around 23 days. Further characterization of these mice demonstrated megalencephaly, enlarged neurons, abnormal neuronal migration, altered progenitor pools, hypomyelination, and an astrogliosis. The similarity of these defects to those of TSC patients establishes this mouse as an excellent model for the study of the neuropathology of TSC and testing novel therapies. We further describe the use of this mouse model to assess the therapeutic potential of the macrolide rapamycin, an inhibitor of mTORC1. We demonstrate that rapamycin administered from postnatal day 10 can extend the life of the mutant animals 5 fold. Since TSC is a neurodevelopmental disorder, we also assessed in utero and/or immediate postnatal treatment of the animals with rapamycin. Amazingly, combined in utero and postnatal rapamycin effected a histologic rescue that was almost indistinguishable from control animals, indicating that dysregulation of mTORC1 plays a large role in TSC neuropathology. In spite of the almost complete histologic rescue, behavioral studies demonstrated that combined treatment resulted in poorer learning and memory than postnatal treatment alone. Postnatally-treated animals behaved similarly to treated controls, suggesting that immediate human treatment in the newborn period might provide the most opportune developmental timepoint for rapamycin administration.

Safety, pharmacokinetics, and antiretroviral activity of multiple doses of ibalizumab (formerly TNX-355), an anti-CD4 monoclonal antibody in HIV-1 infected adults

Context: Despite tremendous strides in HIV treatment over the past decade, resistance remains a major problem. A growing number of patients develop resistance and require new therapies to suppress viral replication. ^ Objective: To assess the safety of multiple administrations of the anti-CD4 receptor (anti-CD4) monoclonal antibody ibalizumab given as intravenous (IV) infusions, in three dosage regimens, in subjects infected with human immunodeficiency virus (HIV-1). ^ Design: Phase 1, multi-center, open-label, randomized clinical trial comparing the safety, pharmacokinetics and antiviral activity of three dosages of ibalizumab. ^ Setting: Six clinical trial sites in the United States. ^ Participants: A total of twenty-two HIV-positive patients on no anti-retroviral therapy or a stable failing regimen. ^ Intervention: Randomized to one of two treatment groups in Arms A and B followed by non-randomized enrollment in Arm C. Patients randomized to Arm A received 10 mg/kg of ibalizumab every 7 days, for a total of 10 doses; patients randomized to Arm B received a total of six doses of ibalizumab; a single loading dose of 10 mg/kg on Day 1 followed by five maintenance doses of 6 mg/kg every 14 days, starting at Week 1. Patients assigned to Arm C received 25 mg/kg of ibalizumab every 14 days for a total of 5 doses. All patients were followed for safety for an additional 7 to 8 weeks. ^ Main Outcome Measures: Clinical and laboratory assessments of safety and tolerability of multiple administrations of ibalizumab in HIV-infected patients. Secondary measures of efficacy include HIV-1 RNA (viral load) measurements. ^ Results: 21 patients were treatment-experienced and 1 was naïve to HIV therapy. Six patients were failing despite therapy and 15 were on no current HIV treatment. Mean baseline viral load (4.78 log 10; range 3.7-5.9) and CD4+ cell counts (332/μL; range 89-494) were similar across cohorts. Mean peak decreases in viral load from baseline of 0.99 log10(1.11 log10, and 0.96 log 10 occurred by Wk 2 in Cohorts A, B and C, respectively. Viral loads decreased by >1.0 log10 in 64%; 4 patients viral loads were suppressed to < 400 copies/mL. Viral loads returned towards baseline by Week 9 with reduced susceptibility to ibalizumab. CD4+ cell counts rose transiently and returned toward baseline. Maximum median elevations above BL in CD4+ cell counts for Cohorts A, B and C were +257, +198 and +103 cells/μL, respectively and occurred within 3 Wks in 16 of 22 subjects. The half-life of ibalizumab was 3-3.5 days and elimination was characteristic of capacity-limited kinetics. Administration of ibalizumab was well tolerated. Four serious adverse events were reported during the study. None of these events were related to study drug. Headache, nausea and cough were the most frequently reported treatment emergent adverse events and there were no laboratory abnormalities related to study drug. ^ Conclusions: Ibalizumab administered either weekly or bi-weekly was safe, well tolerated, and demonstrated antiviral activity. Further studies with ibalizumab in combination with standard antiretroviral treatments are warranted.^

Single nucleotide polymorphisms (SNPs) associated with TGF-beta pathway and their significance in systemic sclerosis - A multilevel analysis

Systemic sclerosis (SSc) or Scleroderma is a complex disease and its etiopathogenesis remains unelucidated. Fibrosis in multiple organs is a key feature of SSc and studies have shown that transforming growth factor-β (TGF-β) pathway has a crucial role in fibrotic responses. For a complex disease such as SSc, expression quantitative trait loci (eQTL) analysis is a powerful tool for identifying genetic variations that affect expression of genes involved in this disease. In this study, a multilevel model is described to perform a multivariate eQTL for identifying genetic variation (SNPs) specifically associated with the expression of three members of TGF-β pathway, CTGF, SPARC and COL3A1. The uniqueness of this model is that all three genes were included in one model, rather than one gene being examined at a time. A protein might contribute to multiple pathways and this approach allows the identification of important genetic variations linked to multiple genes belonging to the same pathway. In this study, 29 SNPs were identified and 16 of them located in known genes. Exploring the roles of these genes in TGF-β regulation will help elucidate the etiology of SSc, which will in turn help to better manage this complex disease. ^