We examined the potential gene x gene interactions and gene x smoking interactions in rheumatoid arthritis (RA) using the candidate gene data sets provided by Genetic Analysis Workshop 15 Problem 2. The multifactor dimensionality reduction (MDR) method was used to test gene x gene interactions among candidate genes. The case-only sample was used to test gene x smoking interactions. The best predictive model was the single-locus model with single-nucleotide polymorphism (SNP) rs2476601 in gene PTPN22. However, no clear gene x gene interaction was identified. Substantial departure from multiplicativity was observed between smoking and SNPs in genes CTLA4, PADI4, MIF, and SNPs on chromosome 5 and one haplotype of PTPN22. The strongest evidence of association was identified between the PTPN22 gene and RA status, which was consistently detected in single SNP association, gene x gene interaction and gene x smoking interaction analyses.

The genetic basis of supraventricular and ventricular ectopy (SVE, VE) remains largely uncharacterized, despite established genetic mechanisms of arrhythmogenesis. To identify novel genetic variants associated with SVE/VE in ancestrally diverse human populations, we conducted a genome-wide association study of electrocardiographically identified SVE and VE in five cohorts including approximately 43,000 participants of African, European and Hispanic/Latino ancestry. In thirteen ancestry-stratified subgroups, we tested multivariable-adjusted associations of SVE and VE with single nucleotide polymorphism (SNP) dosage. We combined subgroup-specific association estimates in inverse variance-weighted, fixed-effects and Bayesian meta-analyses. We also combined fixed-effects meta-analytic t-test statistics for SVE and VE in multi-trait SNP association analyses. No loci reached genome-wide significance in trans-ethnic meta-analyses. However, we found genome-wide significant SNPs intronic to an apoptosis-enhancing gene previously associated with QRS interval duration (FAF1; lead SNP rs7545860; effect allele frequency = 0.02; P = 2.0 × 10-8) in multi-trait analysis among European ancestry participants and near a locus encoding calcium-dependent glycoproteins (DSC3; lead SNP rs8086068; effect allele frequency = 0.17) in meta-analysis of SVE (P = 4.0 × 10-8) and multi-trait analysis (P = 2.9 × 10-9) among African ancestry participants. The novel findings suggest several mechanisms by which genetic variation may predispose to ectopy in humans and highlight the potential value of leveraging pleiotropy in future studies of ectopy-related phenotypes.

Primary biliary cholangitis (PBC) is an autoimmune liver disease with a strong hereditary component. Here, we report a genome-wide association study that included 1,122 PBC cases and 4,036 controls of Han Chinese descent, with subsequent replication in a separate cohort of 907 PBC cases and 2,127 controls. Our results show genome-wide association of 14 PBC risk loci including previously identified 6p21 (HLA-DRA and DPB1), 17q12 (ORMDL3), 3q13.33 (CD80), 2q32.3 (STAT1/STAT4), 3q25.33 (IL12A), 4q24 (NF-κB) and 22q13.1 (RPL3/SYNGR1). We also identified variants in IL21, IL21R, CD28/CTLA4/ICOS, CD58, ARID3A and IL16 as novel PBC risk loci. These new findings and histochemical studies showing enhanced expression of IL21 and IL21R in PBC livers (particularly in the hepatic portal tracks) support a disease mechanism in which the deregulation of the IL21 signalling pathway, in addition to CD4 T-cell activation and T-cell co-stimulation are critical components in the development of PBC.

To profile the differential gene expression of the KEGG Adipocytokine Signaling pathway in omental compared to subcutaneous tissue in normal pregnancy.

Genome-wide association studies (GWASs) have identified loci for erythrocyte traits in primarily European ancestry populations. We conducted GWAS meta-analyses of six erythrocyte traits in 71,638 individuals from European, East Asian, and African ancestries using a Bayesian approach to account for heterogeneity in allelic effects and variation in the structure of linkage disequilibrium between ethnicities. We identified seven loci for erythrocyte traits including a locus (RBPMS/GTF2E2) associated with mean corpuscular hemoglobin and mean corpuscular volume. Statistical fine-mapping at this locus pointed to RBPMS at this locus and excluded nearby GTF2E2. Using zebrafish morpholino to evaluate loss of function, we observed a strong in vivo erythropoietic effect for RBPMS but not for GTF2E2, supporting the statistical fine-mapping at this locus and demonstrating that RBPMS is a regulator of erythropoiesis. Our findings show the utility of trans-ethnic GWASs for discovery and characterization of genetic loci influencing hematologic traits.

