Accounting for baseline differences and measurement error in the analysis of change over time

If change over time is compared in several groups, it is important to take into account baseline values so that the comparison is carried out under the same preconditions. As the observed baseline measurements are distorted by measurement error, it may not be sufficient to include them as covariate. By fitting a longitudinal mixed-effects model to all data including the baseline observations and subsequently calculating the expected change conditional on the underlying baseline value, a solution to this problem has been provided recently so that groups with the same baseline characteristics can be compared. In this article, we present an extended approach where a broader set of models can be used. Specifically, it is possible to include any desired set of interactions between the time variable and the other covariates, and also, time-dependent covariates can be included. Additionally, we extend the method to adjust for baseline measurement error of other time-varying covariates. We apply the methodology to data from the Swiss HIV Cohort Study to address the question if a joint infection with HIV-1 and hepatitis C virus leads to a slower increase of CD4 lymphocyte counts over time after the start of antiretroviral therapy.

Clinical and haematological predictors of antibiotic prescribing for acute cough in adults in Swiss practices - an observational study.

BackgroundAcute cough is a common problem in general practice and is often caused by a self-limiting, viral infection. Nonetheless, antibiotics are often prescribed in this situation, which may lead to unnecessary side effects and, even worse, the development of antibiotic resistant microorganisms worldwide. This study assessed the role of point-of-care C-reactive protein (CRP) testing and other predictors of antibiotic prescription in patients who present with acute cough in general practice.MethodsPatient characteristics, symptoms, signs, and laboratory and X-ray findings from 348 patients presenting to 39 general practitioners with acute cough, as well as the GPs themselves, were recorded by fourth-year medical students during their three-week clerkships in general practice. Patient and clinician characteristics of those prescribed and not-prescribed antibiotics were compared using a mixed-effects model.ResultsOf 315 patients included in the study, 22% were prescribed antibiotics. The two groups of patients, those prescribed antibiotics and those treated symptomatically, differed significantly in age, demand for antibiotics, days of cough, rhinitis, lung auscultation, haemoglobin level, white blood cell count, CRP level and the GP¿s license to self-dispense antibiotics. After regression analysis, only the CRP level, the white blood cell count and the duration of the symptoms were statistically significant predictors of antibiotic prescription.ConclusionsThe antibiotic prescription rate of 22% in adult patients with acute cough in the Swiss primary care setting is low compared to other countries. GPs appear to use point-of-care CRP testing in addition to the duration of clinical symptoms to help them decide whether or not to prescribe antibiotics.

CYP17A1 Enzyme Activity Is Linked to Ambulatory Blood Pressure in a Family-Based Population Study

BACKGROUND Genome-wide association studies have linked CYP17A1 coding for the steroid hormone synthesizing enzyme 17α-hydroxylase (CYP17A1) to blood pressure (BP). We hypothesized that the genetic signal may translate into a correlation of ambulatory BP (ABP) with apparent CYP17A1 activity in a family-based population study and estimated the heritability of CYP17A1 activity. METHODS In the Swiss Kidney Project on Genes in Hypertension, day and night urinary excretions of steroid hormone metabolites were measured in 518 participants (220 men, 298 women), randomly selected from the general population. CYP17A1 activity was assessed by 2 ratios of urinary steroid metabolites: one estimating the combined 17α-hydroxylase/17,20-lyase activity (ratio 1) and the other predominantly 17α-hydroxylase activity (ratio 2). A mixed linear model was used to investigate the association of ABP with log-transformed CYP17A1 activities exploring effect modification by urinary sodium excretion. RESULTS Daytime ABP was positively associated with ratio 1 under conditions of high, but not low urinary sodium excretion (P interaction <0.05). Ratio 2 was not associated with ABP. Heritability estimates (SE) for day and night CYP17A1 activities were 0.39 (0.10) and 0.40 (0.09) for ratio 1, and 0.71 (0.09) and 0.55 (0.09) for ratio 2 (P values <0.001). CYP17A1 activities, assessed with ratio 1, were lower in older participants. CONCLUSIONS Low apparent CYP17A1 activity (assessed with ratio 1) is associated with elevated daytime ABP when salt intake is high. CYP17A1 activity is heritable and diminished in the elderly. These observations highlight the modifying effect of salt intake on the association of CYP17A1 with BP.

