Making genetic biodiversity measurable : a review of statistical multivariate methods to study variability at gene level

Measures of agro-ecosystems genetic variability are essential to sustain scientific-based actions and policies tending to protect the ecosystem services they provide. To build the genetic variability datum it is necessary to deal with a large number and different types of variables. Molecular marker data is highly dimensional by nature, and frequently additional types of information are obtained, as morphological and physiological traits. This way, genetic variability studies are usually associated with the measurement of several traits on each entity. Multivariate methods are aimed at finding proximities between entities characterized by multiple traits by summarizing information in few synthetic variables. In this work we discuss and illustrate several multivariate methods used for different purposes to build the datum of genetic variability. We include methods applied in studies for exploring the spatial structure of genetic variability and the association of genetic data to other sources of information. Multivariate techniques allow the pursuit of the genetic variability datum, as a unifying notion that merges concepts of type, abundance and distribution of variability at gene level.

Antiretroviral therapy (ART) adherence among people living with HIV/AIDS (PLHIV) in the north of Vietnam : a multi-method approach

The antiretroviral therapy (ART) program for People Living with HIV/AIDS (PLHIV) in Vietnam has been scaled up rapidly in recent years (from 50 clients in 2003 to almost 38,000 in 2009). ART success is highly dependent on the ability of the patients to fully adhere to the prescribed treatment regimen. Despite the remarkable extension of ART programs in Vietnam, HIV/AIDS program managers still have little reliable data on levels of ART adherence and factors that might promote or reduce adherence. Several previous studies in Vietnam estimated extremely high levels of ART adherence among their samples, although there are reasons to question the veracity of the conclusion that adherence is nearly perfect. Further, no study has quantitatively assessed the factors influencing ART adherence. In order to reduce these gaps, this study was designed to include several phases and used a multi-method approach to examine levels of ART non-adherence and its relationship to a range of demographic, clinical, social and psychological factors. The study began with an exploratory qualitative phase employing four focus group discussions and 30 in-depth interviews with PLHIV, peer educators, carers and health care providers (HCPs). Survey interviews were completed with 615 PLHIV in five rural and urban out-patient clinics in northern Vietnam using an Audio Computer Assisted Self-Interview (ACASI) and clinical records extraction. The survey instrument was carefully developed through a systematic procedure to ensure its reliability and validity. Cultural appropriateness was considered in the design and implementation of both the qualitative study and the cross sectional survey. The qualitative study uncovered several contrary perceptions between health care providers and HIV/AIDS patients regarding the true levels of ART adherence. Health care providers often stated that most of their patients closely adhered to their regimens, while PLHIV and their peers reported that “it is not easy” to do so. The quantitative survey findings supported the PLHIV and their peers’ point of view in the qualitative study, because non-adherence to ART was relatively common among the study sample. Using the ACASI technique, the estimated prevalence of onemonth non-adherence measured by the Visual Analogue Scale (VAS) was 24.9% and the prevalence of four-day not-on-time-adherence using the modified Adult AIDS Clinical Trials Group (AACTG) instrument was 29%. Observed agreement between the two measures was 84% and kappa coefficient was 0.60 (SE=0.04 and p<0.0001). The good agreement between the two measures in the current study is consistent with those found in previous research and provides evidence of cross-validation of the estimated adherence levels. The qualitative study was also valuable in suggesting important variables for the survey conceptual framework and instrument development. The survey confirmed significant correlations between two measures of ART adherence (i.e. dose adherence and time adherence) and many factors identified in the qualitative study, but failed to find evidence of significant correlations of some other factors and ART adherence. Non-adherence to ART was significantly associated with untreated depression, heavy alcohol use, illicit drug use, experiences with medication side-effects, chance health locus of control, low quality of information from HCPs, low satisfaction with received support and poor social connectedness. No multivariate association was observed between ART adherence and age, gender, education, duration of ART, the use of adherence aids, disclosure of ART, patients’ ability to initiate communication with HCPs or distance between clinic and patients’ residence. This is the largest study yet reported in Asia to examine non-adherence to ART and its possible determinants. The evidence strongly supports recent calls from other developing nations for HIV/AIDS services to provide screening, counseling and treatment for patients with depressive symptoms, heavy use of alcohol and substance use. Counseling should also address fatalistic beliefs about chance or luck determining health outcomes. The data suggest that adherence could be enhanced by regularly providing information on ART and assisting patients to maintain social connectedness with their family and the community. This study highlights the benefits of using a multi-method approach in examining complex barriers and facilitators of medication adherence. It also demonstrated the utility of the ACASI interview method to enhance open disclosure by people living with HIV/AIDS and thus, increase the veracity of self-reported data.

