Nearby Genes Research Articles

Identify novel, shared and disorder-specific genetic architecture of major depressive disorder, insomnia and chronic pain

Major depressive disorder (MDD), insomnia (INS) and chronic pain (CP) often have high comorbidity and show high genetic correlation. Here we aimed to better characterize their novel, shared and disorder-specific genetic architecture. Based on genome-wide association study (GWAS) summary data, we applied the conditional false discovery rate (condFDR) and conjunctional FDR (conjFDR) approach to investigate the novel and overlapped genetic loci for MDD, INS and CP. In addition, putative disorder-specific SNP associations were analyzed by conditioning the other two traits. The functions of the identified genomic loci were explored by performing gene set enrichment analysis (GSEA) for the loci mapped genes. We identified 22, 43 and 91 novel risk loci for MDD, INS and CP. GSEA for the loci mapped genes highlighted olfactory signaling pathway for MDD novel loci, breast cancer related gene set for both INS and CP novel loci, and nervous system related development, structure and activity for CP. Furthermore, we identified three loci jointly associated with the three disorders, including 13q14.3 locus with nearby gene OLFM4, 14q21.1 locus with nearby gene LRFN5 and 5q21.2 locus located in intergenic region. In addition, we identified one specific loci for MDD, 7 for INS and 11 for CP respectively by conditioning the other two traits, which were mapped to 68 genes for MDD, 85 for INS and 100 for CP. The MDD specific genes are enriched in immune system related pathways. This study increases understanding of the genetic architectures underlying MDD, INS and CP. The shared underlying genetic risk may help to explain the high comorbidity rates of the disorders.

MEF2C-AS1 regulates its nearby gene MEF2C to mediate cervical cancer cell malignant phenotypes in vitro

Cervical cancer (CC) is the second most common malignancy among women. GEPIA demonstrated that MEF2C-AS1 and its nearby gene MEF2C present downregulation in CC tissues. We attempted to clarify molecular mechanism between MEF2C-AS1 and MEF2C underlying CC progression. RT-qPCR was used to measure expression levels and subcellular distribution of MEF2C-AS1 and MEF2C in CC cell lines. Gain-of-function assays were conducted to reveal roles of MEF2C-AS1 and MEF2C in CC cell behaviors. Bioinformatics, RNA pull down, and RIP assays were performed to assess association of MEF2C-AS1 or MEF2C with miR-20 b-5p in CC cells. Rescue assays were done to assess regulatory function of the MEF2C-AS1-miR-20 b-5p-MEF2C axis in CC cellular processes. MEF2C-AS1 and its nearby gene MEF2C showed downregulation and had a positive expression correlation in CC tissues. MEF2C-AS1 and MEF2C presented downregulation in CC cells, and they majorly distributed in CC cell cytoplasm. MEF2C-AS1 and MEF2C upregulation repressed CC cell proliferative, migratory, and angiogenic abilities. MEF2C-AS1 competitively bound with miR-20 b-5p to upregulate MEF2C in CC cells. The impacts of MEF2C-AS1 elevation on CC cell proliferative, migratory, and angiogenic capabilities were countervailed by miR-20 b-5p overexpression. The impacts of miR-20 b-5p inhibitor on CC cell proliferative, migratory and angiogenic capabilities were countervailed by MEF2C depletion. To sum up, MEF2C-AS1 and its nearby gene MEF2C present downregulation and serve as tumor suppressors in CC cells. MEF2C-AS1 suppresses CC cell malignancy in vitro through sponging miR-20 b-5p to upregulate MEF2C, which may provide a potential new direction for seeking therapeutic plans of CC.

Dynamics of transcriptome and chromatin accessibility revealed sequential regulation of potential transcription factors during the brown adipose tissue whitening in rabbits.

