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.

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.

Reversible modification of proteins with linkage-specific ubiquitin chains is critical for intracellular signaling. Information on physiological roles and underlying mechanisms of particular ubiquitin linkages during human development are limited. Here, relying on genomic constraint scores, we identify 10 patients with multiple congenital anomalies caused by hemizygous variants in OTUD5, encoding a K48/K63 linkage-specific deubiquitylase. By studying these mutations, we find that OTUD5 controls neuroectodermal differentiation through cleaving K48-linked ubiquitin chains to counteract degradation of select chromatin regulators (e.g., ARID1A/B, histone deacetylase 2, and HCF1), mutations of which underlie diseases that exhibit phenotypic overlap with OTUD5 patients. Loss of OTUD5 during differentiation leads to less accessible chromatin at neuroectodermal enhancers and aberrant gene expression. Our study describes a previously unidentified disorder we name LINKED (LINKage-specific deubiquitylation deficiency-induced Embryonic Defects) syndrome and reveals linkage-specific ubiquitin cleavage from chromatin remodelers as an essential signaling mode that coordinates chromatin remodeling during embryogenesis.

Neo-tetraploid rice (NTR) is a useful new germplasm that developed from the descendants of the autotetraploid rice (ATR) hybrids. NTR showed improved fertility and yield potential, and produced high yield heterosis when crossed with indica ATR for commercial utilization. However, their classification, population structure and genomic feature remain elusive.

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%.

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.

The adaptive immune system uses several strategies to generate a repertoire of T cell receptors (TCR) with sufficient diversity to recognize the universe of potential pathogens. However, it remains unclear how differences in the T cell receptor (TCR) contribute to heterogeneity in T cell state. In this study, we used polychromatic flow cytometry to isolate highly pure CD4+/CD8+ naive and memory T cells, and applied deep sequencing to characterize corresponding TCR β-chain (TCRβ) complementary-determining region 3 (CDR3) repertoires. We find that shorter TCRβ CDR3s with fewer insertions were highly enriched during thymic selection. Antigen-experienced T cells (memory T cells) harbor shorter CDR3s vs. naive T cells. Moreover, the public TCRβ CDR3 clonotypes within cell subsets or interindividual tend to have shorter CDR3 length and a significantly larger size compared with "private" clonotypes. Taken together, shorter CDR3s highly enriched during thymic selection and antigen-driven selection, and further enriched in public T-cell responses. These results indicated that it may be evolutionary pressures drive short CDR3s to recognize most of antigen in nature.

Neo-tetraploid rice (NTR) is a new tetraploid rice germplasm that developed from the crossing and directional selection of different autotetraploid rice lines, which showed high fertility and promising yield potential. However, systematic yield assessment, genome composition and functional variations associated with fertility and yield remain elusive.

Gastric cancer (GC) remains an important cancer worldwide. Further understanding of the molecular mechanisms of gastric carcinogenesis will enhance the diagnosis and treatment of GC.

Autotetraploid rice exhibited hybrid vigor and greater genetic variation compared to diploid rice, but low pollen fertility is a major hindrance for its utilization. Our previous analysis revealed that large number of pollen fertility genes were exhibited down-regulation in autotetraploid rice. Hence, it is of utmost importance to reveal the expression patterns of pollen fertility genes with high accuracy. To find stable reference genes for autotetraploid rice, we compared the pollen development stages between diploid and autotetraploid rice, and 14 candidate genes were selected based on transcriptome analysis to evaluate their expression levels. Autotetraploid rice (i.e. Taichung65-4x) displayed lower seed set (40.40%) and higher percentage of abnormalities during the pollen development process than its diploid counterpart. To detect the candidate reference genes for pollen development of autotetraploid and diploid rice, we used five different algorithms, including NormFinder, BestKeeper, ΔCt method, geNorm and Re-Finder to evaluate their expression patterns stability. Consequently, we identified two genes, Cytochrome b5 and CPI, as the best candidate reference genes for qRT-PCR normalization in autotetraploid and diploid rice during pre-meiosis, meiosis, single microspore and bicellular pollen development stages. However, Cytochrome b5 was found to be the most stably expressed gene during different pollen development stages in autotetraploid rice. The results of our study provide a platform for subsequent gene expression analyses in autotetraploid rice, which could also be used in other polyploid plants.

