Odyssey: a semi-automated pipeline for phasing, imputation, and analysis of genome-wide genetic data

BACKGROUND: Genome imputation, admixture resolution and genome-wide association analyses are timely and computationally intensive processes with many composite and requisite steps. Analysis time increases further when building and installing the run programs required for these analyses. For scientists that may not be as versed in programing language, but want to perform these operations hands on, there is a lengthy learning curve to utilize the vast number of programs available for these analyses. RESULTS: In an effort to streamline the entire process with easy-to-use steps for scientists working with big data, the Odyssey pipeline was developed. Odyssey is a simplified, efficient, semi-automated genome-wide imputation and analysis pipeline, which prepares raw genetic data, performs pre-imputation quality control, phasing, imputation, post-imputation quality control, population stratification analysis, and genome-wide association with statistical data analysis, including result visualization. Odyssey is a pipeline that integrates programs such as PLINK, SHAPEIT, Eagle, IMPUTE, Minimac, and several R packages, to create a seamless, easy-to-use, and modular workflow controlled via a single user-friendly configuration file. Odyssey was built with compatibility in mind, and thus utilizes the Singularity container solution, which can be run on Linux, MacOS, and Windows platforms. It is also easily scalable from a simple desktop to a High-Performance System (HPS). CONCLUSION: Odyssey facilitates efficient and fast genome-wide association analysis automation and can go from raw genetic data to genome: phenome association visualization and analyses results in 3-8 h on average, depending on the input data, choice of programs within the pipeline and available computer resources. Odyssey was built to be flexible, portable, compatible, scalable, and easy to setup. Biologists less familiar with programing can now work hands on with their own big data using this easy-to-use pipeline

3D facial phenotyping by biometric sibling matching used in contemporary genomic methodologies

The analysis of contemporary genomic data typically operates on one-dimensional phenotypic measurements (e.g. standing height). Here we report on a data-driven, family-informed strategy to facial phenotyping that searches for biologically relevant traits and reduces multivariate 3D facial shape variability into amendable univariate measurements, while preserving its structurally complex nature. We performed a biometric identification of siblings in a sample of 424 children, defining 1,048 sib-shared facial traits. Subsequent quantification and analyses in an independent European cohort (n = 8,246) demonstrated significant heritability for a subset of traits (0.17–0.53) and highlighted 218 genome-wide significant loci (38 also study-wide) associated with facial variation shared by siblings. These loci showed preferential enrichment for active chromatin marks in cranial neural crest cells and embryonic craniofacial tissues and several regions harbor putative craniofacial genes, thereby enhancing our knowledge on the genetic architecture of normal-range facial variation

Treatment of allergic rhinitis during and outside the pollen season using mobile technology. A MASK study

Background: The analysis of mobile health (mHealth) data has generated innovative insights into improving allergic rhinitis control, but additive information is needed. A cross-sectional real-world observational study was undertaken in 17 European countries during and outside the estimated pollen season. The aim was to collect novel information including the phenotypic characteristics of the users. Methods: The Allergy Diary–MASK-air–mobile phone app, freely available via Google Play and App, was used to collect the data of daily visual analogue scales (VASs) for overall allergic symptoms and medication use. Fluticasone Furoate (FF), Mometasone Furoate (MF), Azelastine Fluticasone Proprionate combination (MPAzeFlu) and eight oral H1-antihistamines were studied. Phenotypic characteristics were recorded at entry. The ARIA severity score was derived from entry data. This was an a priori planned analysis. Results: 9037 users filled in 70,286 days of VAS in 2016, 2017 and 2018. The ARIA severity score was lower outside than during the pollen season. Severity was similar for all treatment groups during the pollen season, and lower in the MPAzeFlu group outside the pollen season. Days with MPAzeFlu had lower VAS levels and a higher frequency of monotherapy than the other treatments during the season. Outside the season, days with MPAzeFlu also had a higher frequency of monotherapy. The number of reported days was significantly higher with MPAzeFlu during and outside the season than with MF, FF or oral H1-antihistamines. Conclusions: This study shows that the overall efficacy of treatments is similar during and outside the pollen season and indicates that medications are similarly effective during the year

