Search

Background: Genome-wide association studies have identified numerous disease-associated variants, but a vast majority are located in non-coding regions, making it challenging to understand their functional impact. This complexity necessitates new techniques to identify causal variants in non-coding regions and elucidate their specific cellular contexts and mechanisms of action. Here we present work i) examining mutations that create nullomers in the human genome to explore its potential utility in identifying pathogenic mutations and ii) a single-cell multi-omic study identifying the transcriptome and regulome of the human and mouse hypothalamus to identify regulatory regions of obesity-associated variants.Methods: (i) We generated all possible mutations of the human genome that can lead to emergence of a nullomer, and examine where in the genome they emerge. (ii) We apply single-cell RNA and ATAC sequencing to adult hypothalamus samples.Results and Conclusions: (i) Our findings highlight CpG hypermutability and methylated cytosines as key elements leading to resurfacing of nullomers in individuals. We also showcase that nullomers can have applications in disease annotation and pathogenic variant identification. (ii) We identified regulatory elements of hypothalamus cell types and mapped obesity-associated variants to cell-type specific peaks. We validated these regions to be enhancers using CRISPR editing and CRISPRi.

Whole genome sequences allow for new, powerful inferences of evolutionary history. In addition to increasing the scale of application of traditional population genetic inferences by using them on an enormous amount of loci, genomic data also provide valuable information about the location of those sites and correlations among them, which can be leveraged for inference of evolutionary parameters. In this dissertation, I explore three types of evolutionary inferences that are thriving with the increasing availability of genomic data and new methods. In Chapter 1, I use simulations to evaluate methods for inference of ancestral recombination graphs. In Chapter 2, I apply a new method based on the pairwise sequential Markovian coalescent to infer population split times and migration rates from a pair of diploid genomes. In Chapter 3, I use population-level genomic data from the population of a rural village in Ecuador to infer signatures of selection possibly related to the beneficial effects of diet on their cardiovascular health.

The use of next-generation sequencing technology has been extremely useful for genome-wide analysis in many different organisms. For example, transcriptomic analysis based on RNA-seq can identify all changes in gene expression in a given tissue in a single experiment. Other types of experiments, such as ChIP-seq, allow mapping of transcription factor binding sites across the genome. Much of the work in this dissertation focuses on the transcription factor hepatocyte nuclear factor 4a (HNF4a) which is a member of the nuclear receptor superfamily of ligand-dependent transcription factors and abundantly expressed in the liver and intestines. HNF4a has been linked to several diseases including diabetes, liver and colon cancer, inflammatory bowel disease and others. The expression of HNF4a is driven by two promoters (P1 and P2) which result in the expression of 12 different isoforms. Additionally, HNF4a has as its endogenous ligand an essential fatty acid (linoleic acid, LA) that must be obtained from the diet and is known to play a role in gluconeogenesis during fasting. Therefore, we examine fasted vs fed conditions as well as high fat diets with different amounts of LA on gene expression in the liver and intestines, respectively. Finally, HNF4a target genes have been extensively studied and mutations in HNF4a binding sites in regulatory regions of select target genes associated with disease have been identified. However, to date, there has been no systematic approach to identify variants in HNF4a binding sites of target genes on a genome-wide scale.In Chapter 1 of this dissertation, we provide an introduction to the biological and research problems addressed in the remainder of the dissertation and review the relevant literature. In Chapter 2, we analyze transcriptomic (RNA-seq) data from male mice fed three different high fat diets (HFDs) with different amounts of LA from four different tissues across the mouse intestinal tract. We found that different portions of the intestines have unique gene profiles, especially in terms of nuclear receptor signaling, xenobiotic and drug metabolism, and intestinal epithelial barrier function. We found that different types of HFDs impact gene expression in the intestinal tract in different ways, including genes linked to diseases such as inflammatory bowel disease and colon cancer. A network analysis revealed that mice fed the HFDs had altered expression of intestinal genes involved in crucial body functions such as the immune system and the intestinal microbiome which could facilitate bacterial and viral infections such as SARS-CoV-2. To our knowledge, this is the most comprehensive dataset including multiple diets and multiple parts of the intestines and should provide an excellent resource for the scientific community for some time to come.In Chapter 3, we compared the impact of fasting and HNF4a isoforms on liver gene expression and alternative splicing, comparing and contrasting two RNA-seq datasets from male and female mice. We also compared wild type mice (express predominantly P1-HNF4a isoforms) to exon swap mice (a7HMZ) which express only the P2 isoforms of HNF4a which are typically not expressed in the adult liver. We found that a 12-hour fast has a significant effect on alternative splicing and that the P2-HNF4a isoforms seem to play a role. Interestingly, there was a different pattern of alternative splicing in female livers.In Chapter 4, we analyze the impact of single nucleotide polymorphisms (SNPs) on HNF4a binding sites using publicly available datasets of HNF4a ChIP-seq (chromatin immunoprecipitation) from human liver and the liver cancer cell line HepG2 which expresses HNF4a and exhibits a hepatocyte phenotype. For this purpose, we trained a support vector machine (SVM) to predict binding affinity scores based on DNA sequences known to bind HNF4a from protein binding microarray (PBM) data. This model allowed us to identify 10 putative affinity altering SNPs (aaSNPs) in the human liver HNF4a ChIP-seq and six in the HepG2 ChIP-seq that impact the binding of HNF4a to the chromatin. Additionally, some of these identified SNPs were found in genes related to cancer, suggesting a potential role in personalized medicine. These results present a proof concept of affinity disruption in binding between HNF4a and some target genes in a genome-wide analysis.Finally, in Chapter 5, we provide an example of how transcriptomic data could be used to generate new biological hypotheses and test them. Additionally, we provide future directions for Chapter 3 and Chapter 4.

