Shared activity patterns arising at genetic susceptibility loci reveal underlying genomic and cellular architecture of human disease

Temporal dynamics of genetic architecture governing leaf development in Populus

Leaf development is a multifaceted and dynamic process orchestrated by a myriad of genes to shape the proper size and morphology. The dynamic genetic network underlying leaf development remains largely unknown. Utilizing a synergistic genetic approach encompassing dynamic genome-wide association study (GWAS), time-ordered gene co-expression network (TO-GCN) analyses and gene manipulation, we explored the temporal genetic architecture and regulatory network governing leaf development in Populus. We identified 42 time-specific and 18 consecutive genes that displayed different patterns of expression at various time points. We then constructed eight TO-GCNs that covered the cell proliferation, transition, and cell expansion stages of leaf development. Integrating GWAS and TO-GCN, we postulated the functions of 27 causative genes for GWAS and identified PtoGRF9 as a key player in leaf development. Genetic manipulation via overexpression and suppression of PtoGRF9 revealed its primary influence on leaf development by modulating cell proliferation. Furthermore, we elucidated that PtoGRF9 governs leaf development by activating PtoHB21 during the cell proliferation stage and attenuating PtoLD during the transition stage. Our study provides insights into the dynamic genetic underpinnings of leaf development and understanding the regulatory mechanism of PtoGRF9 in this dynamic process.

The relationship between extreme inter-individual variation in macrophage gene expression and genetic susceptibility to inflammatory bowel disease

The differentiation of resident intestinal macrophages from blood monocytes depends upon signals from the macrophage colony-stimulating factor receptor (CSF1R). Analysis of genome-wide association studies (GWAS) indicates that dysregulation of macrophage differentiation and response to microorganisms contributes to susceptibility to chronic inflammatory bowel disease (IBD). Here, we analyzed transcriptomic variation in monocyte-derived macrophages (MDM) from affected and unaffected sib pairs/trios from 22 IBD families and 6 healthy controls. Transcriptional network analysis of the data revealed no overall or inter-sib distinction between affected and unaffected individuals in basal gene expression or the temporal response to lipopolysaccharide (LPS). However, the basal or LPS-inducible expression of individual genes varied independently by as much as 100-fold between subjects. Extreme independent variation in the expression of pairs of HLA-associated transcripts (HLA-B/C, HLA-A/F and HLA-DRB1/DRB5) in macrophages was associated with HLA genotype. Correlation analysis indicated the downstream impacts of variation in the immediate early response to LPS. For example, variation in early expression of IL1B was significantly associated with local SNV genotype and with subsequent peak expression of target genes including IL23A, CXCL1, CXCL3, CXCL8 and NLRP3. Similarly, variation in early IFNB1 expression was correlated with subsequent expression of IFN target genes. Our results support the view that gene-specific dysregulation in macrophage adaptation to the intestinal milieu is associated with genetic susceptibility to IBD.

Mapping the epigenomic landscape of human monocytes following innate immune activation reveals context-specific mechanisms driving endotoxin tolerance

BACKGROUND

Monocytes are key mediators of innate immunity to infection, undergoing profound and dynamic changes in epigenetic state and immune function which are broadly protective but may be dysregulated in disease. Here, we aimed to advance understanding of epigenetic regulation following innate immune activation, acutely and in endotoxin tolerant states.

METHODS

We exposed human primary monocytes from healthy donors (n = 6) to interferon-γ or differing combinations of endotoxin (lipopolysaccharide), including acute response (2 h) and two models of endotoxin tolerance: repeated stimulations (6 + 6 h) and prolonged exposure to endotoxin (24 h). Another subset of monocytes was left untreated (naïve). We identified context-specific regulatory elements based on epigenetic signatures for chromatin accessibility (ATAC-seq) and regulatory non-coding RNAs from total RNA sequencing.

RESULTS

