Sparse hierarchical Bayesian models for detecting relevant antigenic sites in virus evolution

Understanding how virus strains offer protection against closely related emerging strains is vital for creating effective vaccines. For many viruses, including Foot-and-Mouth Disease Virus (FMDV) and the Influenza virus where multiple serotypes often co-circulate, in vitro testing of large numbers of vaccines can be infeasible. Therefore the development of an in silico predictor of cross-protection between strains is important to help optimise vaccine choice. Vaccines will offer cross-protection against closely related strains, but not against those that are antigenically distinct. To be able to predict cross-protection we must understand the antigenic variability within a virus serotype, distinct lineages of a virus, and identify the antigenic residues and evolutionary changes that cause the variability. In this thesis we present a family of sparse hierarchical Bayesian models for detecting relevant antigenic sites in virus evolution (SABRE), as well as an extended version of the method, the extended SABRE (eSABRE) method, which better takes into account the data collection process. The SABRE methods are a family of sparse Bayesian hierarchical models that use spike and slab priors to identify sites in the viral protein which are important for the neutralisation of the virus. In this thesis we demonstrate how the SABRE methods can be used to identify antigenic residues within different serotypes and show how the SABRE method outperforms established methods, mixed-effects models based on forward variable selection or l1 regularisation, on both synthetic and viral datasets. In addition we also test a number of different versions of the SABRE method, compare conjugate and semi-conjugate prior specifications and an alternative to the spike and slab prior; the binary mask model. We also propose novel proposal mechanisms for the Markov chain Monte Carlo (MCMC) simulations, which improve mixing and convergence over that of the established component-wise Gibbs sampler. The SABRE method is then applied to datasets from FMDV and the Influenza virus in order to identify a number of known antigenic residue and to provide hypotheses of other potentially antigenic residues. We also demonstrate how the SABRE methods can be used to create accurate predictions of the important evolutionary changes of the FMDV serotypes. In this thesis we provide an extended version of the SABRE method, the eSABRE method, based on a latent variable model. The eSABRE method takes further into account the structure of the datasets for FMDV and the Influenza virus through the latent variable model and gives an improvement in the modelling of the error. We show how the eSABRE method outperforms the SABRE methods in simulation studies and propose a new information criterion for selecting the random effects factors that should be included in the eSABRE method; block integrated Widely Applicable Information Criterion (biWAIC). We demonstrate how biWAIC performs equally to two other methods for selecting the random effects factors and combine it with the eSABRE method to apply it to two large Influenza datasets. Inference in these large datasets is computationally infeasible with the SABRE methods, but as a result of the improved structure of the likelihood, we are able to show how the eSABRE method offers a computational improvement, leading it to be used on these datasets. The results of the eSABRE method show that we can use the method in a fully automatic manner to identify a large number of antigenic residues on a variety of the antigenic sites of two Influenza serotypes, as well as making predictions of a number of nearby sites that may also be antigenic and are worthy of further experiment investigation.

Modelling the distribution of Bonelli's eagle in Spain: Implications for conservation planning

Bonelli’s eagle, Hieraaetus fasciatus , has recently suffered a severe population decline and is currently endangered. Spain supports about 70% of the European population. We used stepwise logistic regression on a set of environmental, spatial and human variables to model Bonelli’s eagle distribution in the 5167 UTM 10 × 10 km quadrats of peninsular Spain. We obtained a model based on 16 variables, which allowed us to identify favourable and unfavourable areas for this species in Spain, as well as intermediate favourability areas. We assessed the stepwise progression of the model by comparing the model’s predictions in each step with those of the final model, and selected a parsimonious explanatory model based on three variables — slope, July temperature and precipitation — comprising 76% of the predictive capacity of the

An efficient Jeans modelling of axisymmetric galaxies with multiple stellar components