The genomic sequence of the horse has been available since 2007, providing critical resources for discovering important genomic variants regarding both animal health and population structure. However, to fully understand the functional implications of these variants, detailed annotation of the horse genome is required. Currently, the horse genome is annotated using the limited available RNA-seq data, as well as through comparative genomics by translating human and mouse genome annotation. While this approach has served the equine researchers well and led to a number of discoveries improving the care and management of horses, many important questions remain unanswered. The limitation of the current annotation is two pronged. First, a comparative genomics approach is insufficient to identify many genes that are less evolutionarily conserved, especially those that are noncoding. The sole reliance on short-read RNA-seq data also meant that alternate isoforms could not be accurately resolved. Second, epigenomic regulatory elements are crucial to detailed understanding of gene expression network but are yet to be systemically identified in the horse. Many regulatory elements, including enhancers, promoters, and insulators, are not transcribed or transcribed at a very low level, necessitating alternate approaches to identify them. To solve these problems, the Functional Annotation of the Animal Genomes (FAANG) project proposed a systemic approach to tissue collection, phenotyping, and data generation, adopting the blueprint laid out by the Encyclopedia of DNA Elements (ENCODE) project. This thesis describes the equine FAANG team’s effort to map tissue-specific gene expression and regulation in the horse genome. Chapter 1 provides an overview of the equine FAANG project’s approach to functional annotation. Chapter 2 describes an improved transcriptome that includes novel genes and alternate isoforms compared to the current annotation. In Chapter 3, we use ATAC-seq to create a catalog of tissue-specific open chromatin regions, which can serve as proxies to active regulatory elements. Chapter 4 provides a complete annotation of chromatin states across nine tissues. The Addendum detailed our effort to validate assays for transposase accessible chromatin using sequencing (ATAC-seq) in both frozen tissues and cryopreserved nuclei from fresh tissues. This thesis presents the first comprehensive overview of gene expression and their regulation in the horse, enabling interrogation of complex gene regulatory network and further studies of complex traits in horses. Future work should focus on both widening the scope of the equine FAANG project by including more tissue types and developmental stages, as well as refining gene network at single-cell resolution.

Elucidating the genetic changes underlying the evolution of human traits remains an unfinished puzzle. Structural variants (SVs) account for more genetic differences than single-nucleotide polymorphisms between humans and our closest living relatives, chimpanzees, and are a hallmark of great ape evolution. The genomes of great apes are enriched in large interspersed segmental duplications (SDs), defined as duplications larger than 1 kbp with over 90% sequence identity, that sensitize the genome to further genomic rearrangements, including copy-number variation, via non- allelic homologous recombination. Despite their relevance, the identification and characterization of these SVs has been hindered by short reads lengths as they lack enough sequence context to discover breakpoints and cannot unequivocally be mapped to highly identical duplicates. Long-read sequencing technologies overcome these limitations by providing reads thousands of bases long, but the availability of population cohorts remains limited.This thesis studies primate SVs and SDs characterized using diverse sequencing technologies and assesses their representation in reference genomes, variation across modern populations, their putative molecular impacts, and their roles in evolution and adaptation. We found novel SVs, including 88 deletions and 36 inversions, in two chimpanzee individuals sequenced with nanopore and optical mapping. Deletions and inversion breakpoints were depleted within topologically associated domains but enriched in differentially expressed genes between the two species. Focusing on human SDs, we identified eight Mbp of erroneously collapsed duplications in the human reference genome, impacting 48 protein coding and ten medically relevant genes, that are corrected in the first complete sequence of a human genome, T2T-CHM13. Leveraging this new reference, we identified 417 genes embedded in SDs with over 98% sequence identity (SD-98) that are near copy-number (CN) fixed in modern humans (1000 Genomes Project; 1KGP), 205 genes showing stratification between diverse modern populations (VST>95th percentile), and 22 protein-encoding genes showing consistent Tajima’s D outlier values across all humans examined. Our approach highlighted potential relevant human gene duplications, which are priority candidates for functional studies. Finally, we provide a compendium of tools and practices that we recommend be adopted by computational biologists to increase reproducibility in the field.

