Search

Ecological communities can be stable over multiple generations, or rapidly shift into structurally and functionally different configurations. In kelp forest ecosystems, overgrazing by sea urchins can abruptly shift forests into alternative states that are void of macroalgae and primarily dominated by actively grazing sea urchins. Beginning in 2014, a sea urchin outbreak along the central coast of California resulted in a patchy mosaic of remnant forests interspersed with sea urchin barrens. In this study, we used a 14-year subtidal monitoring dataset of invertebrates, algae, and fishes to explore changes in community structure associated with the loss of forests. We found that the spatial mosaic of barrens and forests resulted in a region-wide shift in community structure. However, the magnitude of kelp forest loss and taxonomic-level consequences were spatially heterogeneous. Taxonomic diversity declined across the region, but there were no declines in richness for any group, suggesting compositional redistribution. Baseline ecological and environmental conditions, and sea urchin behaviour, explained the persistence of forests through multiple stressors. These results indicate that spatial heterogeneity in preexisting ecological and environmental conditions can explain patterns of community change.

Disparities in data underlying clinical genomic interpretation is an acknowledged problem, but there is a paucity of data demonstrating it. The All of Us Research Program is collecting data including whole-genome sequences, health records, and surveys for at least a million participants with diverse ancestry and access to healthcare, representing one of the largest biomedical research repositories of its kind. Here, we examine pathogenic and likely pathogenic variants that were identified in the All of Us cohort. The European ancestry subgroup showed the highest overall rate of pathogenic variation, with 2.26% of participants having a pathogenic variant. Other ancestry groups had lower rates of pathogenic variation, including 1.62% for the African ancestry group and 1.32% in the Latino/Admixed American ancestry group. Pathogenic variants were most frequently observed in genes related to Breast/Ovarian Cancer or Hypercholesterolemia. Variant frequencies in many genes were consistent with the data from the public gnomAD database, with some notable exceptions resolved using gnomAD subsets. Differences in pathogenic variant frequency observed between ancestral groups generally indicate biases of ascertainment of knowledge about those variants, but some deviations may be indicative of differences in disease prevalence. This work will allow targeted precision medicine efforts at revealed disparities.

Coronary artery calcification (CAC) is a measure of atherosclerosis and a well-established predictor of coronary artery disease (CAD) events. Here we describe a genome-wide association study (GWAS) of CAC in 22,400 participants from multiple ancestral groups. We confirmed associations with four known loci and identified two additional loci associated with CAC (ARSE and MMP16), with evidence of significant associations in replication analyses for both novel loci. Functional assays of ARSE and MMP16 in human vascular smooth muscle cells (VSMCs) demonstrate that ARSE is a promoter of VSMC calcification and VSMC phenotype switching from a contractile to a calcifying or osteogenic phenotype. Furthermore, we show that the association of variants near ARSE with reduced CAC is likely explained by reduced ARSE expression with the G allele of enhancer variant rs5982944. Our study highlights ARSE as an important contributor to atherosclerotic vascular calcification, and a potential drug target for vascular calcific disease.

Multi-omics datasets are becoming more common, necessitating better integration methods to realize their revolutionary potential. Here, we introduce multi-set correlation and factor analysis (MCFA), an unsupervised integration method tailored to the unique challenges of high-dimensional genomics data that enables fast inference of shared and private factors. We used MCFA to integrate methylation markers, protein expression, RNA expression, and metabolite levels in 614 diverse samples from the Trans-Omics for Precision Medicine/Multi-Ethnic Study of Atherosclerosis multi-omics pilot. Samples cluster strongly by ancestry in the shared space, even in the absence of genetic information, while private spaces frequently capture dataset-specific technical variation. Finally, we integrated genetic data by conducting a genome-wide association study (GWAS) of our inferred factors, observing that several factors are enriched for GWAS hits and trans-expression quantitative trait loci. Two of these factors appear to be related to metabolic disease. Our study provides a foundation and framework for further integrative analysis of ever larger multi-modal genomic datasets.

BackgroundAflatoxin (AF) exposure is associated with child growth faltering in cross-sectional studies, with limited findings from longitudinal studies.ObjectivesTo evaluate the relationship between maternal AF B1-lysine adduct concentration, child AF B1-lysine adduct concentration, and child growth in the first 30 mo of life.MethodsAF B1-lysine adduct was measured in mother-child dyad plasma samples using isotope dilution mass spectrometry. Using linear regression, we assessed the relationship between AF B1-lysine adduct concentration and child weight, height, and head and mid-upper arm circumferences at 1 wk, 6, 12, 18, 24, and 30 mo of age.ResultsIn adjusted models, maternal prenatal AF B1-lysine adduct (pg/μL) was positively associated with newborn anthropometric outcomes; largest beta coefficients for associations between standardized values were for newborn weight-for-age z-score [β = 0.13; 95% confidence interval (CI): 0.02, 0.24; P < 0.05 and β = 0.11; 95% CI: 0.00, 0.22; P < 0.05 for second and third trimester AF, respectively]. Child AF B1-lysine adduct (pg/μL) at 6 mo was negatively associated with head circumference-for-age z-score at 6, 18, 24, and 30 mo, with beta coefficients ranging from β = -0.15; 95% CI: -0.28, -0.02 to β = -0.17; 95% CI: -0.31, -0.03; P < 0.05); 18-mo AF was negatively associated with anthropometric outcomes at 18, 24, and 30 mo, most consistently with length-for-age z-score (β = -0.18; 95% CI: -0.32, -0.04, β = -0.21; 95% CI: -0.35, -0.07, β = -0.18; 95% CI: -0.32, -0.03 at 18, 24 and 30 mo, respectively).ConclusionsChild AF exposure was associated with impaired child growth, but maternal AF exposure was not. Exposure during infancy was linked to persistent deficit in head circumference, implying reduced brain size lasting beyond the age of 2 years. Exposure at 18 mo was linked to persistent linear growth deficit. Further research should elucidate mechanisms through which AF affects child growth.