A simple method to localise pleiotropic susceptibility loci using univariate linkage analyses of correlated traits

Univariate linkage analysis is used routinely to localise genes for human complex traits. Often, many traits are analysed but the significance of linkage for each trait is not corrected for multiple trait testing, which increases the experiment-wise type-I error rate. In addition, univariate analyses do not realise the full power provided by multivariate data sets. Multivariate linkage is the ideal solution but it is computationally intensive, so genome-wide analysis and evaluation of empirical significance are often prohibitive. We describe two simple methods that efficiently alleviate these caveats by combining P-values from multiple univariate linkage analyses. The first method estimates empirical pointwise and genome-wide significance between one trait and one marker when multiple traits have been tested. It is as robust as an appropriate Bonferroni adjustment, with the advantage that no assumptions are required about the number of independent tests performed. The second method estimates the significance of linkage between multiple traits and one marker and, therefore, it can be used to localise regions that harbour pleiotropic quantitative trait loci (QTL). We show that this method has greater power than individual univariate analyses to detect a pleiotropic QTL across different situations. In addition, when traits are moderately correlated and the QTL influences all traits, it can outperform formal multivariate VC analysis. This approach is computationally feasible for any number of traits and was not affected by the residual correlation between traits. We illustrate the utility of our approach with a genome scan of three asthma traits measured in families with a twin proband.

Multivariate analysis methods in the search for the Higgs boson produced in association with top pairs at the ATLAS experiment at LHC.

The Large Hadron Collider, located at the CERN laboratories in Geneva, is the largest particle accelerator in the world. One of the main research fields at LHC is the study of the Higgs boson, the latest particle discovered at the ATLAS and CMS experiments. Due to the small production cross section for the Higgs boson, only a substantial statistics can offer the chance to study this particle properties. In order to perform these searches it is desirable to avoid the contamination of the signal signature by the number and variety of the background processes produced in pp collisions at LHC. Much account assumes the study of multivariate methods which, compared to the standard cut-based analysis, can enhance the signal selection of a Higgs boson produced in association with a top quark pair through a dileptonic final state (ttH channel). The statistics collected up to 2012 is not sufficient to supply a significant number of ttH events; however, the methods applied in this thesis will provide a powerful tool for the increasing statistics that will be collected during the next LHC data taking.

Sediment physicochemical properties condition pollutant association : a multivariate study of polybrominated diphenyl ethers on the Mendoza River irrigation areas

This work is a multidisciplinary environmental study that provides new insights into the relationships between sediment-organic matter characteristics and polybrominated diphenyl ethers (PBDEs) concentration. The aim of the present multivariate study was to correlate factors influencing PBDEs accumulation in sediment by using principal component analysis (PCA). Organic matter studies by Fourier Transform-Infrared spectroscopy and physicochemical analyses (Total Organic Carbon, pH, electrical conductivity) of sediment samples were considered for PCA. Samples were collected from an artificial irrigation network on the Mendoza River irrigation areas. PCA provided a comprehensive analysis of the studied variables, identifying two components that explained 63% of the data variance. Those factors were mainly associated to organic matter degradation degree, which represent a new insight into the relationships between organic matter in sediments and PBDEs fate. In this sense it was possible to determine that not only the content but also the type of organic matter (chemical structure) could be relevant when evaluating PBDEs accumulation and transport in the environment. Typification of organic matter may be a useful tool to predict more feasible areas where PBDE, may accumulate, as well as sediment transportation capability.

Effects of dopamine receptor D4 variation on alcohol and tobacco use and on novelty seeking: Multivariate linkage and association analysis