Brown adipose tissue (BAT) represents a valuable target for treating obesity in humans. BAT losses of thermogenic capacity and gains a “white adipose tissue-like (WAT-like)” phenotype (BAT whitening) under thermoneutral environments, which could lead to potential low therapy responsiveness in BAT-based obesity treatments. However, the epigenetic mechanisms of BAT whitening remain largely unknown. In this study, BATs were collected from rabbits at day0 (D0), D15, D85, and 2 years (Y2). RNA-sequencing (RNA-seq) and the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) were performed to investigate transcriptome and chromatin accessibility of BATs at the four whitening stages, respectively. Our data showed that many genes and chromatin accessible regions (refer to as “peaks”) were identified as significantly changed during BAT whitening in rabbits. The BAT-selective genes downregulated while WAT-selective genes upregulated from D0 to Y2, and the de novo lipogenesis-related genes reached the highest expression levels at D85. Both the highly expressed genes and accessible regions in Y2 were significantly enriched in immune response-related signal pathways. Analysis of different relationships between peaks and their nearby genes found an increased proportion of the synchronous changes between chromatin accessibility and gene expression during BAT whitening. The synergistic changes between the chromatin accessibility of promoter and the gene expression were found in the key adipose genes. The upregulated genes which contained increased peaks were significantly enriched in the PI3K-Akt signaling pathway, steroid biosynthesis, TGF-beta signaling pathway, osteoclast differentiation, and dilated cardiomyopathy. Moreover, the footprinting analysis suggested that sequential regulation of potential transcription factors (TFs) mediated the loss of thermogenic phenotype and the gain of a WAT-like phenotype of BAT. In conclusion, our study provided the transcriptional and epigenetic frameworks for understanding BAT whitening in rabbits for the first time and might facilitate potential insights into BAT-based obesity treatments.

Genome-wide association study of traits in sacred lotus uncovers MITE-associated variants underlying stamen petaloid and petal number variations.

Understanding the genetic variants responsible for floral trait diversity is important for the molecular breeding of ornamental flowers. Widely used in water gardening for thousands of years, the sacred lotus exhibits a wide range of diversity in floral organs. Nevertheless, the genetic variations underlying various morphological characteristics in lotus remain largely unclear. Here, we performed a genome-wide association study of sacred lotus for 12 well-recorded ornamental traits. Given a moderate linkage disequilibrium level of 32.9 kb, we successfully identified 149 candidate genes responsible for seven flower traits and plant size variations, including many pleiotropic genes affecting multiple floral-organ-related traits, such as NnKUP2. Notably, we found a 2.75-kb presence-and-absence genomic fragment significantly associated with stamen petaloid and petal number variations, which was further confirmed by re-examining another independent population dataset with petal number records. Intriguingly, this fragment carries MITE transposons bound by siRNAs and is related to the expression differentiation of a nearby candidate gene between few-petalled and double-petalled lotuses. Overall, these genetic variations and candidate genes responsible for diverse lotus traits revealed by our GWAS highlight the role of transposon variations, particularly MITEs, in shaping floral trait diversity.

The Heterogeneity in the Landscape of Gene Dominance in Maize is Accompanied by Unique Chromatin Environments.

Subgenome dominance after whole-genome duplication (WGD) has been observed in many plant species. However, the degree to which the chromatin environment affects this bias has not been explored. Here, we compared the dominant subgenome (maize1) and the recessive subgenome (maize2) with respect to patterns of sequence substitutions, genes expression, transposable element accumulation, small interfering RNAs, DNA methylation, histone modifications, and accessible chromatin regions (ACRs). Our data show that the degree of bias between subgenomes for all the measured variables does not vary significantly when both of the WGD genes are located in pericentromeric regions. Our data further indicate that the location of maize1 genes in chromosomal arms is pivotal for maize1 to maintain its dominance, but location has a less effect on maize2 homoeologs. In addition to homoeologous genes, we compared ACRs, which often harbor cis-regulatory elements, between the two subgenomes and demonstrate that maize1 ACRs have a higher level of chromatin accessibility, a lower level of sequence substitution, and are enriched in chromosomal arms. Furthermore, we find that a loss of maize1 ACRs near their nearby genes is associated with a reduction in purifying selection and expression of maize1 genes relative to their maize2 homoeologs. Taken together, our data suggest that chromatin environment and cis-regulatory elements are important determinants shaping the divergence and evolution of duplicated genes.

Genomics and epigenetics guided identification of tissue-specific genomic safe harbors

