Discovery and characterization of human exonic transcriptional regulatory elements

Angular limb deformity associated with TSPAN18, NRG3 and NOVA2 in Rambouillet rams

Angular limb deformity (ALD) affects many species of livestock and companion animals. The mechanisms of ALD development are not well understood, but previous research suggests the involvement of genetic risk factors. A case-control genome-wide association study (GWAS) was conducted with 40 ALD-affected and 302 unaffected Rambouillet rams and 40,945 single nucleotide polymorphisms (SNPs). Forelimbs of 6 ALD-affected rams were examined and diagnosed with osteochondrosis. Genome-wide or chromosome-wide significant SNPs were positioned exonic, intronic or within the 3'UTR of genes TSPAN18, NRG3 and NOVA2, respectively. These genes have previously described roles related to angiogenesis and osteoblast, osteoclast and chondrocyte proliferation and differentiation, which suggests the possibility for their involvement in the pathogenesis of osteochondrosis. Functional consequences of SNPs were evaluated through transcription factor binding site analysis, which predicted binding sites for transcription factors of known importance to bone growth, including SOX6, SOX9 and RUNX2. The identification of genetic risk factors for ALD may help to improve animal welfare and production in Rambouillet, a breed known to be at risk for ALD development. This study proposes genes TSPAN18, NRG3 and NOVA2 as targets for further research towards understanding the etiology of ALD in Rambouillet sheep.

Tubulin Alpha-1b as a Potential Biomarker for Lung Adenocarcinoma Diagnosis and Prognosis

Background: Because lung cancer is the main cause of cancer deaths and lung adenocarcinoma (LUAD) accounts for more than 40% of all lung malignancies, it is essential to develop clinically useful biomarkers for the disease. The aim of this investigation is to assess the potential application of tubulin alpha-1b (TUBA1B) as a biomarker for diagnosing and monitoring the outcome of LUAD. Methods: The clinical data of the LUAD patients was retrospectively analyzed. Immunohistochemistry (IHC) analysis of a tissue microarray containing 90 LUAD cases was implemented to examine the expression of TUBA1B. The protein and mRNA levels of TUBA1B in serum were detected by enzyme-linked immunosorbent assay (ELISA) and quantitative real-time PCR (qRT-PCR) analysis respectively. UALCAN was employed to confirm the expression levels and survival probability of TUBA1B in LUAD patients. Results: Compared to adjacent non-cancerous tissues in the microarray, the expression of TUBA1B in LUAD tissues was much higher. The expression of TUBA1B in LUAD was statistically correlated with lymph node status (P = .031). Moreover, patients with higher TUBA1B expression had shorter overall survival (P < .0001). Furthermore, cox multi-factor analysis also suggested that TUBA1B may be an independent predictor for LUAD prognosis (P = .030). The results of TCGA data analysis by UALCAN were consistent with the microarray results, except for that TUBA1B was also significantly correlated with clinical tumor stages. Protein levels of TUBA1B in serum were obviously elevated in LUAD patients than control (P < .0001), and the area under the ROC curve was 0.99. TUBA1B also showed better sensitivity of 92.9% for LUAD than common clinical biomarkers. Conclusion: TUBA1B may be a non-invasive prognostic and diagnostic biomarker for LUAD patients.

Congenital nephrotic syndrome in a Hispanic Guatemalan newborn associated with a NPHS1 variant: A case report

Congenital nephrotic syndrome (CNS) is an autosomal recessive disorder usually detected in the first 3 months of life when the syndromes effects manifest, including edema and a failure to gain weight. A baby boy was admitted to the Neonatal Intensive Care Unit for prematurity (35 weeks) with unremarkable maternal prenatal laboratory tests. The patient had persistent systemic hypertension, hypoproteinemia, hypoalbuminemia and nephrotic range proteinuria. CNS was diagnosed, and genetic testing showed a homozygous variant, c.3024A>G (AGA>AGG) in exon 22 of the nephrin locus. Bioinformatics analysis suggested the genetic condition was likely a result of malfunctional DNA binding sites of transcription factors FOXL1 and FOXC1.