Dynamical models of stellar systems represent a powerful tool to study their internal structure and dynamics, to interpret the observed morphological and kinematical fields, and also to support numerical simulations of their evolution. We present a method especially designed to build axisymmetric Jeans models of galaxies, assumed as stationary and collisionless stellar systems. The aim is the development of a rigorous and flexible modelling procedure of multicomponent galaxies, composed of different stellar and dark matter distributions, and a central supermassive black hole. The stellar components, in particular, are intended to represent different galaxy structures, such as discs, bulges, halos, and can then have different structural (density profile, flattening, mass, scale-length), dynamical (rotation, velocity dispersion anisotropy), and population (age, metallicity, initial mass function, mass-to-light ratio) properties. The theoretical framework supporting the modelling procedure is presented, with the introduction of a suitable nomenclature, and its numerical implementation is discussed, with particular reference to the numerical code JASMINE2, developed for this purpose. We propose an approach for efficiently scaling the contributions in mass, luminosity, and rotational support, of the different matter components, allowing for fast and flexible explorations of the model parameter space. We also offer different methods of the computation of the gravitational potentials associated of the density components, especially convenient for their easier numerical tractability. A few galaxy models are studied, showing internal, and projected, structural and dynamical properties of multicomponent galaxies, with a focus on axisymmetric early-type galaxies with complex kinematical morphologies. The application of galaxy models to the study of initial conditions for hydro-dynamical and $N$-body simulations of galaxy evolution is also addressed, allowing in particular to investigate the large number of interesting combinations of the parameters which determine the structure and dynamics of complex multicomponent stellar systems.

Modelling and classification with quantile-based distributions

In this work, we explore and demonstrate the potential for modeling and classification using quantile-based distributions, which are random variables defined by their quantile function. In the first part we formalize a least squares estimation framework for the class of linear quantile functions, leading to unbiased and asymptotically normal estimators. Among the distributions with a linear quantile function, we focus on the flattened generalized logistic distribution (fgld), which offers a wide range of distributional shapes. A novel naïve-Bayes classifier is proposed that utilizes the fgld estimated via least squares, and through simulations and applications, we demonstrate its competitiveness against state-of-the-art alternatives. In the second part we consider the Bayesian estimation of quantile-based distributions. We introduce a factor model with independent latent variables, which are distributed according to the fgld. Similar to the independent factor analysis model, this approach accommodates flexible factor distributions while using fewer parameters. The model is presented within a Bayesian framework, an MCMC algorithm for its estimation is developed, and its effectiveness is illustrated with data coming from the European Social Survey. The third part focuses on depth functions, which extend the concept of quantiles to multivariate data by imposing a center-outward ordering in the multivariate space. We investigate the recently introduced integrated rank-weighted (IRW) depth function, which is based on the distribution of random spherical projections of the multivariate data. This depth function proves to be computationally efficient and to increase its flexibility we propose different methods to explicitly model the projected univariate distributions. Its usefulness is shown in classification tasks: the maximum depth classifier based on the IRW depth is proven to be asymptotically optimal under certain conditions, and classifiers based on the IRW depth are shown to perform well in simulated and real data experiments.

Modelling, simulations and monitoring for sustainable plant and animal production systems through precision technique

Protected crop production is a modern and innovative approach to cultivating plants in a controlled environment to optimize growth, yield, and quality. This method involves using structures such as greenhouses or tunnels to create a sheltered environment. These productive solutions are characterized by a careful regulation of variables like temperature, humidity, light, and ventilation, which collectively contribute to creating an optimal microclimate for plant growth. Heating, cooling, and ventilation systems are used to maintain optimal conditions for plant growth, regardless of external weather fluctuations. Protected crop production plays a crucial role in addressing challenges posed by climate variability, population growth, and food security. Similarly, animal husbandry involves providing adequate nutrition, housing, medical care and environmental conditions to ensure animal welfare. Then, sustainability is a critical consideration in all forms of agriculture, including protected crop and animal production. Sustainability in animal production refers to the practice of producing animal products in a way that minimizes negative impacts on the environment, promotes animal welfare, and ensures the long-term viability of the industry. Then, the research activities performed during the PhD can be inserted exactly in the field of Precision Agriculture and Livestock farming. Here the focus is on the computational fluid dynamic (CFD) approach and environmental assessment applied to improve yield, resource efficiency, environmental sustainability, and cost savings. It represents a significant shift from traditional farming methods to a more technology-driven, data-driven, and environmentally conscious approach to crop and animal production. On one side, CFD is powerful and precise techniques of computer modeling and simulation of airflows and thermo-hygrometric parameters, that has been applied to optimize the growth environment of crops and the efficiency of ventilation in pig barns. On the other side, the sustainability aspect has been investigated and researched in terms of Life Cycle Assessment analyses.

