European and multi-ancestry genome-wide association meta-analysis of atopic dermatitis highlights importance of systemic immune regulation.

Atopic dermatitis (AD) is a common inflammatory skin condition and prior genome-wide association studies (GWAS) have identified 71 associated loci. In the current study we conducted the largest AD GWAS to date (discovery N = 1,086,394, replication N = 3,604,027), combining previously reported cohorts with additional available data. We identified 81 loci (29 novel) in the European-only analysis (which all replicated in a separate European analysis) and 10 additional loci in the multi-ancestry analysis (3 novel). Eight variants from the multi-ancestry analysis replicated in at least one of the populations tested (European, Latino or African), while two may be specific to individuals of Japanese ancestry. AD loci showed enrichment for DNAse I hypersensitivity and eQTL associations in blood. At each locus we prioritised candidate genes by integrating multi-omic data. The implicated genes are predominantly in immune pathways of relevance to atopic inflammation and some offer drug repurposing opportunities.

Multi-ancestry genome-wide association analyses improve resolution of genes and pathways influencing lung function and chronic obstructive pulmonary disease risk.

Lung-function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry genome-wide association meta-analysis of lung function to date, comprising 580,869 participants, we identified 1,020 independent association signals implicating 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies for selected associated variants as well as trait and pathway-specific genetic risk scores to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

DeepPheWAS: an R package for phenotype generation and association analysis for phenome-wide association studies.

SUMMARY: DeepPheWAS is an R package for phenome-wide association studies that creates clinically curated composite phenotypes and integrates quantitative phenotypes from primary care data, longitudinal trajectories of quantitative measures, disease progression and drug response phenotypes. Tools are provided for efficient analysis of association with any genetic input, under any genetic model, with optional sex-stratified analysis, and for developing novel phenotypes. AVAILABILITY AND IMPLEMENTATION: The DeepPheWAS R package is freely available under GNU general public licence v3.0 from at https://github.com/Richard-Packer/DeepPheWAS. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Exome-wide analysis of copy number variation shows association of the human leukocyte antigen region with asthma in UK Biobank.

BACKGROUND: The role of copy number variants (CNVs) in susceptibility to asthma is not well understood. This is, in part, due to the difficulty of accurately measuring CNVs in large enough sample sizes to detect associations. The recent availability of whole-exome sequencing (WES) in large biobank studies provides an unprecedented opportunity to study the role of CNVs in asthma. METHODS: We called common CNVs in 49,953 individuals in the first release of UK Biobank WES using ClinCNV software. CNVs were tested for association with asthma in a stage 1 analysis comprising 7098 asthma cases and 36,578 controls from the first release of sequencing data. Nominally-associated CNVs were then meta-analysed in stage 2 with an additional 17,280 asthma cases and 115,562 controls from the second release of UK Biobank exome sequencing, followed by validation and fine-mapping. RESULTS: Five of 189 CNVs were associated with asthma in stage 2, including a deletion overlapping the HLA-DQA1 and HLA-DQB1 genes, a duplication of CHROMR/PRKRA, deletions within MUC22 and TAP2, and a duplication in FBRSL1. The HLA-DQA1, HLA-DQB1, MUC22 and TAP2 genes all reside within the human leukocyte antigen (HLA) region on chromosome 6. In silico analyses demonstrated that the deletion overlapping HLA-DQA1 and HLA-DQB1 is likely to be an artefact arising from under-mapping of reads from non-reference HLA haplotypes, and that the CHROMR/PRKRA and FBRSL1 duplications represent presence/absence of pseudogenes within the HLA region. Bayesian fine-mapping of the HLA region suggested that there are two independent asthma association signals. The variants with the largest posterior inclusion probability in the two credible sets were an amino acid change in HLA-DQB1 (glutamine to histidine at residue 253) and a multi-allelic amino acid change in HLA-DRB1 (presence/absence of serine, glycine or leucine at residue 11). CONCLUSIONS: At least two independent loci characterised by amino acid changes in the HLA-DQA1, HLA-DQB1 and HLA-DRB1 genes are likely to account for association of SNPs and CNVs in this region with asthma. The high divergence of haplotypes in the HLA can give rise to spurious CNVs, providing an important, cautionary tale for future large-scale analyses of sequencing data.

Pleiotropic associations of heterozygosity for the SERPINA1 Z allele in the UK Biobank.

