An Organismal CNV Mutator Phenotype Restricted to Early Human Development.

De novo copy number variants (dnCNVs) arising at multiple loci in a personal genome have usually been considered to reflect cancer somatic genomic instabilities. We describe a multiple dnCNV (MdnCNV) phenomenon in which individuals with genomic disorders carry five to ten constitutional dnCNVs. These CNVs originate from independent formation incidences, are predominantly tandem duplications or complex gains, exhibit breakpoint junction features reminiscent of replicative repair, and show increased de novo point mutations flanking the rearrangement junctions. The active CNV mutation shower appears to be restricted to a transient perizygotic period. We propose that a defect in the CNV formation process is responsible for the "CNV-mutator state," and this state is dampened after early embryogenesis. The constitutional MdnCNV phenomenon resembles chromosomal instability in various cancers. Investigations of this phenomenon may provide unique access to understanding genomic disorders, structural variant mutagenesis, human evolution, and cancer biology.

Regulation of T Cell Receptor Signaling by DENND1B in TH2 Cells and Allergic Disease.

The DENN domain is an evolutionary conserved protein module found in all eukaryotes and serves as an exchange factor for Rab-GTPases to regulate diverse cellular functions. Variants in DENND1B are associated with development of childhood asthma and other immune disorders. To understand how DENND1B may contribute to human disease, Dennd1b(-/-) mice were generated and exhibit hyper-allergic responses following antigen challenge. Dennd1b(-/-) TH2, but not other TH cells, exhibit delayed receptor-induced T cell receptor (TCR) downmodulation, enhanced TCR signaling, and increased production of effector cytokines. As DENND1B interacts with AP-2 and Rab35, TH2 cells deficient in AP-2 or Rab35 also exhibit enhanced TCR-mediated effector functions. Moreover, human TH2 cells carrying asthma-associated DENND1B variants express less DENND1B and phenocopy Dennd1b(-/-) TH2 cells. These results provide a molecular basis for how DENND1B, a previously unrecognized regulator of TCR downmodulation in TH2 cells, contributes to asthma pathogenesis and how DENN-domain-containing proteins may contribute to other human disorders.

Haploinsufficiency underlies the neurodevelopmental consequences of SLC6A1 variants.

Heterozygous variants in SLC6A1, encoding the GAT-1 GABA transporter, are associated with seizures, developmental delay, and autism. The majority of affected individuals carry missense variants, many of which are recurrent germline de novo mutations, raising the possibility of gain-of-function or dominant-negative effects. To understand the functional consequences, we performed an in vitro GABA uptake assay for 213 unique variants, including 24 control variants. De novo variants consistently resulted in a decrease in GABA uptake, in keeping with haploinsufficiency underlying all neurodevelopmental phenotypes. Where present, ClinVar pathogenicity reports correlated well with GABA uptake data; the functional data can inform future reports for the remaining 72% of unscored variants. Surface localization was assessed for 86 variants; two-thirds of loss-of-function missense variants prevented GAT-1 from being present on the membrane while GAT-1 was on the surface but with reduced activity for the remaining third. Surprisingly, recurrent de novo missense variants showed moderate loss-of-function effects that reduced GABA uptake with no evidence for dominant-negative or gain-of-function effects. Using linear regression across multiple missense severity scores to extrapolate the functional data to all potential SLC6A1 missense variants, we observe an abundance of GAT-1 residues that are sensitive to substitution. The extent of this missense vulnerability accounts for the clinically observed missense enrichment; overlap with hypermutable CpG sites accounts for the recurrent missense variants. Strategies to increase the expression of the wild-type SLC6A1 allele are likely to be beneficial across neurodevelopmental disorders, though the developmental stage and extent of required rescue remain unknown.

SLC6A1 variant pathogenicity, molecular function and phenotype: a genetic and clinical analysis.