[diet Quality And Associated Factors Among The Elderly: A Population-based Study In Campinas, São Paulo State, Brazil].

The aim of this study was to assess the quality of diet among the elderly and associations with socio-demographic variables, health-related behaviors, and diseases. A population-based cross-sectional study was conducted in a representative sample of 1,509 elderly participants in a health survey in Campinas, São Paulo State, Brazil. Food quality was assessed using the Revised Diet Quality Index (DQI-R). Mean index scores were estimated and a multiple regression model was employed for the adjusted analyses. The highest diet quality scores were associated with age 80 years or older, Evangelical religion, diabetes mellitus, and physical activity, while the lowest scores were associated with home environments shared with three or more people, smoking, and consumption of soft drinks and alcoholic beverages. The findings emphasize a general need for diet quality improvements in the elderly, specifically in subgroups with unhealthy behaviors, who should be targeted with comprehensive strategies.

Falls Among The Elderly: Risk Factors In A Population-based Study.

The aim of the present study was to identify factors associated with the occurrence of falls among elderly adults in a population-based study (ISACamp 2008). A population-based cross-sectional study was carried out with two-stage cluster sampling. The sample was composed of 1,520 elderly adults living in the urban area of the city of Campinas, São Paulo, Brazil. The occurrence of falls was analyzed based on reports of the main accident occurred in the previous 12 months. Data on socioeconomic/demographic factors and adverse health conditions were tested for possible associations with the outcome. Prevalence ratios (PR) were estimated and adjusted for gender and age using the Poisson multiple regression analysis. Falls were more frequent, after adjustment for gender and age, among female elderly participants (PR = 2.39; 95% confidence interval (95% CI) 1.47 - 3.87), elderly adults (80 years old and older) (PR = 2.50; 95% CI 1.61 - 3.88), widowed (PR = 1.74; 95% CI 1.04 - 2.89) and among elderly adults who had rheumatism/arthritis/arthrosis (PR = 1.58; 95% CI 1.00 - 2.48), osteoporosis (PR = 1.71; 95% CI 1.18 - 2.49), asthma/bronchitis/emphysema (PR = 1,73; 95% CI 1.09 - 2.74), headache (PR = 1.59; 95% CI 1.07 - 2.38), mental common disorder (PR = 1.72; 95% CI 1.12 - 2.64), dizziness (PR = 2.82; 95% CI 1.98 - 4.02), insomnia (PR = 1.75; 95% CI 1.16 - 2.65), use of multiple medications (five or more) (PR = 2.50; 95% CI 1.12 - 5.56) and use of cane/walker (PR = 2.16; 95% CI 1.19 - 3,93). The present study shows segments of the elderly population who are more prone to falls through the identification of factors associated with this outcome. The findings can contribute to the planning of public health policies and programs addressed to the prevention of falls.

The Basic Reproduction Number Obtained From Jacobian And Next Generation Matrices - A Case Study Of Dengue Transmission Modelling.

The basic reproduction number is a key parameter in mathematical modelling of transmissible diseases. From the stability analysis of the disease free equilibrium, by applying Routh-Hurwitz criteria, a threshold is obtained, which is called the basic reproduction number. However, the application of spectral radius theory on the next generation matrix provides a different expression for the basic reproduction number, that is, the square root of the previously found formula. If the spectral radius of the next generation matrix is defined as the geometric mean of partial reproduction numbers, however the product of these partial numbers is the basic reproduction number, then both methods provide the same expression. In order to show this statement, dengue transmission modelling incorporating or not the transovarian transmission is considered as a case study. Also tuberculosis transmission and sexually transmitted infection modellings are taken as further examples.

Clinicopathological And Immunohistochemical Evaluation Of Oral And Oropharyngeal Squamous Cell Carcinoma In Chilean Population.