BackgroundGenomic safe harbors are regions of the genome that can maintain transgene expression without disrupting the function of host cells. Genomic safe harbors play an increasingly important role in improving the efficiency and safety of genome engineering. However, limited safe harbors have been identified.ResultsHere, we develop a framework to facilitate searches for genomic safe harbors by integrating information from polymorphic mobile element insertions that naturally occur in human populations, epigenomic signatures, and 3D chromatin organization. By applying our framework to polymorphic mobile element insertions identified in the 1000 Genomes project and the Genotype-Tissue Expression (GTEx) project, we identify 19 candidate safe harbors in blood cells and 5 in brain cells. For three candidate sites in blood, we demonstrate the stable expression of transgene without disrupting nearby genes in host erythroid cells. We also develop a computer program, Genomics and Epigenetic Guided Safe Harbor mapper (GEG-SH mapper), for knowledge-based tissue-specific genomic safe harbor selection.ConclusionsOur study provides a new knowledge-based framework to identify tissue-specific genomic safe harbors. In combination with the fast-growing genome engineering technologies, our approach has the potential to improve the overall safety and efficiency of gene and cell-based therapy in the near future.

Complex chromosomal rearrangements induced by transposons in maize.

Eukaryotic genomes are large and complex, and gene expression can be affected by multiple regulatory elements and their positions within the dynamic chromatin architecture. Transposable elements are known to play important roles in genome evolution, yet questions remain as to how transposable elements alter genome structure and affect gene expression. Previous studies have shown that genome rearrangements can be induced by Reversed Ends Transposition involving termini of Activator and related transposable elements in maize and other plants. Here, we show that complex alleles can be formed by the rapid and progressive accumulation of Activator-induced duplications and rearrangements. The p1 gene enhancer in maize can induce ectopic expression of the nearby p2 gene in pericarp tissue when placed near it via different structural rearrangements. By screening for p2 expression, we identified and studied 5 cases in which multiple sequential transposition events occurred and increased the p1 enhancer copy number. We see active p2 expression due to multiple copies of the p1 enhancer present near p2 in all 5 cases. The p1 enhancer effects are confirmed by the observation that loss of p2 expression is correlated with transposition-induced excision of the p1 enhancers. We also performed a targeted Chromosome Conformation Capture experiment to test the physical interaction between the p1 enhancer and p2 promoter region. Together, our results show that transposon-induced rearrangements can accumulate rapidly and progressively increase genetic variation important for genomic evolution.

Integrated Analysis of Transcriptome, microRNAs, and Chromatin Accessibility Revealed Potential Early B-Cell Factor1-Regulated Transcriptional Networks during the Early Development of Fetal Brown Adipose Tissues in Rabbits.

In domestic mammals, cold stress decreases the survival rate of newborns and increases the cost of management. Brown adipose tissue (BAT) is the main thermogenic organ for cubs, and well-developed fetal BAT (FBAT) is beneficial for newborns to maintain core temperatures during the first several days of life. However, our knowledge of the epigenetic mechanisms during the early development of FBAT remains largely unknown. Rabbits (Oryctolagus cuniculus) are economically important domestic animals. In this study, a histological analysis shows that the tissue content, thermogenic capacity, and lipid content of FBAT dramatically increases from gestational day 21 (G21) to gestational day 24 (G24) in rabbits. RNA-seq, microRNA-seq (miRNA-seq), and the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) show that many genes, miRNAs, and chromatin-accessible regions (referred to as peaks) were identified as significantly changed from G21 to G24, respectively. The upregulated genes from G21 to G24 were significantly enriched in the mitochondrial metabolism and thermogenesis-related signal pathways. The integrated analysis of transcriptome and chromatin accessibility reveals that the peaks in the promoters have a more regulatory effect than peaks in other genomic elements on the expression of their nearby genes in FBATs. The upregulated genes that are associated with increased chromatin accessibility in the promoter regions are involved in the energy metabolism-related signaling pathways. The genes that have a greater tendency to be regulated by miRNAs than the chromatin accessibility in gene promoters are involved in the apelin, insulin, and endocytosis signaling pathways. Furthermore, genome-wide transcription factor (TF) footprinting analysis identifies early B-cell factor1 (EBF1) as playing a key role during early FBAT development. The carbon metabolism, citrate cycle, and PPAR signaling pathways are significantly enriched by the predicted EBF1-regulated cascade TF-network. In conclusion, our work provides a framework for understanding epigenetics regulatory mechanisms underlying the early development of FBAT and identifies potential TF involved in the early development of FBAT in rabbits.