Educational attainment is widely used as a surrogate for socioeconomic status (SES). Low SES is a risk factor for hypertension and high blood pressure (BP). To identify novel BP loci, we performed multi-ancestry meta-analyses accounting for gene-educational attainment interactions using two variables, "Some College" (yes/no) and "Graduated College" (yes/no). Interactions were evaluated using both a 1 degree of freedom (DF) interaction term and a 2DF joint test of genetic and interaction effects. Analyses were performed for systolic BP, diastolic BP, mean arterial pressure, and pulse pressure. We pursued genome-wide interrogation in Stage 1 studies (N = 117 438) and follow-up on promising variants in Stage 2 studies (N = 293 787) in five ancestry groups. Through combined meta-analyses of Stages 1 and 2, we identified 84 known and 18 novel BP loci at genome-wide significance level (P < 5 × 10-8). Two novel loci were identified based on the 1DF test of interaction with educational attainment, while the remaining 16 loci were identified through the 2DF joint test of genetic and interaction effects. Ten novel loci were identified in individuals of African ancestry. Several novel loci show strong biological plausibility since they involve physiologic systems implicated in BP regulation. They include genes involved in the central nervous system-adrenal signaling axis (ZDHHC17, CADPS, PIK3C2G), vascular structure and function (GNB3, CDON), and renal function (HAS2 and HAS2-AS1, SLIT3). Collectively, these findings suggest a role of educational attainment or SES in further dissection of the genetic architecture of BP.

Allele expression (AE) analysis robustly measures cis-regulatory effects. Here, we present and demonstrate the utility of a vast AE resource generated from the GTEx v8 release, containing 15,253 samples spanning 54 human tissues for a total of 431 million measurements of AE at the SNP level and 153 million measurements at the haplotype level. In addition, we develop an extension of our tool phASER that allows effect sizes of cis-regulatory variants to be estimated using haplotype-level AE data. This AE resource is the largest to date, and we are able to make haplotype-level data publicly available. We anticipate that the availability of this resource will enable future studies of regulatory variation across human tissues.

A common genetic variant near MBOAT7 (rs641738C>T) has been previously associated with hepatic fat and advanced histology in NAFLD; however, these findings have not been consistently replicated in the literature. We aimed to establish whether rs641738C>T is a risk factor across the spectrum of NAFLD and to characterise its role in the regulation of related metabolic phenotypes through a meta-analysis.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation.

Autotetraploid rice hybrids have great potential to increase the production, but hybrid sterility is a major hindrance in the utilization of hybrid vigor in polyploid rice, which is mainly caused by pollen abortion. Our previous study showed that double pollen fertility neutral genes, Sa-n and Sb-n, can overcome hybrid sterility in autotetraploid rice. Here, we used an autotetraploid rice line harboring double neutral genes to develop hybrids by crossing with auto- and neo-tetraploid rice, and evaluated heterosis and its underlying molecular mechanism.

Human longevity is heritable, but genome-wide association (GWA) studies have had limited success. Here, we perform two meta-analyses of GWA studies of a rigorous longevity phenotype definition including 11,262/3484 cases surviving at or beyond the age corresponding to the 90th/99th survival percentile, respectively, and 25,483 controls whose age at death or at last contact was at or below the age corresponding to the 60th survival percentile. Consistent with previous reports, rs429358 (apolipoprotein E (ApoE) ε4) is associated with lower odds of surviving to the 90th and 99th percentile age, while rs7412 (ApoE ε2) shows the opposite. Moreover, rs7676745, located near GPR78, associates with lower odds of surviving to the 90th percentile age. Gene-level association analysis reveals a role for tissue-specific expression of multiple genes in longevity. Finally, genetic correlation of the longevity GWA results with that of several disease-related phenotypes points to a shared genetic architecture between health and longevity.

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11-34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count <10 (including singletons) in the imputation target samples, average information content rescued was >86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations.

Aortic calcification is an important independent predictor of future cardiovascular events. We performed a genome-wide association meta-analysis to determine SNPs associated with the extent of abdominal aortic calcification (n = 9,417) or descending thoracic aortic calcification (n = 8,422). Two genetic loci, HDAC9 and RAP1GAP, were associated with abdominal aortic calcification at a genome-wide level (P < 5.0 × 10-8). No SNPs were associated with thoracic aortic calcification at the genome-wide threshold. Increased expression of HDAC9 in human aortic smooth muscle cells promoted calcification and reduced contractility, while inhibition of HDAC9 in human aortic smooth muscle cells inhibited calcification and enhanced cell contractility. In matrix Gla protein-deficient mice, a model of human vascular calcification, mice lacking HDAC9 had a 40% reduction in aortic calcification and improved survival. This translational genomic study identifies the first genetic risk locus associated with calcification of the abdominal aorta and describes a previously unknown role for HDAC9 in the development of vascular calcification.