Bayesian joint modeling of longitudinal and survival data

The joint modeling of longitudinal and survival data is a new approach to many applications such as HIV, cancer vaccine trials and quality of life studies. There are recent developments of the methodologies with respect to each of the components of the joint model as well as statistical processes that link them together. Among these, second order polynomial random effect models and linear mixed effects models are the most commonly used for the longitudinal trajectory function. In this study, we first relax the parametric constraints for polynomial random effect models by using Dirichlet process priors, then three longitudinal markers rather than only one marker are considered in one joint model. Second, we use a linear mixed effect model for the longitudinal process in a joint model analyzing the three markers. In this research these methods were applied to the Primary Biliary Cirrhosis sequential data, which were collected from a clinical trial of primary biliary cirrhosis (PBC) of the liver. This trial was conducted between 1974 and 1984 at the Mayo Clinic. The effects of three longitudinal markers (1) Total Serum Bilirubin, (2) Serum Albumin and (3) Serum Glutamic-Oxaloacetic transaminase (SGOT) on patients' survival were investigated. Proportion of treatment effect will also be studied using the proposed joint modeling approaches. ^ Based on the results, we conclude that the proposed modeling approaches yield better fit to the data and give less biased parameter estimates for these trajectory functions than previous methods. Model fit is also improved after considering three longitudinal markers instead of one marker only. The results from analysis of proportion of treatment effects from these joint models indicate same conclusion as that from the final model of Fleming and Harrington (1991), which is Bilirubin and Albumin together has stronger impact in predicting patients' survival and as a surrogate endpoints for treatment. ^

Preclinical modeling of multi-drug cancer therapies

Anticancer drugs typically are administered in the clinic in the form of mixtures, sometimes called combinations. Only in rare cases, however, are mixtures approved as drugs. Rather, research on mixtures tends to occur after single drugs have been approved. The goal of this research project was to develop modeling approaches that would encourage rational preclinical mixture design. To this end, a series of models were developed. First, several QSAR classification models were constructed to predict the cytotoxicity, oral clearance, and acute systemic toxicity of drugs. The QSAR models were applied to a set of over 115,000 natural compounds in order to identify promising ones for testing in mixtures. Second, an improved method was developed to assess synergistic, antagonistic, and additive effects between drugs in a mixture. This method, dubbed the MixLow method, is similar to the Median-Effect method, the de facto standard for assessing drug interactions. The primary difference between the two is that the MixLow method uses a nonlinear mixed-effects model to estimate parameters of concentration-effect curves, rather than an ordinary least squares procedure. Parameter estimators produced by the MixLow method were more precise than those produced by the Median-Effect Method, and coverage of Loewe index confidence intervals was superior. Third, a model was developed to predict drug interactions based on scores obtained from virtual docking experiments. This represents a novel approach for modeling drug mixtures and was more useful for the data modeled here than competing approaches. The model was applied to cytotoxicity data for 45 mixtures, each composed of up to 10 selected drugs. One drug, doxorubicin, was a standard chemotherapy agent and the others were well-known natural compounds including curcumin, EGCG, quercetin, and rhein. Predictions of synergism/antagonism were made for all possible fixed-ratio mixtures, cytotoxicities of the 10 best-scoring mixtures were tested, and drug interactions were assessed. Predicted and observed responses were highly correlated (r2 = 0.83). Results suggested that some mixtures allowed up to an 11-fold reduction of doxorubicin concentrations without sacrificing efficacy. Taken together, the models developed in this project present a general approach to rational design of mixtures during preclinical drug development. ^

The impact of a child's unexpected hospitalization upon Taiwanese parents' situational anxiety