Homozygosity for the SERPINA1 Z allele causes α1-antitrypsin deficiency, a rare condition that can cause lung and liver disease. However, the effects of Z allele heterozygosity on nonrespiratory phenotypes, and on lung function in the general population, remain unclear. We conducted a large, population-based study to determine Z allele effects on >2400 phenotypes in the UK Biobank (N=303 353). Z allele heterozygosity was strongly associated with increased height (ß=1.02 cm, p=3.91×10-68), and with other nonrespiratory phenotypes including increased risk of gall bladder disease, reduced risk of heart disease and lower blood pressure, reduced risk of osteoarthritis and reduced bone mineral density, increased risk of headache and enlarged prostate, as well as with blood biomarkers of liver function. Heterozygosity was associated with higher height-adjusted forced expiratory volume in 1 s (FEV1) (ß=19.36 mL, p=9.21×10-4) and FEV1/forced vital capacity (ß=0.0031, p=1.22×10-5) in nonsmokers, whereas in smokers, this protective effect was abolished. Furthermore, we show for the first time that sex modifies the association of the Z allele on lung function. We conclude that Z allele heterozygosity and homozygosity exhibit opposing effects on lung function in the UK population, and that these associations are modified by smoking and sex. In exploratory analyses, heterozygosity for the Z allele also showed pleiotropic associations with nonrespiratory health-related traits and disease risk.

SCN3A-Related Neurodevelopmental Disorder: A Spectrum of Epilepsy and Brain Malformation.

OBJECTIVE: Pathogenic variants in SCN3A, encoding the voltage-gated sodium channel subunit Nav1.3, cause severe childhood onset epilepsy and malformation of cortical development. Here, we define the spectrum of clinical, genetic, and neuroimaging features of SCN3A-related neurodevelopmental disorder. METHODS: Patients were ascertained via an international collaborative network. We compared sodium channels containing wild-type versus variant Nav1.3 subunits coexpressed with ß1 and ß2 subunits using whole-cell voltage clamp electrophysiological recordings in a heterologous mammalian system (HEK-293T cells). RESULTS: Of 22 patients with pathogenic SCN3A variants, most had treatment-resistant epilepsy beginning in the first year of life (16/21, 76%; median onset, 2 weeks), with severe or profound developmental delay (15/20, 75%). Many, but not all (15/19, 79%), exhibited malformations of cortical development. Pathogenic variants clustered in transmembrane segments 4 to 6 of domains II to IV. Most pathogenic missense variants tested (10/11, 91%) displayed gain of channel function, with increased persistent current and/or a leftward shift in the voltage dependence of activation, and all variants associated with malformation of cortical development exhibited gain of channel function. One variant (p.Ile1468Arg) exhibited mixed effects, with gain and partial loss of function. Two variants demonstrated loss of channel function. INTERPRETATION: Our study defines SCN3A-related neurodevelopmental disorder along a spectrum of severity, but typically including epilepsy and severe or profound developmental delay/intellectual disability. Malformations of cortical development are a characteristic feature of this unusual channelopathy syndrome, present in >75% of affected individuals. Gain of function at the channel level in developing neurons is likely an important mechanism of disease pathogenesis. ANN NEUROL 2020;88:348-362.

Variants associated with HHIP expression have sex-differential effects on lung function.

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P<5x10 -8) interactions with sex on lung function, and 21 showed suggestive interactions (P<1x10 -6). The strongest signal, from rs7697189 (chr4:145436894) on forced expiratory volume in 1 second (FEV 1) (P=3.15x10 -15), was replicated (P=0.016) in SpiroMeta. The C allele increased FEV 1 more in males (untransformed FEV 1 ß=0.028 [SE 0.0022] litres) than females (ß=0.009 [SE 0.0014] litres), and this effect was not accounted for by differential effects on height, smoking or pubertal age. rs7697189 resides upstream of the hedgehog-interacting protein ( HHIP) gene and was previously associated with lung function and HHIP lung expression. We found HHIP expression was significantly different between the sexes (P=6.90x10 -6), but we could not detect sex differential effects of rs7697189 on expression. Conclusions: We identified a novel genotype-by-sex interaction at a putative enhancer region upstream of the HHIP gene. Establishing the mechanism by which HHIP SNPs have different effects on lung function in males and females will be important for our understanding of lung health and diseases in both sexes.

TMX2 Is a Crucial Regulator of Cellular Redox State, and Its Dysfunction Causes Severe Brain Developmental Abnormalities.