Sleep disordered breathing (SDB) is a common disorder that results in oxidative stress and inflammation and is associated with multiple age-related health outcomes. Epigenetic age acceleration is a DNA methylation (DNAm)-based marker of fast biological aging. We examined the associations of SDB traits with epigenetic age acceleration.

Populations of Mexican American ancestry are at an increased risk for nonalcoholic fatty liver disease. The objective of this study was to determine whether loci in known and novel genes were associated with variation in aspartate aminotransferase (AST) (n = 3,644), alanine aminotransferase (ALT) (n = 3,595), and gamma-glutamyl transferase (GGT) (n = 1,577) levels by conducting the first genome-wide association study (GWAS) of liver enzymes, which commonly measure liver function, in individuals of Mexican American ancestry.

The management of diabetic ulcer (DU) to rescue stalled wound healing remains a paramount clinical challenge due to the spatially and temporally coupled pathological wound microenvironment that features hyperglycemia, biofilm infection, hypoxia and excessive oxidative stress. Here we present a pH-switchable nanozyme cascade catalysis (PNCC) strategy for spatial-temporal modulation of pathological wound microenvironment to rescue stalled healing in DU. The PNCC is demonstrated by employing the nanozyme of clinically approved iron oxide nanoparticles coated with a shell of glucose oxidase (Fe3O4-GOx). The Fe3O4-GOx possesses intrinsic glucose oxidase (GOx), catalase (CAT) and peroxidase (POD)-like activities, and can catalyze pH-switchable glucose-initiated GOx/POD and GOx/CAT cascade reaction in acidic and neutral environment, respectively. Specifically, the GOx/POD cascade reaction generating consecutive fluxes of toxic hydroxyl radical spatially targets the acidic biofilm (pH ~ 5.5), and eradicates biofilm to shorten the inflammatory phase and initiate normal wound healing processes. Furthermore, the GOx/CAT cascade reaction producing consecutive fluxes of oxygen spatially targets the neutral wound tissue, and accelerates the proliferation and remodeling phases of wound healing by addressing the issues of hyperglycemia, hypoxia, and excessive oxidative stress. The shortened inflammatory phase temporally coupled with accelerated proliferation and remodeling phases significantly speed up the normal orchestrated wound-healing cascades. Remarkably, this Fe3O4-GOx-instructed spatial-temporal remodeling of DU microenvironment enables complete re-epithelialization of biofilm-infected wound in diabetic mice within 15 days while minimizing toxicity to normal tissues, exerting great transformation potential in clinical DU management. The proposed PNCC concept offers a new perspective for complex pathological microenvironment remodeling, and may provide a powerful modality for the treatment of microenvironment-associated diseases.

Cultivated rice varieties are all diploid, and polyploidization of rice has long been desired because of its advantages in genome buffering, vigorousness, and environmental robustness. However, a workable route remains elusive. Here, we describe a practical strategy, namely de novo domestication of wild allotetraploid rice. By screening allotetraploid wild rice inventory, we identified one genotype of Oryza alta (CCDD), polyploid rice 1 (PPR1), and established two important resources for its de novo domestication: (1) an efficient tissue culture, transformation, and genome editing system and (2) a high-quality genome assembly discriminated into two subgenomes of 12 chromosomes apiece. With these resources, we show that six agronomically important traits could be rapidly improved by editing O. alta homologs of the genes controlling these traits in diploid rice. Our results demonstrate the possibility that de novo domesticated allotetraploid rice can be developed into a new staple cereal to strengthen world food security.

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Mitochondrial DNA (mtDNA) is present in multiple copies in human cells. We evaluated cross-sectional associations of whole blood mtDNA copy number (CN) with several cardiometabolic disease traits in 408,361 participants of multiple ancestries in TOPMed and UK Biobank. Age showed a threshold association with mtDNA CN: among younger participants (<65 years of age), each additional 10 years of age was associated with 0.03 standard deviation (s.d.) higher level of mtDNA CN (P = 0.0014) versus a 0.14 s.d. lower level of mtDNA CN (P = 1.82 × 10-13) among older participants (≥65 years). At lower mtDNA CN levels, we found age-independent associations with increased odds of obesity (P = 5.6 × 10-238), hypertension (P = 2.8 × 10-50), diabetes (P = 3.6 × 10-7), and hyperlipidemia (P = 6.3 × 10-5). The observed decline in mtDNA CN after 65 years of age may be a key to understanding age-related diseases.