Novel Statistical Methods in Quantitative Genetics : Modeling Genetic Variance for Quantitative Trait Loci Mapping and Genomic Evaluation

This thesis develops and evaluates statistical methods for different types of genetic analyses, including quantitative trait loci (QTL) analysis, genome-wide association study (GWAS), and genomic evaluation. The main contribution of the thesis is to provide novel insights in modeling genetic variance, especially via random effects models. In variance component QTL analysis, a full likelihood model accounting for uncertainty in the identity-by-descent (IBD) matrix was developed. It was found to be able to correctly adjust the bias in genetic variance component estimation and gain power in QTL mapping in terms of precision.  Double hierarchical generalized linear models, and a non-iterative simplified version, were implemented and applied to fit data of an entire genome. These whole genome models were shown to have good performance in both QTL mapping and genomic prediction. A re-analysis of a publicly available GWAS data set identified significant loci in Arabidopsis that control phenotypic variance instead of mean, which validated the idea of variance-controlling genes.  The works in the thesis are accompanied by R packages available online, including a general statistical tool for fitting random effects models (hglm), an efficient generalized ridge regression for high-dimensional data (bigRR), a double-layer mixed model for genomic data analysis (iQTL), a stochastic IBD matrix calculator (MCIBD), a computational interface for QTL mapping (qtl.outbred), and a GWAS analysis tool for mapping variance-controlling loci (vGWAS).

Estimation of Parameters in Random Effect Models with Incidence Matrix Uncertainty

Random effect models have been widely applied in many fields of research. However, models with uncertain design matrices for random effects have been little investigated before. In some applications with such problems, an expectation method has been used for simplicity. This method does not include the extra information of uncertainty in the design matrix is not included. The closed solution for this problem is generally difficult to attain. We therefore propose an two-step algorithm for estimating the parameters, especially the variance components in the model. The implementation is based on Monte Carlo approximation and a Newton-Raphson-based EM algorithm. As an example, a simulated genetics dataset was analyzed. The results showed that the proportion of the total variance explained by the random effects was accurately estimated, which was highly underestimated by the expectation method. By introducing heuristic search and optimization methods, the algorithm can possibly be developed to infer the 'model-based' best design matrix and the corresponding best estimates.

Assessing the physical environment of older people’s residential care facilities: development of the Swedish version of the Sheffield Care Environment Assessment Matrix (S-SCEAM)

Background There is emerging evidence that the physical environment is important for health, quality of life and care, but there is a lack of valid instruments to assess health care environments. The Sheffield Care Environment Assessment Matrix (SCEAM), developed in the United Kingdom, provides a comprehensive assessment of the physical environment of residential care facilities for older people. This paper reports on the translation and adaptation of SCEAM for use in Swedish residential care facilities for older people, including information on its validity and reliability. Methods SCEAM was translated into Swedish and back-translated into English, and assessed for its relevance by experts using content validity index (CVI) together with qualitative data. After modification, the validity assessments were repeated and followed by test-retest and inter-rater reliability tests in six units within a Swedish residential care facility that varied in terms of their environmental characteristics. Results Translation and back translation identified linguistic and semantic related issues. The results of the first content validity analysis showed that more than one third of the items had item-CVI (I-CVI) values less than the critical value of 0.78.  After modifying the instrument, the second content validation analysis resulted in I-CVI scores above 0.78, the suggested criteria for excellent content validity. Test-retest reliability showed high stability (96% and 95% for two independent raters respectively), and inter-rater reliability demonstrated high levels of agreement (95% and 94% on two separate rating occasions). Kappa values were very good for test-retest (κ= 0.903 and 0.869) and inter-rater reliability (κ= 0.851 and 0.832). Conclusions Adapting an instrument to a domestic context is a complex and time-consuming process, requiring an understanding of the culture where the instrument was developed and where it is to be used. A team, including the instrument’s developers, translators, and researchers is necessary to ensure a valid translation and adaption. This study showed preliminary validity and reliability evidence for the Swedish version (S-SCEAM) when used in a Swedish context. Further, we believe that the S-SCEAM has improved compared to the original instrument and suggest that it can be used as a foundation for future developments of the SCEAM model.