Shared regulation and functional relevance of local gene co-expression revealed by single cell analysis

Most human genes are co-expressed with a nearby gene. Previous studies have revealed this local gene co-expression to be widespread across chromosomes and across dozens of tissues. Yet, so far these studies used bulk RNA-seq, averaging gene expression measurements across millions of cells, thus being unclear if this co-expression stems from transcription events in single cells. Here, we leverage single cell datasets in >85 individuals to identify gene co-expression across cells, unbiased by cell-type heterogeneity and benefiting from the co-occurrence of transcription events in single cells. We discover >3800 co-expressed gene pairs in two human cell types, induced pluripotent stem cells (iPSCs) and lymphoblastoid cell lines (LCLs) and (i) compare single cell to bulk RNA-seq in identifying local gene co-expression, (ii) show that many co-expressed genes – but not the majority – are composed of functionally related genes and (iii) using proteomics data, provide evidence that their co-expression is maintained up to the protein level. Finally, using single cell RNA-sequencing (scRNA-seq) and single cell ATAC-sequencing (scATAC-seq) data for the same single cells, we identify gene-enhancer associations and reveal that >95% of co-expressed gene pairs share regulatory elements. These results elucidate the potential reasons for co-expression in single cell gene regulatory networks and warrant a deeper study of shared regulatory elements, in view of explaining disease comorbidity due to affecting several genes. Our in-depth view of local gene co-expression and regulatory element co-activity advances our understanding of the shared regulatory architecture between genes.

Signatures of Convergent Evolution and Natural Selection at the Alcohol Dehydrogenase Gene Region are Correlated with Agriculture in Ethnically Diverse Africans.

The alcohol dehydrogenase (ADH) family of genes encodes enzymes that catalyze the metabolism of ethanol into acetaldehyde. Nucleotide variation in ADH genes can affect the catalytic properties of these enzymes and is associated with a variety of traits, including alcoholism and cancer. Some ADH variants, including the ADH1B*48His (rs1229984) mutation in the ADH1B gene, reduce the risk of alcoholism and are under positive selection in multiple human populations. The advent of Neolithic agriculture and associated increase in fermented foods and beverages is hypothesized to have been a selective force acting on such variants. However, this hypothesis has not been tested in populations outside of Asia. Here, we use genome-wide selection scans to show that the ADH gene region is enriched for variants showing strong signals of positive selection in multiple Afroasiatic-speaking, agriculturalist populations from Ethiopia, and that this signal is unique among sub-Saharan Africans. We also observe strong selection signals at putatively functional variants in nearby lipid metabolism genes, which may influence evolutionary dynamics at the ADH region. Finally, we show that haplotypes carrying these selected variants were introduced into Northeast Africa from a West-Eurasian source within the last ∼2,000 years and experienced positive selection following admixture. These selection signals are not evident in nearby, genetically similar populations that practice hunting/gathering or pastoralist subsistence lifestyles, supporting the hypothesis that the emergence of agriculture shapes patterns of selection at ADH genes. Together, these results enhance our understanding of how adaptations to diverse environments and diets have influenced the African genomic landscape.

The Gallus gallus RJF reference genome reveals an MHCY haplotype organized in gene blocks that contain 107 loci including 45 specialized, polymorphic MHC class I loci, 41 C-type lectin-like loci, and other loci amid hundreds of transposable elements.