Daily allergic multimorbidity in rhinitis using mobile technology:a novel concept of the MASK study

Background: Multimorbidity in allergic airway diseases is well known, but no data exist about the daily dynamics of symptoms and their impact on work. To better understand this, we aimed to assess the presence and control of daily allergic multimorbidity (asthma, conjunctivitis, rhinitis) and its impact on work productivity using a mobile technology, the Allergy Diary. Methods: We undertook a 1-year prospective observational study in which 4 210 users and 32 585 days were monitored in 19 countries. Five visual analogue scales (VAS) assessed the daily burden of the disease (i.e., global evaluation, nose, eyes, asthma and work). Visual analogue scale levels <20/100 were categorized as "Low" burden and VAS levels ≥50/100 as "High" burden. Results: Visual analogue scales global measured levels assessing the global control of the allergic disease were significantly associated with allergic multimorbidity. Eight hypothesis-driven patterns were defined based on "Low" and "High" VAS levels. There were <0.2% days of Rhinitis Low and Asthma High or Conjunctivitis High patterns. There were 5.9% days with a Rhinitis High-Asthma Low pattern. There were 1.7% days with a Rhinitis High-Asthma High-Conjunctivitis Low pattern. A novel Rhinitis High-Asthma High-Conjunctivitis High pattern was identified in 2.9% days and had the greatest impact on uncontrolled VAS global measured and impaired work productivity. Work productivity was significantly correlated with VAS global measured levels. Conclusions: In a novel approach examining daily symptoms with mobile technology, we found considerable intra-individual variability of allergic multimorbidity including a previously unrecognized extreme pattern of uncontrolled multimorbidity

Geolocation with respect to persona privacy for the Allergy Diary app - a MASK study

Background: Collecting data on the localization of users is a key issue for the MASK (Mobile Airways Sentinel network: the Allergy Diary) App. Data anonymization is a method of sanitization for privacy. The European Commission's Article 29 Working Party stated that geolocation information is personal data. To assess geolocation using the MASK method and to compare two anonymization methods in the MASK database to find an optimal privacy method. Methods: Geolocation was studied for all people who used the Allergy Diary App from December 2015 to November 2017 and who reported medical outcomes. Two different anonymization methods have been evaluated: Noise addition (randomization) and k-anonymity (generalization). Results: Ninety-three thousand one hundred and sixteen days of VAS were collected from 8535 users and 54,500 (58. 5%) were geolocalized, corresponding to 5428 users. Noise addition was found to be less accurate than k-anonymity using MASK data to protect the users' life privacy. Discussion: k-anonymity is an acceptable method for the anonymization of MASK data and results can be used for other databases.Peer reviewe

Genome scans of facial features in East Africans and cross-population comparisons reveal novel associations

Facial morphology is highly variable, both within and among human populations, and a sizable portion of this variation is attributable to genetics. Previous genome scans have revealed more than 100 genetic loci associated with different aspects of normal-range facial variation. Most of these loci have been detected in Europeans, with few studies focusing on other ancestral groups. Consequently, the degree to which facial traits share a common genetic basis across diverse sets of humans remains largely unknown. We therefore investigated the genetic basis of facial morphology in an East African cohort. We applied an open-ended data-driven phenotyping approach to a sample of 2,595 3D facial images collected on Tanzanian children. This approach segments the face into hierarchically arranged, multivariate features that capture the shape variation after adjusting for age, sex, height, weight, facial size and population stratification. Genome scans of these multivariate shape phenotypes revealed significant (p < 2.5 × 10−8) signals at 20 loci, which were enriched for active chromatin elements in human cranial neural crest cells and embryonic craniofacial tissue, consistent with an early developmental origin of the facial variation. Two of these associations were in highly conserved regions showing craniofacial-specific enhancer activity during embryological development (5q31.1 and 12q21.31). Six of the 20 loci surpassed a stricter threshold accounting for multiple phenotypes with study-wide significance (p < 6.25 × 10−10). Cross-population comparisons indicated 10 association signals were shared with Europeans (seven sharing the same associated SNP), and facilitated fine-mapping of causal variants at previously reported loci. Taken together, these results may point to both shared and population-specific components to the genetic architecture of facial variation