In oral and oropharyngeal squamous cell carcinoma (OCSCC and OPSCC) exist an association between clinical and histopathological parameters with cell proliferation, basal lamina, connective tissue degradation and surrounding stroma markers. We evaluated these associations in Chilean patients. A convenience sample of 37 cases of OCSCC (n=16) and OPSCC (n=21) was analyzed clinically (TNM, clinical stage) and histologically (WHO grade of differentiation, pattern of tumor invasion). We assessed the expression of p53, Ki67, HOXA1, HOXB7, type IV collagen (ColIV) and carcinoma-associated fibroblast (α-SMA-positive cells). Additionally we conducted a univariate/bivariate analysis to assess the relationship of these variables with survival rates. Males were mostly affected (56.2% OCSCC, 76.2% OPSCC). Patients were mainly diagnosed at III/IV clinical stages (68.8% OCSCC, 90.5% OPSCC) with a predominantly infiltrative pattern invasion (62.9% OCSCC, 57.1% OPSCC). Significant association between regional lymph nodes (N) and clinical stage with OCSCC-HOXB7 expression (Chi-Square test P < 0.05) was observed. In OPSCC a statistically significant association exists between p53, Ki67 with gender (Chi-Square test P < 0.05). In OCSCC and OPSCC was statistically significant association between ki67 with HOXA1, HOXB7, and between these last two antigens (Pearson's Correlation test P < 0.05). Furthermore OPSCC-p53 showed significant correlation when it was compared with α-SMA (Kendall's Tau-c test P < 0.05). Only OCSCC-pattern invasion and OPSCC-primary tumor (T) pattern resulted associated with survival at the end of the follow up period (Chi-Square Likelihood Ratio, P < 0.05). Clinical, histological and immunohistochemical features are similar to seen in other countries. Cancer proliferation markers were associated strongly from each other. Our sample highlights prognostic value of T and pattern of invasion, but the conclusions may be limited and should be considered with caution (small sample). Many cases were diagnosed in the advanced stages of the disease, which suggests that the diagnosis of OCSCC and OPSCC is made late.

[access To Basic Health Services And Associated Factors: A Population-based Study].

This study sought to identify factors involved in access to the services of a basic health unit. It is a cross-sectional, population-based study involving 101 randomly-selected families residing in the area covered by the health unit. An adult resident of each household was interviewed. The response variable was whether or not the resident frequented the health unit if he/she or anyone in the family required assistance to resolve a health issue. The independent variables investigated were service provision aspects, demographic and socio-economic characteristics, individual habits, morbidities and use of the health unit. In addition to descriptive and univariate analysis, logistic regression was applied in the multivariate analysis. The results show that access to the basic health unit is associated with the treatment received previously (OR = 3,224) with accessibility (OR = 0,146) and micro-area of residence (OR = 10,918). These findings suggest that access is related to the impressions created by the care received at the health unit and is based on experiences with the service, but can also be strongly modulated by individual aspects and factors related to the territory.

Prevalence Of Urinary Incontinence And Its Association With Multimorbidity In Women Aged 50 Years Or Older: A Population-based Study.

To evaluate the prevalence and associated risk factors for urinary incontinence, as well as its association with multimorbidity among Brazilian women aged 50 or over. This was a secondary analysis of a cross-sectional population-based study including 622 women 50 years or older, conducted in the city of Campinas-SP-Brazil. The dependent variable was Urinary Incontinence (UI), defined as any complaint of urine loss. The independent variables were sociodemographic data, health-related habits, self-perception of health and functional capacity evaluation. Statistical analysis was carried out using the Chi-square test and Poisson regression. The mean age of the women was 64. UI was prevalent in 52.3% of these women: Mixed UI (26.6%), Urge UI (13.2%) and Stress UI (12.4%). Factors associated with a higher prevalence of UI were hypertension (OR 1.21, CI 1:01-1:47, P = 0.004), osteoarthritis (OR 1.24, CI 1:03-1:50, P = 0.022), physical activity ≥3 days/week (OR 1.21, CI 1:01-1:44, P = 0.039), BMI ≥ 25 at the time of the interview (OR 1.25, CI 1:04-1:49, P = 0.018), negative self-perception of health (OR 1.23, CI 1:06-1:44 P = 0.007) and limitations in daily living activities (PR 1:56 CI 1:16-2:10, P = 0.004). The prevalence of UI was high. Mixed incontinence was the most frequent type of UI. Many associated factors can be prevented or improved. Thus, health policies targeted at these combined factors could reduce their prevalence rate and possibly decrease the prevalence of UI. Neurourol. Urodynam. © 2014 Wiley Periodicals, Inc.

Equity Of Access To Pap Smears: Population-based Study In Campinas, São Paulo, Brazil.