Expression of matrix metalloproteinases, type IV collagen, and interleukin-10 in rabbits treated with morphine after lamellar keratectomy

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

Evaluation of an average numerator relationship matrix model and a Bayesian hierarchical model for growth traits in Nellore cattle with uncertain paternity

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

Parameterized fast decoupled power flow methods for obtaining the maximum loading point of power systems. Part I. Mathematical modeling

The conventional Newton and fast decoupled power flow (FDPF) methods have been considered inadequate to obtain the maximum loading point of power systems due to ill-conditioning problems at and near this critical point. It is well known that the PV and Q-theta decoupling assumptions of the fast decoupled power flow formulation no longer hold in the vicinity of the critical point. Moreover, the Jacobian matrix of the Newton method becomes singular at this point. However, the maximum loading point can be efficiently computed through parameterization techniques of continuation methods. In this paper it is shown that by using either theta or V as a parameter, the new fast decoupled power flow versions (XB and BX) become adequate for the computation of the maximum loading point only with a few small modifications. The possible use of reactive power injection in a selected PV bus (Q(PV)) as continuation parameter (mu) for the computation of the maximum loading point is also shown. A trivial secant predictor, the modified zero-order polynomial which uses the current solution and a fixed increment in the parameter (V, theta, or mu) as an estimate for the next solution, is used in predictor step. These new versions are compared to each other with the purpose of pointing out their features, as well as the influence of reactive power and transformer tap limits. The results obtained with the new approach for the IEEE test systems (14, 30, 57 and 118 buses) are presented and discussed in the companion paper. The results show that the characteristics of the conventional method are enhanced and the region of convergence around the singular solution is enlarged. In addition, it is shown that parameters can be switched during the tracing process in order to efficiently determine all the PV curve points with few iterations. (C) 2003 Elsevier B.V. All rights reserved.

Matrix metalloproteinase-9 expression in pterygium

OBJETIVO: Avaliar a expressão da metaloprotease de matriz (MMP)-9 nos pterígios. MÉTODOS: Foi realizado na Faculdade de Medicina de Botucatu estudo prospectivo, aleatório, com o intuito de avaliar a expressão da metaloprotease de matriz na cápsula de Tenon normal e de pterígios primários e recidivados, usando o método da imuno-histoquímica e o sistema computadorizado de análise de imagem. Os resultados foram avaliados estatisticamente. RESULTADOS: A expressão da metaloprotease de matriz foi semelhante na cápsula de Tenon normal e nos pterígios primários e recidivados. CONCLUSÃO: A expressão da metaloprotease de matriz na cápsula de Tenon normal e nos pterígios primários ou recidivados é semelhante, o que nos leva a concluir que esta metaloprotease de matriz não esteja envolvida na gênese ou na recidiva do pterígio.

Influence of Bioactive Glass and/or Acellular Dermal Matrix on Bone Healing of Surgically Created Defects in Rat Tibiae: A Histological and Histometric Study