PKR kinase directly regulates tau expression and Alzheimer's disease-related tau phosphorylation

Deposition of extensively hyperphosphorylated tau in specific brain cells is a clear pathological hallmark in Alzheimer's disease and a number of other neurodegenerative disorders, collectively termed the tauopathies. Furthermore, hyperphosphorylation of tau prevents it from fulfilling its physiological role as a microtubule‐stabilizing protein and leaves it increasingly vulnerable to self‐assembly, suggestive of a central underlying role of hyperphosphorylation as a contributing factor in the etiology of these diseases. Via in vitro phosphorylation and regulation of kinase activity within cells and acute brain tissue, we reveal that the inflammation associated kinase, protein kinase R (PKR), directly phosphorylates numerous abnormal and disease‐modifying residues within tau including Thr181, Ser199/202, Thr231, Ser262, Ser396, Ser404 and Ser409. Similar to disease processes, these PKR‐mediated phosphorylations actively displace tau from microtubules in cells. In addition, PKR overexpression and knockdown, respectively, increase and decrease tau protein and mRNA levels in cells. This regulation occurs independent of noncoding transcriptional elements, suggesting an underlying mechanism involving intra‐exonic regulation of the tau‐encoding microtubule‐associated protein tau (MAPT) gene. Finally, acute encephalopathy in wild type mice, induced by intracranial Langat virus infection, results in robust inflammation and PKR upregulation accompanied by abnormally phosphorylated full‐length‐ and truncated tau. These findings indicate that PKR, independent of other kinases and upon acute brain inflammation, is capable of triggering pathological modulation of tau, which, in turn, might form the initial pathologic seed in several tauopathies such as Alzheimer's disease and Chronic traumatic encephalopathy where inflammation is severe.

Predicting Functional Effects of Synonymous Variants: A Systematic Review and Perspectives

Recent advances in high-throughput experimentation have put the exploration of genome sequences at the forefront of precision medicine. In an effort to interpret the sequencing data, numerous computational methods have been developed for evaluating the effects of genome variants. Interestingly, despite the fact that every person has as many synonymous (sSNV) as non-synonymous single nucleotide variants, our ability to predict their effects is limited. The paucity of experimentally tested sSNV effects appears to be the limiting factor in development of such methods. Here, we summarize the details and evaluate the performance of nine existing computational methods capable of predicting sSNV effects. We used a set of observed and artificially generated variants to approximate large scale performance expectations of these tools. We note that the distribution of these variants across amino acid and codon types suggests purifying evolutionary selection retaining generated variants out of the observed set; i.e., we expect the generated set to be enriched for deleterious variants. Closer inspection of the relationship between the observed variant frequencies and the associated prediction scores identifies predictor-specific scoring thresholds of reliable effect predictions. Notably, across all predictors, the variants scoring above these thresholds were significantly more often generated than observed. which confirms our assumption that the generated set is enriched for deleterious variants. Finally, we find that while the methods differ in their ability to identify severe sSNV effects, no predictor appears capable of definitively recognizing subtle effects of such variants on a large scale.

Precise genome-wide base editing by the CRISPR Nickase system in yeast

The CRISPR/Cas9 system has been applied to efficient genome editing in many eukaryotic cells. However, the bases that can be edited by this system have been limited to those within the protospacer adjacent motif (PAM) and guide RNA-targeting sequences. In this study, we developed a genome-wide base editing technology, "CRISPR Nickase system" that utilizes a single Cas9 nickase. This system was free from the limitation of editable bases that was observed in the CRISPR/Cas9 system, and was able to precisely edit bases up to 53 bp from the nicking site. In addition, this system showed no off-target editing, in contrast to the CRISPR/Cas9 system. Coupling the CRISPR Nickase system with yeast gap repair cloning enabled the construction of yeast mutants within only five days. The CRISPR Nickase system provides a versatile and powerful technology for rapid, site-specific, and precise base editing in yeast.