MHCY is a second major histocompatibility complex-like gene region in chickens originally identified by the presence of major histocompatibility complex class I-like and class II-like gene sequences. Up to now, the MHCY gene region has been poorly represented in genomic sequence data. A high density of repetitive sequence and multiple members of several gene families prevented the accurate assembly of short-read sequence data for MHCY. Identified here by single-molecule real-time sequencing sequencing of BAC clones for the Gallus gallus Red Jungle Fowl reference genome are 107 MHCY region genes (45 major histocompatibility complex class I-like, 41 c-type-lectin-like, 8 major histocompatibility complex class IIβ, 8 LENG9-like, 4 zinc finger protein loci, and a single only zinc finger-like locus) located amid hundreds of retroelements within 4 contigs representing the region. Sequences obtained for nearby ribosomal RNA genes have allowed MHCY to be precisely mapped with respect to the nucleolar organizer region. Gene sequences provide insights into the unusual structure of the MHCY class I molecules. The MHCY class I loci are polymorphic and group into 22 types based on predicted amino acid sequences. Some MHCY class I loci are full-length major histocompatibility complex class I genes. Others with altered gene structure are considered gene candidates. The amino acid side chains at many of the polymorphic positions in MHCY class I are directed away rather than into the antigen-binding groove as is typical of peptide-binding major histocompatibility complex class I molecules. Identical and nearly identical blocks of genomic sequence contribute to the observed multiplicity of identical MHCY genes and the large size (>639 kb) of the Red Jungle Fowl MHCY haplotype. Multiple points of hybridization observed in fluorescence in situ hybridization suggest that the Red Jungle Fowl MHCY haplotype is made up of linked, but physically separated genomic segments. The unusual gene content, the evidence of highly similar duplicated segments, and additional evidence of variation in haplotype size distinguish polymorphic MHCY from classical polymorphic major histocompatibility complex regions.

Development of SSR Markers Linked to Stress Responsive Genes along Tomato Chromosome 3 (Solanum lycopersicum L.).

This study aimed to develop novel SSR markers in tomato. Several BAC clones along chromosome 3 in tomato were selected based on their content. The criteria was the availability of genes, either directly or indirectly related to stress response (drought, salinity, and heat) in tomato. A total of 20 novel in silico SSR markers were developed and 96 important nearby genes were identified. The identified nearby genes represent different tomato genes involved in plant growth and development and biotic and abiotic stress tolerance. The developed SSR markers were assessed using tomato landraces. A total of 29 determinate and semi-determinate local tomato landraces collected from diverse environments were utilized. A total of 33 alleles with mean of 1.65 alleles per locus were scored, showing 100% polymorphic patterns, with a mean of 0.18 polymorphism information content (PIC) values. The mean of observed and expected heterozygosity were 0.19 and 0.24, respectively. The mean value of the Jaccard similarity index was used for clustering the landraces. The developed microsatellite markers showed potential to assess genetic variability among tomato landraces. The genetic distance information reported in this study can be used by breeders in future genetic improvement of tomato for tolerance against diverse stresses.

Fine-mapping and replication of EWAS loci harboring putative epigenetic alterations associated with AD neuropathology in a large collection of human brain tissue samples.

Our previous epigenome-wide association study (EWAS) of Alzheimer's disease (AD) in human brain identified 71 CpGs associated with AD pathology. However, due to low coverage of the Illumina platform, many important CpGs might have been missed. In a large collection of human brain tissue samples (N = 864), we fine-mapped previous EWAS loci by targeted bisulfite sequencing and examined their associations with AD neuropathology. DNA methylation was also linked to gene expression of the same brain cortex. Our targeted sequencing captured 130 CpGs (∼1.2 kb), 93 of which are novel. Of the 130 CpGs, 57 sites (only 17 included in previous EWAS) and 12 gene regions (e.g., ANK1, BIN1, RHBDF2, SPG7, PODXL) were significantly associated with amyloid load. DNA methylation in some regions was associated with expression of nearby genes. Targeted methylation sequencing can validate previous EWAS loci and discover novel CpGs associated with AD pathology.

Characterization of Transposon-Derived Accessible Chromatin Regions in Rice (Oryza Sativa).

Growing evidence indicates that transposons or transposable elements (TEs)-derived accessible chromatin regions (ACRs) play essential roles in multiple biological processes by interacting with trans-acting factors. However, the function of TE-derived ACRs in the regulation of gene expression in the rice genome has not been well characterized. In this study, we examined the chromatin dynamics in six types of rice tissues and found that ~8% of ACRs were derived from TEs and exhibited distinct levels of accessibility and conservation as compared to those without TEs. TEs exhibited a TE subtype-dependent impact on ACR formation, which can be mediated by changes in the underlying DNA methylation levels. Moreover, we found that tissue-specific TE-derived ACRs might function in the tissue development through the modulation of nearby gene expression. Interestingly, many genes in domestication sweeps were found to overlap with TE-derived ACRs, suggesting their potential functions in the rice domestication. In addition, we found that the expression divergence of 1070 duplicate gene pairs were associated with TE-derived ACRs and had distinct distributions of TEs and ACRs around the transcription start sites (TSSs), which may experience different selection pressures. Thus, our study provides some insights into the biological implications of TE-derived ACRs in the rice genome. Our results imply that these ACRs are likely involved in the regulation of tissue development, rice domestication and functional divergence of duplicated genes.