The purpose of this study was to histologically analyze the influence of bioactive glass and/or acellular dermal matrix on bone healing in surgically created defects in the tibiae of 64 rats (Rattus norvegicus, albinus, Wistar). Materials and Methods: A 4-mm X 3-mm unicortical defect was created on the anterolateral surface of the tibia. Animals were divided into 4 groups: C, control; BG, the defect was filled with bioactive glass; ADM, the defect was covered with acellular dermal matrix; and BG/ADM, the defect was filled with bioactive glass and covered with acellular dermal matrix. Animals were sacrificed at 10 or 30 days postoperatively, and the specimens were removed for histologic processing. The formation of new bone in the cortical area of the defect was evaluated histomorphometrically. Results: At 10 and 30 days postoperatively, groups C (39.65% +/- 5.63% / 63.34% +/- 5.22%) and ADM (38.12% +/- 5.53 / 58.96% +/- 7.05%) presented a larger amount of bone formation compared to the other groups (P<.05). In the same periods, groups BG (13.10% +/- 6.29% / 29.5% +/- 5.56%) and BG/ADM (20.72% +/- 8.31% / 24.19% +/- 6.69%) exhibited statistically similar new bone formation. However, unlike the other groups, group BG/ADM did not present a significant increase in bone formation between the 2 time points. Conclusion: Based on these results, it can be concluded that all of the materials used in this study delayed bone healing in non-critical-size defects. INT J ORAL MAXILLOFAC IMPLANTS 2008;23:811-817

Long-term effect of carbodiimide on dentin matrix and resin-dentin bonds

Objectives: To characterize the interaction of 1-Ethyl-3-[3-dimethylaminopropyl] carbodiimide Hydrochloride (EDC) with dentin matrix and its effect on the resin-dentin bond. Methods: Changes to the stiffness of demineralized dentin fragments treated with EDC/N-hydroxysuccinimide (NHS) in different solutions were evaluated at different time points. The resistance against enzymatic degradation was indirectly evaluated by ultimate tensile strength (UTS) test of demineralized dentin treated or not with EDC/NHS and subjected to collagenase digestion. Short- and long-term evaluations of the strength of resin-dentin interfaces treated with EDC/NHS for 1 h were performed using microtensile bond strength (mu TBS) test. All data (MPa) were individually analyzed using ANOVA and Tukey HSD tests (alpha = 0.05). Results: The different exposure times significantly increased the stiffness of dentin (p < 0.0001, control-5.15 and EDC/NHS-29.50), while no differences were observed among the different solutions of EDC/NHS (p = 0.063). Collagenase challenge did not affect the UTS values of EDC/NHS group (6.08) (p > 0.05), while complete degradation was observed for the control group (p = 0.0008, control-20.84 and EDC/NHS-43.15). EDC/NHS treatment did not significantly increase resin-dentin mu TBS, but the values remained stable after 12 months water storage (p < 0.05). Conclusions: Biomimetic use of EDC/NHS to induce exogenous collagen cross-links resulted in increased mechanical properties and stability of dentin matrix and dentin-resin interfaces. (C) 2010 Wiley Periodicals, Inc. J Biomed Mater Res Part B: Appl Biomater 94B: 250-255, 2010.

Staining methods applied to glycol methacrylate embedded tissue sections