The redox state of the neural progenitors regulates physiological processes such as neuronal differentiation and dendritic and axonal growth. The relevance of endoplasmic reticulum (ER)-associated oxidoreductases in these processes is largely unexplored. We describe a severe neurological disorder caused by bi-allelic loss-of-function variants in thioredoxin (TRX)-related transmembrane-2 (TMX2); these variants were detected by exome sequencing in 14 affected individuals from ten unrelated families presenting with congenital microcephaly, cortical polymicrogyria, and other migration disorders. TMX2 encodes one of the five TMX proteins of the protein disulfide isomerase family, hitherto not linked to human developmental brain disease. Our mechanistic studies on protein function show that TMX2 localizes to the ER mitochondria-associated membranes (MAMs), is involved in posttranslational modification and protein folding, and undergoes physical interaction with the MAM-associated and ER folding chaperone calnexin and ER calcium pump SERCA2. These interactions are functionally relevant because TMX2-deficient fibroblasts show decreased mitochondrial respiratory reserve capacity and compensatory increased glycolytic activity. Intriguingly, under basal conditions TMX2 occurs in both reduced and oxidized monomeric form, while it forms a stable dimer under treatment with hydrogen peroxide, recently recognized as a signaling molecule in neural morphogenesis and axonal pathfinding. Exogenous expression of the pathogenic TMX2 variants or of variants with an in vitro mutagenized TRX domain induces a constitutive TMX2 polymerization, mimicking an increased oxidative state. Altogether these data uncover TMX2 as a sensor in the MAM-regulated redox signaling pathway and identify it as a key adaptive regulator of neuronal proliferation, migration, and organization in the developing brain.

Author Correction: New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries.

Reduced lung function predicts mortality and is key to the diagnosis of chronic obstructive pulmonary disease (COPD). In a genome-wide association study in 400,102 individuals of European ancestry, we define 279 lung function signals, 139 of which are new. In combination, these variants strongly predict COPD in independent populations. Furthermore, the combined effect of these variants showed generalizability across smokers and never smokers, and across ancestral groups. We highlight biological pathways, known and potential drug targets for COPD and, in phenome-wide association studies, autoimmune-related and other pleiotropic effects of lung function-associated variants. This new genetic evidence has potential to improve future preventive and therapeutic strategies for COPD.

Missense Mutations of the Pro65 Residue of PCGF2 Cause a Recognizable Syndrome Associated with Craniofacial, Neurological, Cardiovascular, and Skeletal Features.

PCGF2 encodes the polycomb group ring finger 2 protein, a transcriptional repressor involved in cell proliferation, differentiation, and embryogenesis. PCGF2 is a component of the polycomb repressive complex 1 (PRC1), a multiprotein complex which controls gene silencing through histone modification and chromatin remodelling. We report the phenotypic characterization of 13 patients (11 unrelated individuals and a pair of monozygotic twins) with missense mutations in PCGF2. All the mutations affected the same highly conserved proline in PCGF2 and were de novo, excepting maternal mosaicism in one. The patients demonstrated a recognizable facial gestalt, intellectual disability, feeding problems, impaired growth, and a range of brain, cardiovascular, and skeletal abnormalities. Computer structural modeling suggests the substitutions alter an N-terminal loop of PCGF2 critical for histone biding. Mutant PCGF2 may have dominant-negative effects, sequestering PRC1 components into complexes that lack the ability to interact efficiently with histones. These findings demonstrate the important role of PCGF2 in human development and confirm that heterozygous substitutions of the Pro65 residue of PCGF2 cause a recognizable syndrome characterized by distinctive craniofacial, neurological, cardiovascular, and skeletal features.

De novo mutations in GRIN1 cause extensive bilateral polymicrogyria.