Background. The parents of a sick child likely experience situational anxiety due to their young child being unexpectedly hospitalized. The emotional upheaval may be great enough that their anxiety inhibits them in providing positive support to their hospitalized child. Because anxiety affects psychological distress as well as behavioral distress, identifying parental distress helps parents improving their coping mechanisms. ^ Purpose. The study compared situational anxiety levels between Taiwanese fathers and mothers and focused on differences between parental anxiety levels at the beginning of the child's unplanned hospitalization and at time of discharge. The study also identified factors related to the parents' distress and use of coping mechanisms. ^ Methods. A descriptive, comparative research design was used to determine the difference between the anxiety levels of 62 Taiwanese father-mother dyads during the situational crisis of their child's unexpected hospitalization. The Mandarin version (M) of Visual Analog Scale (VAS-M), State-Trait Anxiety Inventory (STAI-M), and the Index of Parent Participation/Hospitalized Child (IPP/HC-M) were used to differentiate maternal and paternal anxiety levels and identify factors related to the parents' distress. Questionnaires were completed by parents within 24-36 hours of the child's hospital admission and within 24 hours prior to discharge. A paired t-test, two sample t-test, and linear mixed regression model were used to test and support the study hypothesis. ^ Results. The findings reveal that the mothers' anxiety levels did not significantly differ from the fathers' anxiety level when their child had a sudden admission to the hospital. In particular, parental state anxiety levels did not decrease during the child's hospital stay and subsequent discharge. Moreover, anxiety levels did not differ between parents regardless of whether the child's disease was acute or chronic. The most effective factor related to parental situational anxiety was parental perception of the severity of the child's illness. ^ Conclusions. Parental anxiety was found to be significantly related to changes in their perception of the severity of their child's illness. However, the study was not able to illustrate how parental involvement in the child's hospital care was related to parental perception of the severity of their child's illness. Future studies, using a qualitative approach to gamer more information as to what variables influence parental anxiety during a situational crisis, may provide a richer database from which to modify key variables as well as the instruments used to improve the quality of the data obtained. ^

Motivational interviewing and feedback for college drinkers: Handling missing values with complete case analysis and multiple imputation

Objective: In this secondary data analysis, three statistical methodologies were implemented to handle cases with missing data in a motivational interviewing and feedback study. The aim was to evaluate the impact that these methodologies have on the data analysis. ^ Methods: We first evaluated whether the assumption of missing completely at random held for this study. We then proceeded to conduct a secondary data analysis using a mixed linear model to handle missing data with three methodologies (a) complete case analysis, (b) multiple imputation with explicit model containing outcome variables, time, and the interaction of time and treatment, and (c) multiple imputation with explicit model containing outcome variables, time, the interaction of time and treatment, and additional covariates (e.g., age, gender, smoke, years in school, marital status, housing, race/ethnicity, and if participants play on athletic team). Several comparisons were conducted including the following ones: 1) the motivation interviewing with feedback group (MIF) vs. the assessment only group (AO), the motivation interviewing group (MIO) vs. AO, and the intervention of the feedback only group (FBO) vs. AO, 2) MIF vs. FBO, and 3) MIF vs. MIO.^ Results: We first evaluated the patterns of missingness in this study, which indicated that about 13% of participants showed monotone missing patterns, and about 3.5% showed non-monotone missing patterns. Then we evaluated the assumption of missing completely at random by Little's missing completely at random (MCAR) test, in which the Chi-Square test statistic was 167.8 with 125 degrees of freedom, and its associated p-value was p=0.006, which indicated that the data could not be assumed to be missing completely at random. After that, we compared if the three different strategies reached the same results. For the comparison between MIF and AO as well as the comparison between MIF and FBO, only the multiple imputation with additional covariates by uncongenial and congenial models reached different results. For the comparison between MIF and MIO, all the methodologies for handling missing values obtained different results. ^ Discussions: The study indicated that, first, missingness was crucial in this study. Second, to understand the assumptions of the model was important since we could not identify if the data were missing at random or missing not at random. Therefore, future researches should focus on exploring more sensitivity analyses under missing not at random assumption.^

(Table 1) Water column biogenic silica concentrations and Si isotopic composition during Marion Dufresne cruise MD166 off South Africa