Small RNAs are short, non-coding RNA molecules that have been identified in a wide range of species across all three domains of life. In mammals, canonical small RNAs such as microRNAs (miRNAs) and Piwi-interacting RNAs (piRNAs) are tissue-specifically distributed and can regulate gene expression at both transcriptional and post-transcriptional levels, associating with various fundamental functions such as gene silencing and retrotransposon control. While miRNAs and piRNAs have been extensively investigated, the existence and mechanism of other non-canonical mammalian small RNAs remain underexplored.With the extensive use of high-throughput sequencing technologies in the past decade, small RNA diversity is rapidly growing. However, the conventional small RNA library construction method lacks the detection capability for non-canonical small RNAs that carry specific terminal and internal modifications, especially for tRNA-derived small RNAs (tsRNAs) and rRNA-derived small RNAs (rsRNAs). In addition, existing downstream bioinformatics tools are mostly focused on analyzing canonical small RNAs such as miRNAs and piRNAs, while the annotation of other small RNAs remains rudimentary or ignored.To reveal a panoramic view of small RNAs, a small RNA annotation pipeline (that is, SPORTS) is developed to simultaneously and comparatively annotate both canonical and non-canonical small RNAs, along with RNA modification prediction capacity. Moreover, a small RNA library preparation procedure (that is, PANDORA-seq) is optimized to comprehensively capture modified small RNAs such as tsRNAs and rsRNAs.The improved RNA-seq and bioinformatics strategy leads to a new and surprising small RNA landscape that tsRNAs and rsRNAs have a higher abundance than canonical small RNAs in a majority of mouse and human tissues/cells that have been examined. Those mammalian tsRNAs and rsRNAs also exhibit tissue- and cell-specific patterns and the expression level of those small RNAs are dynamically altered during the generation of induced pluripotent stem cells (iPSCs) based on PANDORA-seq. Those newly identified small RNAs also display translational regulation during embryonic stem cell differentiation and have a role in regulating embryonic stem cell lineage fate based on the transcriptomic changes after their transfection.

Xanthomonas arboricola pv. corylina is the causal agent of bacterial blight of hazelnut. The bacterium has been listed as an A2 quarantine pathogen in Europe since 1978 and on the regulated non-quarantine pest list since 2019. Three isolates from various geographic regions and isolated at different times were sequenced using a hybrid approach with short- and long-read technologies to generate closed genome and plasmid sequences in order to better understand the biology of this pathogen.

Insertions and deletions (indels) enable evolution and cause disease. Due to technical challenges, indels are left out of most mutational scans, limiting our understanding of them in disease, biology, and evolution. We develop a low cost and bias method, DIMPLE, for systematically generating deletions, insertions, and missense mutations in genes, which we test on a range of targets, including Kir2.1. We use DIMPLE to study how indels impact potassium channel structure, disease, and evolution. We find deletions are most disruptive overall, beta sheets are most sensitive to indels, and flexible loops are sensitive to deletions yet tolerate insertions.

Only 40 years ago, just three individuals made the journey from Finland/Russia to found the current wolf population in the southwest of Sweden. This population, that to this date descends from less than 10 founders, has a substantial increased extinction risk due to inbreeding. Several previous studies have used SNPs to monitor the level of inbreeding and homozygosity in the population, as well as measure immigration and the inflow of new genetic material. This study uses both short- and long-read data to discover structural variants (SVs) and small indels in the population, so that they may be used to extend the analyses and provide more insight into the current state of the Scandinavian wolf population. After the calling of the SVs, strict filtering and manual curation were applied to the data, thereby removing many false positive calls and increasing confidence in the remaining SVs. Short-read and long-read SV-callers found 31,800 and 57,821 SVs respectively, with relatively little overlap between the two sets. By far, the most common SV-types were deletions and insertions, at about 30,000 each with varying length ranging from a 50 base pairs to several tens of Mbp. Analyses on the data, such as PCAs and parent-offspring trio analyses, reveal high-confidence calls and consistent results between SV-types and SV-callers, as well as a low estimated genotyping error rate. PCAs performed on the SVs resembled those performed on SNPs, which strengthens the credibility of the identified variants. Finally, this study suggests several alternative steps for possible improvement to the dataset, along with some proposals for subsequent research topics that may use the variants discovered in this study.

Genome-wide association studies (GWAS) have suggested that the underlying genetic basis of complex traits and disease is driven by large numbers of non-coding variants with modest effects that likely act by modifying gene regulation. Towards understanding the regulatory impact of genetic variation, expression quantitative trait loci (eQTL) analyses have been performed across dozens of human tissues to link the influence of genetic variants on gene expression levels. While these eQTL studies have provided important biological insights, they are still limited by not considering the contexts in which these variants function, including the stage of development and cell type. Specifically, others have shown increased disease risk in adulthood has links to fetal origins, suggesting that characterizing gene expression in fetal-like cells could identify genetic variants that are associated with adult traits, but function primarily or solely during development. Additionally, as eQTL studies are typically performed across bulk tissues, unaccounted for cellular heterogeneity present in bulk gene expression measurements can affect genotype-gene expression associations. Thus, it is important to identify regulatory variants that alter gene expression in both primitive and adult cell types and to characterize cellular heterogeneity across tissues to comprehensively understand the genetic basis of complex traits and disease.Here, I present two studies, which utilize gene expression data from fetal-like and adult tissues to characterize cellular heterogeneity at distinct stages of human development. I have examined the cellular heterogeneity in fetal-like induced pluripotent stem cell (iPSC)-derived cardiovascular progenitor cells (CVPCs) using single cell (sc)RNA-seq data to identify cell populations that emerge as a result of the cardiac differentiation. Further, I deconvoluted 180 iPSC-CVPCs and identified factors innate to iPSCs that impacted cardiac fate. Next, I showed that mouse scRNA-seq can be used as an alternative to human scRNA-seq for the deconvolution of adult GTEx bulk tissues and considering cell composition eQTL studies powered the discovery of novel eQTLs, some of which were cell-type-associated and colocalized with GWAS disease loci.