We present an atlas of differential gene expression for endotoxin and interferon response, identifying widespread context specific changes. Across assayed states, only 24-29% of genes showing differential exon usage are also differential at the gene level. Overall, 19.9% (6,884 of 34,616) of repeatedly observed ATAC peaks were differential in at least one condition, the majority upregulated on stimulation and located in distal regions (64.1% vs 45.9% of non-differential peaks) within which sequences were less conserved than non-differential peaks. We identified enhancer-derived RNA signatures specific to different monocyte states that correlated with chromatin accessibility changes. The endotoxin tolerance models showed distinct chromatin accessibility and transcriptomic signatures, with integrated analysis identifying genes and pathways involved in the inflammatory response, detoxification, metabolism and wound healing. We leveraged eQTL mapping for the same monocyte activation states to link potential enhancers with specific genes, identifying 1,946 unique differential ATAC peaks with 1,340 expression associated genes. We further use this to inform understanding of reported GWAS, for example involving FCHO1 and coronary artery disease.

CONCLUSION

This study reports context-specific regulatory elements based on transcriptomic profiling and epigenetic signatures for enhancer-derived RNAs and chromatin accessibility in immune tolerant monocyte states, and demonstrates the informativeness of linking such elements and eQTL to inform future mechanistic studies aimed at defining therapeutic targets of immunosuppression and diseases.

The contribution of evolutionarily volatile promoters to molecular phenotypes and human trait variation

BACKGROUND

Promoters are sites of transcription initiation that harbour a high concentration of phenotype-associated genetic variation. The evolutionary gain and loss of promoters between species (collectively, termed turnover) is pervasive across mammalian genomes and may play a prominent role in driving human phenotypic diversity.

RESULTS

We classified human promoters by their evolutionary history during the divergence of mouse and human lineages from a common ancestor. This defined conserved, human-inserted and mouse-deleted promoters, and a class of functional-turnover promoters that align between species but are only active in humans. We show that promoters of all evolutionary categories are hotspots for substitution and often, insertion mutations. Loci with a history of insertion and deletion continue that mode of evolution within contemporary humans. The presence of an evolutionary volatile promoter within a gene is associated with increased expression variance between individuals, but only in the case of human-inserted and mouse-deleted promoters does that correspond to an enrichment of promoter-proximal genetic effects. Despite the enrichment of these molecular quantitative trait loci (QTL) at evolutionarily volatile promoters, this does not translate into a corresponding enrichment of phenotypic traits mapping to these loci.

CONCLUSIONS

Promoter turnover is pervasive in the human genome, and these promoters are rich in molecularly quantifiable but phenotypically inconsequential variation in gene expression. However, since evolutionarily volatile promoters show evidence of selection, coupled with high mutation rates and enrichment of QTLs, this implicates them as a source of evolutionary innovation and phenotypic variation, albeit with a high background of selectively neutral expression variation.

Macrophage control of Crohn's disease

The intestinal tract is the body's largest mucosal surface and permanently exposed to microbial and environmental signals. Maintaining a healthy intestine requires the presence of sentinel grounds keeper cells, capable of controlling immunity and tissue homeostasis through specialized functions. Intestinal macrophages are such cells and important players in steady-state functions and during acute and chronic inflammation. Crohn's disease, a chronic inflammatory condition of the intestinal tract is proposed to be the consequence of an altered immune system through microbial and environmental stimulation. This hypothesis suggests an involvement of macrophages in the regulation of this pathology. Within this chapter, we will discuss intestinal macrophage development and highlight data suggesting their implication in chronic intestinal pathologies like Crohn's disease.

Genetic Determinants of Plasma Low-Density Lipoprotein Cholesterol Levels: Monogenicity, Polygenicity, and "Missing" Heritability

Changes in plasma low-density lipoprotein cholesterol (LDL-c) levels relate to a high risk of developing some common and complex diseases. LDL-c, as a quantitative trait, is multifactorial and depends on both genetic and environmental factors. In the pregenomic age, targeted genes were used to detect genetic factors in both hyper- and hypolipidemias, but this approach only explained extreme cases in the population distribution. Subsequently, the genetic basis of the less severe and most common dyslipidemias remained unknown. In the genomic age, performing whole-exome sequencing in families with extreme plasma LDL-c values identified some new candidate genes, but it is unlikely that such genes can explain the majority of inexplicable cases. Genome-wide association studies (GWASs) have identified several single-nucleotide variants (SNVs) associated with plasma LDL-c, introducing the idea of a polygenic origin. Polygenic risk scores (PRSs), including LDL-c-raising alleles, were developed to measure the contribution of the accumulation of small-effect variants to plasma LDL-c. This paper discusses other possibilities for unexplained dyslipidemias associated with LDL-c, such as mosaicism, maternal effect, and induced epigenetic changes. Future studies should consider gene-gene and gene-environment interactions and the development of integrated information about disease-driving networks, including phenotypes, genotypes, transcription, proteins, metabolites, and epigenetics.