Southern Ocean biogeochemical processes have an impact on global marine primary production and global elemental cycling, e.g. by likely controlling glacial-interglacial pCO2 variation. In this context, the natural silicon isotopic composition (d30Si) of sedimentary biogenic silica has been used to reconstruct past Si-consumption:supply ratios in the surface waters. We present a new dataset in the Southern Ocean from a IPY-GEOTRACES transect (Bonus-GoodHope) which includes for the first time summer d30Si signatures of suspended biogenic silica (i) for the whole water column at three stations and (ii) in the mixed layer at seven stations from the subtropical zone up to the Weddell Gyre. In general, the isotopic composition of biogenic opal exported to depth was comparable to the opal leaving the mixed layer and did not seem to be affected by any diagenetic processes during settling, even if an effect of biogenic silica dissolution cannot be ruled out in the northern part of the Weddell Gyre. We develop a mechanistic understanding of the processes involved in the modern Si-isotopic balance, by implementing a mixed layer model. We observe that the accumulated biogenic silica (sensu Rayleigh distillation) should satisfactorily describe the d30Si composition of biogenic silica exported out of the mixed layer, within the limit of the current analytical precision on the d30Si. The failures of previous models (Rayleigh and steady state) become apparent especially at the end of the productive period in the mixed layer, when biogenic silica production and export are low. This results from (1) a higher biogenic silica dissolution:production ratio imposing a lower net fractionation factor and (2) a higher Si-supply:Si-uptake ratio supplying light Si-isotopes into the mixed layer. The latter effect is especially expressed when the summer mixed layer becomes strongly Si-depleted, together with a large vertical silicic acid gradient, e.g. in the Polar Front Zone and at the Polar Front.

Modelling seed germination in forest tree species through survival analysis. The Pinus pinea L. case study

The direct application of existing models for seed germination may often be inadequate in the context of ecology and forestry germination experiments. This is because basic model assumptions are violated and variables available to forest managers are rarely used. In this paper, we present a method which addresses the aforementioned shortcomings. The approach is illustrated through a case study of Pinus pinea L. Our findings will also shed light on the role of germination in the general failure of natural regeneration in managed forests of this species. The presented technique consists of a mixed regression model based on survival analysis. Climate and stand covariates were tested. Data for fitting the model were gathered from a 5-year germination experiment in a mature, managed P. pinea stand in the Northern Plateau of Spain in which two different stand densities can be found. The model predictions proved to be unbiased and highly accurate when compared with the training data. Germination in P. pinea was controlled through thermal variables at stand level. At microsite level, low densities negatively affected the probability of germination. A time-lag in the response was also detected. Overall, the proposed technique provides a reliable alternative to germination modelling in ecology/forestry studies by using accessible/ suitable variables. The P. pinea case study highlights the importance of producing unbiased predictions. In this species, the occurrence and timing of germination suggest a very different regeneration strategy from that understood by forest managers until now, which may explain the high failure rate of natural regeneration in managed stands. In addition, these findings provide valuable information for the management of P. pinea under climate-change conditions.

Las nuevas medidas agroenergéticas mexicanas: un modelo microeconométrico para la evaluación de su impacto sobre los ecosistemas naturales en los estados de Hidalgo, Querétaro y Tamaulipas (México)

La ley para la Promoción y Desarrollo de Biocombustibles aprobada en México en 2007 permite la producción de bioetanol y biodiesel. Esta producción puede entrar en conflicto con la producción de alimentos y con los ecosistemas naturales y en esta tesis se desarrolla un modelo microeconométrico que puede servir de base para anticiparse a esos conflictos y para diseñar medidas de política agraria orientadas a potenciar la compatibilidad de la producción de biocombustibles con la de alimentos y con la conservación de los ecosistemas naturales. A partir de una muestra de explotaciones de tres Estados de México – Hidalgo, Querétaro y Tamaulipas- y de un modelo logit multinomial mixto, se estima la elasticidad de la superficie destinada a cultivos alimentarios respecto a cambios en los márgenes económicos de los cultivos agroenergéticos. Esa elasticidad resulta ser significativa. Mostramos que su estimación es útil para anticipar cambios en la superficie destinada a los cultivos alimentarios y a los forestales. Se evalúa el impacto de varios escenarios relativos a los márgenes brutos de los cultivos sobre las decisiones de los agricultores y se muestra la utilidad del modelo para detectar tendencias de cambio a largo plazo en la alternativa de cultivos, incluyendo los forestales. ABSTRACT The Law for the Promotion and Development of Biofuels in Mexico adopted in 2007 allows for the production of bioethanol and biodiesel. This production may conflict with food production and natural ecosystems and this thesis develops a microeconometric model that can serve as a basis to anticipate such conflicts and to implement agricultural policy measures designed to enhance the compatibility of biofuels with production food and natural ecosystems conservation. We estimate the elasticity of the area devoted to food crops with respect to changes in economic margins of energy crops, using a sample of farms in three states of Mexico - Hidalgo, Queretaro and Tamaulipas - , and a multinomial mixed logit model. We found that this elasticity is significant. And we show how it can be useful to anticipate changes in area under food crops and forests. The impact of various scenarios about gross margins on farmers' decisions is assessed and it is shown the usefulness of the model to detect trends of long-term change in the crops area, including forests.