The intersection of the genetic architectures of orofacial clefts and normal facial variation

Unaffected relatives of individuals with non-syndromic cleft lip with or without cleft palate (NSCL/P) show distinctive facial features. The presence of this facial endophenotype is potentially an expression of underlying genetic susceptibility to NSCL/P in the larger unselected population. To explore this hypothesis, we first partitioned the face into 63 partially overlapping regions representing global-to-local facial morphology and then defined endophenotypic traits by contrasting the 3D facial images from 264 unaffected parents of individuals with NSCL/P versus 3,171 controls. We observed distinct facial features between parents and controls across 59 global-to-local facial segments at nominal significance (p ≤ 0.05) and 52 segments at Bonferroni corrected significance (p < 1.2 × 10–3), respectively. Next, we quantified these distinct facial features as univariate traits in another dataset of 8,246 unaffected European individuals and performed a genome-wide association study. We identified 29 independent genetic loci that were associated (p < 5 × 10–8) with at least one of the tested endophenotypic traits, and nine genetic loci also passed the study-wide threshold (p < 8.47 × 10–10). Of the 29 loci, 22 were in proximity of loci previously associated with normal facial variation, 18 were near genes that show strong evidence in orofacial clefting (OFC), and another 10 showed some evidence in OFC. Additionally, polygenic risk scores for NSCL/P showed associations with the endophenotypic traits. This study thus supports the hypothesis of a shared genetic architecture of normal facial development and OFC

Cumulative Prognostic Score Predicting Mortality in Patients Older Than 80 Years Admitted to the ICU.

OBJECTIVES: To develop a scoring system model that predicts mortality within 30 days of admission of patients older than 80 years admitted to intensive care units (ICUs). DESIGN: Prospective cohort study. SETTING: A total of 306 ICUs from 24 European countries. PARTICIPANTS: Older adults admitted to European ICUs (N = 3730; median age = 84 years [interquartile range = 81-87 y]; 51.8% male). MEASUREMENTS: Overall, 24 variables available during ICU admission were included as potential predictive variables. Multivariable logistic regression was used to identify independent predictors of 30-day mortality. Model sensitivity, specificity, and accuracy were evaluated with receiver operating characteristic curves. RESULTS: The 30-day-mortality was 1562 (41.9%). In multivariable analysis, these variables were selected as independent predictors of mortality: age, sex, ICU admission diagnosis, Clinical Frailty Scale, Sequential Organ Failure Score, invasive mechanical ventilation, and renal replacement therapy. The discrimination, accuracy, and calibration of the model were good: the area under the curve for a score of 10 or higher was .80, and the Brier score was .18. At a cut point of 10 or higher (75% of all patients), the model predicts 30-day mortality in 91.1% of all patients who die. CONCLUSION: A predictive model of cumulative events predicts 30-day mortality in patients older than 80 years admitted to ICUs. Future studies should include other potential predictor variables including functional status, presence of advance care plans, and assessment of each patient's decision-making capacity

Multimessenger NuEM Alerts with AMON

The Astrophysical Multimessenger Observatory Network (AMON), has developed a real-time multi-messenger alert system. The system performs coincidence analyses of datasets from gamma-ray and neutrino detectors, making the Neutrino-Electromagnetic (NuEM) alert channel. For these analyses, AMON takes advantage of sub-threshold events, i.e., events that by themselves are not significant in the individual detectors. The main purpose of this channel is to search for gamma-ray counterparts of neutrino events. We will describe the different analyses that make-up this channel and present a selection of recent results