This work was funded by the United States Department of Agriculture–Agricultural Research Service appropriated project 2034-22000-012-00D and the Xylella fastidiosa Active Containment Through a Multidisciplinary-Oriented Research Strategy (XF-ACTORS) project. Illumina sequencing for BB08-1 was carried out at the University of California, Davis, Genome Center DNA Technologies and Expression Analysis Cores, supported by National Institute of Health Shared Instrumentation Grant 1S10OD010786-01.

Stroke is one of the leading causes of patients' death and long‐term disability worldwide, and ischaemic stroke (IS) accounts for nearly 80% of all strokes. Differential genes and weighted gene co‐expression network analysis (WGCNA) in male and female patients with IS were compared. The authors used cell type identification by estimating relative subsets of RNA transcripts (CIBERSORT) to analyse the distribution pattern of immune subtypes between male and female patients. In this study, 141 up‐regulated and 61 down‐regulated genes were gathered and distributed into five modules in response to their correlation degree to clinical traits. The criterion for Gene Ontology (GO) term and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway indicated that detailed analysis had the potential to enhance clinical prediction and to identify the gender‐related mechanism. After that, the expression levels of hub genes were measured via the quantitative real‐time PCR (qRT‐PCR) method. Finally, CCL20, ICAM1 and PTGS2 were identified and these may be some promising targets for sex differences in IS. Besides, the hub genes were further verified by rat experiments. Furthermore, these CIBERSORT results showed that T cells CD8 and Monocytes may be the target for the treatment of male and female patients, respectively.

Hypertension is a leading risk factor for stroke, heart disease, chronic kidney disease, vascular cognitive impairment, and Alzheimer’s disease. Previous genetic studies have nominated hundreds of genes linked to hypertension, and renal and cognitive diseases. Some have been advanced as candidate genes by showing that they can alter blood pressure or renal and cerebral vascular function in knockout animals; however, final validation of the causal variants and underlying mechanisms has remained elusive. This review chronicles 40 years of work, from the initial identification of adducin (ADD) as an ACTIN-binding protein suggested to increase blood pressure in Milan hypertensive rats, to the discovery of a mutation in ADD1 as a candidate gene for hypertension in rats that were subsequently linked to hypertension in man. More recently, a recessive K572Q mutation in ADD3 was identified in Fawn-Hooded Hypertensive (FHH) and Milan Normotensive (MNS) rats that develop renal disease, which is absent in resistant strains. ADD3 dimerizes with ADD1 to form functional ADD protein. The mutation in ADD3 disrupts a critical ACTIN-binding site necessary for its interactions with actin and spectrin to regulate the cytoskeleton. Studies using Add3 KO and transgenic strains, as well as a genetic complementation study in FHH and MNS rats, confirmed that the K572Q mutation in ADD3 plays a causal role in altering the myogenic response and autoregulation of renal and cerebral blood flow, resulting in increased susceptibility to hypertension-induced renal disease and cerebral vascular and cognitive dysfunction.

KiSS-1, first identified as an anti-metastasis gene in melanoma, encodes C-terminally amidated peptide products, including kisspeptin-145, kisspeptin-54, kisspeptin-14, kisspeptin-13 and kisspeptin-10. These products are endogenous ligands coupled to G protein-coupled receptor 54 (GPR54)/hOT7T175/AXOR12. To date, the regulatory activities of the KiSS-1/GPR54 system, such as puberty initiation, antitumor metastasis, fertility in adulthood, hypothalamic-pituitary-gonadal axis (HPG axis) feedback, and trophoblast invasion, have been investigated intensively. Accumulating evidence has demonstrated that KiSS-1 played a key role in reproduction and served as a promising biomarker relative to the diagnosis, identification of therapeutic targets and prognosis in various carcinomas, while few studies have systematically summarized its subjective factors and concluded the functions of KiSS-1/GPR54 signaling in physiology homeostasis and cancer biology. In this review, we retrospectively summarized the regulators of the KiSS-1/GPR54 system in different animal models and reviewed its functions according to physiological homeostasis regulations and above all, cancer biology, which provided us with a profound understanding of applying the KiSS-1/GPR54 system into medical applications.

Neurodegenerative diseases, characterized by progressive neural loss, have been some of the most challenging medical problems in aging societies. Treatment strategies such as symptom management have little impact on disease progression, while intervention with specific disease mechanisms may only slow down disease progression. One therapeutic strategy that has the potential to reverse the disease phenotype is to replenish neurons and rebuild the pathway lost to degeneration. Although it is generally believed that the central nervous system has lost the capability to regenerate, increasing evidence indicates that the brain is more plastic than previously thought, containing perhaps the biggest repertoire of cells with latent neurogenic programs in the body. This review focuses on key advances in generating new neurons through in situ neuronal reprogramming, which is tied to fundamental questions regarding adult neurogenesis, cell source, and mechanisms for neuronal reprogramming, as well as the ability of new neurons to integrate into the existing circuitry.