Exploring the role of social capital influence variables on travel behaviour

This paper explores the potential role of individual trip characteristics and social capital network variables in the choice of transport mode. A sample of around 100 individuals living or working in one suburb of Madrid (i.e. Las Rosas district of Madrid) participated in a smartphone short panel survey, entering travel data for an entire working week. A Mixed Logit model was estimated with this data to analyze shifts to metro as a consequence of the opening of two new stations in the area. Apart from classical explanatory variables, such as travel time and cost, gender, license and car ownership, the model incorporated two “social capital network” variables: participation in voluntary activities and receiving help for various tasks (i.e. child care, housekeeping, etc.). Both variables improved the capacity of the model to explain transport mode shifts. Further, our results confirm that the shift towards metro was higher in the case of people “helped” and lower for those participating in some voluntary activities.

A CDKN2-like polymorphism in Xiphophorus LG V is associated with UV-B-induced melanoma formation in platyfish–swordtail hybrids

The genetic basis of spontaneous melanoma formation in spotted dorsal (Sd) Xiphophorus platyfish–swordtail hybrids has been studied for decades, and is adequately explained by a two-gene inheritance model involving a sex-linked oncogene, Xmrk, and an autosomal tumor suppressor, DIFF. The Xmrk oncogene encodes a receptor tyrosine kinase related to EGFR; the nature of the DIFF tumor suppressor gene is unknown. We analyzed the genetic basis of UV-B-induced melanoma formation in closely related, spotted side platyfish–swordtail hybrids, which carry a different sex-linked pigment pattern locus, Sp. We UV-irradiated spotted side Xiphophorus platyfish–swordtail backcross hybrids to induce melanomas at frequencies 6-fold higher than occur spontaneously in unirradiated control animals. To identify genetic determinants of melanoma susceptibility in this UV-inducible Xiphophorus model, we genotyped individual animals from control and UV-irradiated experimental regimes using allozyme and DNA restriction fragment length polymorphisms and tested for joint segregation of genetic markers with pigmentation phenotype and UV-induced melanoma formation. Joint segregation results show linkage of a CDKN2-like DNA polymorphism with UV-B-induced melanoma formation in these hybrids. The CDKN2-like polymorphism maps to Xiphophorus linkage group V and exhibits recombination fractions with ES1 and MDH2 allozyme markers consistent with previous localization of the DIFF tumor suppressor locus. Our results indicate that the CDKN2-like sequence we have cloned and mapped is a candidate for the DIFF tumor suppressor gene.

Avaliação in vivo da efetividade e do pH de géis clareadores no clareamento em consultório em 12 meses de acompanhamento