Polymicrogyria is a malformation of cortical development. The aetiology of polymicrogyria remains poorly understood. Using whole-exome sequencing we found de novo heterozygous missense GRIN1 mutations in 2 of 57 parent-offspring trios with polymicrogyria. We found nine further de novo missense GRIN1 mutations in additional cortical malformation patients. Shared features in the patients were extensive bilateral polymicrogyria associated with severe developmental delay, postnatal microcephaly, cortical visual impairment and intractable epilepsy. GRIN1 encodes GluN1, the essential subunit of the N-methyl-d-aspartate receptor. The polymicrogyria-associated GRIN1 mutations tended to cluster in the S2 region (part of the ligand-binding domain of GluN1) or the adjacent M3 helix. These regions are rarely mutated in the normal population or in GRIN1 patients without polymicrogyria. Using two-electrode and whole-cell voltage-clamp analysis, we showed that the polymicrogyria-associated GRIN1 mutations significantly alter the in vitro activity of the receptor. Three of the mutations increased agonist potency while one reduced proton inhibition of the receptor. These results are striking because previous GRIN1 mutations have generally caused loss of function, and because N-methyl-d-aspartate receptor agonists have been used for many years to generate animal models of polymicrogyria. Overall, our results expand the phenotypic spectrum associated with GRIN1 mutations and highlight the important role of N-methyl-d-aspartate receptor signalling in the pathogenesis of polymicrogyria.

Deep Brain Stimulation in Three Related Cases of North Sea Progressive Myoclonic Epilepsy from South Africa.

We report on a white Afrikaans family from eastern South Africa with three members affected with North Sea progressive myoclonus epilepsy, resulting from a homozygous founder GOSR2 mutation (c.430G>T, p.Gly144Trp). The mutation was identified by exomic sequencing in a research study investigating childhood onset ataxias. All three subjects presented with ataxia, tremor, early gait difficulties, and myoclonic and generalized tonic clonic (GTC) epilepsy. Each patient underwent deep brain stimulation of the caudal Zona Incerta before coming to the attention of the authors. In each case there was a reduction in GTC seizures, and two patients exhibited a reduction in involuntary movements, as evaluated during long-term follow-up. In one case there was an improvement in gait and stance when assessed while the stimulation was on.

Behavioral Variation in Gorillas: Evidence of Potential Cultural Traits.

The question of whether any species except humans exhibits culture has generated much debate, partially due to the difficulty of providing conclusive evidence from observational studies in the wild. A starting point for demonstrating the existence of culture that has been used for many species including chimpanzees and orangutans is to show that there is geographic variation in the occurrence of particular behavioral traits inferred to be a result of social learning and not ecological or genetic influences. Gorillas live in a wide variety of habitats across Africa and they exhibit flexibility in diet, behavior, and social structure. Here we apply the 'method of exclusion' to look for the presence/absence of behaviors that could be considered potential cultural traits in well-habituated groups from five study sites of the two species of gorillas. Of the 41 behaviors considered, 23 met the criteria of potential cultural traits, of which one was foraging related, nine were environment related, seven involved social interactions, five were gestures, and one was communication related. There was a strong positive correlation between behavioral dissimilarity and geographic distance among gorilla study sites. Roughly half of all variation in potential cultural traits was intraspecific differences (i.e. variability among sites within a species) and the other 50% of potential cultural traits were differences between western and eastern gorillas. Further research is needed to investigate if the occurrence of these traits is influenced by social learning. These findings emphasize the importance of investigating cultural traits in African apes and other species to shed light on the origin of human culture.

De novo point mutations in patients diagnosed with ataxic cerebral palsy.

Cerebral palsy is a sporadic disorder with multiple likely aetiologies, but frequently considered to be caused by birth asphyxia. Genetic investigations are rarely performed in patients with cerebral palsy and there is little proven evidence of genetic causes. As part of a large project investigating children with ataxia, we identified four patients in our cohort with a diagnosis of ataxic cerebral palsy. They were investigated using either targeted next generation sequencing or trio-based exome sequencing and were found to have mutations in three different genes, KCNC3, ITPR1 and SPTBN2. All the mutations were de novo and associated with increased paternal age. The mutations were shown to be pathogenic using a combination of bioinformatics analysis and in vitro model systems. This work is the first to report that the ataxic subtype of cerebral palsy can be caused by de novo dominant point mutations, which explains the sporadic nature of these cases. We conclude that at least some subtypes of cerebral palsy may be caused by de novo genetic mutations and patients with a clinical diagnosis of cerebral palsy should be genetically investigated before causation is ascribed to perinatal asphyxia or other aetiologies.

CGAT: computational genomics analysis toolkit.

Computational genomics seeks to draw biological inferences from genomic datasets, often by integrating and contextualizing next-generation sequencing data. CGAT provides an extensive suite of tools designed to assist in the analysis of genome scale data from a range of standard file formats. The toolkit enables filtering, comparison, conversion, summarization and annotation of genomic intervals, gene sets and sequences. The tools can both be run from the Unix command line and installed into visual workflow builders, such as Galaxy.