The Yin and Yang of codon usage

The genetic code is degenerate. With the exception of two amino acids (Met and Trp), all other amino acid residues are each encoded by multiple, so-called synonymous codons. Synonymous codons were initially presumed to have entirely equivalent functions, however, the finding that synonymous codons are not present at equal frequencies in genes/genomes suggested that codon choice might have functional implications beyond amino acid coding. The pattern of non-uniform codon use (known as codon usage bias) varies between organisms and represents a unique feature of an organism. Organism-specific codon choice is related to organism-specific differences in populations of cognate tRNAs. This implies that, in a given organism, frequently used codons will be translated more rapidly than infrequently used ones and vice versa A theory of codon-tRNA co-evolution (necessary to balance accurate and efficient protein production) was put forward to explain the existence of codon usage bias. This model suggests that selection favours preferred (frequent) over un-preferred (rare) codons in order to sustain efficient protein production in cells and that a given un-preferred codon will have the same effect on an organism's fitness regardless of its position within an mRNA's open reading frame. However, many recent studies refute this prediction. Un-preferred codons have been found to have important functional roles and their effects appeared to be position-dependent. Synonymous codon usage affects the efficiency/stringency of mRNA decoding, mRNA biogenesis/stability, and protein secretion and folding. This review summarizes recent developments in the field that have identified novel functions of synonymous codons and their usage.

A systematic, large-scale comparison of transcription factor binding site models

Background

The modelling of gene regulation is a major challenge in biomedical research. This process is dominated by transcription factors (TFs) and mutations in their binding sites (TFBSs) may cause the misregulation of genes, eventually leading to disease. The consequences of DNA variants on TF binding are modelled in silico using binding matrices, but it remains unclear whether these are capable of accurately representing in vivo binding.

In this study, we present a systematic comparison of binding models for 82 human TFs from three freely available sources: JASPAR matrices, HT-SELEX-generated models and matrices derived from protein binding microarrays (PBMs). We determined their ability to detect experimentally verified “real” in vivo TFBSs derived from ENCODE ChIP-seq data. As negative controls we chose random downstream exonic sequences, which are unlikely to harbour TFBS. All models were assessed by receiver operating characteristics (ROC) analysis.

Results

While the area-under-curve was low for most of the tested models with only 47 % reaching a score of 0.7 or higher, we noticed strong differences between the various position-specific scoring matrices with JASPAR and HT-SELEX models showing higher success rates than PBM-derived models. In addition, we found that while TFBS sequences showed a higher degree of conservation than randomly chosen sequences, there was a high variability between individual TFBSs.

Conclusions

Our results show that only few of the matrix-based models used to predict potential TFBS are able to reliably detect experimentally confirmed TFBS.

We compiled our findings in a freely accessible web application called ePOSSUM (http:/mutationtaster.charite.de/ePOSSUM/) which uses a Bayes classifier to assess the impact of genetic alterations on TF binding in user-defined sequences. Additionally, ePOSSUM provides information on the reliability of the prediction using our test set of experimentally confirmed binding sites.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2729-8) contains supplementary material, which is available to authorized users.

Long noncoding RNA expression profiles in gut tissues constitute molecular signatures that reflect the types of microbes

The gut microbiota is commonly referred to as a hidden organ due to its pivotal effects on host physiology, metabolism, nutrition and immunity. The gut microbes may be shaped by environmental and host genetic factors, and previous studies have focused on the roles of protein-coding genes. Here we show a link between long non-coding RNA (lncRNA) expression and gut microbes. By repurposing exon microarrays and comparing the lncRNA expression profiles between germ-free, conventional and different gnotobiotic mice, we revealed subgroups of lncRNAs that were specifically enriched in each condition. A nearest shrunken centroid methodology was applied to obtain lncRNA-based signatures to identify mice in different conditions. The lncRNA-based prediction model successfully identified different gnotobiotic mice from conventional and germ-free mice, and also discriminated mice harboring transplanted microbes from fecal samples of mice or zebra fishes. To achieve optimal prediction accuracy, fewer lncRNAs were required in the prediction model than protein-coding genes. Taken together, our study demonstrated the effecacy of lncRNA expression profiles in discriminating the types of microbes in the gut. These results also provide a resource of gut microbe-associated lncRNAs for the development of lncRNA biomarkers and the identification of functional lncRNAs in host-microbes interactions.