Genetic variants in the SLC6A1 gene can cause a broad phenotypic disease spectrum by altering the protein function. Thus, systematically curated clinically relevant genotype-phenotype associations are needed to understand the disease mechanism and improve therapeutic decision-making. We aggregated genetic and clinical data from 172 individuals with likely pathogenic/pathogenic (lp/p) SLC6A1 variants and functional data for 184 variants (14.1% lp/p). Clinical and functional data were available for a subset of 126 individuals. We explored the potential associations of variant positions on the GAT1 3D structure with variant pathogenicity, altered molecular function and phenotype severity using bioinformatic approaches. The GAT1 transmembrane domains 1, 6 and extracellular loop 4 (EL4) were enriched for patient over population variants. Across functionally tested missense variants (n = 156), the spatial proximity from the ligand was associated with loss-of-function in the GAT1 transporter activity. For variants with complete loss of in vitro GABA uptake, we found a 4.6-fold enrichment in patients having severe disease versus non-severe disease (P = 2.9 × 10-3, 95% confidence interval: 1.5-15.3). In summary, we delineated associations between the 3D structure and variant pathogenicity, variant function and phenotype in SLC6A1-related disorders. This knowledge supports biology-informed variant interpretation and research on GAT1 function. All our data can be interactively explored in the SLC6A1 portal (https://slc6a1-portal.broadinstitute.org/).

LACC1 Regulates TNF and IL-17 in Mouse Models of Arthritis and Inflammation.

Both common and rare genetic variants of laccase domain-containing 1 (LACC1, previously C13orf31) are associated with inflammatory bowel disease, leprosy, Behcet disease, and systemic juvenile idiopathic arthritis. However, the functional relevance of these variants is unclear. In this study, we use LACC1-deficient mice to gain insight into the role of LACC1 in regulating inflammation. Following oral administration of Citrobacter rodentium, LACC1 knockout (KO) mice had more severe colon lesions compared with wildtype (WT) controls. Immunization with collagen II, a collagen-induced arthritis (CIA) model, resulted in an accelerated onset of arthritis and significantly worse arthritis and inflammation in LACC1 KO mice. Similar results were obtained in a mannan-induced arthritis model. Serum and local TNF in CIA paws and C. rodentium colons were significantly increased in LACC1 KO mice compared with WT controls. The percentage of IL-17A-producing CD4+ T cells was elevated in LACC1 KO mice undergoing CIA as well as aged mice compared with WT controls. Neutralization of IL-17, but not TNF, prevented enhanced mannan-induced arthritis in LACC1 KO mice. These data provide new mechanistic insight into the function of LACC1 in regulating TNF and IL-17 during inflammatory responses. We hypothesize that these effects contribute to immune-driven pathologies observed in individuals carrying LACC1 variants.

Germline genetic polymorphisms influence tumor gene expression and immune cell infiltration.

Cancer immunotherapy has emerged as an effective therapy in a variety of cancers. However, a key challenge in the field is that only a subset of patients who receive immunotherapy exhibit durable response. It has been hypothesized that host genetics influences the inherent immune profiles of patients and may underlie their differential response to immunotherapy. Herein, we systematically determined the association of common germline genetic variants with gene expression and immune cell infiltration of the tumor. We identified 64,094 expression quantitative trait loci (eQTLs) that associated with 18,210 genes (eGenes) across 24 human cancers. Overall, eGenes were enriched for their being involved in immune processes, suggesting that expression of immune genes can be shaped by hereditary genetic variants. We identified the endoplasmic reticulum aminopeptidase 2 (ERAP2) gene as a pan-cancer type eGene whose expression levels stratified overall survival in a subset of patients with bladder cancer receiving anti-PD-L1 (atezolizumab) therapy. Finally, we identified 103 gene signature QTLs (gsQTLs) that were associated with predicted immune cell abundance within the tumor microenvironment. Our findings highlight the impact of germline SNPs on cancer-immune phenotypes and response to therapy; and these analyses provide a resource for integration of germline genetics as a component of personalized cancer immunotherapy.

Paired Immunoglobulin-like Type 2 Receptor Alpha G78R variant alters ligand binding and confers protection to Alzheimer's disease.

Paired Immunoglobulin-like Type 2 Receptor Alpha (PILRA) is a cell surface inhibitory receptor that recognizes specific O-glycosylated proteins and is expressed on various innate immune cell types including microglia. We show here that a common missense variant (G78R, rs1859788) of PILRA is the likely causal allele for the confirmed Alzheimer's disease risk locus at 7q21 (rs1476679). The G78R variant alters the interaction of residues essential for sialic acid engagement, resulting in >50% reduced binding for several PILRA ligands including a novel ligand, complement component 4A, and herpes simplex virus 1 (HSV-1) glycoprotein B. PILRA is an entry receptor for HSV-1 via glycoprotein B, and macrophages derived from R78 homozygous donors showed significantly decreased levels of HSV-1 infection at several multiplicities of infection compared to homozygous G78 macrophages. We propose that PILRA G78R protects individuals from Alzheimer's disease risk via reduced inhibitory signaling in microglia and reduced microglial infection during HSV-1 recurrence.

Previously reported placebo-response-associated variants do not predict patient outcomes in inflammatory disease Phase III trial placebo arms.

In clinical trials, a placebo response refers to improvement in disease symptoms arising from the psychological effect of receiving a treatment rather than the actual treatment under investigation. Previous research has reported genomic variation associated with the likelihood of observing a placebo response, but these studies have been limited in scope and have not been validated. Here, we analyzed whole-genome sequencing data from 784 patients undergoing placebo treatment in Phase III Asthma or Rheumatoid Arthritis trials to assess the impact of previously reported variation on patient outcomes in the placebo arms and to identify novel variants associated with the placebo response. Contrary to expectations based on previous reports, we did not observe any statistically significant associations between genomic variants and placebo treatment outcome. Our findings suggest that the biological origin of the placebo response is complex and likely to be variable between disease areas.

The Ox40/Ox40 Ligand Pathway Promotes Pathogenic Th Cell Responses, Plasmablast Accumulation, and Lupus Nephritis in NZB/W F1 Mice.

Ox40 ligand (Ox40L) locus genetic variants are associated with the risk for systemic lupus erythematosus (SLE); however, it is unclear how Ox40L contributes to SLE pathogenesis. In this study, we evaluated the contribution of Ox40L and its cognate receptor, Ox40, using in vivo agonist and antagonist approaches in the NZB × NZW (NZB/W) F1 mouse model of SLE. Ox40 was highly expressed on several CD4 Th cell subsets in the spleen and kidney of diseased mice, and expression correlated with disease severity. Treatment of aged NZB/W F1 mice with agonist anti-Ox40 mAbs potently exacerbated renal disease, which was accompanied by activation of kidney-infiltrating T cells and cytokine production. The agonist mAbs also induced activation and inflammatory gene expression in splenic CD4 T cells, including IFN-regulated genes, increased the number of follicular helper T cells and plasmablasts in the spleen, and led to elevated levels of serum IgM and enhanced renal glomerular IgM deposition. In a type I IFN-accelerated lupus model, treatment with an antagonist Ox40:Fc fusion protein significantly delayed the onset of severe proteinuria and improved survival. These data support the hypothesis that the Ox40/Ox40L pathway drives cellular and humoral autoimmune responses during lupus nephritis in NZB/W F1 mice and emphasize the potential clinical value of targeting this pathway in human lupus.

Mutational History of a Human Cell Lineage from Somatic to Induced Pluripotent Stem Cells.

The accuracy of replicating the genetic code is fundamental. DNA repair mechanisms protect the fidelity of the genome ensuring a low error rate between generations. This sustains the similarity of individuals whilst providing a repertoire of variants for evolution. The mutation rate in the human genome has recently been measured to be 50-70 de novo single nucleotide variants (SNVs) between generations. During development mutations accumulate in somatic cells so that an organism is a mosaic. However, variation within a tissue and between tissues has not been analysed. By reprogramming somatic cells into induced pluripotent stem cells (iPSCs), their genomes and the associated mutational history are captured. By sequencing the genomes of polyclonal and monoclonal somatic cells and derived iPSCs we have determined the mutation rates and show how the patterns change from a somatic lineage in vivo through to iPSCs. Somatic cells have a mutation rate of 14 SNVs per cell per generation while iPSCs exhibited a ten-fold lower rate. Analyses of mutational signatures suggested that deamination of methylated cytosine may be the major mutagenic source in vivo, whilst oxidative DNA damage becomes dominant in vitro. Our results provide insights for better understanding of mutational processes and lineage relationships between human somatic cells. Furthermore it provides a foundation for interpretation of elevated mutation rates and patterns in cancer.

geneAttribution: trait agnostic identification of candidate genes associated with noncoding variation.

Motivation: We have developed geneAttribution, an R package that assigns candidate causal gene(s) to a risk variant identified by a genetic association study such as a GWAS. The method combines user-supplied functional annotation such as expression quantitative trait loci (eQTL) or Hi-C genome conformation data and reports the most likely candidate genes. In the absence of annotation data, geneAttribution relies on the distances between the genes and the input variant. Availability and Implementation: The package is freely available from http://www.bioconductor.org/ . A quick-start vignette is included with the package. Contact: wustera@gene.com.

DeNovoGear: de novo indel and point mutation discovery and phasing.

We present DeNovoGear software for analyzing de novo mutations from familial and somatic tissue sequencing data. DeNovoGear uses likelihood-based error modeling to reduce the false positive rate of mutation discovery in exome analysis and fragment information to identify the parental origin of germ-line mutations. We used DeNovoGear on human whole-genome sequencing data to produce a set of predicted de novo insertion and/or deletion (indel) mutations with a 95% validation rate.

Comorbidities confound metabolomics studies of human disease.

The co-occurrence of multiple chronic conditions, termed multimorbidity, presents an expanding global health challenge, demanding effective diagnostics and treatment strategies. Chronic ailments such as obesity, diabetes, and cardiovascular diseases have been linked to metabolites interacting between the host and microbiota. In this study, we investigated the impact of co-existing conditions on risk estimations for 1375 plasma metabolites in 919 individuals from population-based Estonian Biobank cohort using liquid chromatography mass spectrometry (LC-MS) method. We leveraged annually linked national electronic health records (EHRs) data to delineate comorbidities in incident cases and controls for the 14 common chronic conditions. Among the 254 associations observed across 13 chronic conditions, we primarily identified disease-specific risk factors (92%, 217/235), with most predictors (93%, 219/235) found to be related to the gut microbiome upon cross-referencing recent literature data. Accounting for comorbidities led to a reduction of common metabolite predictors across various conditions. In conclusion, our study underscores the potential of utilizing biobank-linked retrospective and prospective EHRs for the disease-specific profiling of diverse multifactorial chronic conditions.

Rare CIDEC coding variants enriched in age-related macular degeneration patients with small low-luminance deficit cause lipid droplet and fat storage defects.

BACKGROUND: The basis of Age-related macular degeneration (AMD) genetic risk has been well documented; however, few studies have looked at genetic biomarkers of disease progression or treatment response within advanced AMD patients. Here we report the first genome-wide analysis of genetic determinants of low-luminance vision deficit (LLD), which is seen as predictive of visual acuity loss and anti-VEGF treatment response in neovascular AMD patients. METHODS: AMD patients were separated into small- and large-LLD groups for comparison and whole genome sequencing was performed. Genetic determinants of LLD were assessed by common and rare variant genetic analysis. Follow-up functional analysis of rare coding variants identified by the burden test was then performed in vitro. RESULTS: We identified four coding variants in the CIDEC gene. These rare variants were only present in patients with a small LLD, which has been previously shown to indicate better prognosis and better anti-VEGF treatment response. Our in vitro functional characterization of these CIDEC alleles revealed that all decrease the binding affinity between CIDEC and the lipid droplet fusion effectors PLIN1, RAB8A and AS160. The rare CIDEC alleles all cause a hypomorphic defect in lipid droplet fusion and enlargement, resulting in a decreased fat storage capability in adipocytes. CONCLUSIONS: As we did not detect CIDEC expression in the ocular tissue affected by AMD, our results suggest that the CIDEC variants do not play a direct role in the eye and influence low-luminance vision deficit via an indirect and systemic effect related to fat storage capacity.

Spial: analysis of subtype-specific features in multiple sequence alignments of proteins.

MOTIVATION: Spial (Specificity in alignments) is a tool for the comparative analysis of two alignments of evolutionarily related sequences that differ in their function, such as two receptor subtypes. It highlights functionally important residues that are either specific to one of the two alignments or conserved across both alignments. It permits visualization of this information in three complementary ways: by colour-coding alignment positions, by sequence logos and optionally by colour-coding the residues of a protein structure provided by the user. This can aid in the detection of residues that are involved in the subtype-specific interaction with a ligand, other proteins or nucleic acids. Spial may also be used to detect residues that may be post-translationally modified in one of the two sets of sequences. AVAILABILITY: http://www.mrc-lmb.cam.ac.uk/genomes/spial/; supplementary information is available at http://www.mrc-lmb.cam.ac.uk/genomes/spial/help.html.

Cross-species chemogenomic profiling reveals evolutionarily conserved drug mode of action.

We present a cross-species chemogenomic screening platform using libraries of haploid deletion mutants from two yeast species, Saccharomyces cerevisiae and Schizosaccharomyces pombe. We screened a set of compounds of known and unknown mode of action (MoA) and derived quantitative drug scores (or D-scores), identifying mutants that are either sensitive or resistant to particular compounds. We found that compound-functional module relationships are more conserved than individual compound-gene interactions between these two species. Furthermore, we observed that combining data from both species allows for more accurate prediction of MoA. Finally, using this platform, we identified a novel small molecule that acts as a DNA damaging agent and demonstrate that its MoA is conserved in human cells.

Identifying therapeutic drug targets using bidirectional effect genes.

Prioritizing genes for translation to therapeutics for common diseases has been challenging. Here, we propose an approach to identify drug targets with high probability of success by focusing on genes with both gain of function (GoF) and loss of function (LoF) mutations associated with opposing effects on phenotype (Bidirectional Effect Selected Targets, BEST). We find 98 BEST genes for a variety of indications. Drugs targeting those genes are 3.8-fold more likely to be approved than non-BEST genes. We focus on five genes (IGF1R, NPPC, NPR2, FGFR3, and SHOX) with evidence for bidirectional effects on stature. Rare protein-altering variants in those genes result in significantly increased risk for idiopathic short stature (ISS) (OR = 2.75, p = 3.99 × 10-8). Finally, using functional experiments, we demonstrate that adding an exogenous CNP analog (encoded by NPPC) rescues the phenotype, thus validating its potential as a therapeutic treatment for ISS. Our results show the value of looking for bidirectional effects to identify and validate drug targets.

Genomic analysis reveals a tight link between transcription factor dynamics and regulatory network architecture.

Although several studies have provided important insights into the general principles of biological networks, the link between network organization and the genome-scale dynamics of the underlying entities (genes, mRNAs, and proteins) and its role in systems behavior remain unclear. Here we show that transcription factor (TF) dynamics and regulatory network organization are tightly linked. By classifying TFs in the yeast regulatory network into three hierarchical layers (top, core, and bottom) and integrating diverse genome-scale datasets, we find that the TFs have static and dynamic properties that are similar within a layer and different across layers. At the protein level, the top-layer TFs are relatively abundant, long-lived, and noisy compared with the core- and bottom-layer TFs. Although variability in expression of top-layer TFs might confer a selective advantage, as this permits at least some members in a clonal cell population to initiate a response to changing conditions, tight regulation of the core- and bottom-layer TFs may minimize noise propagation and ensure fidelity in regulation. We propose that the interplay between network organization and TF dynamics could permit differential utilization of the same underlying network by distinct members of a clonal cell population.

Influence of genetic copy number variants of the human GLUT3 glucose transporter gene SLC2A3 on protein expression, glycolysis and rheumatoid arthritis risk: A genetic replication study.

OBJECTIVES: The gene encoding glucose transporter 3 (GLUT3, SLC2A3) is present in the human population at variable copy number. An overt disease phenotype of SLC2A3 copy number variants has not been reported; however, deletion of SLC2A3 has been previously reported to protect carriers from rheumatoid arthritis, implicating GLUT3 as a therapeutic target in rheumatoid arthritis. Here we aim to perform functional analysis of GLUT3 copy number variants in immune cells, and test the reported protective association of the GLUT3 copy number variants for rheumatoid arthritis in a genetic replication study. METHODS: Cells from genotyped healthy controls were analyzed for SLC2A3/GLUT3 expression and glycolysis capacity. We genotyped the SLC2A3 copy number variant in four independent cohorts of rheumatoid arthritis and controls and one cohort of multiple sclerosis and controls. RESULTS: Heterozygous deletion of SLC2A3 correlates directly with expression levels of GLUT3 and influences glycolysis rates in the human immune system. The frequency of the SLC2A3 copy number variant is not different between rheumatoid arthritis, multiple sclerosis and control groups. CONCLUSIONS: Despite a robust SLC2A3 gene copy number dependent phenotype, our study of large groups of rheumatoid arthritis cases and controls provides no evidence for rheumatoid arthritis disease protection in deletion carriers. These data emphasize the importance of well powered replication studies to confirm or refute genetic associations, particularly for relatively rare variants.

Conservation and evolutionary dynamics of the agr cell-to-cell communication system across firmicutes.

We present evidence that the agr cell-to-cell communication system is present across firmicutes, including the human pathogen Clostridium perfringens. Although we find that the agr system is evolutionarily conserved and that the general functions which it regulates are similar in different species, the individual regulated genes are not the same. This suggests that the regulatory network controlled by agr is dynamic and evolves rapidly.