The dopamine D4 receptor gene contains a polymorphic sequence consisting of a variable number of 48-base-pair (bp) repeats, and there have been a number of reports that this polymorphism is associated with variation in novelty seeking or in substance abuse and addictive behaviors. In this study we have assessed the linkage and association of DRD4 genotype with novelty seeking, alcohol use, and smoking in a sample of 377 dizygotic twin pairs and 15 single twins recruited from the Australian Twin Registry (ATR). We found no evidence of linkage or association of the DRD4 locus with any of the phenotypes. We made use of repeated measures for some phenotypes to increase power by multivariate genetic analysis, but allelic effects were still non-significant. Specifically, it has been suggested that the DRD4 7-repeat allele is associated with increased novelty seeking in males but we found no evidence for this, despite considerable power to do so. We conclude that DRD4 variation does not have an effect on use of alcohol and the problems that arise from it, on smoking, or on novelty seeking behavior. (C) 2003 Wiley-Liss, Inc.

Characterisation of soils with the use of instrumental techniques : a multivariate forensic study

The value of soil evidence in the forensic discipline is well known. However, it would be advantageous if an in-situ method was available that could record responses from tyre or shoe impressions in ground soil at the crime scene. The development of optical fibres and emerging portable NIR instruments has unveiled a potential methodology which could permit such a proposal. The NIR spectral region contains rich chemical information in the form of overtone and combination bands of the fundamental infrared absorptions and low-energy electronic transitions. This region has in the past, been perceived as being too complex for interpretation and consequently was scarcely utilized. The application of NIR in the forensic discipline is virtually non-existent creating a vacancy for research in this area. NIR spectroscopy has great potential in the forensic discipline as it is simple, nondestructive and capable of rapidly providing information relating to chemical composition. The objective of this study is to investigate the ability of NIR spectroscopy combined with Chemometrics to discriminate between individual soils. A further objective is to apply the NIR process to a simulated forensic scenario where soil transfer occurs. NIR spectra were recorded from twenty-seven soils sampled from the Logan region in South-East Queensland, Australia. A series of three high quartz soils were mixed with three different kaolinites in varying ratios and NIR spectra collected. Spectra were also collected from six soils as the temperature of the soils was ramped from room temperature up to 6000C. Finally, a forensic scenario was simulated where the transferral of ground soil to shoe soles was investigated. Chemometrics methods such as the commonly known Principal Component Analysis (PCA), the less well known fuzzy clustering (FC) and ranking by means of multicriteria decision making (MCDM) methodology were employed to interpret the spectral results. All soils were characterised using Inductively Coupled Plasma Optical Emission Spectroscopy and X-Ray Diffractometry. Results were promising revealing NIR combined with Chemometrics is capable of discriminating between the various soils. Peak assignments were established by comparing the spectra of known minerals with the spectra collected from the soil samples. The temperature dependent NIR analysis confirmed the assignments of the absorptions due to adsorbed and molecular bound water. The relative intensities of the identified NIR absorptions reflected the quantitative XRD and ICP characterisation results. PCA and FC analysis of the raw soils in the initial NIR investigation revealed that the soils were primarily distinguished on the basis of their relative quartz and kaolinte contents, and to a lesser extent on the horizon from which they originated. Furthermore, PCA could distinguish between the three kaolinites used in the study, suggesting that the NIR spectral region was sensitive enough to contain information describing variation within kaolinite itself. The forensic scenario simulation PCA successfully discriminated between the ‘Backyard Soil’ and ‘Melcann® Sand’, as well as the two sampling methods employed. Further PCA exploration revealed that it was possible to distinguish between the various shoes used in the simulation. In addition, it was possible to establish association between specific sampling sites on the shoe with the corresponding site remaining in the impression. The forensic application revealed some limitations of the process relating to moisture content and homogeneity of the soil. These limitations can both be overcome by simple sampling practices and maintaining the original integrity of the soil. The results from the forensic scenario simulation proved that the concept shows great promise in the forensic discipline.

Pollutant concentrations within households in Lao PDR and association with housing characteristics and occupants' activities