Reassessing the "duon" hypothesis of protein evolution

There are two distinct types of DNA sequences, namely coding sequences and regulatory sequences, in a genome. A recent study of the occupancy of transcription factors (TFs) in human cells suggested that protein-coding sequences also serve as the codes of TF occupancy, and proposed a "duon" hypothesis in which up to 15% of codons of human protein genes are constrained by the additional coding requirements that regulate gene expression. This hypothesis challenges our basic understanding on the human genome. We reanalyzed the data and found that the previous study was confounded by ascertainment bias related to base composition. Using an unbiased comparison in which G/C and A/T sites are considered separately, we reveal a similar level of conservation between TF-bound codons and TF-depleted codons, suggesting largely no extra purifying selection provided by the TF occupancy on the codons of human genes. Given the generally short binding motifs of TFs and the open chromatin structure during transcription, we argue that the occupancy of TFs on protein-coding sequences is mostly passive and evolutionarily neutral, with to-be-determined functions in the regulation of gene expression.

Exonic transcription factor binding directs codon choice and affects protein evolution

Genomes contain both a genetic code specifying amino acids and a regulatory code specifying transcription factor (TF) recognition sequences. We used genomic deoxyribonuclease I footprinting to map nucleotide resolution TF occupancy across the human exome in 81 diverse cell types. We found that ~15% of human codons are dual-use codons ("duons") that simultaneously specify both amino acids and TF recognition sites. Duons are highly conserved and have shaped protein evolution, and TF-imposed constraint appears to be a major driver of codon usage bias. Conversely, the regulatory code has been selectively depleted of TFs that recognize stop codons. More than 17% of single-nucleotide variants within duons directly alter TF binding. Pervasive dual encoding of amino acid and regulatory information appears to be a fundamental feature of genome evolution.

Osmotic stress changes the expression and subcellular localization of the Batten disease protein CLN3

Juvenile CLN3 disease (formerly known as juvenile neuronal ceroid lipofuscinosis) is a fatal childhood neurodegenerative disorder caused by mutations in the CLN3 gene. CLN3 encodes a putative lysosomal transmembrane protein with unknown function. Previous cell culture studies using CLN3-overexpressing vectors and/or anti-CLN3 antibodies with questionable specificity have also localized CLN3 in cellular structures other than lysosomes. Osmoregulation of the mouse Cln3 mRNA level in kidney cells was recently reported. To clarify the subcellular localization of the CLN3 protein and to investigate if human CLN3 expression and localization is affected by osmotic changes we generated a stably transfected BHK (baby hamster kidney) cell line that expresses a moderate level of myc-tagged human CLN3 under the control of the human ubiquitin C promoter. Hyperosmolarity (800 mOsm), achieved by either NaCl/urea or sucrose, dramatically increased the mRNA and protein levels of CLN3 as determined by quantitative real-time PCR and Western blotting. Under isotonic conditions (300 mOsm), human CLN3 was found in a punctate vesicular pattern surrounding the nucleus with prominent Golgi and lysosomal localizations. CLN3-positive early endosomes, late endosomes and cholesterol/sphingolipid-enriched plasma membrane microdomain caveolae were also observed. Increasing the osmolarity of the culture medium to 800 mOsm extended CLN3 distribution away from the perinuclear region and enhanced the lysosomal localization of CLN3. Our results reveal that CLN3 has multiple subcellular localizations within the cell, which, together with its expression, prominently change following osmotic stress. These data suggest that CLN3 is involved in the response and adaptation to cellular stress.