Search

Abstract The utilization of optimal orthodontic force is crucial to prevent undesirable side effects and ensure efficient tooth movement during orthodontic treatment. However, the sensitivity of existing detection techniques is not sufficient, and the criteria for evaluating optimal force have not been yet established. Here, by employing 3D finite element analysis methodology, we found that the apical distal region (A-D region) of mesial roots is particularly sensitive to orthodontic force in rats. Tartrate-resistant acidic phosphatase (TRAP)-positive osteoclasts began accumulating in the A-D region under the force of 40 grams (g), leading to alveolar bone resorption and tooth movement. When the force reached 80 g, TRAP-positive osteoclasts started appearing on the root surface in the A-D region. Additionally, micro-computed tomography revealed a significant root resorption at 80 g. Notably, the A-D region was identified as a major contributor to whole root resorption. It was determined that 40 g is the minimum effective force for tooth movement with minimal side effects according to the analysis of tooth movement, inclination, and hyalinization. These findings suggest that the A-D region with its changes on the root surface is an important consideration and sensitive indicator when evaluating orthodontic forces for a rat model. Collectively, our investigations into this region would aid in offering valuable implications for preventing and minimizing root resorption during patients’ orthodontic treatment.

Summary: Background: Understanding formation of the human tissue resident memory T cell (TRM) repertoire requires longitudinal access to human non-lymphoid tissues. Methods: By applying flow cytometry and next generation sequencing to serial blood, lymphoid tissue, and gut samples from 16 intestinal transplantation (ITx) patients, we assessed the origin, distribution, and specificity of human TRMs at phenotypic and clonal levels. Findings: Donor age ≥1 year and blood T cell macrochimerism (peak level ≥4%) were associated with delayed establishment of stable recipient TRM repertoires in the transplanted ileum. T cell receptor (TCR) overlap between paired gut and blood repertoires from ITx patients was significantly greater than that in healthy controls, demonstrating increased gut-blood crosstalk after ITx. Crosstalk with the circulating pool remained high for years of follow-up. TCR sequences identifiable in pre-Tx recipient gut but not those in lymphoid tissues alone were more likely to populate post-Tx ileal allografts. Clones detected in both pre-Tx gut and lymphoid tissue had distinct transcriptional profiles from those identifiable in only one tissue. Recipient T cells were distributed widely throughout the gut, including allograft and native colon, which had substantial repertoire overlap. Both alloreactive and microbe-reactive recipient T cells persisted in transplanted ileum, contributing to the TRM repertoire. Interpretation: Our studies reveal human intestinal TRM repertoire establishment from the circulation, preferentially involving lymphoid tissue counterparts of recipient intestinal T cell clones, including TRMs. We have described the temporal and spatial dynamics of this active crosstalk between the circulating pool and the intestinal TRM pool. Funding: This study was funded by the National Institute of Allergy and Infectious Diseases (NIAID) P01 grant AI106697.

Abstract Bassoon (BSN) is a component of a hetero-dimeric presynaptic cytomatrix protein that orchestrates neurotransmitter release with Piccolo (PCLO) from glutamatergic neurons throughout the brain. Heterozygous missense variants in BSN have previously been associated with neurodegenerative disorders in humans. We performed an exome-wide association analysis of ultra-rare variants in about 140,000 unrelated individuals from the UK Biobank to search for new genes associated with obesity. We found that rare heterozygous predicted loss of function (pLoF) variants in BSN are associated with higher BMI with p-value of 3.6e-12 in the UK biobank cohort. Additionally, we identified two individuals (one of whom has a de novo variant) with a heterozygous pLoF variant in a cohort of early onset or extreme obesity and report the clinical histories of these individuals with non-syndromic obesity with no history of neurobehavioral or cognitive disability. The BMI association was replicated in the All of Us whole genome sequencing data. Heterozygous pLoF BSN variants constitute a new etiology for obesity.

Abstract Rare or de novo variants have substantial contribution to human diseases, but the statistical power to identify risk genes by rare variants is generally low due to rarity of genotype data. Previous studies have shown that risk genes usually have high expression in relevant cell types, although for many conditions the identity of these cell types are largely unknown. Recent efforts in single cell atlas in human and model organisms produced large amount of gene expression data. Here we present VBASS, a Bayesian method that integrates single-cell expression and de novo variant (DNV) data to improve power of disease risk gene discovery. VBASS models disease risk prior as a function of expression profiles, approximated by deep neural networks. It learns the weights of neural networks and parameters of Gamma-Poisson likelihood models of DNV counts jointly from expression and genetics data. On simulated data, VBASS shows proper error rate control and better power than state-of-the-art methods. We applied VBASS to published datasets and identified more candidate risk genes with supports from literature or data from independent cohorts. VBASS can be generalized to integrate other types of functional genomics data in statistical genetics analysis.