The paper presents the results of a study conducted to investigate indoor air quality within residential dwellings in Lao PDR. Results from PM 10, CO, and NO2 measurements inside 167 dwellings in Lao PDR over a five month period (December 2005-April 2006) are discussed as a function of household characteristics and occupant activities. Extremely high PM10 and NO2 concentrations (12 h mean PM10 concentrations 1275 ± 98 μg m-3 and 1183 ± 99 μg m-3 in Vientiane and Bolikhamxay provinces, respectively; 12 h mean NO2 concentrations 1210 ± 94 μg m-3 and 561 ± 45 μg m-3 in Vientiane and Bolikhamxay, respectively) were measured within the dwellings. Correlations, ANOVA analysis (univariate and multivariate), and linear regression results suggest a substantial contribution from cookingandsmoking.The PM10 concentrations were significantly higher in houses without a chimney compared to houses in which cooking occurred on a stove with a chimney. However, no significant differences in pollutantconcentrations were observed as a function of cooking location. Furthermore, PM10 and NO2 concentrations were higher in houses in which smoking occurred, suggestive of a relationship between increased indoor concentrations and smoking (0.05 < p < 0.10). Resuspension of dust from soil floors was another significant source of PM10 inside the house (634 μg m-3, p < 0.05).

Regression and hypothesis tests for multivariate GNSS state time series

A satellite based observation system can continuously or repeatedly generate a user state vector time series that may contain useful information. One typical example is the collection of International GNSS Services (IGS) station daily and weekly combined solutions. Another example is the epoch-by-epoch kinematic position time series of a receiver derived by a GPS real time kinematic (RTK) technique. Although some multivariate analysis techniques have been adopted to assess the noise characteristics of multivariate state time series, statistic testings are limited to univariate time series. After review of frequently used hypotheses test statistics in univariate analysis of GNSS state time series, the paper presents a number of T-squared multivariate analysis statistics for use in the analysis of multivariate GNSS state time series. These T-squared test statistics have taken the correlation between coordinate components into account, which is neglected in univariate analysis. Numerical analysis was conducted with the multi-year time series of an IGS station to schematically demonstrate the results from the multivariate hypothesis testing in comparison with the univariate hypothesis testing results. The results have demonstrated that, in general, the testing for multivariate mean shifts and outliers tends to reject less data samples than the testing for univariate mean shifts and outliers under the same confidence level. It is noted that neither univariate nor multivariate data analysis methods are intended to replace physical analysis. Instead, these should be treated as complementary statistical methods for a prior or posteriori investigations. Physical analysis is necessary subsequently to refine and interpret the results.

The association between pterygium and conjunctival ultraviolet autofluorescence : The Norfolk Island Eye Study

Purpose: To investigate the association between conjunctival ultraviolet autofluorescence (UVAF), a biomarker of ocular ultraviolet radiation (UVR) exposure, and prevalent pterygium. Methods: We conducted a cross-sectional study on Norfolk Island, South Pacific. All permanent residents aged ‡15 were invited to participate. Participants completed a sun exposure questionnaire and underwent autorefraction and slit lamp biomicroscope examination. Area of conjunctival UVAF (sum of temporal ⁄ nasal area in right and left eyes) was determined using computerized methods. Multivariate logistic and linear regression models were used to estimate the associations with pterygia and UVAF, respectively. Results: Of 641 participants, 70 people (10.9%) had pterygium in one or both eyes, and prevalence was higher in males (15.0% versus 7.7%, p = 0.003). Significant independent associations with pterygium in any eye were UVAF (per 10 mm2) [odds ratio (OR) 1.16, 95% confidence interval (CI) 1.16–1.28, p = 0.002], tanning skin phenotype (OR 2.17,1.20–3.92, p = 0.010) and spending more than three-quarters of the day outside (OR 2.22, 1.20–4.09, p = 0.011). Increasing quartile of UVAF was associated with increased risk of pterygium following adjustment of age, sex and time outdoors (pTrend = 0.002). Independent associations with increasing UVAF (per 10 mm2) were decreasing age, time outdoors, skin type and male gender (all p < 0.001). UVAF area correlated well with the duration of outdoor activity (pTrend < 0.001). Conclusion: Pterygium occurs in approximately one-tenth of Norfolk Islanders. Increasing conjunctival UVAF is associated with prevalent pterygia, confirming earlier epidemiological, laboratory and ray-tracing studies that pterygia are associated with UVR. Protection from the sun should be encouraged to reduce the prevalence of pterygium in the community.