Barth Syndrome is an X-linked disorder of mitochondrial cardiolipin metabolism caused by pathogenic variants in TAFAZZIN with pleiotropic effects including cardiomyopathy, neutropenia, growth delay, and skeletal myopathy. Management requires a multidisciplinary approach to the organ-specific manifestations including specialists from cardiology, hematology, nutrition, physical therapy, genetics, and metabolism. Currently, treatment is centered on management of specific clinical features, and is not targeted towards remediating the underlying biochemical defect. However, two clinical trials have been recently undertaken which target the mitochondrial pathology of this disease: a study to examine the effects of elamipretide, a cardiolipin targeted agent, and a study to examine the effects of bezafibrate, a peroxisome proliferator-activated receptor (PPAR) agonist. Treatments to directly target the defective TAFAZZIN pathway are under development, including enzyme and gene therapies. This article is protected by copyright. All rights reserved.

Long non-coding RNAs (lncRNAs) have been gaining importance in the field of cancer research in recent years. PRNCR1 (prostate cancer-associated non-coding RNA1) is a 12.7 kb, intron-less lncRNA found to play an oncogenic role in malignancy of diverse organs including prostate, breast, lung, oral cavity, colon and rectum. Single-nucleotide polymorphisms (SNPs) of PRNCR1 locus have been found to be associated with cancer susceptibility in different populations. In this review, an attempt has been made for the first time to summarize all sorts of available data on PRNCR1 to date from relevant databases (GeneCard, LncExpDB, Ensembl genome browser, and PubMed). As functional roles of PRNCR1, miRNA (microRNA) sponging was mostly highlighted in the pathogenesis of different cancer; in addition, an association of the lncRNA with chromatin-modifying complex to enhance androgen receptor-mediated gene transcription was reported in prostate cancer. Diagnostic and prognostic importance of PRNCR1 was found in some malignancies suggesting potency of the lncRNA to serve as a clinical biomarker. For PRNCR1 SNPs, although cancer susceptibility of the risk alleles/genotypes was reported in different populations, majorities of the findings were not replicated and underlying molecular mechanisms remained unexplored. Therapeutic implication of PRNCR1 was not studied well and future research may come up in this direction for intervening novel strategies to fight against cancer.

The European Bioinformatics Institute (EMBL-EBI) maintains a comprehensive range of freely available and up-to-date molecular data resources, which includes over 40 resources covering every major data type in the life sciences. This year's service update for EMBL-EBI includes new resources, PGS Catalog and AlphaFold DB, and updates on existing resources, including the COVID-19 Data Platform, trRosetta and RoseTTAfold models introduced in Pfam and InterPro, and the launch of Genome Integrations with Function and Sequence by UniProt and Ensembl. Furthermore, we highlight projects through which EMBL-EBI has contributed to the development of community-driven data standards and guidelines, including the Recommended Metadata for Biological Images (REMBI), and the BioModels Reproducibility Scorecard. Training is one of EMBL-EBI’s core missions and a key component of the provision of bioinformatics services to users: this year's update includes many of the improvements that have been developed to EMBL-EBI’s online training offering.

Funder: Kennedy Trust Rheumatology Research Prize Studentship, Funder: DFG Cluster of Excellence ���Precision Medicine in Chronic In-flammation��� (PMI; ID: EXC2167), Funder: EC | EC Seventh Framework Programm | FP7 Ideas: European Research Council (FP7-IDEAS-ERC - Specific Programme: ���Ideas��� Implementing the Seventh Framework Programme of the European Community for Research, Technological Development and Demonstration Activities (2007 to 2013)); doi: https://doi.org/10.13039/100011199; Grant(s): 715772, Funder: NWO-VIDI grant 016.178.056, the Netherlands Heart Foundation CVON grant 2018-27, and NWO Gravitation grant ExposomeNL, Funder: Li Ka Shing Foundation (Li Ka Shing Foundation Limited); doi: https://doi.org/10.13039/100007421, Irritable bowel syndrome (IBS) results from disordered brain���gut interactions. Identifying susceptibility genes could highlight the underlying pathophysiological mechanisms. We designed a digestive health questionnaire for UK Biobank and combined identified cases with IBS with independent cohorts. We conducted a genome-wide association study with 53,400 cases and 433,201 controls and replicated significant associations in a 23andMe panel (205,252 cases and 1,384,055 controls). Our study identified and confirmed six genetic susceptibility loci for IBS. Implicated genes included NCAM1, CADM2, PHF2/FAM120A, DOCK9, CKAP2/TPTE2P3 and BAG6. The first four are associated with mood and anxiety disorders, expressed in the nervous system, or both. Mirroring this, we also found strong genome-wide correlation between the risk of IBS and anxiety, neuroticism and depression (rg > 0.5). Additional analyses suggested this arises due to shared pathogenic pathways rather than, for example, anxiety causing abdominal symptoms. Implicated mechanisms require further exploration to help understand the altered brain���gut interactions underlying IBS.