Comprehensive Characterization of Transcriptional Activity during Influenza A Virus Infection Reveals Biases in Cap-Snatching of Host RNA Sequences

Infection with influenza A virus (IAV) infection is responsible for an estimated 500,000 deaths and up to 5 million cases of severe respiratory illness each year. In this study, we looked at human primary immune cells (macrophages) infected with IAV. Our method allows us to look at both the host and the virus in parallel. We used these data to explore a process known as “cap-snatching,” where IAV snatches a short nucleotide sequence from capped host RNA. This process was believed to be random. We demonstrate biased snatching of numerous host RNAs, including those associated with snRNA transcription, and avoidance of host transcripts encoding host ribosomal proteins, which are required by IAV for replication. We then describe the transcriptional landscape of the host response to IAV, observing new features, including a failure of IAV-treated MDMs to induce feedback inhibitors of inflammation, seen in response to other treatments.

Macrophages in the lung detect and respond to influenza A virus (IAV), determining the nature of the immune response. Using terminal-depth cap analysis of gene expression (CAGE), we quantified transcriptional activity of both host and pathogen over a 24-h time course of IAV infection in primary human monocyte-derived macrophages (MDMs). This method allowed us to observe heterogenous host sequences incorporated into IAV mRNA, “snatched” 5′ RNA caps, and corresponding RNA sequences from host RNAs. In order to determine whether cap-snatching is random or exhibits a bias, we systematically compared host sequences incorporated into viral mRNA (“snatched”) against a complete survey of all background host RNA in the same cells, at the same time. Using a computational strategy designed to eliminate sources of bias due to read length, sequencing depth, and multimapping, we were able to quantify overrepresentation of host RNA features among the sequences that were snatched by IAV. We demonstrate biased snatching of numerous host RNAs, particularly small nuclear RNAs (snRNAs), and avoidance of host transcripts encoding host ribosomal proteins, which are required by IAV for replication. We then used a systems approach to describe the transcriptional landscape of the host response to IAV, observing many new features, including a failure of IAV-treated MDMs to induce feedback inhibitors of inflammation, seen in response to other treatments.

IMPORTANCE Infection with influenza A virus (IAV) infection is responsible for an estimated 500,000 deaths and up to 5 million cases of severe respiratory illness each year. In this study, we looked at human primary immune cells (macrophages) infected with IAV. Our method allows us to look at both the host and the virus in parallel. We used these data to explore a process known as “cap-snatching,” where IAV snatches a short nucleotide sequence from capped host RNA. This process was believed to be random. We demonstrate biased snatching of numerous host RNAs, including those associated with snRNA transcription, and avoidance of host transcripts encoding host ribosomal proteins, which are required by IAV for replication. We then describe the transcriptional landscape of the host response to IAV, observing new features, including a failure of IAV-treated MDMs to induce feedback inhibitors of inflammation, seen in response to other treatments.

Host susceptibility to severe influenza A virus infection

Most people exposed to a new flu virus do not notice any symptoms. A small minority develops critical illness. Some of this extremely broad variation in susceptibility is explained by the size of the initial inoculum or the influenza exposure history of the individual; some is explained by generic host factors, such as frailty, that decrease resilience following any systemic insult. Some demographic factors (pregnancy, obesity, and advanced age) appear to confer a more specific susceptibility to severe illness following infection with influenza viruses. As with other infectious diseases, a substantial component of susceptibility is determined by host genetics. Several genetic susceptibility variants have now been reported with varying levels of evidence. Susceptible hosts may have impaired intracellular controls of viral replication (e.g. IFITM3, TMPRS22 variants), defective interferon responses (e.g. GLDC, IRF7/9 variants), or defects in cell-mediated immunity with increased baseline levels of systemic inflammation (obesity, pregnancy, advanced age). These mechanisms may explain the prolonged viral replication reported in critically ill patients with influenza: patients with life-threatening disease are, by definition, abnormal hosts. Understanding these molecular mechanisms of susceptibility may in the future enable the design of host-directed therapies to promote resilience.