Background/objectivesObesity, defined as excessive fat accumulation that represents a health risk, is increasing in adults and children, reaching global epidemic proportions. Body mass index (BMI) correlates with body fat and future health risk, yet differs in prediction by fat distribution, across populations and by age. Nonetheless, few genetic studies of BMI have been conducted in ancestrally diverse populations. Gene expression association with BMI was assessed in the Multi-Ethnic Study of Atherosclerosis (MESA) in four self-identified race and ethnicity (SIRE) groups to identify genes associated with obesity.Subjects/methodsRNA-sequencing was performed on 1096 MESA participants (37.8% white, 24.3% Hispanic, 28.4% African American, and 9.5% Chinese American) and linear models were used to assess the association of expression from each gene for its effect on BMI, adjusting for age, sex, sequencing center, study site, five expression and four genetic principal components in each self-identified race group. Sample-size-weighted meta-analysis was performed to identify genes with BMI-associated expression across ancestry groups.ResultsWithin individual SIRE groups, there were zero to three genes whose expression is significantly (p

Sleep Disordered Breathing (SDB) is a common disease associated with increased risk for cardiometabolic, cardiovascular, and cognitive diseases. How SDB affects the molecular environment is still poorly understood. We study the association of three SDB measures with gene expression measured using RNA-seq in multiple blood tissues from the Multi-Ethnic Study of Atherosclerosis. We develop genetic instrumental variables for the associated transcripts as polygenic risk scores (tPRS), then generalize and validate the tPRS in the Women's Health Initiative. We measure the associations of the validated tPRS with SDB and serum metabolites in Hispanic Community Health Study/Study of Latinos. Here we find differential gene expression by blood cell type in relation to SDB traits and link P2XR4 expression to average oxyhemoglobin saturation during sleep and butyrylcarnitine (C4) levels. These findings can be used to develop interventions to alleviate the effect of SDB on the human molecular environment.

OBJECTIVES: In an evolving era of immunotherapeutic options for persistent or recurrent laryngeal squamous cell carcinoma (LSCC), there is a need for improved biomarkers of treatment response and survival to inform optimal treatment selection and prognostication. Herein, our primary objective was to explore correlations between tumor infiltrating lymphocytes (TILs) and PD-L1 Combined Positive Score (CPS). Secondarily, we sought to explore their combined association with survival outcomes in patients with persistent or recurrent LSCC treated with salvage surgery. MATERIALS AND METHODS: This was a retrospective cohort study at a single academic medical center. Immunohistochemistry staining for TILs and PD-L1 was performed on a tissue microarray of persistent or recurrent LSCC pathologic specimens. Correlations between TIL subsets and PD-L1 CPS were examined using Pearsons correlation coefficient and survival outcomes were analyzed with the Kaplan-Meier method and log-rank tests. RESULTS: Only CD103+ TILs showed a statistically significant, weakly-positive correlation with PD-L1 CPS (r2 = 0.264, p

ObjectiveBioselection to assess tumor response after induction chemotherapy has been introduced as an alternative treatment strategy to total laryngectomy for patients with advanced larynx squamous cell carcinoma (LSCC). Tumor-infiltrating lymphocytes (TILs) have proven to serve as prognostic biomarkers in head and neck cancer but have not been evaluated as a way to select patients for treatment paradigms. The aim of this study is to evaluate the role of pretreatment TILs in patients with advanced LSCC undergoing the bioselection paradigm.Study designRetrospective study.SettingTertiary care hospital.MethodsPatients with advanced LSCC treated with bioselection and available tissue were included (N = 76). Patients were stratified into CD8-low and CD8-high cohorts by using the median TIL count. Kaplan-Meier survival analysis and multivariate cox regression were performed with SPSS version 26 (IBM).ResultsAfter controlling for tobacco use, tumor site, and stage, a high CD8 TIL count was an independent predictor of improved 5-year disease-specific survival (hazard ratio, 0.17 [95% CI, 0.03-0.84]; P = .03). CD8 TIL counts did not predict response to induction chemotherapy; however, subgroup analysis of patients treated with chemoradiation therapy revealed that CD8 TIL count was significantly associated with degree of response (P = .012).ConclusionThese findings support prior data published by our group showing that TILs are predictive of disease-specific survival in patients with head and neck cancer. CD8 TIL counts were significantly associated with degree of clinical response after induction chemotherapy. These results suggest that pretreatment assessment of tumor-infiltrating CD8 cells could be useful in selecting patients.

BackgroundIn utero or early-life exposure to aflatoxin, which contaminates staple crops in disadvantaged settings, may compromise pregnancy and infant outcomes, but investigations into the extent, persistence, and determinants of aflatoxin exposure at these life stages have lacked longitudinal data collection and broad geographic representation.ObjectivesAflatoxin exposure and selected determinants thereof were characterized in mother-child dyads with serial plasma/serum samples in prenatal, perinatal, and early life in Malawi and Bangladesh.MethodsCirculating aflatoxin B1 (AFB1)-lysine albumin adducts were measured in dyads from Bangladesh (n = 573; maternal first and third trimester, 3 mo postpartum, cord blood, infant 24 mo) and Malawi (n = 255; maternal second and third trimester, 6 mo postpartum, infant 6 and 18 mo) with isotope dilution mass spectrometry. We examined AFB1-lysine adduct magnitude, persistence, seasonality, and associations with infant feeding, and estimated daily AFB1 intake.ResultsMaternal AFB1-lysine was higher in Malawi (98% detectable; median: 0.469, IQR: 0.225-1.027 pg/µL) than in Bangladesh (59%; 0.030, nondetectable [nd]-0.077 pg/µL). Although estimated dietary exposure in Malawi was temporally stable (648 ng AFB1/day), estimated intake in Bangladesh was reduced by 94% between rainy and winter seasons (98 to 6 ng/day). AFB1-lysine was low in cord blood from Bangladesh (15% detectable; 0.045, 0.031-0.088 pg/µL among detectable) and in Malawian infants at 6 mo of age (0.072, nd-0.236 pg/µL), but reached maternal concentrations by 18 or 24 mo (Bangladesh: 0.034, nd-0.063 pg/µL; Malawi: 0.370, 0.195-0.964 pg/µL). In Malawian infants, exclusive breastfeeding at 3 mo was associated with 58% lower AFB1-lysine concentrations at 6 mo compared with other feeding modes (P = 0.010).ConclusionsAmong pregnant women, aflatoxin exposure was persistently high in Malawi, while lower and seasonal in Bangladesh. Infants were partially protected from exposure in utero and with exclusive breastfeeding, but exposures reached adult levels by 18-24 mo of age. The Bangladesh and Malawi trials are registered at clinicaltrials.gov as NCT00860470 and NCT01239693.

In vitro evolution and whole genome analysis were used to comprehensively identify the genetic determinants of chemical resistance in Saccharomyces cerevisiae. Sequence analysis identified many genes contributing to the resistance phenotype as well as numerous amino acids in potential targets that may play a role in compound binding. Our work shows that compound-target pairs can be conserved across multiple species. The set of 25 most frequently mutated genes was enriched for transcription factors, and for almost 25 percent of the compounds, resistance was mediated by one of 100 independently derived, gain-of-function SNVs found in a 170 amino acid domain in the two Zn2C6 transcription factors YRR1 and YRM1 (p

Obstructive sleep apnea (OSA) is a common disorder associated with increased risk of cardiovascular disease and mortality. Iron and heme metabolism, implicated in ventilatory control and OSA comorbidities, was associated with OSA phenotypes in recent admixture mapping and gene enrichment analyses. However, its causal contribution was unclear. In this study, we performed pathway-level transcriptional Mendelian randomization (MR) analysis to investigate the causal relationships between iron and heme related pathways and OSA. In primary analysis, we examined the expression level of four iron/heme Reactome pathways as exposures and four OSA traits as outcomes using cross-tissue cis-eQTLs from the Genotype-Tissue Expression portal and published genome-wide summary statistics of OSA. We identify a significant putative causal association between up-regulated heme biosynthesis pathway with higher sleep time percentage of hypoxemia (p = 6.14 × 10-3). This association is supported by consistency of point estimates in one-sample MR in the Multi-Ethnic Study of Atherosclerosis using high coverage DNA and RNA sequencing data generated by the Trans-Omics for Precision Medicine project. Secondary analysis for 37 additional iron/heme Gene Ontology pathways did not reveal any significant causal associations. This study suggests a causal association between increased heme biosynthesis and OSA severity.

Large datasets of hundreds to thousands of individuals measuring RNA-seq in observational studies are becoming available. Many popular software packages for analysis of RNA-seq data were constructed to study differences in expression signatures in an experimental design with well-defined conditions (exposures). In contrast, observational studies may have varying levels of confounding transcript-exposure associations; further, exposure measures may vary from discrete (exposed, yes/no) to continuous (levels of exposure), with non-normal distributions of exposure. We compare popular software for gene expression-DESeq2, edgeR and limma-as well as linear regression-based analyses for studying the association of continuous exposures with RNA-seq. We developed a computation pipeline that includes transformation, filtering and generation of empirical null distribution of association P-values, and we apply the pipeline to compute empirical P-values with multiple testing correction. We employ a resampling approach that allows for assessment of false positive detection across methods, power comparison and the computation of quantile empirical P-values. The results suggest that linear regression methods are substantially faster with better control of false detections than other methods, even with the resampling method to compute empirical P-values. We provide the proposed pipeline with fast algorithms in an R package Olivia, and implemented it to study the associations of measures of sleep disordered breathing with RNA-seq in peripheral blood mononuclear cells in participants from the Multi-Ethnic Study of Atherosclerosis.

BACKGROUND:Indications for and efficacy of paratracheal nodal dissection (PTND) in patients undergoing laryngectomy (salvage) for persistent or recurrent laryngeal squamous cell carcinoma are not well-defined. METHODS:A retrospective cohort study was performed for patients undergoing salvage laryngectomy with clinically and radiographically negative neck disease between 1998 and 2015 (n = 210). Univariate and multivariate Cox regression analyses were performed. RESULTS:PTND was performed on 77/210 patients (36%). The PTND cohort had a greater proportion of advanced T classification (rT3/rT4) tumors (78%) than subjects without PTND (55%; p = 0.001). There was a 14% rate of occult nodal metastases in the paratracheal basin; of these, 55% did not have pathologic lateral neck disease. Multivariate analysis controlling for tumor site, tumor stage, and pathologic lateral neck disease demonstrated that PTND was associated with improved overall survival [OS] (p = 0.03; hazard ratio [HR] 0.60, 95% confidence interval [CI] 0.38-0.96), disease-free survival [DFS] (p = 0.03; HR 0.55, 95% CI 0.31-0.96), and distant DFS survival (p = 0.01; HR 0.29, 95% CI 0.11-0.77). The rate of hypocalcemia did not differ between subjects who underwent bilateral PTND, unilateral PTND, or no PTND (p = 0.19 at discharge, p = 0.17 at last follow-up). CONCLUSIONS:PTND at the time of salvage laryngectomy was more common in patients with rT3/rT4 tumors and was associated with improved OS and DFS, with no effect on hypocalcemia. In patients undergoing PTND, the finding of occult paratracheal metastases was often independent of lateral neck metastases.

ObjectivesTo report our institutional experience, management, and outcomes of cutaneous periauricular squamous cell carcinoma (SCC).Study designRetrospective chart review.SettingTertiary academic center.SubjectsPatients undergoing treatment of cutaneous periauricular SCC from 2000 to 2016.ResultsA total of 112 patients had a median follow-up of 24.5 months, a mean ± SD age of 75.7 ± 10.6 years, and a strong male predominance (93.8%). Site distribution shows 87 (77.7%) auricular, 26 (23.2%) preauricular, and 10 (8.8%) postauricular lesions. Of auricular lesions, tumors involved the tragus (n = 3, 3.4%), helix/antihelix (n = 47, 54.0%), conchal bowl (n = 31, 35.6%), external auditory canal (n = 18, 16.1%), and lobule (n = 3, 3.4%). Most patients presented at stage I (52.7%) versus stages II (28.6%), III (6.3%), and IV (12.5%). Patients were largely treated surgically with primary tumor resection ranging from wide local excision to lateral temporal bone resection (± parotidectomy and neck dissection), with 17.0% and 5.4% receiving adjuvant radiation and chemoradiation, respectively. Metastatic spread was seen to the parotid (25.9%) and neck (26.8%), with most common cervical spread to level II. Overall survival, disease-specific survival, and disease-free survival at 3 years were 62%, 89%, and 56%, respectively. Nodal disease was associated with worse disease-specific survival (P < .001) and disease-free survival (P = .042). Pre- and postauricular sites were associated with worse overall survival (P = .007) relative to auricular sites.ConclusionAmong cutaneous SCC, periauricular subsites pose treatment challenges related to surrounding anatomy and represent a unique tumor population. The reported propensity toward recurrence and patterns of metastasis may better guide treatment of aggressive tumors to include regional nodal dissection.

To develop a map of cell-cell communication mediated by extracellular RNA (exRNA), the NIH Extracellular RNA Communication Consortium created the exRNA Atlas resource (https://exrna-atlas.org). The Atlas version 4P1 hosts 5,309 exRNA-seq and exRNA qPCR profiles from 19 studies and a suite of analysis and visualization tools. To analyze variation between profiles, we apply computational deconvolution. The analysis leads to a model with six exRNA cargo types (CT1, CT2, CT3A, CT3B, CT3C, CT4), each detectable in multiple biofluids (serum, plasma, CSF, saliva, urine). Five of the cargo types associate with known vesicular and non-vesicular (lipoprotein and ribonucleoprotein) exRNA carriers. To validate utility of this model, we re-analyze an exercise response study by deconvolution to identify physiologically relevant response pathways that were not detected previously. To enable wide application of this model, as part of the exRNA Atlas resource, we provide tools for deconvolution and analysis of user-provided case-control studies.

BackgroundRecurrent laryngeal squamous cell carcinomas (LSCCs) are associated with poor outcomes, without reliable biomarkers to identify patients who may benefit from adjuvant therapies. Given the emergence of tumor-infiltrating lymphocytes (TIL) as a biomarker in head and neck squamous cell carcinoma, we generated predictive models to understand the utility of CD4+, CD8+ and/or CD103+ TIL status in patients with advanced LSCC.MethodsTissue microarrays were constructed from salvage laryngectomy specimens of 183 patients with recurrent/persistent LSCC and independently stained for CD4+, CD8+, and CD103+ TIL content. Cox proportional hazards regression analysis was employed to assess combinations of CD4+, CD8+, and CD103+ TIL levels for prediction of overall survival (OS), disease-specific survival (DSS), and disease-free survival (DFS) in patients with recurrent/persistent LSCC.ResultsHigh tumor CD103+ TIL content was associated with significantly improved OS, DSS, and DFS and was a stronger predictor of survival in recurrent/persistent LSCC than either high CD8+ or CD4+ TIL content. On multivariate analysis, an "immune-rich" phenotype, in which tumors were enriched for both CD103+ and CD4+ TILs, conferred a survival benefit (OS hazard ratio: 0.28, p = 0.0014; DSS hazard ratio: 0.09, p = 0.0015; DFS hazard ratio: 0.18, p = 0.0018) in recurrent/persistent LSCC.ConclusionsAn immune profile driven by CD103+ TIL content, alone and in combination with CD4+ TIL content, is a prognostic biomarker of survival in patients with recurrent/persistent LSCC. Predictive models described herein may thus prove valuable in prognostic stratification and lead to personalized treatment paradigms for this patient population.

BACKGROUND:We sought to describe targeted DNA sequencing data of persistent/recurrent laryngeal squamous cell carcinoma (LSCC) and to compare gene-specific alteration frequencies with that of primary, untreated LSCC specimens from The Cancer Genome Atlas (TCGA). METHODS:The tumors of 21 patients with persistent/recurrent LSCC were subjected to targeted DNA sequencing using the Ion AmpliSeq Comprehensive Cancer Panel. Gene-specific alteration frequencies were compared (Chi-Square test) to primary, untreated LSCC sequencing data from TCGA using the cBioPortal platform. RESULTS:Persistent/recurrent LSCC was characterized by a high rate of inactivating alterations in TP53 (38.1%) and CDKN2A (33%), amplification events of CCND1 (19.1%), and ERBB2 (14.3%), and NOTCH1 (19.1%) mutations. Comparison of primary vs persistent/recurrent LSCC revealed significant differences in alteration frequencies of eight critical genes: BAP1, CDKN2A, DCUN1D1, MSH2, MTOR, PIK3CA, TET2, and TP53. CONCLUSIONS:Our results provide preliminary support for a distinct mutational profile of persistent/recurrent LSCC that requires validation in larger cohorts.

Individuals with African ancestry have extensive genomic diversity but have been underrepresented in genomic research. There is also extensive global diversity in the exposome (the totality of human environmental exposures from conception onwards) which should be considered for integrative genomic and environmental health research in Africa. To address current research gaps, we organized a workshop on environmental health research in Africa in conjunction with the H3Africa Consortium and the African Society of Human Genetics meetings in Kigali, Rwanda. The workshop was open to all researchers with an interest in environmental health in Africa and involved presentations from experts within and outside of the Consortium. This workshop highlighted innovative research occurring on the African continent related to environmental health and the interplay between the environment and the human genome. Stories of success, challenges, and collaborative opportunities were discussed through presentations, breakout sessions, poster presentations, and a panel discussion. The workshop informed participants about environmental risk factors that can be incorporated into current or future epidemiology studies and addressed research design considerations, biospecimen collection and storage, biomarkers for measuring chemical exposures, laboratory strategies, and statistical methodologies. Inclusion of environmental exposure measurements with genomic data, including but not limited to H3Africa projects, can offer a strong platform for building gene-environment (G x E) research in Africa. Opportunities to leverage existing resources and add environmental exposure data for ongoing and planned studies were discussed. Future directions include expanding the measurement of both genomic and exposomic risk factors and incorporating sophisticated statistical approaches for analyzing high dimensional G x E data. A better understanding of how environmental and genomic factors interact with nutrition and infection is also needed. Considering that the environment represents many modifiable risk factors, these research findings can inform intervention and prevention efforts towards improving global health.

ObjectiveThe pathogenetic underpinnings of extrapulmonary small cell carcinomas (EPSCCs) of the head and neck are poorly understood. We sought to describe the clinical case and whole-exome DNA sequencing data of a patient with sinonasal Schneiderian carcinoma in situ whose tumor progressed to small cell carcinoma (SCC).Study designCase report and whole-exome sequencing of tumor DNA.SettingAcademic medical center.Subjects and methodsA 52-year-old man with sinonasal Schneiderian carcinoma in situ whose tumor progressed to small cell carcinoma. We performed whole-exome genetic sequencing and copy-number variation (CNV) analysis of tumor and normal DNA extracted from flash-frozen, paraffin-embedded (FFPE) samples.ResultsA total of 93 high-confidence, nonsynonymous somatic mutation events were identified in sinonasal SCC, including loss-of-function mutations in TP53, MAML3, a transcriptional coactivator of the Notch pathway, and GAS6, an activating ligand of the TAM family of tyrosine kinase receptors. Focal amplifications of chromosomal regions 6p25-11.1, containing SOX4 and VEGFA, and 14q32.1-32.3, containing AKT1 and the Notch inhibitory ligand DLK1, were also seen. Further CNV analysis revealed deletions in the critical cell cycle regulators CDKN2A, RB1, RBL1, and RBL2 and the chromatin modifier EP300.ConclusionsSmall cell carcinoma may rarely arise from sinonasal Schneiderian carcinoma in situ and exhibits similar genomic aberrations (eg, SOX amplification, Notch pathway inactivation) to pulmonary small cell carcinoma.

Protein engineering by directed evolution can alter proteins' structures, properties, and functions. However, membrane proteins, despite their importance to living organisms, remain relatively unexplored as targets for protein engineering and directed evolution. This gap in capabilities likely results from the tendency of membrane proteins to aggregate and fail to overexpress in bacteria cells. For example, the membrane protein caveolin-1 has been implicated in many cell signaling pathways and diseases, yet the full-length protein is too aggregation-prone for detailed mutagenesis, directed evolution, and biophysical characterization. Using a phage-displayed library of full-length caveolin-1 variants, directed evolution with alternating subtractive and functional selections isolated a full-length, soluble variant, termed cavsol, for expression in E. coli. Cavsol folds correctly and binds to its known protein ligands HIV gp41, the catalytic domain of cAMP-dependent protein kinase A, and the polymerase I and transcript release factor. As expected, cavsol does not bind off-target proteins. Cellular studies show that cavsol retains the parent protein's ability to localize at the cellular membrane. Unlike truncated versions of caveolin, cavsol forms large, oligomeric complexes consisting of approximately >50 monomeric units without requiring additional cellular components. Cavsol's secondary structure is a mixture of α-helices and β-strands. Isothermal titration calorimetry experiments reveal that cavsol binds to gp41 and PKA with low micromolar binding affinity (KD). In addition to the insights into caveolin structure and function, the approach applied here could be generalized to other membrane proteins.

This corrects the article DOI: 10.1038/sdata.2017.179.

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1-5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.

ObjectivePlasma levels of the fibrinogen degradation product D-dimer are higher among African Americans (AAs) compared with those of European ancestry and higher among women compared with men. Among AAs, little is known of the genetic architecture of D-dimer or the relationship of D-dimer to incident cardiovascular disease.Approach and resultsWe measured baseline D-dimer in 4163 AAs aged 21 to 93 years from the prospective JHS (Jackson Heart Study) cohort and assessed association with incident cardiovascular disease events. In participants with whole genome sequencing data (n=2980), we evaluated common and rare genetic variants for association with D-dimer. Each standard deviation higher baseline D-dimer was associated with a 20% to 30% increased hazard for incident coronary heart disease, stroke, and all-cause mortality. Genetic variation near F3 was associated with higher D-dimer (rs2022030, β=0.284, P=3.24×10-11). The rs2022030 effect size was nearly 3× larger among women (β=0.373, P=9.06×10-13) than among men (β=0.135, P=0.06; P interaction =0.009). The sex by rs2022030 interaction was replicated in an independent sample of 10 808 multiethnic men and women (P interaction =0.001). Finally, the African ancestral sickle cell variant (HBB rs334) was significantly associated with higher D-dimer in JHS (β=0.507, P=1.41×10-14), and this association was successfully replicated in 1933 AAs (P=2.3×10-5).ConclusionsThese results highlight D-dimer as an important predictor of cardiovascular disease risk in AAs and suggest that sex-specific and African ancestral genetic effects of the F3 and HBB loci contribute to the higher levels of D-dimer among women and AAs.

To identify novel coding association signals and facilitate characterization of mechanisms influencing glycemic traits and type 2 diabetes risk, we analyzed 109,215 variants derived from exome array genotyping together with an additional 390,225 variants from exome sequence in up to 39,339 normoglycemic individuals from five ancestry groups. We identified a novel association between the coding variant (p.Pro50Thr) in AKT2 and fasting plasma insulin (FI), a gene in which rare fully penetrant mutations are causal for monogenic glycemic disorders. The low-frequency allele is associated with a 12% increase in FI levels. This variant is present at 1.1% frequency in Finns but virtually absent in individuals from other ancestries. Carriers of the FI-increasing allele had increased 2-h insulin values, decreased insulin sensitivity, and increased risk of type 2 diabetes (odds ratio 1.05). In cellular studies, the AKT2-Thr50 protein exhibited a partial loss of function. We extend the allelic spectrum for coding variants in AKT2 associated with disorders of glucose homeostasis and demonstrate bidirectional effects of variants within the pleckstrin homology domain of AKT2.

Attention-Deficit Hyperactivity Disorder (ADHD) has high heritability; however, studies of common variation account for

Essentially all biochemistry and most molecular biology experiments require recombinant proteins. However, large, hydrophobic proteins typically aggregate into insoluble and misfolded species, and are directed into inclusion bodies. Current techniques to fold proteins recovered from inclusion bodies rely on denaturation followed by dialysis or rapid dilution. Such approaches can be time consuming, wasteful, and inefficient. Here, we describe rapid protein folding using a vortex fluidic device (VFD). This process uses mechanical energy introduced into thin films to rapidly and efficiently fold proteins. With the VFD in continuous flow mode, large volumes of protein solution can be processed per day with 100-fold reductions in both folding times and buffer volumes.

Background & aimsNeonatal sclerosing cholangitis (NSC) is a severe neonatal-onset cholangiopathy commonly leading to liver transplantation (LT) for end-stage liver disease in childhood. Liver biopsy findings histopathologically resemble those in biliary atresia (BA); however, in NSC extrahepatic bile ducts are patent, whilst in BA their lumina are obliterated. NSC is commonly seen in consanguineous kindreds, suggesting autosomal recessive inheritance.MethodsFrom 29 NSC patients (24 families) identified, DNA was available in 24 (21 families). Thirteen (7 male) patients (12 families) of consanguineous parentage were selected for whole exome sequencing. Sequence variants were filtered for homozygosity, pathogenicity, minor allele frequency, quality score, and encoded protein expression pattern.ResultsFour of 13 patients were homozygous and two were compound heterozygous for mutations in the doublecortin domain containing 2 gene (DCDC2), which encodes DCDC2 protein and is expressed in cholangiocyte cilia. Another 11 patients were sequenced: one (with one sibling pair) was compound heterozygous for DCDC2 mutations. All mutations were protein-truncating. In available liver tissue from patients with DCDC2 mutations, immunostaining for human DCDC2 and the ciliary protein acetylated alpha-tubulin (ACALT) showed no expression (n=6) and transmission electron microscopy found that cholangiocytes lacked primary cilia (n=5). DCDC2 and ACALT were expressed in NSC patients without DCDC2 mutations (n=22). Of the patients carrying DCDC2 mutations, one died awaiting LT; five came to LT, of whom one died 2years later. The other 4 are well.ConclusionAmong 24 NSC patients with available DNA, 7 had mutations in DCDC2 (6 of 19 families). NSC patients in substantial proportion harbour mutations in DCDC2. Their disease represents a novel liver-based ciliopathy.Lay summaryNeonatal sclerosing cholangitis (NSC) is a rare genetic form of liver disease presenting in infancy. Through next generation sequencing we identified mutations in the gene encoding for doublecortin domain containing 2 (DCDC2) protein in a group of NSC children. DCDC2 is a signalling and structural protein found in primary cilia of cholangiocytes. Cholangiocytes are the cells forming the biliary system which is the draining system of the liver.

BackgroundIn addition to lowering low density lipoprotein cholesterol (LDL-C), statin therapy also raises high density lipoprotein cholesterol (HDL-C) levels. Inter-individual variation in HDL-C response to statins may be partially explained by genetic variation.Methods and resultsWe performed a meta-analysis of genome-wide association studies (GWAS) to identify variants with an effect on statin-induced high density lipoprotein cholesterol (HDL-C) changes. The 123 most promising signals with p

Neutrophils are a key component of innate immunity. Individuals with low neutrophil count are susceptible to frequent infections. Linkage and association between congenital neutropenia and a single rare missense variant in TCIRG1 have been reported in a single family. Here, we report on nine rare missense variants at evolutionarily conserved sites in TCIRG1 that are associated with lower absolute neutrophil count (ANC; p = 0.005) in 1,058 participants from three cohorts: Atherosclerosis Risk in Communities (ARIC), Coronary Artery Risk Development in Young Adults (CARDIA), and Jackson Heart Study (JHS) of the NHLBI Grand Opportunity Exome Sequencing Project (GO ESP). These results validate the effects of TCIRG1 coding variation on ANC and suggest that this gene may be associated with a spectrum of mild to severe effects on ANC.

Waist-to-hip ratio (WHR), a relative comparison of waist and hip circumferences, is an easily accessible measurement of body fat distribution, in particular central abdominal fat. A high WHR indicates more intra-abdominal fat deposition and is an established risk factor for cardiovascular disease and type 2 diabetes. Recent genome-wide association studies have identified numerous common genetic loci influencing WHR, but the contributions of rare variants have not been previously reported. We investigated rare variant associations with WHR in 1510 European-American and 1186 African-American women from the National Heart, Lung, and Blood Institute-Exome Sequencing Project. Association analysis was performed on the gene level using several rare variant association methods. The strongest association was observed for rare variants in IKBKB (P=4.0 × 10(-8)) in European-Americans, where rare variants in this gene are predicted to decrease WHRs. The activation of the IKBKB gene is involved in inflammatory processes and insulin resistance, which may affect normal food intake and body weight and shape. Meanwhile, aggregation of rare variants in COBLL1, previously found to harbor common variants associated with WHR and fasting insulin, were nominally associated (P=2.23 × 10(-4)) with higher WHR in European-Americans. However, these significant results are not shared between African-Americans and European-Americans that may be due to differences in the allelic architecture of the two populations and the small sample sizes. Our study indicates that the combined effect of rare variants contribute to the inter-individual variation in fat distribution through the regulation of insulin response.

Caveolin-1 is a target for academic and pharmaceutical research due to its many cellular roles and associated diseases. We report peptide WL47 (1), a small, high-affinity, selective disrupter of caveolin-1 oligomers. Developed and optimized through screening and analysis of synthetic peptide libraries, ligand 1 has 7500-fold improved affinity compared to its T20 parent ligand and an 80% decrease in sequence length. Ligand 1 will permit targeted study of caveolin-1 function.

We performed whole-genome sequencing (WGS) of 208 genomes from 53 families affected by simplex autism. For the majority of these families, no copy-number variant (CNV) or candidate de novo gene-disruptive single-nucleotide variant (SNV) had been detected by microarray or whole-exome sequencing (WES). We integrated multiple CNV and SNV analyses and extensive experimental validation to identify additional candidate mutations in eight families. We report that compared to control individuals, probands showed a significant (p = 0.03) enrichment of de novo and private disruptive mutations within fetal CNS DNase I hypersensitive sites (i.e., putative regulatory regions). This effect was only observed within 50 kb of genes that have been previously associated with autism risk, including genes where dosage sensitivity has already been established by recurrent disruptive de novo protein-coding mutations (ARID1B, SCN2A, NR3C2, PRKCA, and DSCAM). In addition, we provide evidence of gene-disruptive CNVs (in DISC1, WNT7A, RBFOX1, and MBD5), as well as smaller de novo CNVs and exon-specific SNVs missed by exome sequencing in neurodevelopmental genes (e.g., CANX, SAE1, and PIK3CA). Our results suggest that the detection of smaller, often multiple CNVs affecting putative regulatory elements might help explain additional risk of simplex autism.

Multiple pterygium syndrome (MPS) is a phenotypically and genetically heterogeneous group of rare Mendelian conditions characterized by multiple pterygia, scoliosis, and congenital contractures of the limbs. MPS typically segregates as an autosomal-recessive disorder, but rare instances of autosomal-dominant transmission have been reported. Whereas several mutations causing recessive MPS have been identified, the genetic basis of dominant MPS remains unknown. We identified four families affected by dominantly transmitted MPS characterized by pterygia, camptodactyly of the hands, vertebral fusions, and scoliosis. Exome sequencing identified predicted protein-altering mutations in embryonic myosin heavy chain (MYH3) in three families. MYH3 mutations underlie distal arthrogryposis types 1, 2A, and 2B, but all mutations reported to date occur in the head and neck domains. In contrast, two of the mutations found to cause MPS in this study occurred in the tail domain. The phenotypic overlap among persons with MPS, coupled with physical findings distinct from other conditions caused by mutations in MYH3, suggests that the developmental mechanism underlying MPS differs from that of other conditions and/or that certain functions of embryonic myosin might be perturbed by disruption of specific residues and/or domains. Moreover, the vertebral fusions in persons with MPS, coupled with evidence of MYH3 expression in bone, suggest that embryonic myosin plays a role in skeletal development.

Freeman-Sheldon syndrome, or distal arthrogryposis type 2A (DA2A), is an autosomal-dominant condition caused by mutations in MYH3 and characterized by multiple congenital contractures of the face and limbs and normal cognitive development. We identified a subset of five individuals who had been putatively diagnosed with "DA2A with severe neurological abnormalities" and for whom congenital contractures of the limbs and face, hypotonia, and global developmental delay had resulted in early death in three cases; this is a unique condition that we now refer to as CLIFAHDD syndrome. Exome sequencing identified missense mutations in the sodium leak channel, non-selective (NALCN) in four families affected by CLIFAHDD syndrome. We used molecular-inversion probes to screen for NALCN in a cohort of 202 distal arthrogryposis (DA)-affected individuals as well as concurrent exome sequencing of six other DA-affected individuals, thus revealing NALCN mutations in ten additional families with "atypical" forms of DA. All 14 mutations were missense variants predicted to alter amino acid residues in or near the S5 and S6 pore-forming segments of NALCN, highlighting the functional importance of these segments. In vitro functional studies demonstrated that NALCN alterations nearly abolished the expression of wild-type NALCN, suggesting that alterations that cause CLIFAHDD syndrome have a dominant-negative effect. In contrast, homozygosity for mutations in other regions of NALCN has been reported in three families affected by an autosomal-recessive condition characterized mainly by hypotonia and severe intellectual disability. Accordingly, mutations in NALCN can cause either a recessive or dominant condition characterized by varied though overlapping phenotypic features, perhaps based on the type of mutation and affected protein domain(s).

Recombinant protein overexpression of large proteins in bacteria often results in insoluble and misfolded proteins directed to inclusion bodies. We report the application of shear stress in micrometer-wide, thin fluid films to refold boiled hen egg white lysozyme, recombinant hen egg white lysozyme, and recombinant caveolin-1. Furthermore, the approach allowed refolding of a much larger protein, cAMP-dependent protein kinase A (PKA). The reported methods require only minutes, which is more than 100 times faster than conventional overnight dialysis. This rapid refolding technique could significantly shorten times, lower costs, and reduce waste streams associated with protein expression for a wide range of industrial and research applications.

We conducted a genome-wide association analysis of 7 subfractions of low density lipoproteins (LDLs) and 3 subfractions of intermediate density lipoproteins (IDLs) measured by gradient gel electrophoresis, and their response to statin treatment, in 1868 individuals of European ancestry from the Pharmacogenomics and Risk of Cardiovascular Disease study. Our analyses identified four previously-implicated loci (SORT1, APOE, LPA, and CETP) as containing variants that are very strongly associated with lipoprotein subfractions (log(10)Bayes Factor > 15). Subsequent conditional analyses suggest that three of these (APOE, LPA and CETP) likely harbor multiple independently associated SNPs. Further, while different variants typically showed different characteristic patterns of association with combinations of subfractions, the two SNPs in CETP show strikingly similar patterns--both in our original data and in a replication cohort--consistent with a common underlying molecular mechanism. Notably, the CETP variants are very strongly associated with LDL subfractions, despite showing no association with total LDLs in our study, illustrating the potential value of the more detailed phenotypic measurements. In contrast with these strong subfraction associations, genetic association analysis of subfraction response to statins showed much weaker signals (none exceeding log(10)Bayes Factor of 6). However, two SNPs (in APOE and LPA) previously-reported to be associated with LDL statin response do show some modest evidence for association in our data, and the subfraction response proles at the LPA SNP are consistent with the LPA association, with response likely being due primarily to resistance of Lp(a) particles to statin therapy. An additional important feature of our analysis is that, unlike most previous analyses of multiple related phenotypes, we analyzed the subfractions jointly, rather than one at a time. Comparisons of our multivariate analyses with standard univariate analyses demonstrate that multivariate analyses can substantially increase power to detect associations. Software implementing our multivariate analysis methods is available at http://stephenslab.uchicago.edu/software.html.

ImportanceThe association between sickle cell trait (SCT) and chronic kidney disease (CKD) is uncertain.ObjectiveTo describe the relationship between SCT and CKD and albuminuria in self-identified African Americans.Design, setting, and participantsUsing 5 large, prospective, US population-based studies (the Atherosclerosis Risk in Communities Study [ARIC, 1987-2013; n = 3402], Jackson Heart Study [JHS, 2000-2012; n = 2105], Coronary Artery Risk Development in Young Adults [CARDIA, 1985-2006; n = 848], Multi-Ethnic Study of Atherosclerosis [MESA, 2000-2012; n = 1620], and Women's Health Initiative [WHI, 1993-2012; n = 8000]), we evaluated 15,975 self-identified African Americans (1248 participants with SCT [SCT carriers] and 14,727 participants without SCT [noncarriers]).Main outcomes and measuresPrimary outcomes were CKD (defined as an estimated glomerular filtration rate [eGFR] of 30 mg/g or albumin excretion rate >30 mg/24 hours), and decline in eGFR (defined as a decrease of >3 mL/min/1.73 m2 per year). Effect sizes were calculated separately for each cohort and were subsequently meta-analyzed using a random-effects model.ResultsA total of 2233 individuals (239 of 1247 SCT carriers [19.2%] vs 1994 of 14,722 noncarriers [13.5%]) had CKD, 1298 (140 of 675 SCT carriers [20.7%] vs 1158 of 8481 noncarriers [13.7%]) experienced incident CKD, 1719 (150 of 665 SCT carriers [22.6%] vs 1569 of 8249 noncarriers [19.0%]) experienced decline in eGFR, and 1322 (154 of 485 SCT carriers [31.8%] vs 1168 of 5947 noncarriers [19.6%]) had albuminuria during the study period. Individuals with SCT had an increased risk of CKD (odds ratio [OR], 1.57 [95% CI, 1.34-1.84]; absolute risk difference [ARD], 7.6% [95% CI, 4.7%-10.8%]), incident CKD (OR, 1.79 [95% CI, 1.45-2.20]; ARD, 8.5% [95% CI, 5.1%-12.3%]), and decline in eGFR (OR, 1.32 [95% CI, 1.07-1.61]; ARD, 6.1% [95% CI, 1.4%-13.0%]) compared with noncarriers. Sickle cell trait was also associated with albuminuria (OR, 1.86 [95% CI, 1.49-2.31]; ARD, 12.6% [95% CI, 7.7%-17.7%]).Conclusions and relevanceAmong African Americans in these cohorts, the presence of SCT was associated with an increased risk of CKD, decline in eGFR, and albuminuria, compared with noncarriers. These findings suggest that SCT may be associated with the higher risk of kidney disease in African Americans.

Whole exome sequencing has proven to be a powerful tool for understanding the genetic architecture of human disease. Here we apply it to more than 2,500 simplex families, each having a child with an autistic spectrum disorder. By comparing affected to unaffected siblings, we show that 13% of de novo missense mutations and 43% of de novo likely gene-disrupting (LGD) mutations contribute to 12% and 9% of diagnoses, respectively. Including copy number variants, coding de novo mutations contribute to about 30% of all simplex and 45% of female diagnoses. Almost all LGD mutations occur opposite wild-type alleles. LGD targets in affected females significantly overlap the targets in males of lower intelligence quotient (IQ), but neither overlaps significantly with targets in males of higher IQ. We estimate that LGD mutation in about 400 genes can contribute to the joint class of affected females and males of lower IQ, with an overlapping and similar number of genes vulnerable to contributory missense mutation. LGD targets in the joint class overlap with published targets for intellectual disability and schizophrenia, and are enriched for chromatin modifiers, FMRP-associated genes and embryonically expressed genes. Most of the significance for the latter comes from affected females.

Statins effectively lower LDL cholesterol levels in large studies and the observed interindividual response variability may be partially explained by genetic variation. Here we perform a pharmacogenetic meta-analysis of genome-wide association studies (GWAS) in studies addressing the LDL cholesterol response to statins, including up to 18,596 statin-treated subjects. We validate the most promising signals in a further 22,318 statin recipients and identify two loci, SORT1/CELSR2/PSRC1 and SLCO1B1, not previously identified in GWAS. Moreover, we confirm the previously described associations with APOE and LPA. Our findings advance the understanding of the pharmacogenetic architecture of statin response.

Gordon syndrome (GS), or distal arthrogryposis type 3, is a rare, autosomal-dominant disorder characterized by cleft palate and congenital contractures of the hands and feet. Exome sequencing of five GS-affected families identified mutations in piezo-type mechanosensitive ion channel component 2 (PIEZO2) in each family. Sanger sequencing revealed PIEZO2 mutations in five of seven additional families studied (for a total of 10/12 [83%] individuals), and nine families had an identical c.8057G>A (p.Arg2686His) mutation. The phenotype of GS overlaps with distal arthrogryposis type 5 (DA5) and Marden-Walker syndrome (MWS). Using molecular inversion probes for targeted sequencing to screen PIEZO2, we found mutations in 24/29 (82%) DA5-affected families and one of two MWS-affected families. The presence of cleft palate was significantly associated with c.8057G>A (Fisher's exact test, adjusted p value < 0.0001). Collectively, although GS, DA5, and MWS have traditionally been considered separate disorders, our findings indicate that they are etiologically related and perhaps represent variable expressivity of the same condition.

Elucidating genetic causes of cholestasis has proved to be important in understanding the physiology and pathophysiology of the liver. Here we show that protein-truncating mutations in the tight junction protein 2 gene (TJP2) cause failure of protein localization and disruption of tight-junction structure, leading to severe cholestatic liver disease. These findings contrast with those in the embryonic-lethal knockout mouse, highlighting differences in redundancy in junctional complexes between organs and species.

Elevated low-density lipoprotein cholesterol (LDL-C) is a treatable, heritable risk factor for cardiovascular disease. Genome-wide association studies (GWASs) have identified 157 variants associated with lipid levels but are not well suited to assess the impact of rare and low-frequency variants. To determine whether rare or low-frequency coding variants are associated with LDL-C, we exome sequenced 2,005 individuals, including 554 individuals selected for extreme LDL-C (>98(th) or

Statins are prescribed widely to lower plasma low-density lipoprotein (LDL) concentrations and cardiovascular disease risk and have been shown to have beneficial effects in a broad range of patients. However, statins are associated with an increased risk, albeit small, of clinical myopathy and type 2 diabetes. Despite evidence for substantial genetic influence on LDL concentrations, pharmacogenomic trials have failed to identify genetic variations with large effects on either statin efficacy or toxicity, and have produced little information regarding mechanisms that modulate statin response. Here we identify a downstream target of statin treatment by screening for the effects of in vitro statin exposure on genetic associations with gene expression levels in lymphoblastoid cell lines derived from 480 participants of a clinical trial of simvastatin treatment. This analysis identified six expression quantitative trait loci (eQTLs) that interacted with simvastatin exposure, including rs9806699, a cis-eQTL for the gene glycine amidinotransferase (GATM) that encodes the rate-limiting enzyme in creatine synthesis. We found this locus to be associated with incidence of statin-induced myotoxicity in two separate populations (meta-analysis odds ratio = 0.60). Furthermore, we found that GATM knockdown in hepatocyte-derived cell lines attenuated transcriptional response to sterol depletion, demonstrating that GATM may act as a functional link between statin-mediated lowering of cholesterol and susceptibility to statin-induced myopathy.

Opsismodysplasia is a rare, autosomal-recessive skeletal dysplasia characterized by short stature, characteristic facial features, and in some cases severe renal phosphate wasting. We used linkage analysis and whole-genome sequencing of a consanguineous trio to discover that mutations in inositol polyphosphate phosphatase-like 1 (INPPL1) cause opsismodysplasia with or without renal phosphate wasting. Evaluation of 12 families with opsismodysplasia revealed that INPPL1 mutations explain ~60% of cases overall, including both of the families in our cohort with more than one affected child and 50% of the simplex cases.

Plasma concentrations of total cholesterol, low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides are among the most important risk factors for coronary artery disease (CAD) and are targets for therapeutic intervention. We screened the genome for common variants associated with plasma lipids in >100,000 individuals of European ancestry. Here we report 95 significantly associated loci (P < 5 x 10(-8)), with 59 showing genome-wide significant association with lipid traits for the first time. The newly reported associations include single nucleotide polymorphisms (SNPs) near known lipid regulators (for example, CYP7A1, NPC1L1 and SCARB1) as well as in scores of loci not previously implicated in lipoprotein metabolism. The 95 loci contribute not only to normal variation in lipid traits but also to extreme lipid phenotypes and have an impact on lipid traits in three non-European populations (East Asians, South Asians and African Americans). Our results identify several novel loci associated with plasma lipids that are also associated with CAD. Finally, we validated three of the novel genes-GALNT2, PPP1R3B and TTC39B-with experiments in mouse models. Taken together, our findings provide the foundation to develop a broader biological understanding of lipoprotein metabolism and to identify new therapeutic opportunities for the prevention of CAD.

BackgroundStatins effectively lower total and plasma LDL-cholesterol, but the magnitude of decrease varies among individuals. To identify single nucleotide polymorphisms (SNPs) contributing to this variation, we performed a combined analysis of genome-wide association (GWA) results from three trials of statin efficacy.Methods and principal findingsBayesian and standard frequentist association analyses were performed on untreated and statin-mediated changes in LDL-cholesterol, total cholesterol, HDL-cholesterol, and triglyceride on a total of 3932 subjects using data from three studies: Cholesterol and Pharmacogenetics (40 mg/day simvastatin, 6 weeks), Pravastatin/Inflammation CRP Evaluation (40 mg/day pravastatin, 24 weeks), and Treating to New Targets (10 mg/day atorvastatin, 8 weeks). Genotype imputation was used to maximize genomic coverage and to combine information across studies. Phenotypes were normalized within each study to account for systematic differences among studies, and fixed-effects combined analysis of the combined sample were performed to detect consistent effects across studies. Two SNP associations were assessed as having posterior probability greater than 50%, indicating that they were more likely than not to be genuinely associated with statin-mediated lipid response. SNP rs8014194, located within the CLMN gene on chromosome 14, was strongly associated with statin-mediated change in total cholesterol with an 84% probability by Bayesian analysis, and a p-value exceeding conventional levels of genome-wide significance by frequentist analysis (P = 1.8 x 10(-8)). This SNP was less significantly associated with change in LDL-cholesterol (posterior probability = 0.16, P = 4.0 x 10(-6)). Bayesian analysis also assigned a 51% probability that rs4420638, located in APOC1 and near APOE, was associated with change in LDL-cholesterol.Conclusions and significanceUsing combined GWA analysis from three clinical trials involving nearly 4,000 individuals treated with simvastatin, pravastatin, or atorvastatin, we have identified SNPs that may be associated with variation in the magnitude of statin-mediated reduction in total and LDL-cholesterol, including one in the CLMN gene for which statistical evidence for association exceeds conventional levels of genome-wide significance.Trial registrationPRINCE and TNT are not registered. CAP is registered at Clinicaltrials.gov NCT00451828.

BackgroundTetrahymena thermophila, a widely studied model for cellular and molecular biology, is a binucleated single-celled organism with a germline micronucleus (MIC) and somatic macronucleus (MAC). The recent draft MAC genome assembly revealed low sequence repetitiveness, a result of the epigenetic removal of invasive DNA elements found only in the MIC genome. Such low repetitiveness makes complete closure of the MAC genome a feasible goal, which to achieve would require standard closure methods as well as removal of minor MIC contamination of the MAC genome assembly. Highly accurate preliminary annotation of Tetrahymena's coding potential was hindered by the lack of both comparative genomic sequence information from close relatives and significant amounts of cDNA evidence, thus limiting the value of the genomic information and also leaving unanswered certain questions, such as the frequency of alternative splicing.ResultsWe addressed the problem of MIC contamination using comparative genomic hybridization with purified MIC and MAC DNA probes against a whole genome oligonucleotide microarray, allowing the identification of 763 genome scaffolds likely to contain MIC-limited DNA sequences. We also employed standard genome closure methods to essentially finish over 60% of the MAC genome. For the improvement of annotation, we have sequenced and analyzed over 60,000 verified EST reads from a variety of cellular growth and development conditions. Using this EST evidence, a combination of automated and manual reannotation efforts led to updates that affect 16% of the current protein-coding gene models. By comparing EST abundance, many genes showing apparent differential expression between these conditions were identified. Rare instances of alternative splicing and uses of the non-standard amino acid selenocysteine were also identified.ConclusionWe report here significant progress in genome closure and reannotation of Tetrahymena thermophila. Our experience to date suggests that complete closure of the MAC genome is attainable. Using the new EST evidence, automated and manual curation has resulted in substantial improvements to the over 24,000 gene models, which will be valuable to researchers studying this model organism as well as for comparative genomics purposes.