Background Rhinoviruses (RV) are associated with the development and exacerbations of asthma and chronic obstructive pulmonary disease. They've also been linked to more severe diseases like pneumonia, acute bronchiolitis, croup, and otitis media. Because of the hypervariable sequences in the same serotypes, no effective vaccine against rhinoviruses has been developed to date. With the availability of new full-length genome sequences for all RV-A and RV-B serotyped strains, this study used bioinformatics to find a suitable RV strain with the highest similarity matrices to the other strains. Methods The full genomic sequences of all known different RV-A and -B prototypes were downloaded from the National Centre for Biotechnology Information (NCBI) and divided into minor low-density lipoprotein receptor (LDLR) and major intercellular adhesion molecule groups (ICAM). The sequences were edited using Biological Sequence Alignment Editor, v 7.2.0 (BioEdit software) to study each capsid protein (VP1, VP2, VP3, and VP4) and analyzed using the EMBL-EBI ClustalW server and the more current Clustal Omega tool for the calculation of the identities and similarities. Results We analyzed and predicted immunogenic motifs from capsid proteins that are conserved across distinct RV serotypes using a bioinformatics technique. The amino acid sequences of VP3 were found to be the most varied, while VP4 was the most conserved protein among all RV-A and RV-B strains. Among all strains studied, RV-74 demonstrated the highest degree of homology to other strains and could be a potential genetic source for recombinant protein production. Nine highly conserved regions with a minimum length of 9-mers were identified, which could serve as potential immune targets against rhinoviruses. Conclusion Therefore, bioinformatics analysis conducted in the current study has paved the way for the selection of immunogenic targets. Bioinformatically, the ideal strain’s capsid protein is suggested to contain the most common RVs immunogenic sites.

N6-methyladenosine (m6A) RNA methylation is an emerging area of epigenetics, which is a reversible and dynamic modification mediating by 'writers' (methylase, adding methyl groups, METTL3, METTL14, and WTAP), 'erasers' (demethylase, deleting methyl groups, FTO and ALKBH5), and 'readers' (YTHDF1-3, YTHDC1 and YTHDC2). Recent studies in human, animal models and cell levels have disclosed a critical role of m6A modification in regulating the homeostasis of metabolic processes and cardiovascular function. Evidence from these studies identify m6A as a candidate of biomarker and therapeutic target for metabolic abnormality and cardiovascular diseases (CVD). Comprehensive understanding of the complexity of m6A regulation in metabolic diseases and CVD will be helpful for us to understand the pathogenesis of CVD. In this review, we discuss the regulatory role of m6A in metabolic abnormality and CVD. We will emphasize the clinical relevance of m6A dysregulation in CVD.

Graphical abstract, Farmers apply broiler chicken litter to soils to enrich organic matter and provide crops with nutrients, following varying periods of stockpiling. However, litter frequently harbors fecal-derived microbial pathogens and associated antibiotic resistance genes (ARGs), and may be a source of microbial contamination of produce. We coupled a cutting-edge Loop Genomics long-read 16S rRNA amplicon-sequencing platform with high-throughput qPCR that targeted a suite of ARGs, to assess temporal (five time points over a 60-day period) and spatial (top, middle and bottom layers) microbiome and resistome dynamics in a broiler litter stockpile. We focused on potentially pathogenic species from the Enterobacteriaceae, Enterococcaceae and Staphylococcaceae families associated with food-borne disease. Bacterial diversity was significantly lower in the middle of the stockpile, where targeted pathogens were lowest and Bacillaceae were abundant. E. coli was the most abundant Enterobacteriaceae species, and high levels of the opportunistic pathogen Enterococcus faecium were detected. Correlation analyses revealed that the latter was significantly associated with aminoglycoside (aac(6′)-Ib(aka aacA4), aadA5), tetracycline (tetG), vancomycin (vanC), phenicol (floR) and MLSB (mphB) resistance genes. Staphylococcaceae were primarily non-pathogenic, but extremely low levels of the opportunistic pathogen S. aureus were detected, as was the opportunistic pathogen S. saprophyticus, which was linked to vancomycin (vanSA, vanC1), MLSB (vatE, ermB) and tetracycline (tetK) resistance genes. Collectively, we found that stockpile microbiomes and resistomes are strongly dictated by temporal fluctuations and spatial heterogeneity. Insights from this study can be exploited to improve stockpile management practice to support sustainable antimicrobial resistance mitigation policies in the future.

COVID-19 is a serious infectious disease that has recently swept the world, and research on its causative virus, SARS-CoV-2, remains insufficient. Therefore, this study uses bioinformatics analysis techniques to explore the human digestive tract diseases that may be caused by SARS-CoV-2 infection. The gene expression profile data set, numbered GSE149312, is from the Gene Expression Omnibus (GEO) database and is divided into a 24-h group and a 60-h group. R software is used to analyze and screen out differentially expressed genes (DEGs) and then gene ontology (GO) term and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses are performed. In KEGG, the pathway of non-alcoholic fatty liver disease exists in both the 24-h group and 60-h group. STRING is used to establish a protein–protein interaction (PPI) network, and Cytoscape is then used to visualize the PPI and define the top 12 genes of the node as the hub genes. Through verification, nine statistically significant hub genes are identified: AKT1, TIMP1, NOTCH, CCNA2, RRM2, TTK, BUB1B, KIF20A, and PLK1. In conclusion, the results of this study can provide a certain direction and basis for follow-up studies of SARS-CoV-2 infection of the human digestive tract and provide new insights for the prevention and treatment of diseases caused by SARS-CoV-2. Supplementary Information The online version contains supplementary material available at 10.1007/s10528-021-10144-w.