Association between xenobiotic gene polymorphisms and non-Hodgkin's lymphoma risk

The last four decades have seen a significant increase in the incidence of non-Hodgkin's lymphoma (NHL) as a possible result of increasing environmental carcinogen exposure, particularly pesticides and solvents. Based on the increasing evidence for an association between carcinogen exposure-related cancer risk and xenobiotic gene polymorphisms, we have undertaken a case-control study of xenobiotic gene polymorphisms in individuals with a diagnosis of NHL. Polymorphisms of six xenobiotic genes (CYP1A1, GSTT1, GSTM1, PON1, NAT1, NAT2) were characterized in 169 individuals with NHL and 205 normal controls using polymerase chain reaction-based methods. Polymorphic frequencies were compared using Fisher's exact tests, and odds ratios for NHL risk were calculated. Among the NHL group, the incidence of GSTT1 null and PON1 BB genotypes were significantly increased compared with controls, 34% vs 14%, and 24% vs 11% respectively. Adjusted odds ratios calculated from multivariate analyses demonstrated that GSTT1 null conferred a fourfold increase in NHL risk (OR = 4.27; 95% CI, 2.40-7.61, P < 0.001) and PON1 BB a 2.9-fold increase (OR = 2.92; 95% CI, 1.49-5.72, P = 0.002). Furthermore, GSTT1 null combined with PON1 BB or GSTM1 null conferred an additional risk of NHL. This is the first time that a PON1 gene polymorphism has been shown to be associated with cancer risk. We conclude that the two polymorphisms, GSTT1 null and PON1 BB, are common genetic traits that pose low individual risk but may be important determinants of overall population NHL risk, particularly among groups exposed to NHL-related carcinogens.

Bivariate genome-wide association study of genetically correlated neuroimaging phenotypes from DTI and MRI through a seemingly unrelated regression model

Large multisite efforts (e.g., the ENIGMA Consortium), have shown that neuroimaging traits including tract integrity (from DTI fractional anisotropy, FA) and subcortical volumes (from T1-weighted scans) are highly heritable and promising phenotypes for discovering genetic variants associated with brain structure. However, genetic correlations (rg) among measures from these different modalities for mapping the human genome to the brain remain unknown. Discovering these correlations can help map genetic and neuroanatomical pathways implicated in development and inherited risk for disease. We use structural equation models and a twin design to find rg between pairs of phenotypes extracted from DTI and MRI scans. When controlling for intracranial volume, the caudate as well as related measures from the limbic system - hippocampal volume - showed high rg with the cingulum FA. Using an unrelated sample and a Seemingly Unrelated Regression model for bivariate analysis of this connection, we show that a multivariate GWAS approach may be more promising for genetic discovery than a univariate approach applied to each trait separately.

Psychological resources, their social antecedents, and association with well-being and health behaviour in early adulthood