To determine the prevalence of the Papanicolaou exam among women aged 20 to 59 years in the city of Campinas (state of São Paulo, Brazil) and to analyze associations between this test and affiliation to private health insurance plans as well as socioeconomic/demographic variables and health-related behavior. To do so, a population-based, cross-sectional study was carried out. Statistical analyses took the study design into account. Despite the significant socioeconomic differences between women with and without private health plans, no differences between these groups were found regarding having been submitted to the Papanicolaou test. In fact no differences were found as to socioeconomic and health variables analyzed. Among all variables analyzed, only marital status was significantly associated with having undergone the test. The Brazilian public health system accounted for 55.7% of the exams. The present findings indicate social equity in the city of Campinas regarding the preventive exam for cervical cancer in the age group studied.

The Prevalence Of The Metabolically Healthy Obese Phenotype In An Aging Population And Its Association With Subclinical Cardiovascular Disease: The Brazilian Study On Healthy Aging.

Current literature has elucidated a new phenotype, metabolically healthy obese (MHO), with risks of cardiovascular disease similar to that of normal weight individuals. Few studies have examined the MHO phenotype in an aging population, especially in association with subclinical CVD. This cross sectional study population consisted of 208 octogenarians and older. Anthropometrics, biochemical, and radiological parameters were measured to assess obesity, metabolic health (assessed by the National Cholesterol Education Program -Adult Treatment Panel (NCEP-ATP III) criteria), and subclinical measures of CVD. The prevalence of MHO was 13.5% (N = 28). No significant association with MHO was noted for age, coronary artery calcium score, cIMT, or hs-CRP > 3 mg/dl (p = NS). Our results suggest that the MHO phenotype exists in the elderly; however, subclinical CVD measures were not different in sub-group analysis suggesting traditional metabolic risk factor algorithms may not be accurate in the very elderly.

Survival Extrapolation In The Presence Of Cause Specific Hazards.

Health economic evaluations require estimates of expected survival from patients receiving different interventions, often over a lifetime. However, data on the patients of interest are typically only available for a much shorter follow-up time, from randomised trials or cohorts. Previous work showed how to use general population mortality to improve extrapolations of the short-term data, assuming a constant additive or multiplicative effect on the hazards for all-cause mortality for study patients relative to the general population. A more plausible assumption may be a constant effect on the hazard for the specific cause of death targeted by the treatments. To address this problem, we use independent parametric survival models for cause-specific mortality among the general population. Because causes of death are unobserved for the patients of interest, a polyhazard model is used to express their all-cause mortality as a sum of latent cause-specific hazards. Assuming proportional cause-specific hazards between the general and study populations then allows us to extrapolate mortality of the patients of interest to the long term. A Bayesian framework is used to jointly model all sources of data. By simulation, we show that ignoring cause-specific hazards leads to biased estimates of mean survival when the proportion of deaths due to the cause of interest changes through time. The methods are applied to an evaluation of implantable cardioverter defibrillators for the prevention of sudden cardiac death among patients with cardiac arrhythmia. After accounting for cause-specific mortality, substantial differences are seen in estimates of life years gained from implantable cardioverter defibrillators.

Factors Associated With The Age Of The Onset Of Diabetes In Women Aged 50 years Or More: A Population-based Study.

Investigate factors associated with the onset of diabetes in women aged more than 49 years. Cross-sectional, population-based study using self-reports with 622 women. The dependent variable was the age of occurrence of diabetes using the life table method. Cox multiple regression models were adjusted to analyse the onset of diabetes according to predictor variables. Sociodemographic, clinical and behavioural factors were evaluated. Of the 622 women interviewed, 22.7% had diabetes. The mean age at onset was 56 years. The factors associated with the age of occurrence of diabetes were self-rated health (very good, good) (coefficient=-0.792; SE of the coefficient=0.215; p=0.0001), more than two individuals living in the household (coefficient=0.656, SE of the coefficient=0.223; p=0.003), and body mass index (BMI) (kg/m(2)) at 20-30 years of age (coefficient= 0.056, SE of the coefficient=0.023; p=0.014). Self-rated health considered good or very good was associated with a higher rate of survival without diabetes. Sharing a home with two or more other people and a weight increase at 20-30 years of age was associated with the onset of type 2 diabetes.