The cornerstone of management in patients with acute promyelocytic leukemia (APL) is early diagnosis and prompt initiation of treatment with an all-trans retinoic acid (ATRA)-based regimen. Identification of the t(15;17)(PML-RARA) chromosomal translocation through conventional cytogenetics fluorescence in-situ hybridization (FISH) or detection of the promyelocytic leukemia-retinoic acid receptor alpha (PML-RARα) fusion through RT-PCR represent the current standard of care for diagnosing APL. However, about 1–2% of patients with APL have a variant translocation involving other fusion partners with RARα besides PML. These patients present a unique diagnostic and clinical challenge in that conventional cytogenetics in addition to FISH and/or RT-PCR for PML-RARα may fail to identify these clinically relevant genetic lesions leading to an inappropriate diagnosis and treatment. We present two cases of patients who had APL with variant translocations whose bone marrow specimens were sent to the University of Michigan for enrollment in the MI-ONCOSEQ study (HUM00067928) after standard testing failed to identify PML-RARα or t(15;17) despite a phenotypic concern for this diagnosis. In these two patients, whole exome and transcriptome profiling via the MI-ONCOSEQ platform identified a PRKAR1A-RARα fusion in one patient and ZBTB16-RARα fusion in another patient. These cases illustrate the utility of whole exome and transcriptome profiling in diagnosing variant translocations in patients in whom there is a high clinical suspicion for APL based on hematopathology review.

Obesity is a complex disease with multifactorial causes, and its prevalence is becoming a serious health crisis. For this reason, there is a crucial need to identify novel targets and players. With this aim in mind, we analyzed via RNA-sequencing the subcutaneous adipose tissue of normal weight and obesity-affected women, highlighting the differential expression in the two tissues. We specifically focused on long non-coding RNAs, as 6 of these emerged as dysregulated in the diseased-tissue (COL4A2-AS2, RPS21-AS, PELATON, ITGB2-AS1, ACER2-AS and CTEPHA1). For each of them, we performed both a thorough in silico dissection and in vitro validation, to predict their function during adipogenesis. We report the lncRNAs expression during adipose derived stem cells differentiation to adipocytes as model of adipogenesis and their potential modulation by adipogenesis-related transcription factors (C/EBPs and PPARγ). Moreover, inhibiting CTEPHA1 expression we investigated its impact on adipogenesis-related transcription factors, showing its significative dysregulation of C/EBPα expression. Lastly, we dissected the subcellular localization, pathway involvement and disease-correlation for coding differentially expressed genes. Together, these findings highlight a transcriptional deregulation at the basis of obesity, impacted by both coding and long non-coding RNAs.

Identification of cancer-predisposing germline variants in childhood cancer patients is important for therapeutic decisions, disease surveillance and risk assessment for patients, and potentially, also for family members. We investigated the spectrum and prevalence of pathogenic germline variants in selected childhood cancer patients with features suggestive of genetic predisposition to cancer. Germline DNA was subjected to exome sequencing to filter variants in 1048 genes of interest including 176 known cancer predisposition genes (CPGs). An enrichment burden analysis compared rare deleterious germline CPG variants in the patient cohort with those in a healthy aged control population. A subset of predicted deleterious variants in novel candidate CPGs was investigated further by examining matched tumor samples, and the functional impact of AXIN1 variants was analyzed in cultured cells. Twenty-two pathogenic/likely pathogenic (P/LP) germline variants detected in 13 CPGs were identified in 19 of 76 patients (25.0%). Unclear association with the diagnosed cancer types was observed in 11 of 19 patients carrying P/LP CPG variants. The burden of rare deleterious germline variants in autosomal dominant CPGs was significantly higher in study patients versus healthy aged controls. A novel AXIN1 frameshift variant (Ser321fs) may impact the regulation of β-catenin levels. Selection of childhood cancer patients for germline testing based on features suggestive of an underlying genetic predisposition could help to identify carriers of clinically relevant germline CPG variants, and streamline the integration of germline genomic testing in the pediatric oncology clinic.

SOX10 belongs to a family of 20 SRY (sex-determining region Y)-related high mobility group box-containing (SOX) proteins, most of which contribute to cell type specification and differentiation of various lineages. The first clue that SOX10 is essential for development, especially in the neural crest, came with the discovery that heterozygous mutations occurring within and around SOX10 cause Waardenburg syndrome type 4. Since then, heterozygous mutations have been reported in Waardenburg syndrome type 2 (Waardenburg syndrome type without Hirschsprung disease), PCWH or PCW (peripheral demyelinating neuropathy, central dysmyelination, Waardenburg syndrome, with or without Hirschsprung disease), intestinal manifestations beyond Hirschsprung (ie, chronic intestinal pseudo-obstruction), Kallmann syndrome and cancer. All of these diseases are consistent with the regulatory role of SOX10 in various neural crest derivatives (melanocytes, the enteric nervous system, Schwann cells and olfactory ensheathing cells) and extraneural crest tissues (inner ear, oligodendrocytes). The recent evolution of medical practice in constitutional genetics has led to the identification of SOX10 variants in atypical contexts, such as isolated hearing loss or neurodevelopmental disorders, making them more difficult to classify in the absence of both a typical phenotype and specific expertise. Here, we report novel mutations and review those that have already been published and their functional consequences, along with current understanding of SOX10 function in the affected cell types identified through in vivo and in vitro models. We also discuss research options to increase our understanding of the origin of the observed phenotypic variability and improve the diagnosis and medical care of affected patients.