Hypertrophic cardiomyopathy (HCM) represents one of the most common inherited heart diseases. To identify key molecules involved in the development of HCM, gene expression patterns of the heart tissue samples in HCM patients from multiple microarray and RNA-seq platforms were investigated.

Regulatory B cells (Breg) are considered as immunosuppressive cells. Different subsets of Breg cells have been identified both in human beings and in mice. However, there is a lack of unique markers to identify Breg cells, and the heterogeneity of Breg cells in different organs needs to be further illuminated. In this study, we performed high-throughput single-cell RNA sequencing (scRNA-seq) and single-cell B-cell receptor sequencing (scBCR-seq) of B cells from the murine spleen, liver, mesenteric lymph nodes, bone marrow, and peritoneal cavity to better define the phenotype of these cells. Breg cells were identified based on the expression of immunosuppressive genes and IL-10-producing B (B10) cell-related genes, to define B10 and non-B10 subsets in Breg cells based on the score of the B10 gene signatures. Moreover, we characterized 19 common genes significantly expressed in Breg cells, including Fcrl5, Zbtb20, Ccdc28b, Cd9, and Ptpn22, and further analyzed the transcription factor activity in defined Breg cells. Last, a BCR analysis was used to determine the clonally expanded clusters and the relationship of Breg cells across different organs. We demonstrated that Atf3 may potentially modulate the function of Breg cells as a transcription factor and that seven organ-specific subsets of Breg cells are found. Depending on gene expression and functional modules, non-B10 Breg cells exhibited activated the TGF-β pathway, thus suggesting that non-B10 Breg cells have specific immunosuppressive properties different from conventional B10 cells. In conclusion, our work provides new insights into Breg cells and illustrates their transcriptional profiles and BCR repertoire in different organs under physiological conditions.

Following muscle injury, fibro-adipogenic progenitors (FAPs) are rapidly activated and undergo apoptosis at the resolution stage, which is required for proper muscle regeneration. When excessive FAPs remain, it contributes to fibrotic and fatty infiltration, impairing muscle recovery. Mechanisms controlling FAP apoptosis remain poorly defined. We hypothesized that AMP-activated protein kinase (AMPK) in FAPs mediates their apoptosis during the muscle regeneration.

Hepatocellular carcinoma (HCC) is one of the most common primary hepatic malignancies. E2F transcription factors play an important role in the tumorigenesis and progression of HCC, mainly through the RB/E2F pathway. Prognostic models for HCC based on gene signatures have been developed rapidly in recent years; however, their discriminating ability at the single-cell level remains elusive, which could reflect the underlying mechanisms driving the sample bifurcation. In this study, we constructed and validated a predictive model based on E2F expression, successfully stratifying patients with HCC into two groups with different survival risks. Then we used a single-cell dataset to test the discriminating ability of the predictive model on infiltrating T cells, demonstrating remarkable cellular heterogeneity as well as altered cell fates. We identified distinct cell subpopulations with diverse molecular characteristics. We also found that the distribution of cell subpopulations varied considerably across onset stages among patients, providing a fundamental basis for patient-oriented precision evaluation. Moreover, single-sample gene set enrichment analysis revealed that subsets of CD8+ T cells with significantly different cell adhesion levels could be associated with different patterns of tumor cell dissemination. Therefore, our findings linked the conventional prognostic gene signature to the immune microenvironment and cellular heterogeneity at the single-cell level, thus providing deeper insights into the understanding of HCC tumorigenesis.

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.