The success of entering work life, young people s psychological resources and self-reported well-being were studied in a longitudinal setting from a life-span developmental-contextual perspective in early adulthood. The aim was to analyse how psychosocial characteristics in early childhood and adolescence predict successful entrance into work life, how this is associated with well-being, and to assess the level of psychological resources such as dispositional optimism, personal meaning of work and coping in early adulthood. The role of these and social support, in the relationship between regional factors (such as place of residence and migration), self-reported health and life satisfaction was studied. The association between a specific coping strategy, i.e. eating and drinking in a stressful situation and eating habits, was studied to demonstrate how coping is associated with health behaviour. Multivariate methods, including binary logistic regression analyses and ANOVA, were used for statistical analyses. The subjects were members of the Northern Finland 1966 Birth Cohort, which consists of all women and men born in 1966 in the two northernmost provinces of Finland (n= 12,058). The most recent follow-up, at the age of 31 years when 11,637 subjects were alive, took place in 1997-1998. The results show, first, that social resources in the childhood family and adolescence school achievement predict entrance into the labour market. Secondly, psychosocial resources were found to mediate the relationship between migration from rural to urban areas, and subjective well-being. Thirdly, psychological resources at entrance into the labour market were found to develop from early infancy on. They are, however, influenced later by work history. Fourthly, stress-related eating and drinking, as a way of coping, was found to be directly associated with unhealthy eating habits and alcohol use. Gender differences were found in psychosocial resources predicting, and being associated with success in entering the labour market. For men, the role of attitudinal and psychological factors seems to be especially important in entrance into work life and in the development of psychological resources. For women, academic attainment was more important for successfully entering work life, and lack of emotional social support was a risk factor for stress-related eating only among women. Stress-related eating and drinking habits were predicted by a long history of unemployment as well as a low level of education among both genders, but not excluding an academic degree among men. The results emphasize the role of childhood psychosocial factors in preventing long-term unemployment and in enhancing psychological well-being in early adulthood. Success in entering work life, in terms of continuous work history, plays a crucial role for well-being and the amount of psychological resources in early adulthood. The results emphasize the crucial role of enhancing psychological resources for promoting positive health behaviour and diminishing regional differences in subjective well-being.

A versatile gene-based test for genome-wide association studies

We have derived a versatile gene-based test for genome-wide association studies (GWAS). Our approach, called VEGAS (versatile gene-based association study), is applicable to all GWAS designs, including family-based GWAS, meta-analyses of GWAS on the basis of summary data, and DNA-pooling-based GWAS, where existing approaches based on permutation are not possible, as well as singleton data, where they are. The test incorporates information from a full set of markers (or a defined subset) within a gene and accounts for linkage disequilibrium between markers by using simulations from the multivariate normal distribution. We show that for an association study using singletons, our approach produces results equivalent to those obtained via permutation in a fraction of the computation time. We demonstrate proof-of-principle by using the gene-based test to replicate several genes known to be associated on the basis of results from a family-based GWAS for height in 11,536 individuals and a DNA-pooling-based GWAS for melanoma in approximately 1300 cases and controls. Our method has the potential to identify novel associated genes; provide a basis for selecting SNPs for replication; and be directly used in network (pathway) approaches that require per-gene association test statistics. We have implemented the approach in both an easy-to-use web interface, which only requires the uploading of markers with their association p-values, and a separate downloadable application.

The association of chronic kidney disease and dialysis treatment with foot ulceration and major amputation

Objective The objective of this study was to investigate the risk of chronic kidney disease (CKD) stage 4-5 and dialysis treatment on incidence of foot ulceration and major lower extremity amputation in comparison to CKD stage 3. Methods In this retrospective study, all individuals who visited our hospital between 2006 and 2012 because of CKD stages 3 to 5 or dialysis treatment were included. Medical records were reviewed for incidence of foot ulceration and major amputation. The time from CKD 3, CKD 4-5, and dialysis treatment until first foot ulceration and first major lower extremity amputation was calculated and analyzed by Kaplan-Meier curves and multivariate Cox proportional hazards model. Diabetes mellitus, peripheral arterial disease, peripheral neuropathy, and foot deformities were included for potential confounding. Results A total of 669 individuals were included: 539 in CKD 3, 540 in CKD 4-5, and 259 in dialysis treatment (individuals could progress from one group to the next). Unadjusted foot ulcer incidence rates per 1000 patients per year were 12 for CKD 3, 47 for CKD 4-5, and 104 for dialysis (P < .001). In multivariate analyses, the hazard ratio for incidence of foot ulceration was 4.0 (95% confidence interval [CI], 2.6-6.3) in CKD 4-5 and 7.6 (95% CI, 4.8-12.1) in dialysis treatment compared with CKD 3. Hazard ratios for incidence of major amputation were 9.5 (95% CI, 2.1-43.0) and 15 (95% CI, 3.3-71.0), respectively. Conclusions CKD 4-5 and dialysis treatment are independent risk factors for foot ulceration and major amputation compared with CKD 3. Maximum effort is needed in daily clinical practice to prevent foot ulcers and their devastating consequences in all individuals with CKD 4-5 or dialysis treatment.