Vocal learning, the ability to imitate sounds from conspecifics and the environment, is a key component of human spoken language and learned song in three independently evolved avian groups—oscine songbirds, parrots, and hummingbirds. Humans and each of these three bird clades exhibit specialized behavioral, neuroanatomical, and brain gene expression convergence related to vocal learning, speech, and song. To understand the evolutionary basis of vocal learning gene specializations and convergence, we searched for and identified accelerated genomic regions (ARs), a marker of positive selection, specific to vocal learning birds. We found avian vocal learner-specific ARs, and they were enriched in noncoding regions near genes with known speech functions or brain gene expression specializations in humans and vocal learning birds, including FOXP2, NEUROD6, ZEB2, and MEF2C, and near genes with major neurodevelopmental functions, including NR2F1, NRP2, and BCL11B. We also found enrichment near the SFARI class S genes associated with syndromic vocal communication forms of autism spectrum disorders. These findings reveal strong candidate noncoding regions near genes for the evolutionary adaptations that distinguish vocal learning species from their close vocal nonlearning relatives and provide further evidence of molecular convergence between birdsong and human spoken language.

The leucine-rich repeat receptor-like protein kinase (LRR-RLK) gene family is the largest family of the receptor-like protein kinases (RLKs) superfamily in higher plants, which is involved in regulating the plant growth and development, stress responses, signal transduction and so on. However, no comprehensive analyses of LRR-RLKs have been reported in sugarcane. Here, we performed a comprehensive analysis of the LRR-RLK gene family in sugarcane ancestor species Saccharum spontaneum. A total of 437 LRR-RLK genes were identified and categorized into 14 groups based on a maximum likelihood phylogenetic tree. The chromosome location showed an uneven distribution on all 32 chromosomes in sugarcane. Subsequently, the exon–intron organization structure and conserved motif arrangement were relatively conserved among the same groups or subgroups and between Arabidopsis and S. spontaneum genomes. Furthermore, the promoter sequences analyses showed that sugarcane LRR-RLK genes (SsLRR-RLKs) were strongly regulated by various environmental stimuli, phytohormonal factors and transcription factors (TFs). Eventually, the expression profiles of SsLRR-RLK genes at different stresses were analyzed based on RNA-seq data, suggesting their potential roles in the regulation of sugarcane responses to diverse abiotic and biotic stress. Overall, the findings provide insight into the potential functional roles and lay the foundation for further functional study.

Many efforts have been made in recent years to investigate the alterations in protein-coding genes as well as non-coding RNAs that are playing an emerging role in the development and progression of cancers. These miRNAs are short non-coding functional RNAs that are involved in the regulation of transcriptome. In different studies, it was found that human miRNA-149 is an important microRNA that is functioning either as onco-miRNAs or acting as tumor suppressors, in different conditions. Many of the miRNAs are regulating different SNPs of FOXE1 in different studies which are causing low-to-moderate penetrance of genes that initiates the development of thyroid cancer. The involvement of SNPs in miRNA-149 gene rs2292832 and FOXE1 rs3758249 with PTC for better disease prognosis and management was determined in this study and the relation between these SNPs at the genotypic level was also evaluated. PTC patients with age and gender-matched controls were recruited in the present study. Blood samples were collected in EDTA vacutainer followed by DNA extraction by the organic method. Genotyping of rs2292832 and rs3758249 was done by ARMS-PCR and PCR- RFLP respectively. Statistical analyses were carried out by using SPSS software (version 20). The mutation T>C in miRNA-149 rs2292832 was significantly associated with thyroid cancer (p-value 0.0004, C did not show significant association with the disease (p-value 0.124244, > 0.05). Moreover, no correlation of rs2292832 at the genotype level was observed with rs3758249. miRNA-149 gene SNP rs2292832 was observed in strong association with thyroid cancer. Lack of genetic association of rs3758249 of FOXE1 gene has been ruled for the disease. The statistically significant association of rs2292832 with thyroid cancer depicts its mechanistic involvement at the cellular level in Papillary Thyroid Carcinoma.

We report a 38-year-old Saudi male with Ehlers-Danlos Syndrome (EDS). The patient presented with rare and unusual neurological manifestations, including but not limited to ophthalmoplegia and myopathic pattern on his electromyography. In addition to hand weakness, there was skin hyperextensibility, joint hyperflexibility, and frontal baldness. Next-generation sequencing was performed on target exon sequences, using whole exome sequencing and Burrows-Wheeler Aligner for alignment/base calling. Genome Analysis Toolkit and reference genome Homo sapiens (UCSC hg19) were used for sequence processing and analysis. Variant classification was done according to standard international recommendations. A novel homozygous variant, NM_019105.6: c.8488C>T p.(Gln2830*), was detected in the TNXB gene. This variant is not reported in the literature nor dbSNP or gnomAD databases. Additionally, this variant is predicted to create a premature stop codon and produce a truncated protein or nonsense-mediated mRNA decay. Hence, it is classified as a likely pathogenic variant. The same point variant was found in a heterozygous state in the patient's father and sister. Both presented with milder symptoms associated with Ehlers-Danlos syndromes and heritable connective tissue disorders. Therefore, the patient was diagnosed as a tenascin-X (TNX) deficient type of EDS known as classical-like Ehlers-Danlos syndrome. TNX deficient patients may present with clinical and electrophysiological manifestations that are unusual in EDS like frontal baldness, ophthalmoplegia, and myotonia, which mimic myotonic dystrophy type I. Clinicians should be aware of the potential overlap of symptoms among these two diseases to ensure correct diagnosis is made.