Overhead transmission lines (OTLs) have always been the major means of power delivery. With the significant increase of transmission voltage and transmission capacity, the dimensions of transmission towers are increasing accordingly, resulting in extensive occupation of land resources. Towers with composite cross arms are a promising solution to this problem, considering the remarkable performance of composite line insulators. In this research, a full-scale alternating current (AC) 500 kV model of a transmission tower with composite cross arms is manufactured and applied under a lightning overvoltage of different polarities. The developing process of streamer–leader discharge is recorded with a high-speed camera, and the major path of the flashover is identified. The flashover voltages are measured and corrected to standard conditions while considering the air humidity and air density, and clearly confirm the polarity effect. The tower’s lightning-withstand level is calculated based on the tower structure and the flashover characteristics. Based on the results obtained from full-scale experiments, the feasibility of composite cross arms is confirmed, and a structural optimization is proposed.

Abstract. Background:. Colorectal cancer (CRC) still ranks the top in morbidity and mortality of cancers worldwide, posing a huge threat and burden to the society. We aimed to determine the age-standardized incidence, mortality, and disability-adjusted life years (DALYs) of CRC and explore potential changes in the temporal trends of the CRC burden in Shanghai during 2002 to 2016. Methods:. The cancer statistics and demographics were obtained from the Cancer Registry and the Statistics Bureau of Pudong New Area, respectively. Data from 2002 to 2016 were included and analyzed retrospectively. DALYs were calculated using DisMod and the age-standardized rates (ASRs) were obtained according to Segi world standard population. Joinpoint regression was used to measure the trends in CRC incidence and to estimate the annual percent change. Results:. The increasing trend of CRC ASR incidence halted after 2014, coinciding with the introduction of the Shanghai CRC screening program. The ASRs of mortality and DALYs increased, at 0.42% (P 74 years. The benefits of the screening program have been partially proven by the ASRs of CRC incidence, providing important insights into better and wider application of screening programs.

hPSCs model of Maturity Onset Diabetes of the Young caused by mutations in the transcription factor HNF1A (HNF1A-MODY), regulates the expression of genes required for the formation of dense-core insulin granules and calcium-dependent insulin secretion, demonstrating a basis to treat HNF1A-MODY patients with sulfonylureas.

Abstract Background DNA sequencing is increasingly incorporated into the routine care of cancer patients, many of whom also carry inherited, moderate/high-penetrance variants associated with other diseases. Yet, the prevalence and consequence of such variants remain unclear. Methods We analyzed the germline genomes of 10,389 adult cancer cases in the TCGA cohort, identifying pathogenic/likely pathogenic variants in autosomal-dominant genes, autosomal-recessive genes, and 59 medically actionable genes curated by the American College of Molecular Genetics (i.e., the ACMG 59 genes). We also analyzed variant- and gene-level expression consequences in carriers. Results The affected genes exhibited varying pan-ancestry and population-specific patterns, and overall, the European population showed the highest frequency of pathogenic/likely pathogenic variants. We further identified genes showing expression consequence supporting variant functionality, including altered gene expression, allelic specific expression, and mis-splicing determined by a massively parallel splicing assay. Conclusions Our results demonstrate that expression-altering variants are found in a substantial fraction of cases and illustrate the yield of genomic risk assessments for a wide range of diseases across diverse populations.

Insufficient bone matrix formation caused by diabetic chronic inflammation can result in bone nonunion, which is perceived as a worldwide epidemic, with a substantial socioeconomic and public health burden. Macrophages in microenvironment orchestrate the inflammation and launch the process of bone remodeling and repair, but aberrant activation of macrophages can drive drastic inflammatory responses during diabetic bone regeneration. In diabetes mellitus, the proliferation of resident macrophages in bone microenvironment is limited, while enhanced myeloid differentiation of hematopoietic stem cells (HSCs) leads to increased and constant monocyte recruitment and thus macrophages shift toward the classic pro-inflammatory phenotype, which leads to the deficiency of bone regeneration. In this review, we systematically summarized the anomalous origin of macrophages under diabetic conditions. Moreover, we evaluated the deficit of pro-regeneration macrophages in the diabetic inflammatory microenvironment. Finally, we further discussed the latest developments on strategies based on targeting macrophages to promote diabetic bone regeneration. Briefly, this review aimed to provide a basis for modulating the biological functions of macrophages to accelerate bone regeneration and rescue diabetic fracture healing in the future.

Neuromyelitis optica spectrum disorders (NMOSD) are rare, debilitating autoimmune diseases of the central nervous system. Many NMOSD patients have antibodies to Aquaporin-4 (AQP4). Prior studies show associations of NMOSD with individual Human Leukocyte Antigen (HLA) alleles and with mutations in the complement pathway and potassium channels. HLA allele associations with NMOSD are inconsistent between populations, suggesting complex relationships between the identified alleles and risk of disease. We used a retrospective case-control approach to identify contributing genetic variants in patients who met the diagnostic criteria for NMOSD and their unaffected family members. Potentially deleterious variants identified in NMOSD patients were compared to members of their families who do not have the disease and to existing databases of human genetic variation. HLA sequences from patients from Belgrade, Serbia, were compared to the frequency of HLA haplotypes in the general population in Belgrade. We analyzed exome sequencing on 40 NMOSD patients and identified rare inherited variants in the complement pathway and potassium channel genes. Haplotype analysis further detected two haplotypes, HLA-A*01, B*08, DRB1*03 and HLA-A*01, B*08, C*07, DRB1*03, DQB1*02, which were more prevalent in NMOSD patients than in unaffected individuals. In silico modeling indicates that HLA molecules within these haplotypes are predicted to bind AQP4 at several sites, potentially contributing to the development of autoimmunity. Our results point to possible autoimmune and neurodegenerative mechanisms that cause NMOSD, and can be used to investigate potential NMOSD drug targets.

Abstract Background Pulmonary arterial hypertension (PAH) is a lethal vasculopathy characterized by pathogenic remodeling of pulmonary arterioles leading to increased pulmonary pressures, right ventricular hypertrophy, and heart failure. PAH can be associated with other diseases (APAH: connective tissue diseases, congenital heart disease, and others) but often the etiology is idiopathic (IPAH). Mutations in bone morphogenetic protein receptor 2 (BMPR2) are the cause of most heritable cases but the vast majority of other cases are genetically undefined. Methods To identify new risk genes, we utilized an international consortium of 4241 PAH cases with exome or genome sequencing data from the National Biological Sample and Data Repository for PAH, Columbia University Irving Medical Center, and the UK NIHR BioResource – Rare Diseases Study. The strength of this combined cohort is a doubling of the number of IPAH cases compared to either national cohort alone. We identified protein-coding variants and performed rare variant association analyses in unrelated participants of European ancestry, including 1647 IPAH cases and 18,819 controls. We also analyzed de novo variants in 124 pediatric trios enriched for IPAH and APAH-CHD. Results Seven genes with rare deleterious variants were associated with IPAH with false discovery rate smaller than 0.1: three known genes (BMPR2, GDF2, and TBX4), two recently identified candidate genes (SOX17, KDR), and two new candidate genes (fibulin 2, FBLN2; platelet-derived growth factor D, PDGFD). The new genes were identified based solely on rare deleterious missense variants, a variant type that could not be adequately assessed in either cohort alone. The candidate genes exhibit expression patterns in lung and heart similar to that of known PAH risk genes, and most variants occur in conserved protein domains. For pediatric PAH, predicted deleterious de novo variants exhibited a significant burden compared to the background mutation rate (2.45×, p = 2.5e−5). At least eight novel pediatric candidate genes carrying de novo variants have plausible roles in lung/heart development. Conclusions Rare variant analysis of a large international consortium identified two new candidate genes—FBLN2 and PDGFD. The new genes have known functions in vasculogenesis and remodeling. Trio analysis predicted that ~ 15% of pediatric IPAH may be explained by de novo variants.

Accurate prediction of variant pathogenicity is essential to understanding genetic risks in disease. Here, the authors present a deep neural network method for prediction of missense variant pathogenicity, MVP, and demonstrate its utility in prioritizing de novo variants contributing to developmental disorders.

Abstract Background Accurate prediction of protein structure is fundamentally important to understand biological function of proteins. Template-based modeling, including protein threading and homology modeling, is a popular method for protein tertiary structure prediction. However, accurate template-query alignment and template selection are still very challenging, especially for the proteins with only distant homologs available. Results We propose a new template-based modelling method called ThreaderAI to improve protein tertiary structure prediction. ThreaderAI formulates the task of aligning query sequence with template as the classical pixel classification problem in computer vision and naturally applies deep residual neural network in prediction. ThreaderAI first employs deep learning to predict residue-residue aligning probability matrix by integrating sequence profile, predicted sequential structural features, and predicted residue-residue contacts, and then builds template-query alignment by applying a dynamic programming algorithm on the probability matrix. We evaluated our methods both in generating accurate template-query alignment and protein threading. Experimental results show that ThreaderAI outperforms currently popular template-based modelling methods HHpred, CNFpred, and the latest contact-assisted method CEthreader, especially on the proteins that do not have close homologs with known structures. In particular, in terms of alignment accuracy measured with TM-score, ThreaderAI outperforms HHpred, CNFpred, and CEthreader by 56, 13, and 11%, respectively, on template-query pairs at the similarity of fold level from SCOPe data. And on CASP13’s TBM-hard data, ThreaderAI outperforms HHpred, CNFpred, and CEthreader by 16, 9 and 8% in terms of TM-score, respectively. Conclusions These results demonstrate that with the help of deep learning, ThreaderAI can significantly improve the accuracy of template-based structure prediction, especially for distant-homology proteins.

Abstract Background The contribution of somatic mosaicism, or genetic mutations arising after oocyte fertilization, to congenital heart disease (CHD) is not well understood. Further, the relationship between mosaicism in blood and cardiovascular tissue has not been determined. Methods We developed a new computational method, EM-mosaic (Expectation-Maximization-based detection of mosaicism), to analyze mosaicism in exome sequences derived primarily from blood DNA of 2530 CHD proband-parent trios. To optimize this method, we measured mosaic detection power as a function of sequencing depth. In parallel, we analyzed our cohort using MosaicHunter, a Bayesian genotyping algorithm-based mosaic detection tool, and compared the two methods. The accuracy of these mosaic variant detection algorithms was assessed using an independent resequencing method. We then applied both methods to detect mosaicism in cardiac tissue-derived exome sequences of 66 participants for which matched blood and heart tissue was available. Results EM-mosaic detected 326 mosaic mutations in blood and/or cardiac tissue DNA. Of the 309 detected in blood DNA, 85/97 (88%) tested were independently confirmed, while 7/17 (41%) candidates of 17 detected in cardiac tissue were confirmed. MosaicHunter detected an additional 64 mosaics, of which 23/46 (50%) among 58 candidates from blood and 4/6 (67%) of 6 candidates from cardiac tissue confirmed. Twenty-five mosaic variants altered CHD-risk genes, affecting 1% of our cohort. Of these 25, 22/22 candidates tested were confirmed. Variants predicted as damaging had higher variant allele fraction than benign variants, suggesting a role in CHD. The estimated true frequency of mosaic variants above 10% mosaicism was 0.14/person in blood and 0.21/person in cardiac tissue. Analysis of 66 individuals with matched cardiac tissue available revealed both tissue-specific and shared mosaicism, with shared mosaics generally having higher allele fraction. Conclusions We estimate that ~ 1% of CHD probands have a mosaic variant detectable in blood that could contribute to cardiac malformations, particularly those damaging variants with relatively higher allele fraction. Although blood is a readily available DNA source, cardiac tissues analyzed contributed ~ 5% of somatic mosaic variants identified, indicating the value of tissue mosaicism analyses.

Abstract Background Group 1 pulmonary arterial hypertension (PAH) is a rare disease with high mortality despite recent therapeutic advances. Pathogenic remodeling of pulmonary arterioles leads to increased pulmonary pressures, right ventricular hypertrophy, and heart failure. Mutations in bone morphogenetic protein receptor type 2 and other risk genes predispose to disease, but the vast majority of non-familial cases remain genetically undefined. Methods To identify new risk genes, we performed exome sequencing in a large cohort from the National Biological Sample and Data Repository for PAH (PAH Biobank, n = 2572). We then carried out rare deleterious variant identification followed by case-control gene-based association analyses. To control for population structure, only unrelated European cases (n = 1832) and controls (n = 12,771) were used in association tests. Empirical p values were determined by permutation analyses, and the threshold for significance defined by Bonferroni’s correction for multiple testing. Results Tissue kallikrein 1 (KLK1) and gamma glutamyl carboxylase (GGCX) were identified as new candidate risk genes for idiopathic PAH (IPAH) with genome-wide significance. We note that variant carriers had later mean age of onset and relatively moderate disease phenotypes compared to bone morphogenetic receptor type 2 variant carriers. We also confirmed the genome-wide association of recently reported growth differentiation factor (GDF2) with IPAH and further implicate T-box 4 (TBX4) with child-onset PAH. Conclusions We report robust association of novel genes KLK1 and GGCX with IPAH, accounting for ~ 0.4% and 0.9% of PAH Biobank cases, respectively. Both genes play important roles in vascular hemodynamics and inflammation but have not been implicated in PAH previously. These data suggest new genes, pathogenic mechanisms, and therapeutic targets for this lethal vasculopathy.

Large whole-exome sequencing studies have suggested that the genetic architecture of syndromic congenital heart disease (CHD) is different from sporadic forms. Here, Watkins et al. estimate the relative contribution of damaging recessive and de novo genotypes to CHD in 2391 trios and find them to be associated with different gene functions.

Energy balance is controlled by interconnected brain regions in the hypothalamus, brainstem, cortex, and limbic system. Gene expression signatures of these regions can help elucidate the pathophysiology underlying obesity. RNA sequencing was conducted on P56 C57BL/6NTac male mice and E14.5 C57BL/6NTac embryo punch biopsies in 16 obesity-relevant brain regions. The expression of 190 known obesity-associated genes (monogenic, rare, and low-frequency coding variants; GWAS; syndromic) was analyzed in each anatomical region. Genes associated with these genetic categories of obesity had localized expression patterns across brain regions. Known monogenic obesity causal genes were highly enriched in the arcuate nucleus of the hypothalamus and developing hypothalamus. The obesity-associated genes clustered into distinct “modules” of similar expression profile, and these were distinct from expression modules formed by similar analysis with genes known to be associated with other disease phenotypes (type 1 and type 2 diabetes, autism, breast cancer) in the same energy balance–relevant brain regions.

Summary: KIF1A-associated neurological disorder (KAND) encompasses a group of rare neurodegenerative conditions caused by variants in KIF1A, a gene that encodes an anterograde neuronal microtubule (MT) motor protein. Here we characterize the natural history of KAND in 117 individuals using a combination of caregiver or self-reported medical history, a standardized measure of adaptive behavior, clinical records, and neuropathology. We developed a heuristic severity score using a weighted sum of common symptoms to assess disease severity. Focusing on 100 individuals, we compared the average clinical severity score for each variant with in silico predictions of deleteriousness and location in the protein. We found increased severity is strongly associated with variants occurring in protein regions involved with ATP and MT binding: the P loop, switch I, and switch II. For a subset of variants, we generated recombinant proteins, which we used to assess transport in vivo by assessing neurite tip accumulation and to assess MT binding, motor velocity, and processivity using total internal reflection fluorescence microscopy. We find all modeled variants result in defects in protein transport, and we describe three classes of protein dysfunction: reduced MT binding, reduced velocity and processivity, and increased non-motile rigor MT binding. The rigor phenotype is consistently associated with the most severe clinical phenotype, while reduced MT binding is associated with milder clinical phenotypes. Our findings suggest the clinical phenotypic heterogeneity in KAND likely reflects and parallels diverse molecular phenotypes. We propose a different way to describe KAND subtypes to better capture the breadth of disease severity.

Abstract Background Classic dendritic cells (cDCs) play a central role in the immune system by processing and presenting antigens to activate T cells, and consist of two major subsets: CD141+ cDC (cDC1) and CD1c+ cDC (cDC2). A population of migratory precursor cells, the pre-cDCs, is the immediate precursors to both cDC subsets. Previous studies showed that there were two pre-committed pre-cDC subpopulations. However, the key molecular drivers of pre-commitment in human pre-cDCs were not investigated. Results To identify the key molecular drivers for pre-commitment in human pre-cDCs, we performed single cell RNA sequencing (RNA-Seq) of two cDC subsets and pre-cDCs, and bulk RNA-Seq of pre-cDCs and cDCs from human peripheral blood. We found that pre-DC subpopulations cannot be separated by either variable genes within pre-cDCs or differentially expressed genes between cDC1 and cDC2. In contrast, they were separated by 16 transcription factors that are themselves differentially expressed or have regulated targets enriched in the differentially expressed genes between bulk cDC1 and cDC2, with one subpopulation close to cDC1 and the other close to cDC2. More importantly, these two pre-cDC sub-populations are correlated with ratio of IRF8 to IRF4 expression level more than their individual expression level. We also verified these findings using three recently published datasets. Conclusions In this study, we demonstrate that single cell transcriptome profiling can reveal pre-cDCs differentiation map, and our results suggest the concept that combinatorial dose of transcription factors determines cell differentiation fate.

Abstract Background Recent large-scale cancer sequencing studies have discovered many novel cancer driver genes (CDGs) in human cancers. Some studies also suggest that CDG mutations contribute to cancer-associated epigenomic and transcriptomic alterations across many cancer types. Here we aim to improve our understanding of the connections between CDG mutations and altered cancer cell epigenomes and transcriptomes on pan-cancer level and how these connections contribute to the known association between epigenome and transcriptome. Method Using multi-omics data including somatic mutation, DNA methylation, and gene expression data of 20 cancer types from The Cancer Genome Atlas (TCGA) project, we conducted a pan-cancer analysis to identify CDGs, when mutated, have strong associations with genome-wide methylation or expression changes across cancer types, which we refer as methylation driver genes (MDGs) or expression driver genes (EDGs), respectively. Results We identified 32 MDGs, among which, eight are known chromatin modification or remodeling genes. Many of the remaining 24 MDGs are connected to chromatin regulators through either regulating their transcription or physically interacting with them as potential co-factors. We identified 29 EDGs, 26 of which are also MDGs. Further investigation on target genes’ promoters methylation and expression alteration patterns of these 26 overlapping driver genes shows that hyper-methylation of target genes’ promoters are significantly associated with down-regulation of the same target genes and hypo-methylation of target genes’ promoters are significantly associated with up-regulation of the same target genes. Conclusion This finding suggests a pivotal role for genetically driven changes in chromatin remodeling in shaping DNA methylation and gene expression patterns during tumor development.

Summary: The diaphragm is critical for respiration and separation of the thoracic and abdominal cavities, and defects in diaphragm development are the cause of congenital diaphragmatic hernias (CDH), a common and often lethal birth defect. The genetic etiology of CDH is complex. Single-nucleotide variants (SNVs), insertions/deletions (indels), and structural variants (SVs) in more than 150 genes have been associated with CDH, although few genes are recurrently mutated in multiple individuals and mutated genes are incompletely penetrant. This suggests that multiple genetic variants in combination, other not-yet-investigated classes of variants, and/or nongenetic factors contribute to CDH etiology. However, no studies have comprehensively investigated in affected individuals the contribution of all possible classes of variants throughout the genome to CDH etiology. In our study, we used a unique cohort of four individuals with isolated CDH with samples from blood, skin, and diaphragm connective tissue and parental blood and deep whole-genome sequencing to assess germline and somatic de novo and inherited SNVs, indels, and SVs. In each individual we found a different mutational landscape that included germline de novo and inherited SNVs and indels in multiple genes. We also found in two individuals a 343 bp deletion interrupting an annotated enhancer of the CDH-associated gene GATA4, and we hypothesize that this common SV (found in 1%–2% of the population) acts as a sensitizing allele for CDH. Overall, our comprehensive reconstruction of the genetic architecture of four CDH individuals demonstrates that the etiology of CDH is heterogeneous and multifactorial.

Abstract Background Pulmonary arterial hypertension (PAH) is a rare disease characterized by distinctive changes in pulmonary arterioles that lead to progressive pulmonary arterial pressures, right-sided heart failure, and a high mortality rate. Up to 30% of adult and 75% of pediatric PAH cases are associated with congenital heart disease (PAH-CHD), and the underlying etiology is largely unknown. There are no known major risk genes for PAH-CHD. Methods To identify novel genetic causes of PAH-CHD, we performed whole exome sequencing in 256 PAH-CHD patients. We performed a case-control gene-based association test of rare deleterious variants using 7509 gnomAD whole genome sequencing population controls. We then screened a separate cohort of 413 idiopathic and familial PAH patients without CHD for rare deleterious variants in the top association gene. Results We identified SOX17 as a novel candidate risk gene (p = 5.5e−7). SOX17 is highly constrained and encodes a transcription factor involved in Wnt/β-catenin and Notch signaling during development. We estimate that rare deleterious variants contribute to approximately 3.2% of PAH-CHD cases. The coding variants identified include likely gene-disrupting (LGD) and deleterious missense, with most of the missense variants occurring in a highly conserved HMG-box protein domain. We further observed an enrichment of rare deleterious variants in putative targets of SOX17, many of which are highly expressed in developing heart and pulmonary vasculature. In the cohort of PAH without CHD, rare deleterious variants of SOX17 were observed in 0.7% of cases. Conclusions These data strongly implicate SOX17 as a new risk gene contributing to PAH-CHD as well as idiopathic/familial PAH. Replication in other PAH cohorts and further characterization of the clinical phenotype will be important to confirm the precise role of SOX17 and better estimate the contribution of genes regulated by SOX17.

Predicting haploinsufficient genes helps to understand the genetic risk underlying developmental disorders. Here, the authors develop a Random Forest-based method that uses epigenomic data to predict haploinsufficiency, Episcore, which is complementary to methods based on mutation intolerance scores.

There is a great need of developing highly sensitive mass spectrometry-based proteomics analysis for small cell populations. Here, the authors establish a robotically controlled chip-based nanodroplet processing platform and demonstrate its ability to profile the proteome from 10–100 mammalian cells.

Tissue-resident memory T cells (TRMs) in mice mediate optimal protective immunity to infection and vaccination, while in humans, the existence and properties of TRMs remain unclear. Here, we use a unique human tissue resource to determine whether human tissue memory T cells constitute a distinct subset in diverse mucosal and lymphoid tissues. We identify a core transcriptional profile within the CD69+ subset of memory CD4+ and CD8+ T cells in lung and spleen that is distinct from that of CD69− TEM cells in tissues and circulation and defines human TRMs based on homology to the transcriptional profile of mouse CD8+ TRMs. Human TRMs in diverse sites exhibit increased expression of adhesion and inhibitory molecules, produce both pro-inflammatory and regulatory cytokines, and have reduced turnover compared with circulating TEM, suggesting unique adaptations for in situ immunity. Together, our results provide a unifying signature for human TRM and a blueprint for designing tissue-targeted immunotherapies.

Congenital diaphragmatic hernias (CDHs) and structural anomalies of the diaphragm are a common class of congenital birth defects that are associated with significant morbidity and mortality due to associated pulmonary hypoplasia, pulmonary hypertension and heart failure. In ∼30% of CDH patients, genomic analyses have identified a range of genetic defects, including chromosomal anomalies, copy number variants and sequence variants. The affected genes identified in CDH patients include transcription factors, such as GATA4, ZFPM2, NR2F2 and WT1, and signaling pathway components, including members of the retinoic acid pathway. Mutations in these genes affect diaphragm development and can have pleiotropic effects on pulmonary and cardiac development. New therapies, including fetal endoscopic tracheal occlusion and prenatal transplacental fetal treatments, aim to normalize lung development and pulmonary vascular tone to prevent and treat lung hypoplasia and pulmonary hypertension, respectively. Studies of the association between particular genetic mutations and clinical outcomes should allow us to better understand the origin of this birth defect and to improve our ability to predict and identify patients most likely to benefit from specialized treatment strategies.

Congenital heart disease (CHD) is a disorder that occurs in ∼1% of live births. Here the authors describe a genome-wide allele-specific expression analyses in CHD patients, identifying five new genes involved in CHD and showing that paternally-expressed imprinted genes are monoallelic, while maternally-expressed imprinted genes are biallelic.

Genome interpretation and analysis of allelic activity requires appropriate haplotype phasing. Here the authors present phASER, a fast and accurate method for variant phrasing from RNA-seq and genome sequencing data.

Congenital diaphragmatic hernia (CDH) is a severe birth defect that is often accompanied by other congenital anomalies. Previous exome sequencing studies for CDH have supported a role of de novo damaging variants but did not identify any recurrently mutated genes. To investigate further the genetics of CDH, we analyzed de novo coding variants in 362 proband-parent trios including 271 new trios reported in this study. We identified four unrelated individuals with damaging de novo variants in MYRF (P = 5.3x10(-8)), including one likely gene-disrupting (LGD) and three deleterious missense (D-mis) variants. Eight additional individuals with de novo LGD or missense variants were identified from our other genetic studies or from the literature. Common phenotypes of MYRF de novo variant carriers include CDH, congenital heart disease and genitourinary abnormalities, suggesting that it represents a novel syndrome. MYRF is a membrane associated transcriptional factor highly expressed in developing diaphragm and is depleted of LGD variants in the general population. All de novo missense variants aggregated in two functional protein domains. Analyzing the transcriptome of patient-derived diaphragm fibroblast cells suggest that disease associated variants abolish the transcription factor activity. Furthermore, we showed that the remaining genes with damaging variants in CDH significantly overlap with genes implicated in other developmental disorders. Gene expression patterns and patient phenotypes support pleiotropic effects of damaging variants in these genes on CDH and other developmental disorders. Finally, functional enrichment analysis implicates the disruption of regulation of gene expression, kinase activities, intra-cellular signaling, and cytoskeleton organization as pathogenic mechanisms in CDH.

Short-read sequencing has inherent limitations in the characterisation of long repeat elements. Shi and Guo et al.combine single-molecule real-time sequencing and IrysChip to construct a Chinese reference genome that fills many gaps in the reference genome, and identify novel spliced genes.

The ability to conduct investigation of cellular transcription, signaling, and function at the single-cell level has opened opportunities to examine heterogeneous populations at unprecedented resolutions. Although methods have been developed to evaluate high-dimensional transcriptomic and proteomic data (relating to cellular mRNA and protein), there has not been a method to evaluate corresponding high-dimensional functionomic data (relating to cellular functions) from single cells. Here, we present a protocol to quantitatively measure the differentiation potentials of single human hematopoietic stem and progenitor cells, and then cluster the cells according to these measurements. High dimensional functionomic analysis of cell potential allows cell function to be linked to molecular mechanisms within the same progenitor population.

Stevens-Johnson syndrome (SJS) and Toxic Epidermal Necrolysis (TEN) represent rare but serious adverse drug reactions (ADRs). Both are characterized by distinctive blistering lesions and significant mortality rates. While there is evidence for strong drug-specific genetic predisposition related to HLA alleles, recent genome wide association studies (GWAS) on European and Asian populations have failed to identify genetic susceptibility alleles that are common across multiple drugs. We hypothesize that this is a consequence of the low to moderate effect size of individual genetic risk factors. To test this hypothesis we developed Pointer, a new algorithm that assesses the aggregate effect of multiple low risk variants on a pathway using a gene set enrichment approach. A key advantage of our method is the capability to associate SNPs with genes by exploiting physical proximity as well as by using expression quantitative trait loci (eQTLs) that capture information about both cis- and trans-acting regulatory effects. We control for known bias-inducing aspects of enrichment based analyses, such as: 1) gene length, 2) gene set size, 3) presence of biologically related genes within the same linkage disequilibrium (LD) region, and, 4) genes shared among multiple gene sets. We applied this approach to publicly available SJS/TEN genome-wide genotype data and identified the ABC transporter and Proteasome pathways as potentially implicated in the genetic susceptibility of non-drug-specific SJS/TEN. We demonstrated that the innovative SNP-to-gene mapping phase of the method was essential in detecting the significant enrichment for those pathways. Analysis of an independent gene expression dataset provides supportive functional evidence for the involvement of Proteasome pathways in SJS/TEN cutaneous lesions. These results suggest that Pointer provides a useful framework for the integrative analysis of pharmacogenetic GWAS data, by increasing the power to detect aggregate effects of multiple low risk variants. The software is available for download at https://sourceforge.net/projects/pointergsa/.