Functional analysis of the Vsx2 super-enhancer uncovers distinct cis-regulatory circuits controlling Vsx2 expression during retinogenesis.

Vsx2 is a transcription factor essential for retinal proliferation and bipolar cell differentiation, but the molecular mechanisms underlying its developmental roles are unclear. Here, we have profiled VSX2 genomic occupancy during mouse retinogenesis, revealing extensive retinal genetic programs associated with VSX2 during development. VSX2 binds and transactivates its enhancer in association with the transcription factor PAX6. Mice harboring deletions in the Vsx2 regulatory landscape exhibit specific abnormalities in retinal proliferation and in bipolar cell differentiation. In one of those deletions, a complete loss of bipolar cells is associated with a bias towards photoreceptor production. VSX2 occupies cis-regulatory elements nearby genes associated with photoreceptor differentiation and homeostasis in the adult mouse and human retina, including a conserved region nearby Prdm1, a factor implicated in the specification of rod photoreceptors and suppression of bipolar cell fate. VSX2 interacts with the transcription factor OTX2 and can act to suppress OTX2-dependent enhancer transactivation of the Prdm1 enhancer. Taken together, our analyses indicate that Vsx2 expression can be temporally and spatially uncoupled at the enhancer level, and they illuminate important mechanistic insights into how VSX2 is engaged with gene regulatory networks that are essential for retinal proliferation and cell fate acquisition.

HPV Integration Site Mapping: A Rapid Method of Viral Integration Site (VIS) Analysis and Visualization Using Automated Workflows in CLC Microbial Genomics.

Human papillomavirus (HPV) integration within the host genome may contribute to carcinogenesis through various disruptive mechanisms. With next-generation sequencing (NGS), identification of viral and host genomic breakpoints and chimeric sequences are now possible. However, a simple, streamlined bioinformatics workflow has been non-existent until recently. Here, we tested two new, automated workflows in CLC Microbial Genomics, i.e., Viral Hybrid Capture (VHC) Data Analysis and Viral Integration Site (VIS) Identification for software performance and efficiency. The workflows embedded with HPV and human reference genomes were used to analyze a publicly available NGS dataset derived from pre- and cancerous HPV+ cervical cytology of 21 Gabonese women. The VHC and VIS workflow median runtimes were 19 and 7 min per sample, respectively. The VIS dynamic graphical outputs included read mappings, virus-host genomic breakpoints, and virus-host integration circular plots. Key findings, including disrupted and nearby genes, were summarized in an auto-generated report. Overall, the VHC and VIS workflows proved to be a rapid and accurate means of localizing viral-host integration site(s) and identifying disrupted and neighboring human genes. Applying HPV VIS-mapping to pre- or invasive tumors will advance our understanding of viral oncogenesis and facilitate the discovery of prognostic biomarkers and therapeutic targets.

DArTseq molecular markers associated with the piping leaf margin phenotype in pineapple (Ananas comosus L.)

This study sought to understand the genetic basis of the piping leaf margin phenotype in pineapple. To achieve this aim, a genome-wide association study (GWAS) using mixed linear regression and logistic regression analysis was conducted on three pineapple diversity panels including seedling populations segregating for spiny, spiny-tip and piping leaf margins. This study identified single nucleotide polymorphism (SNP) markers associated with the piping and spiny-tip leaf margin phenotypes. A broad quantitative trait locus (QTL) positioned on chromosome 23 between positions 240,475 and 2,369,197 bp was the most highly associated with piping leaf margin in all analyses. Major candidate genes proposed are a Zinc finger protein 2, a Zinc finger protein 3, a WUSCHEL-related homeobox 2, a WUSCHEL-related homeobox 1 and a Zinc finger protein CONSTANS-like. Some other genes of a lower association, linked or nearby genes of interest, are also considered potentially involved to varying degrees. All candidate genes are known to be involved in aspects of stem cell maintenance, cell proliferation, epidermal cell differentiation, organogenesis, leaf polarity, cell wall modification or hormone signalling. It is possible each plays a role in either differentiation or morphological aspects of the spiny-tip and piping leaf margin phenotypes. It is expected the relative role of each associated gene might vary with genetic background.