O objetivo deste estudo in vivo, internacional, randomizado e duplo cedo foi avaliar comparativamente a efetividade e o pH de diferentes géis clareadores na técnica de clareamento em consultório, com e sem o emprego de fonte de luz híbrida em função do grau de alteração de cor, sensibilidade e manutenção do tratamento ao longo de 12 meses de acompanhamento. Foram selecionados 48 voluntários de acordo com os critérios de inclusão e exclusão. Os pacientes foram divididos, de forma randomizada, em 4 grupos de 12 participantes cada, onde: Grupo EXP10 5 aplicações do gel de peróxido de hidrogênio a 10% (Gel Experimental DMC Equipamentos) e ativação de luz híbrida de LED (violeta)/Laser (Experimental DMC Equipamentos) com 7′ e 30″ por aplicação, com tempo total de 37′30; Grupo LP15 5 aplicações do gel de peróxido de hidrogênio 15% (Lase Peroxide Lite DMC Equipamentos) seguindo mesmo protocolo do grupo EXP10; Grupo TB35LH 3 aplicações do gel de peróxido de hidrogênio a 35% (Total Blanc Office - DFL) e ativação de luz híbrida de LED (azul)/Laser (Whitening Lase II DMC Equipamentos) de 7′ e 30″ por aplicação, com tempo total de 22′30″; Grupo TB35 3 aplicações do gel de peróxido de hidrogênio a 35% (Total Blanc Office - DFL) sem ativação com fonte de luz, totalizando 45″. A determinação dos valores de pH foi realizada com o peagômetro digital (Sentron Model 1001, Sentron) nos tempos inicial e após o término do protocolo clareador. A aferição da cor foi feita com espectofotômetro VITA Easyshade antes do clareamento, após 24 horas, 1 semana, 1, 6 e 12 meses. A sensibilidade dentária e grau de satisfação dos pacientes foram avaliados por meio do questionário VAS e IPS antes, imediatamente após o clareamento, 24 horas e uma semana após. Os resultados da alteração do pH receberam tratamento estatístico pela ANOVA e teste de Bonferroni a 0,05%. Os resultados indicaram que o pH aumentou do momento inicial para o final para todos os protocolos. Não houve diferenças significativas entre os protocolos TB35 e TB35LH em nenhum dos momentos, e o pH médio do grupo EXP10 foi maior em comparação aos outros três grupos nos dois momentos avaliados. Os resultados do ΔE receberam tratamento estatístico pela ANOVA e teste de Bonferroni a 0,05%. Os resultados indicaram que não houve diferença significativa entre os grupos LP15, TB35 e TB35LH. O ΔE médio observado após 24 horas foi estatisticamente maior que para os outros tempos (inicial, 1 semana, 1 mês, 6 e 12 meses). Para análise da sensibilidade foi construído um modelo linear misto e atribuídos postos (ranks) aos valores de Δ e teste de Bonferroni a 0,05% para comparações pareadas. Não houve diferença nos valores da sensibilidade imediatamente e 24 horas após o tratamento, com relação ao momento inicial. Houve diferença significativa entre Δ1 e Δ3 indicando que a sensação de dor após uma semana do tratamento foi menor do que as observadas nos instantes imediato e após 24 horas. Para os resultados de satisfação foi construído um modelo linear misto e atribuídos postos (ranks) e o Método de Bonferroni (0,05%) foi utilizado para as comparações pareadas do efeito de tempo. Os resultados indicam queda nos níveis de satisfação entre os períodos imediato e um ano e entre os períodos 24 horas e um ano. Todos os géis clareadores apresentaram mínima variação do pH nos tempos avaliados, entretanto houve um aumento do pH da primeira para a última aplicação em todos os grupos estudados e o grupo EXP10 apresentou os maiores valores de pH seguido do LP15, TB35LH e TB35 apresentaram os valores mais baixos de pH. Os grupos LP15, TB35 e TB35LH apresentaram menor variação da cor ao longo de 12 meses de acompanhamento. O efeito do protocolo clareador não influenciou a sensibilidade dos pacientes e após uma semana a sensibilidade retornaram aos níveis normais. O nível de satisfação dos pacientes foi significativo em relação ao tempo e não aos protocolos clareadores, os pacientes do grupo TB35 mostraram-se mais insatisfeitos ao longo da pesquisa.

Incorporação de informações genômicas no desenvolvimento de índices econômicos para a seleção de bovinos leiteiros