In human immune system, V(D)J recombination produces an enormously large repertoire of immunoglobulins (Ig) so that they can tackle different antigens from bacteria, viruses and tumor cells. Several studies have demonstrated the utility of next-generation sequencers such as Roche 454 and Illumina Genome Analyzer to characterize the repertoire of immunoglobulins. However, these techniques typically require separation of B cell population from whole blood and require a few weeks for running the sequencers, so it may not be practical to implement them in clinical settings. Recently, the Ion Torrent personal genome sequencer has emerged as a tabletop personal genome sequencer that can be operated in a time-efficient and cost-effective manner. In this study, we explored the technical feasibility to use multiplex PCR for amplifying V(D)J recombination for IgH, directly from whole blood, then sequence the amplicons by the Ion Torrent sequencer. The whole process including data generation and analysis can be completed in one day. We tested the method in a pilot study on patients with benign, atypical and malignant meningiomas. Despite the noisy data, we were able to compare the samples by their usage frequencies of the V segment, as well as their somatic hypermutation rates. In summary, our study suggested that it is technically feasible to perform clinical monitoring of V(D)J recombination within a day by personal genome sequencers.

Marked prolongation of the QT interval on the electrocardiogram associated with the polymorphic ventricular tachycardia Torsades de Pointes is a serious adverse event during treatment with antiarrhythmic drugs and other culprit medications, and is a common cause for drug relabeling and withdrawal. Although clinical risk factors have been identified, the syndrome remains unpredictable in an individual patient. Here we used genome-wide association analysis to search for common predisposing genetic variants. Cases of drug-induced Torsades de Pointes (diTdP), treatment tolerant controls, and general population controls were ascertained across multiple sites using common definitions, and genotyped on the Illumina 610k or 1M-Duo BeadChips. Principal Components Analysis was used to select 216 Northwestern European diTdP cases and 771 ancestry-matched controls, including treatment-tolerant and general population subjects. With these sample sizes, there is 80% power to detect a variant at genome-wide significance with minor allele frequency of 10% and conferring an odds ratio of ≥2.7. Tests of association were carried out for each single nucleotide polymorphism (SNP) by logistic regression adjusting for gender and population structure. No SNP reached genome wide-significance; the variant with the lowest P value was rs2276314, a non-synonymous coding variant in C18orf21 (p = 3×10(-7), odds ratio = 2, 95% confidence intervals: 1.5-2.6). The haplotype formed by rs2276314 and a second SNP, rs767531, was significantly more frequent in controls than cases (p = 3×10(-9)). Expanding the number of controls and a gene-based analysis did not yield significant associations. This study argues that common genomic variants do not contribute importantly to risk for drug-induced Torsades de Pointes across multiple drugs.