Network analysis of canine brain morphometry links tumour risk to oestrogen deficiency and accelerated brain ageing

Structural ‘brain age’ is a valuable but complex biomarker for several brain disorders. The dog is an unrivalled comparator for neurological disease modeling, however canine brain morphometric diversity creates computational and statistical challenges. Using a data-driven approach, we explored complex interactions between patient metadata, brain morphometry, and neurological disease. Twenty-four morphometric parameters measured from 286 canine brain magnetic resonance imaging scans were combined with clinical parameters to generate 9,438 data points. Network analysis was used to cluster patients according to their brain morphometry profiles. An ‘aged-brain’ profile, defined by a small brain width and volume combined with ventriculomegaly, was revealed in the Boxer breed. Key features of this profile were paralleled in neutered female dogs which, relative to un-neutered females, had an 11-fold greater risk of developing brain tumours. Boxer dog and geriatric dog groups were both enriched for brain tumour diagnoses, despite a lack of geriatric Boxers within the cohort. Our findings suggest that advanced brain ageing enhances brain tumour risk in dogs and may be influenced by oestrogen deficiency—a risk factor for dementia and brain tumours in humans. Morphometric features of brain ageing in dogs, like humans, might better predict neurological disease risk than patient chronological age.

Integrating Coexpression Networks with GWAS to Prioritize Causal Genes in Maize

Genome-wide association studies (GWAS) have identified loci linked to hundreds of traits in many different species. Yet, because linkage equilibrium implicates a broad region surrounding each identified locus, the causal genes often remain unknown. This problem is especially pronounced in nonhuman, nonmodel species, where functional annotations are sparse and there is frequently little information available for prioritizing candidate genes. We developed a computational approach, Camoco, that integrates loci identified by GWAS with functional information derived from gene coexpression networks. Using Camoco, we prioritized candidate genes from a large-scale GWAS examining the accumulation of 17 different elements in maize (Zea mays) seeds. Strikingly, we observed a strong dependence in the performance of our approach based on the type of coexpression network used: expression variation across genetically diverse individuals in a relevant tissue context (in our case, roots that are the primary elemental uptake and delivery system) outperformed other alternative networks. Two candidate genes identified by our approach were validated using mutants. Our study demonstrates that coexpression networks provide a powerful basis for prioritizing candidate causal genes from GWAS loci but suggests that the success of such strategies can highly depend on the gene expression data context. Both the software and the lessons on integrating GWAS data with coexpression networks generalize to species beyond maize.

Exploring Coronary Artery Disease GWAs Targets With Functional Links to Immunometabolism

Finding genetic variants that cause functional disruption or regulatory change among the many implicated GWAs variants remains a key challenge to translating the findings from GWAs to therapeutic treatments. Defining the causal mechanisms behind the variants require functional screening experiments that can be complex and costly. Prioritizing variants for functional characterization using techniques that capture important functional and regulatory elements can assist this. The genetic architecture of complex traits such as cardiovascular disease and type II diabetes comprise an enormously large number of variants of small effect contributing to heritability and spread throughout the genome. This makes it difficult to distinguish which variants or core genes are most relevant for prioritization and how they contribute to the regulatory networks that become dysregulated leading to disease. Despite these challenges, recent GWAs for CAD prioritized genes associated with lipid metabolism, coagulation and adhesion along with novel signals related to innate immunity, adipose tissue and, vascular function as important core drivers of risk. We focus on three examples of novel signals associated with CAD which affect risk through missense or UTR mutations indicating their potential for therapeutic modification. These variants play roles in adipose tissue function vascular function and innate immunity which form the cornerstones of immuno-metabolism. In addition we have explored the putative, but potentially important interactions between the environment, specifically food and nutrition, with respect to key processes.