The use of glycol methacrylate (GMA) avoids some technical artifacts, which are usually observed in paraffin-embedded sections, providing good morphological resolution. on the other hand, weak staining have been mentioned during the use of different methods in plastic sections. In the present study, changes in the histological staining procedures have been assayed during the use of staining and histochemical methods in different GMA-embedded tissues.Samples of tongue, submandibular and sublingual glands, cartilage, portions of respiratory tract and nervous ganglion were fixed in 4% formaldehyde and embedded in glycol methacrylate. The sections of tongue and nervous ganglion were stained by H&E. Picrosirius, Toluidine Blue and Sudan Black B methods were applied, respectively, for identification of collagen fibers in submandibular gland, sulfated glycosaminoglycans in cartilage (metachromasia) and myelin lipids in nervous ganglion. Periodic Acid-Schiff (PAS) method was used for detection of glycoconjugates in submandibular gland and cartilage while AB/PAS combined methods were applied for detection of mucins in the respiratory tract. In addition, a combination of Alcian Blue (AB) and Picrosirius methods was also assayed in the sublingual gland sections.The GMA-embedded tissue sections showed an optimal morphological integrity and were favorable to the staining methods employed in the present study. In the sections of tongue and nervous ganglion, a good contrast of basophilic and acidophilic structures was obtained by H&E. An intense eosinophilia was observed either in the striated muscle fibers or in the myelin sheaths in which the lipids were preserved and revealed by Sudan Black B. In the cartilage matrix, a strong metachromasia was revealed by Toluidine Blue in the negatively-charged glycosaminoglycans. In the chondrocytes, glycogen granules were intensely positive to PAS method. Extracellular glycoproteins were also PAS positive in the basal membrane and in the region occupied by the lamina externa and reticular fibers surrounding each smooth muscle cells of the blood vessels. In the epithelial cells of the respiratory tract, acid and neutral mucins were histochemically detected by AB and PAS methods, respectively. Moreover, granules containing acid and neutral mucins were revealed in purple by AB and PAS concomitantly. In the sublingual gland sections, a distinct affinity of acid mucins by AB (in turquoise-blue) and collagen fibers by Picrosirius (in red) was obtained when these methods were combined. Although some routine dyes used in paraffin sections have showed a weak stain in historesin sections, our results showed that different dyes could be applied in GMA sections if modified staining procedures were assayed. Therefore, appropriate staining contrast and, thus, detection of one or different substances in a same section can be acquired in association to the good morphological resolution provided by GMA. (C) 2003 Elsevier Ltd. All rights reserved.

Effect of green Coffea arabica L. seed oil on extracellular matrix components and water-channel expression in in vitro and ex vivo human skin models

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

Bioinformatical and in vitro approaches to essential oil-induced matrix metalloproteinase inhibition

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Autogenous demineralized dentin matrix for tissue engineering applications: Radiographic and histomorphometric studies

Purpose: This work evaluated the osteoconductive properties of autogenous demineralized dentin matrix (ADDM) on surgical bone defects in the parietal bone of rabbits, using the guided bone regeneration technique and polytetrafluoroethylene (PTFE) membrane. Materials and Methods: Surgical bone defects were created in 24 adult rabbits and repaired with either ADDM and PTFE (experimental group) or PTFE alone (control group). The ADDM had been obtained from the central incisors of the experimental rabbits. The rabbits were sacrificed after 15, 30, 60, and 90 days and the defects examined radiographically and histologically. Results: Radiographically, the defects in the experimental animals achieved radiopacity more quickly than the defects in the control group. Discussion: After 15, 30, 60, and 90 days of observation following surgery, the ADDM slices appeared to stimulate new bone formation. The dentin slices were completely incorporated into the new bone tissue and were resorbed during the bone repair. Conclusions: Bone repair was accelerated on the bone defects treated with ADDM when compared to the control group.

Bose-Einstein condensation dynamics in three dimensions by the pseudospectral and finite-difference methods

We suggest a pseudospectral method for solving the three-dimensional time-dependent Gross-Pitaevskii (GP) equation, and use it to study the resonance dynamics of a trapped Bose-Einstein condensate induced by a periodic variation in the atomic scattering length. When the frequency of oscillation of the scattering length is an even multiple of one of the trapping frequencies along the x, y or z direction, the corresponding size of the condensate executes resonant oscillation. Using the concept of the differentiation matrix, the partial-differential GP equation is reduced to a set of coupled ordinary differential equations, which is solved by a fourth-order adaptive step-size control Runge-Kutta method. The pseudospectral method is contrasted with the finite-difference method for the same problem, where the time evolution is performed by the Crank-Nicholson algorithm. The latter method is illustrated to be more suitable for a three-dimensional standing-wave optical-lattice trapping potential.

An introduction to numerical methods in low-dimensional quantum systems

This is an introductory course to the Lanczos Method and Density Matrix Renormalization Group Algorithms (DMRG), two among the leading numerical techniques applied in studies of low-dimensional quantum models. The idea of studying the models on clusters of a finite size in order to extract their physical properties is briefly discussed. The important role played by the model symmetries is also examined. Special emphasis is given to the DMRG.