Despite the growing number of genome-wide association studies (GWASs), it remains unclear to what extent gene-by-gene and gene-by-environment interactions influence complex traits in humans. The magnitude of genetic interactions in complex traits has been difficult to quantify because GWASs are generally underpowered to detect individual interactions of small effect. Here, we develop a method to test for genetic interactions that aggregates information across all trait-associated loci. Specifically, we test whether SNPs in regions of European ancestry shared between European American and admixed African American individuals have the same causal effect sizes. We hypothesize that in African Americans, the presence of genetic interactions will drive the causal effect sizes of SNPs in regions of European ancestry to be more similar to those of SNPs in regions of African ancestry. We apply our method to two traits: gene expression in 296 African Americans and 482 European Americans in the Multi-Ethnic Study of Atherosclerosis (MESA) and low-density lipoprotein cholesterol (LDL-C) in 74K African Americans and 296K European Americans in the Million Veteran Program (MVP). We find significant evidence for genetic interactions in our analysis of gene expression; for LDL-C, we observe a similar point estimate, although this is not significant, most likely due to lower statistical power. These results suggest that gene-by-gene or gene-by-environment interactions modify the effect sizes of causal variants in human complex traits.

We aimed to characterize the ocular phenotype of patients with ROSAH (retinal dystrophy, optic nerve edema, splenomegaly, anhidrosis, and headache) syndrome and their response to therapy.

The pyrin inflammasome acts as a guard of RhoA GTPases and is central to immune defenses against RhoA-manipulating pathogens. Pyrin activation proceeds in two steps. Yet, the second step is still poorly understood. Using cells constitutively activated for the pyrin step 1, a chemical screen identifies etiocholanolone and pregnanolone, two catabolites of testosterone and progesterone, acting at low concentrations as specific step 2 activators. High concentrations of these metabolites fully and rapidly activate pyrin, in a human specific, B30.2 domain-dependent manner and without inhibiting RhoA. Mutations in MEFV, encoding pyrin, cause two distinct autoinflammatory diseases pyrin-associated autoinflammation with neutrophilic dermatosis (PAAND) and familial Mediterranean fever (FMF). Monocytes from PAAND patients, and to a lower extent from FMF patients, display increased responses to these metabolites. This study identifies an unconventional pyrin activation mechanism, indicates that endogenous steroid catabolites can drive autoinflammation, through the pyrin inflammasome, and explains the "steroid fever" described in the late 1950s upon steroid injection in humans.

Current treatment guidelines recommend ursodeoxycholic acid (UDCA) as the first-line treatment for new-diagnosed primary biliary cholangitis (PBC) patients. However, up to 40% patients are insensitive to UDCA monotherapy, and evaluation of UDCA response at 12 months may result in long period of ineffective treatment. We aimed to develop a new criterion to reliably identify non-response patients much earlier.

Ocular choroidal melanoma (OCM) is the most common ocular primary malignant tumor in adults, and there is an increasing emphasis on its early detection and treatment worldwide. The main obstacle in early detection of OCM is its overlapping clinical features with benign choroidal nevus. Thus, we propose ultrasound localization microscopy (ULM) based on the image deconvolution algorithm to assist the diagnosis of small OCM in early stages. Furthermore, we develop ultrasound (US) plane wave imaging based on three-frame difference algorithm to guide the placement of the probe on the field of view. A high-frequency Verasonics Vantage system and an L22-14v linear array transducer were used to perform experiments on both custom-made modules in vitro and a SD rat with ocular choroidal melanoma in vivo. The results demonstrate that our proposed deconvolution method implement more robust microbubble (MB) localization, reconstruction of microvasculature network in a finer grid and more precise flow velocity estimation. The excellent performance of US plane wave imaging was successfully validated on the flow phantom and in an in vivo OCM model. In the future, the super-resolution ULM, a critical complementary imaging modality, can provide doctors with conclusive suggestions for early diagnosis of OCM, which is significant for the treatment and prognosis of patients.