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 588,452 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.

Chronic respiratory diseases, such as chronic obstructive pulmonary disease (COPD), asthma and interstitial lung diseases are frequently occurring disorders with a polygenic basis that account for a large global burden of morbidity and mortality. Recent large-scale genetic epidemiology studies have identified associations between genetic variation and individual respiratory diseases and linked specific genetic variants to quantitative traits related to lung function. These associations have improved our understanding of the genetic basis and mechanisms underlying common lung diseases. Moreover, examining the overlap between genetic associations of different respiratory conditions, along with evidence for gene–environment interactions, has yielded additional biological insights into affected molecular pathways. This genetic information could inform the assessment of respiratory disease risk and contribute to stratified treatment approaches.

Naturally occurring autopolyploid species, such as the autotetraploid potato Solanum tuberosum, face a variety of challenges during meiosis. These include proper pairing, recombination and correct segregation of multiple homologous chromosomes, which can form complex multivalent configurations at metaphase I, and in turn alter allelic segregation ratios through double reduction. Here, we present a reference map of meiotic stages in diploid and tetraploid S. tuberosum using fluorescence in situ hybridisation (FISH) to differentiate individual meiotic chromosomes 1 and 2. A diploid-like behaviour at metaphase I involving bivalent configurations was predominant in all three tetraploid varieties. The crossover frequency per bivalent was significantly reduced in the tetraploids compared with a diploid variety, which likely indicates meiotic adaptation to the autotetraploid state. Nevertheless, bivalents were accompanied by a substantial frequency of multivalents, which varied by variety and by chromosome (7–48%). We identified possible sites of synaptic partner switching, leading to multivalent formation, and found potential defects in the polymerisation and/or maintenance of the synaptonemal complex in tetraploids. These findings demonstrate the rise of S. tuberosum as a model for autotetraploid meiotic recombination research and highlight constraints on meiotic chromosome configurations and chiasma frequencies as an important feature of an evolved autotetraploid meiosis.

Background: While the possible sources underlying the so-called ‘missing heritability’ evident in current genome-wide association studies (GWAS) of complex traits have been actively pursued in recent years, resolving this mystery remains a challenging task. Studying heritability of genome-wide gene expression traits can shed light on the goal of understanding the relationship between phenotype and genotype. Here we used microarray gene expression measurements of lymphoblastoid cell lines and genome-wide SNP genotype data from 210 HapMap individuals to examine the heritability of gene expression traits. Results: Heritability levels for expression of 10,720 genes were estimated by applying variance component model analyses and 1,043 expression quantitative loci (eQTLs) were detected. Our results indicate that gene expression traits display a bimodal distribution of heritability, one peak close to 0% and the other summit approaching 100%. Such a pattern of the within-population variability of gene expression heritability is common among different HapMap populations of unrelated individuals but different from that obtained in the CEU and YRI trio samples. Higher heritability levels are shown by housekeeping genes and genes associated with cis eQTLs. Both cis and trans eQTLs make comparable cumulative contributions to the heritability. Finally, we modelled gene-gene interactions (epistasis) for genes with multiple eQTLs and revealed that epistasis was not prevailing in all genes but made a substantial contribution in explaining total heritability for some genes analysed. Conclusions: We utilised a mixed effect model analysis for estimating genetic components from population based samples. On basis of analyses of genome-wide gene expression from four HapMap populations, we demonstrated detailed exploitation of the distribution of genetic heritabilities for expression traits from different populations, and highlighted the importance of studying interaction at the gene expression level as an important source of variation underlying missing heritability.

Thyroid hormones play a critical role in regulation of multiple physiological functions and thyroid dysfunction is associated with substantial morbidity. Here, we use electronic health records to undertake a genome-wide association study of thyroid-stimulating hormone (TSH) levels, with a total sample size of 247,107. We identify 158 novel genetic associations, more than doubling the number of known associations with TSH, and implicate 112 putative causal genes, of which 76 are not previously implicated. A polygenic score for TSH is associated with TSH levels in African, South Asian, East Asian, Middle Eastern and admixed American ancestries, and associated with hypothyroidism and other thyroid disease in South Asians. In Europeans, the TSH polygenic score is associated with thyroid disease, including thyroid cancer and age-of-onset of hypothyroidism and hyperthyroidism. We develop pathway-specific genetic risk scores for TSH levels and use these in phenome-wide association studies to identify potential consequences of pathway perturbation. Together, these findings demonstrate the potential utility of genetic associations to inform future therapeutics and risk prediction for thyroid diseases.

DNA methylation in the genome plays a fundamental role in the regulation of gene expression and is widespread in the genome of eukaryotic species. For example, in higher vertebrates, there is a “global” methylation pattern involving complete methylation of CpG sites genome-wide, except in promoter regions that are typically enriched for CpG dinucleotides, or so called “CpG islands.” Here, we comprehensively examined and compared the distribution of CpG sites within ten model eukaryotic species and linked the observed patterns to the role of DNA methylation in controlling gene transcription. The analysis revealed two distinct but conserved methylation patterns for gene promoters in human and mouse genomes, involving genes with distinct distributions of promoter CpGs and gene expression patterns. Comparative analysis with four other higher vertebrates revealed that the primary regulatory role of the DNA methylation system is highly conserved in higher vertebrates.