Cancer invasion and metastasis are closely associated with activities within the degradome; however, little is known about whether these activities can be detected in the blood of cancer patients.The peptidome-degradome profiles of pooled blood plasma sampled from 15 breast cancer patients (BCP) and age, race, and menopausal status matched control healthy persons (HP) were globally characterized using advanced comprehensive separations combined with tandem Fourier transform mass spectrometry and new data analysis approaches that facilitated top-down peptidomic analysis. The BCP pool displayed 71 degradome protein substrates that encompassed 839 distinct peptidome peptides. In contrast, the HP 50 degradome substrates found encompassed 425 peptides. We find that the ratios of the peptidome peptide relative abundances can vary as much as >4000 fold between BCP and HP. The experimental results also show differential degradation of substrates in the BCP sample in their functional domains, including the proteolytic and inhibitory sites of the plasmin-antiplasmin and thrombin-antithrombin systems, the main chains of the extracellular matrix protection proteins, the excessive degradation of innate immune system key convertases and membrane attack complex components, as well as several other cancer suppressor proteins.Degradomics-peptidomics profiling of blood plasma is highly sensitive to changes not evidenced by conventional bottom-up proteomics and potentially provides unique signatures of possible diagnostic utility.

Whole-genome sequencing represents a promising approach to pinpoint chemically induced mutations in genetic model organisms, thereby short-cutting time-consuming genetic mapping efforts.We compare here the ability of two leading high-throughput platforms for paired-end deep sequencing, SOLiD (ABI) and Genome Analyzer (Illumina; "Solexa"), to achieve the goal of mutant detection. As a test case we used a mutant C. elegans strain that harbors a mutation in the lsy-12 locus which we compare to the reference wild-type genome sequence. We analyzed the accuracy, sensitivity, and depth-coverage characteristics of the two platforms. Both platforms were able to identify the mutation that causes the phenotype of the mutant C. elegans strain, lsy-12. Based on a 4 MB genomic region in which individual variants were validated by Sanger sequencing, we observe tradeoffs between rates of false positives and false negatives when using both platforms under similar coverage and mapping criteria.In conclusion, whole-genome sequencing conducted by either platform is a viable approach for the identification of single-nucleotide variations in the C. elegans genome.