Epigenetic Reprogramming Leads to Downregulation of CD4 and Functional Changes in African Green Monkey Memory CD4+ T Cells.

African green monkeys (AGMs), Chlorocebus pygerythrus, are a natural host for a lentivirus related to HIV, SIV. SIV-infected AGMs rarely progress to AIDS despite robust viral replication. Though multiple mechanisms are involved, a primary component is the animals' ability to downregulate CD4 expression on mature CD4+ Th cells, rendering these cells resistant to infection by SIV. These CD8αα+ T cells retain functional characteristics of CD4+ Th cells while simultaneously acquiring abilities of cytotoxic CD8αβ+ T cells. To determine mechanisms underlying functional differences between T cell subsets in AGMs, chromatin accessibility in purified populations was determined by assay for transposase-accessible chromatin sequencing. Differences in chromatin accessibility alone were sufficient to cluster cells by subtype, and accessibility at the CD4 locus reflected changes in CD4 expression. DNA methylation at the CD4 locus also correlated with inaccessible chromatin. By associating accessible regions with nearby genes, gene expression was found to correlate with accessibility changes. T cell and immune system activation pathways were identified when comparing regions that changed accessibility from CD4+ T cells to CD8αα+ T cells. Different transcription factor binding sites are revealed as chromatin accessibility changes, and these differences may elicit downstream changes in differentiation. This comprehensive description of the epigenetic landscape of AGM T cells identified genes and pathways that could have translational value in therapeutic approaches recapitulating the protective effects CD4 downregulation.

Identifying genes targeted by disease-associated non-coding SNPs with a protein knowledge graph.

Genome-wide association studies (GWAS) have identified many single nucleotide polymorphisms (SNPs) that play important roles in the genetic heritability of traits and diseases. With most of these SNPs located on the non-coding part of the genome, it is currently assumed that these SNPs influence the expression of nearby genes on the genome. However, identifying which genes are targeted by these disease-associated SNPs remains challenging. In the past, protein knowledge graphs have often been used to identify genes that are associated with disease, also referred to as “disease genes”. Here, we explore whether protein knowledge graphs can be used to identify genes that are targeted by disease-associated non-coding SNPs by testing and comparing the performance of six existing methods for a protein knowledge graph, four of which were developed for disease gene identification. We compare our performance against two baselines: (1) an existing state-of-the-art method that is based on guilt-by-association, and (2) the leading assumption that SNPs target the nearest gene on the genome. We test these methods with four reference sets, three of which were obtained by different means. Furthermore, we combine methods to investigate whether their combination improves performance. We find that protein knowledge graphs that include predicate information perform comparable to the current state of the art, achieving an area under the receiver operating characteristic curve (AUC) of 79.6% on average across all four reference sets. Protein knowledge graphs that lack predicate information perform comparable to our other baseline (genetic distance) which achieved an AUC of 75.7% across all four reference sets. Combining multiple methods improved performance to 84.9% AUC. We conclude that methods for a protein knowledge graph can be used to identify which genes are targeted by disease-associated non-coding SNPs.

A super pan-genomic landscape of rice

Pan-genomes from large natural populations can capture genetic diversity and reveal genomic complexity. Using de novo long-read assembly, we generated a graph-based super pan-genome of rice consisting of a 251-accession panel comprising both cultivated and wild species of Asian and African rice. Our pan-genome reveals extensive structural variations (SVs) and gene presence/absence variations. Additionally, our pan-genome enables the accurate identification of nucleotide-binding leucine-rich repeat genes and characterization of their inter- and intraspecific diversity. Moreover, we uncovered grain weight-associated SVs which specify traits by affecting the expression of their nearby genes. We characterized genetic variants associated with submergence tolerance, seed shattering and plant architecture and found independent selection for a common set of genes that drove adaptation and domestication in Asian and African rice. This super pan-genome facilitates pinpointing of lineage-specific haplotypes for trait-associated genes and provides insights into the evolutionary events that have shaped the genomic architecture of various rice species.