A eficiência econômica da bovinocultura leiteira está relacionada à utilização de animais que apresentem, concomitantemente, bom desempenho quanto à produção, reprodução, saúde e longevidade. Nisto, o índice de seleção configura-se como ferramenta importante ao aumento da lucratividade nesse sistema, visto que permite a seleção de reprodutores para várias características simultaneamente, considerando a relação entre elas bem como a relevância econômica das mesmas. Com a recente disponibilidade de dados genômicos tornou-se ainda possível expandir a abrangência e acurácia dos índices de seleção por meio do aumento do número e qualidade das informações consideradas. Nesse contexto, dois estudos foram desenvolvidos. No primeiro, o objetivo foi estimar parâmetros genéticos e valores genéticos (VG) para características relacionadas à produção e qualidade do leite incluindo-se a informação genômica na avaliação genética. Foram utilizadas medidas de idade ao primeiro parto (IPP), produção de leite (PROD), teor de gordura (GOR), proteína (PROT), lactose, caseína, escore de células somáticas (ECS) e perfil de ácidos graxos de 4.218 vacas bem como os genótipos de 755 vacas para 57.368 polimorfismos de nucleotídeo único (SNP). Os componentes de variância e VG foram obtidos por meio de um modelo misto animal, incluindo-se os efeitos de grupos de contemporâneas, ordem de lactação, dias em lactação e os efeitos aditivo genético, ambiente permanente e residual. Duas abordagens foram desenvolvidas: uma tradicional, na qual a matriz de relacionamentos é baseada no pedigree; e uma genômica, na qual esta matriz é construída combinando-se a informação de pedigree e dos SNP. As herdabilidades variaram de 0,07 a 0,39. As correlações genéticas entre PROD e os componentes do leite variaram entre -0,45 e -0,13 enquanto correlações altas e positivas foram estimadas entre GOR e os ácidos graxos. O uso da abordagem genômica não alterou as estimativas de parâmetros genéticos; contudo, houve aumento entre 1,5% e 6,8% na acurácia dos VG, à exceção de IPP, para a qual houve uma redução de 1,9%. No segundo estudo, o objetivo foi incorporar a informação genômica no desenvolvimento de índices econômicos de seleção. Neste, os VG para PROD, GOR, PROT, teor de ácidos graxos insaturados (INSAT), ECS e vida produtiva foram combinados em índices de seleção ponderados por valores econômicos estimados sob três cenários de pagamento: exclusivamente por volume de leite (PAG1); por volume e por componentes do leite (PAG2); por volume e componentes do leite incluindo INSAT (PAG3). Esses VG foram preditos a partir de fenótipos de 4.293 vacas e genótipos de 755 animais em um modelo multi-característica sob as abordagens tradicional e genômica. O uso da informação genômica influenciou os componentes de variância, VG e a resposta à seleção. Entretanto, as correlações de ranking entre as abordagens foram altas nos três cenários, com valores entre 0,91 e 0,99. Diferenças foram principalmente observadas entre PAG1 e os demais cenários, com correlações entre 0,67 e 0,88. A importância relativa das características e o perfil dos melhores animais foram sensíveis ao cenário de remuneração considerado. Assim, verificou-se como essencial a consideração dos valores econômicos das características na avaliação genética e decisões de seleção.

Expression profiling and sequence diversity of novel DREB genes from common bean (Phaseolus vulgaris L.) and their association with drought-related traits

Common bean is a major dietary component in several countries, but its productivity is negatively affected by abiotic stresses. Dissecting candidate genes involved in abiotic stress tolerance is a paramount step toward the improvement of common bean performance under such constraints. Thereby, this thesis presents a systematic analysis of the DEHYDRATION RESPONSIVE ELEMENT-BINDING (DREB) gene subfamily, which encompasses genes that regulate several processes during stress responses, but with limited information for common bean. First, a series of in silico analyses with sequences retrieved from the P. vulgaris genome on Phytozome supported the categorization of 54 putative PvDREB genes distributed within six phylogenetic subgroups (A-1 to A-6), along the 11 chromosomes. Second, we cloned four novel PvDREB genes and determined their inducibility-factors, including the dehydration-, salinity- and cold-inducible genes PvDREB1F and PvDREB5A, and the dehydration- and cold-inducible genes PvDREB2A and PvDREB6B. Afterwards, nucleotide polymorphisms were searched through Sanger sequencing along those genes, revealing a high number of single nucleotide polymorphisms within PvDREB6B by the comparison of Mesoamerican and Andean genotypes. The nomenclature of PvDREB6B is discussed in details. Furthermore, we used the BARCBean6K_3 SNP platform to identify and genotype the closest SNP to each one of the 54 PvDREB genes. We selected PvDREB6B for a broader study encompassing a collection of wild common bean accessions of Mesoamerican origin. The population structure of the wild beans was accessed using sequence polymorphisms of PvDREB6B. The genetic clusters were partially associated with variation in latitude, altitude, precipitation and temperature throughout the areas such beans are distributed. With an emphasis on drought stress, an adapted tube-screening method in greenhouse conditions enabled the phenotyping of several drought-related traits in the wild collection. Interestingly, our data revealed a correlation between root depth, plant height and biomass and the environmental data of the location of the accessions. Correlation was also observed between the population structure determined through PvDREB6B and the environmental data. An association study combining data from the SNP array and DREB polymorphisms enabled the detection of SNP associated with drought-related traits through a compressed mixed linear model (CMLM) analysis. This thesis highlighted important features of DREB genes in common bean, revealing candidates for further strategies aimed at improvement of abiotic stress tolerance, with emphasis on drought tolerance