Search

The contactin-associated protein-like 2 (CNTNAP2) gene is a member of the neurexin superfamily. CNTNAP2 was first implicated in the cortical dysplasia-focal epilepsy (CDFE) syndrome, a recessive disease characterized by intellectual disability, epilepsy, language impairments and autistic features. Associated SNPs and heterozygous deletions in CNTNAP2 were subsequently reported in autism, schizophrenia and other psychiatric or neurological disorders. We aimed to comprehensively examine evidence for the role of CNTNAP2 in susceptibility to psychiatric disorders, by the analysis of multiple classes of genetic variation in large genomic datasets. In this study we used: i) summary statistics from the Psychiatric Genomics Consortium (PGC) GWAS for seven psychiatric disorders; ii) examined all reported CNTNAP2 structural variants in patients and controls; iii) performed cross-disorder analysis of functional or previously associated SNPs; and iv) conducted burden tests for pathogenic rare variants using sequencing data (4,483 ASD and 6,135 schizophrenia cases, and 13,042 controls). The distribution of CNVs across CNTNAP2 in psychiatric cases from previous reports was no different from controls of the database of genomic variants. Gene-based association testing did not implicate common variants in autism, schizophrenia or other psychiatric phenotypes. The association of proposed functional SNPs rs7794745 and rs2710102, reported to influence brain connectivity, was not replicated; nor did predicted functional SNPs yield significant results in meta-analysis across psychiatric disorders at either SNP-level or gene-level. Disrupting CNTNAP2 rare variant burden was not higher in autism or schizophrenia compared to controls. Finally, in a CNV mircroarray study of an extended bipolar disorder family with 5 affected relatives we previously identified a 131kb deletion in CNTNAP2 intron 1, removing a FOXP2 transcription factor binding site. Quantitative-PCR validation and segregation analysis of this CNV revealed imperfect segregation with BD. This large comprehensive study indicates that CNTNAP2 may not be a robust risk gene for psychiatric phenotypes. [ABSTRACT FROM AUTHOR]

Transcriptomic profiling of complex tissues by single-nucleus RNA-sequencing (snRNA-seq) affords some advantages over single-cell RNA-sequencing (scRNA-seq). snRNA-seq provides less biased cellular coverage, does not appear to suffer cell isolation-based transcriptional artifacts, and can be applied to archived frozen specimens. We used well-matched snRNA-seq and scRNA-seq datasets from mouse visual cortex to compare cell type detection. Although more transcripts are detected in individual whole cells (~11,000 genes) than nuclei (~7,000 genes), we demonstrate that closely related neuronal cell types can be similarly discriminated with both methods if intronic sequences are included in snRNA-seq analysis. We estimate that the nuclear proportion of total cellular mRNA varies from 20% to over 50% for large and small pyramidal neurons, respectively. Together, these results illustrate the high information content of nuclear RNA for characterization of cellular diversity in brain tissues. [ABSTRACT FROM AUTHOR]

The species-rich diatom family Chaetocerotaceae is common in the coastal marine phytoplankton worldwide where it is responsible for a substantial part of the primary production. Despite its relevance for the global cycling of carbon and silica, many species are still described only morphologically, and numerous specimens do not fit any described taxa. Nowadays, studies to assess plankton biodiversity deploy high throughput sequencing metabarcoding of the 18S rDNA V4 region, but to translate the gathered metabarcodes into biologically meaningful taxa, there is a need for reference barcodes. However, 18S reference barcodes for this important family are still relatively scarce. We provide 18S rDNA and partial 28S rDNA reference sequences of 443 morphologically characterized chaetocerotacean strains. We gathered 164 of the 216 18S sequences and 244 of the 413 28S sequences of strains from the Gulf of Naples, Atlantic France, and Chile. Inferred phylogenies showed 84 terminal taxa in seven principal clades. Two of these clades included terminal taxa whose rDNA sequences contained spliceosomal and Group IC1 introns. Regarding the commonly used metabarcode markers in planktonic diversity studies, all terminal taxa can be discriminated with the 18S V4 hypervariable region; its primers fit their targets in all but two species, and the V4-tree topology is similar to that of the 18S. Hence V4-metabarcodes of unknown Chaetocerotaceae are assignable to the family. Regarding the V9 hypervariable region, most terminal taxa can be discriminated, but several contain introns in their primer targets. Moreover, poor phylogenetic resolution of the V9 region affects placement of metabarcodes of putative but unknown chaetocerotacean taxa, and hence, uncertainty in taxonomic assignment, even of higher taxa. [ABSTRACT FROM AUTHOR]

Glutathione S-Transferases (GSTs) are phase II detoxification enzymes that may have evolved in response to changes of environmental substrates. GST genes formed a multigene family and in mammals, there are six classes known as Alpha, Mu, Omega, Pi, Theta, and Zeta. Recent studies in phase I detoxification system specifically the cytochrome P450s provided a general explanation on why genes from a common origin such as those in a multigene family have both phylogenetically stable and unstable genes. Genes that participate in core functions of organisms such as development and physiology are stable whereas genes that play a role in detoxification are unstable and evolve in a process known as birth-death evolution, which is characterised by frequent gene gains and losses. The generality of the birth-death model at explaining the evolution of detoxification enzymes beyond the phase I enzyme has not been comprehensively explored. This work utilized 383 Gst genes and 300 pseudogenes across 22 mammalian species to study gene gains and losses. GSTs vary greatly in their phylogenetic stability despite their overall sequence similarity. Stable Gst genes from Omega and Zeta classes do not show fluctuation in gene numbers from human to opossum. These genes play a role in biosynthesis related functions. Unstable genes that include Alpha, Mu, Pi and Theta undergo frequent gene gain and loss in a process known as birth-death evolution. Gene members of these four classes are well known for their roles in detoxification. Our positive selection screen identified five positively selected sites in mouse GSTA3. Previous studies showed two of these sites (108H and 208E) were biochemically tested as important residues that conferred catalytic activity against the toxic aflatoxin B1-8,9-epoxide. The functional significance against aflatoxin of the remaining three positively selected sites warrant further investigation. [ABSTRACT FROM AUTHOR]

Copy number variants (CNVs) are a pervasive source of genetic variation and evolutionary potential, but the dynamics and diversity of CNVs within evolving populations remain unclear. Long-term evolution experiments in chemostats provide an ideal system for studying the molecular processes underlying CNV formation and the temporal dynamics with which they are generated, selected, and maintained. Here, we developed a fluorescent CNV reporter to detect de novo gene amplifications and deletions in individual cells. We used the CNV reporter in Saccharomyces cerevisiae to study CNV formation at the GAP1 locus, which encodes the general amino acid permease, in different nutrient-limited chemostat conditions. We find that under strong selection, GAP1 CNVs are repeatedly generated and selected during the early stages of adaptive evolution, resulting in predictable dynamics. Molecular characterization of CNV-containing lineages shows that the CNV reporter detects different classes of CNVs, including aneuploidies, nonreciprocal translocations, tandem duplications, and complex CNVs. Despite GAP1’s proximity to repeat sequences that facilitate intrachromosomal recombination, breakpoint analysis revealed that short inverted repeat sequences mediate formation of at least 50% of GAP1 CNVs. Inverted repeat sequences are also found at breakpoints at the DUR3 locus, where CNVs are selected in urea-limited chemostats. Analysis of 28 CNV breakpoints indicates that inverted repeats are typically 8 nucleotides in length and separated by 40 bases. The features of these CNVs are consistent with origin-dependent inverted-repeat amplification (ODIRA), suggesting that replication-based mechanisms of CNV formation may be a common source of gene amplification. We combined the CNV reporter with barcode lineage tracking and found that 102–104 independent CNV-containing lineages initially compete within populations, resulting in extreme clonal interference. However, only a small number (18–21) of CNV lineages ever constitute more than 1% of the CNV subpopulation, and as selection progresses, the diversity of CNV lineages declines. Our study introduces a novel means of studying CNVs in heterogeneous cell populations and provides insight into their dynamics, diversity, and formation mechanisms in the context of adaptive evolution. [ABSTRACT FROM AUTHOR]

As whole exome sequencing (WES) becomes more widely used in the clinical realm, a wealth of unanalyzed information will be routinely generated. Using WES read depth data to predict copy number variation (CNV) could extend the diagnostic utility of this previously underutilized data by providing clinically important information such as previously unsuspected deletions or duplications. We evaluated ExomeDepth, a free R package, in addition to an aneuploidy prediction method, to detect CNVs in WES data. First, in a blinded pilot study, five out of five genomic alterations were correctly identified from clinical samples with previously defined chromosomal gains or losses, including submicroscopic deletions, duplications, and chromosomal trisomy. We then examined CNV calls among 53 patients participating in the NCGENES research study and undergoing WES, who had existing clinical chromosomal microarray (CMA) data that could be used for validation. For unique CNVs that overlap well with WES coverage regions, sensitivity was 89% for deletions and 65% for duplications. While specificity of the algorithm calls remains a concern, this is less of an issue at high threshold filtering levels. When applied to all 672 patients from the exome sequencing study, ExomeDepth identified eleven diagnostically relevant CNVs ranging in size from a two exon deletion to whole chromosome duplications, as well as numerous other CNVs with varying clinical significance. This opportunistic analysis of WES data yields an additional 1.6% of patients in this study with pathogenic or likely pathogenic CNVs that are clinically relevant to their phenotype as well as clinically relevant secondary findings. Finally, we demonstrate the potential value of copy number analysis in cases where a single heterozygous likely or known pathogenic single nucleotide alteration is identified in a gene associated with an autosomal recessive condition. [ABSTRACT FROM AUTHOR]

KSHV is endemic in Uganda and the HIV epidemic has dramatically increased the incidence of Kaposi sarcoma (KS). To investigate the role of KSHV in the development of KS, we obtained KS biopsies from ART-naïve, HIV-positive individuals in Uganda and analyzed the tumors using RNAseq to globally characterize the KSHV transcriptome. Phylogenetic analysis of ORF75 sequences from 23 tumors revealed 6 distinct genetic clusters with KSHV strains exhibiting M, N or P alleles. RNA reads mapping to specific unique coding sequence (UCDS) features were quantitated using a gene feature file previously developed to globally analyze and quantitate KSHV transcription in infected endothelial cells. A pattern of high level expression was detected in the KSHV latency region that was common to all KS tumors. The clear majority of transcription was derived from the downstream latency transcript promoter P3(LTd) flanking ORF72, with little evidence of transcription from the P1(LTc) latency promoter, which is constitutive in KSHV-infected lymphomas and tissue-culture cells. RNAseq data provided evidence of alternate P3(LTd) transcript editing, splicing and termination resulting in multiple gene products, with 90% of the P3(LTd) transcripts spliced to release the intronic source of the microRNAs K1-9 and 11. The spliced transcripts encode a regulatory uORF upstream of Kaposin A with alterations in intervening repeat sequences yielding novel or deleted Kaposin B/C-like sequences. Hierarchical clustering and PCA analysis of KSHV transcripts revealed three clusters of tumors with different latent and lytic gene expression profiles. Paradoxically, tumors with a latent phenotype had high levels of total KSHV transcription, while tumors with a lytic phenotype had low levels of total KSHV transcription. Morphologically distinct KS tumors from the same individual showed similar KSHV gene expression profiles suggesting that the tumor microenvironment and host response play important roles in the activation level of KSHV within the infected tumor cells. [ABSTRACT FROM AUTHOR]

Mitochondrial DNA (mtDNA) haplogroups (hgs) are evolutionarily conserved sets of mtDNA SNP-haplotypes with characteristic geographical distribution. Associations of hgs with disease and physiological characteristics have been reported, but have frequently not been reproducible. Using 418 mtDNA SNPs on the PsychChip (Illumina), we assessed the spatio-temporal distribution of mtDNA hgs in Denmark from DNA isolated from 24,642 geographically un-biased dried blood spots (DBS), collected from 1981 to 2005 through the Danish National Neonatal Screening program. ADMIXTURE was used to establish the genomic ancestry of all samples using a reference of 100K+ autosomal SNPs in 2,248 individuals from nine populations. Median-joining analysis determined that the hgs were highly variable, despite being typically Northern European in origin, suggesting multiple founder events. Furthermore, considerable heterogeneity and variation in nuclear genomic ancestry was observed. Thus, individuals with hg H exhibited 95%, and U hgs 38.2% - 92.5%, Danish ancestry. Significant clines between geographical regions and rural and metropolitan populations were found. Over 25 years, macro-hg L increased from 0.2% to 1.2% (p = 1.1*E-10), and M from 1% to 2.4% (p = 3.7*E-8). Hg U increased among the R macro-hg from 14.1% to 16.5% (p = 1.9*E-3). Genomic ancestry, geographical skewedness, and sub-hg distribution suggested that the L, M and U increases are due to immigration. The complex spatio-temporal dynamics and genomic ancestry of mtDNA in the Danish population reflect repeated migratory events and, in later years, net immigration. Such complexity may explain the often contradictory and population-specific reports of mito-genomic association with disease. [ABSTRACT FROM AUTHOR]

Increasing β-carotene (a vitamin A precursor) content in Triticum turgidum L. ssp. durum (durum wheat) grains is important to improve pasta nutritional quality. Studies in other species show that altering the expression of LCYE genes increases the flux towards the β-β branch, accumulating higher β-carotene levels. Durum wheat is a tetraploid species that has two LCYE genes (LCYE-A and LCYE-B) associated to the A and B genomes. The objective of this work was to produce durum wheat LCYE mutants through EMS to potentially increase β-carotene content. The LCYE point mutations created with EMS were identified using a Kronos TILLING (Targeting Induced Local Lesion IN Genomes) mutant population. Specific primers that amplified exons 3 through 10 of the LCYE genes were designed and validated. To simplify the TILLING procedure, fragments were digested with CJE (Celery Juice Extract) and visualized on 2% agarose gels. 6X mutant pools were identified, which showed cleavage products and then made into 2X pools to identify mutant individuals. LCYE mutants were then sequenced and evaluated with BLOSUM62, SIFT and PSSM algorithms. Mutants with substitutions W437*, P334L and G368R in LCYE-A and P405L, G352R and T393I in LCYE-B predicted to affect protein function were selected. Substitution W437* increased β-carotene in 75% and overall total carotenoids content in leaves of the mutant 2426 (A1 mutant line), but no significant differences relative to the control were found in grains through HPLC. Finally, the increased levels of β-carotene on leaves have potential applications to improving plant resistance under contaminated environmental conditions. [ABSTRACT FROM AUTHOR]

Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong genetic component and etiology characterized by chronic inflammation and autoantibody production. The purpose of this study was to ascertain copy number variation (CNV) in SLE using a case-control design in an admixed Brazilian population. The whole-genome detection of CNV was performed using Cytoscan HD array in SLE patients and healthy controls. The best CNV candidates were then evaluated by quantitative real-time PCR in a larger cohort or validated using droplet digital PCR. Logistic regression models adjusted for sex and ancestry covariates was applied to evaluate the association between CNV with SLE susceptibility. The data showed a synergistic effect between the FCGR3B and ADAM3A loci with the presence of deletions in both loci significantly increasing the risk to SLE (5.9-fold) compared to the deletion in the single FCGR3B locus (3.6-fold). In addition, duplications in these genes were indeed more frequent in healthy subjects, suggesting that high FCGR3B/ADAM3A gene copy numbers are protective factors against to disease development. Overall, 21 rare CNVs were identified in SLE patients using a four-step pipeline created for identification of rare variants. Furthermore, heterozygous deletions overlapping the CFHR4, CFHR5 and HLA-DPB2 genes were described for the first time in SLE patients. Here we present the first genome-wide CNV study of SLE patients in a tri-hybrid population. The results show that novel susceptibility loci to SLE can be found once the distribution of structural variants is analyzed throughout the whole genome. [ABSTRACT FROM AUTHOR]

Group II introns are ancient retroelements that significantly shaped the origin and evolution of contemporary eukaryotic genomes. These self-splicing ribozymes share a common ancestor with the telomerase enzyme, the spliceosome machinery as well as the highly abundant spliceosomal introns and non-LTR retroelements. More than half of the human genome thus consists of various elements that evolved from ancient group II introns, which altogether significantly contribute to key functions and genetic diversity in eukaryotes. Similarly, group II intron-related elements in bacteria such as abortive phage infection (Abi) retroelements, diversity generating retroelements (DGRs) and some CRISPR-Cas systems have evolved to confer important functions to their hosts. In sharp contrast, since bacterial group II introns are scarce, irregularly distributed and frequently spread by lateral transfer, they have mainly been considered as selfish retromobile elements with no beneficial function to their host. Here we unveil a new group II intron function that generates genetic diversity at the RNA level in bacterial cells. We demonstrate that Ll.LtrB, the model group II intron from Lactococcus lactis, recognizes specific sequence motifs within cellular mRNAs by base pairing, and invades them by reverse splicing. Subsequent splicing of ectopically inserted Ll.LtrB, through circularization, induces a novel trans-splicing pathway that generates exon 1-mRNA and mRNA-mRNA intergenic chimeras. Our data also show that recognition of upstream alternative circularization sites on intron-interrupted mRNAs release Ll.LtrB circles harboring mRNA fragments of various lengths at their splice junction. Intergenic trans-splicing and alternative circularization both produce novel group II intron splicing products with potential new functions. Overall, this work describes new splicing pathways in bacteria that generate, similarly to the spliceosome in eukaryotes, genetic diversity at the RNA level while providing additional functional and evolutionary links between group II introns, spliceosomal introns and the spliceosome. [ABSTRACT FROM AUTHOR]

Previous studies have identified multiple loci for inherited susceptibility to glioma development, including the regulator of telomere elongation helicase 1 (RTEL1). However, the association between RTEL1 variants and risk of glioma has not been well understood. Therefore, we sought to comprehensively examine the genetic interaction between RTEL1 variants and risk of glioma with respect to defined histological and molecular subtypes. We employed a case-control study involving 250 adult glioma patients with previous molecular alterations and 375 population–based controls within Korean populations. Statistical analyses on the association between RTEL1 single nucleotide polymorphisms (SNPs) and glioma risk were conducted using unconditional logistic regression. Additional conditional and stepwise analyses were performed on significant RTEL1 SNPs. We detected significant associations (Bonferroni P < .05) between six SNPs (rs6089953, rs3848669, rs6010620, rs3787089, rs6062302, and rs115303435) and risk of glioma in the Korean subjects. The two coding variants, rs6062302 (D664D) and rs115303435 (A1059T), were plausibly causal variants and were independent among the significantly associated RTEL1 variants. The glioma subgroup analyses showed that the causal variants (rs6062302 and rs115303435) may be associated with increased risk of glioma regardless of histological grades and molecular alterations. This study provides a deeper understanding of relationships between RTEL1 variants and risk of glioma. Further studies are required to ascertain the impact of those variants on glioma susceptibility. [ABSTRACT FROM AUTHOR]

CFTR modulators have revolutionized the treatment of individuals with cystic fibrosis (CF) by improving the function of existing protein. Unfortunately, almost half of the disease-causing variants in CFTR are predicted to introduce premature termination codons (PTC) thereby causing absence of full-length CFTR protein. We hypothesized that a subset of nonsense and frameshift variants in CFTR allow expression of truncated protein that might respond to FDA-approved CFTR modulators. To address this concept, we selected 26 PTC-generating variants from four regions of CFTR and determined their consequences on CFTR mRNA, protein and function using intron-containing minigenes expressed in 3 cell lines (HEK293, MDCK and CFBE41o-) and patient-derived conditionally reprogrammed primary nasal epithelial cells. The PTC-generating variants fell into five groups based on RNA and protein effects. Group A (reduced mRNA, immature (core glycosylated) protein, function <1% (n = 5)) and Group B (normal mRNA, immature protein, function <1% (n = 10)) variants were unresponsive to modulator treatment. However, Group C (normal mRNA, mature (fully glycosylated) protein, function >1% (n = 5)), Group D (reduced mRNA, mature protein, function >1% (n = 5)) and Group E (aberrant RNA splicing, mature protein, function > 1% (n = 1)) variants responded to modulators. Increasing mRNA level by inhibition of NMD led to a significant amplification of modulator effect upon a Group D variant while response of a Group A variant was unaltered. Our work shows that PTC-generating variants should not be generalized as genetic ‘nulls’ as some may allow generation of protein that can be targeted to achieve clinical benefit. [ABSTRACT FROM AUTHOR]

Holy basil (Ocimum sanctum L.) and sweet basil (Ocimum basilicum L.) are the most commonly grown basil species in India for essential oil production and biosynthesis of potentially volatile and non-volatile phytomolecules with commercial significance. The aroma, flavor and pharmaceutical value of Ocimum species is a significance of its essential oil, which contains most of the monoterpenes and sesquiterpenes. A large number of plants have been studied for characterization and identification of terpene synthase genes, involved in terpenoids biosynthesis. The goal of this study is to discover and identify the putative functional terpene synthase genes in O. sanctum. HMMER search was performed by using a set of 13 well sequenced and annotated plant genomes including the newly sequenced genome of O. sanctum with Pfam-A database locally, using HMMER 3.0 hmmsearch for the two Pfam domains (PF01397 and PF03936). Using this search method 81 putative terpene synthases genes (OsaTPS) were identified in O. sanctum; the study further reveals 47 OsaTPS were putatively functional genes, 19 partial OsaTPS, and 15 OsaTPS as probably pseudogenes. All these identified OsaTPS genes were compared with other plant species, and phylogenetic analysis reveals the subfamily classification of OsaTPS in TPS-a, -b, -c, -e, -f and TPS-g subfamilies clusters. This genome-wide identification of OsaTPS genes, their phylogenetic analysis and secondary metabolite pathway mapping predictions together provide a comprehensive understanding of the TPS gene family in Ocimum sanctum and offer opportunities for the characterization and functional validation of numbers of terpene synthase genes. [ABSTRACT FROM AUTHOR]

Many eukaryotic genes play essential roles in multiple biological processes in several different tissues. Conditional mutants are needed to analyze genes with such pleiotropic functions. In vertebrates, conditional gene inactivation has only been feasible in the mouse, leaving other model systems to rely on surrogate experimental approaches such as overexpression of dominant negative proteins and antisense-based tools. Here, we have developed a simple and straightforward method to integrate loxP sequences at specific sites in the zebrafish genome using the CRISPR/Cas9 technology and oligonucleotide templates for homology directed repair. We engineered conditional (floxed) mutants of tbx20 and fleer, and demonstrate excision of exons flanked by loxP sites using tamoxifen-inducible CreERT2 recombinase. To demonstrate broad applicability of our method, we also integrated loxP sites into two additional genes, aldh1a2 and tcf21. The ease of this approach will further expand the use of zebrafish to study various aspects of vertebrate biology, especially post-embryonic processes such as regeneration. [ABSTRACT FROM AUTHOR]

Root-knot nematodes (RKNs, Meloidogyne spp.) are destructive plant parasites with a wide host range. They severely reduce crop quality and yield worldwide. Tobacco is a versatile model plant organism for studying RKNs-host interactions and a key plant material for molecular research. Long noncoding RNAs (lncRNAs) play critical roles in post transcriptional and transcriptional regulation in a wide range of biological pathways, especially plant development and stress response. In the present study, we obtained 5,206 high-confidence lncRNAs based on RNA sequencing data. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses revealed that the target genes of these lncRNAs are mainly involved in plant biotic and abiotic stresses, plant hormone signal transduction, induced systemic resistance, plant-type hypersensitive response, plant-type cell wall organization or biogenesis. The 565 differentially expressed lncRNAs found to be involved in nematode stress response were validated by quantitative PCR using 15 randomly-selected lncRNA genes. Our study provides insights into the molecular mechanisms of RKNs-plant interactions that might help preventing nematode damages to crops. [ABSTRACT FROM AUTHOR]

Clustered copy number variants (CNVs) as detected by chromosomal microarray analysis (CMA) are often reported as germline chromothripsis. However, such cases might need further investigations by massive parallel whole genome sequencing (WGS) in order to accurately define the underlying complex rearrangement, predict the occurrence mechanisms and identify additional complexities. Here, we utilized WGS to delineate the rearrangement structure of 21 clustered CNV carriers first investigated by CMA and identified a total of 83 breakpoint junctions (BPJs). The rearrangements were further sub-classified depending on the patterns observed: I) Cases with only deletions (n = 8) often had additional structural rearrangements, such as insertions and inversions typical to chromothripsis; II) cases with only duplications (n = 7) or III) combinations of deletions and duplications (n = 6) demonstrated mostly interspersed duplications and BPJs enriched with microhomology. In two cases the rearrangement mutational signatures indicated both a breakage-fusion-bridge cycle process and haltered formation of a ring chromosome. Finally, we observed two cases with Alu- and LINE-mediated rearrangements as well as two unrelated individuals with seemingly identical clustered CNVs on 2p25.3, possibly a rare European founder rearrangement. In conclusion, through detailed characterization of the derivative chromosomes we show that multiple mechanisms are likely involved in the formation of clustered CNVs and add further evidence for chromoanagenesis mechanisms in both “simple” and highly complex chromosomal rearrangements. Finally, WGS characterization adds positional information, important for a correct clinical interpretation and deciphering mechanisms involved in the formation of these rearrangements. [ABSTRACT FROM AUTHOR]

Coniothyrium glycines, the causal agent of soybean red leaf blotch, is a USDA APHIS-listed Plant Pathogen Select Agent and potential threat to US agriculture. Sequencing of the C. glycines mt genome revealed a circular 98,533-bp molecule with a mean GC content of 29.01%. It contains twelve of the mitochondrial genes typically involved in oxidative phosphorylation (atp6, cob, cox1-3, nad1-6, and nad4L), one for a ribosomal protein (rps3), four for hypothetical proteins, one for each of the small and large subunit ribosomal RNAs (rns and rnl) and a set of 30 tRNAs. Genes were encoded on both DNA strands with cox1 and cox2 occurring as adjacent genes having no intergenic spacers. Likewise, nad2 and nad3 are adjacent with no intergenic spacers and nad5 is immediately followed by nad4L with an overlap of one base. Thirty-two introns, comprising 54.1% of the total mt genome, were identified within eight protein-coding genes and the rnl. Eighteen of the introns contained putative intronic ORFs with either LAGLIDADG or GIY-YIG homing endonuclease motifs, and an additional eleven introns showed evidence of truncated or degenerate endonuclease motifs. One intron possessed a degenerate N-acetyl-transferase domain. C. glycines shares some conservation of gene order with other members of the Pleosporales, most notably nad6-rnl-atp6 and associated conserved tRNA clusters. Phylogenetic analysis of the twelve shared protein coding genes agrees with commonly accepted fungal taxonomy. C. glycines represents the second largest mt genome from a member of the Pleosporales sequenced to date. This research provides the first genomic information on C. glycines, which may provide targets for rapid diagnostic assays and population studies. [ABSTRACT FROM AUTHOR]

A nonsense variant (p.W358X) of human pancreatic lipase related protein 2 (PNLIPRP2) is present in different ethnic populations with a high allele frequency. In cell culture experiments, the truncated protein mainly accumulates inside the cells and causes endoplasmic reticulum stress. Here, we tested the hypothesis that variant p.W358X might increase risk for chronic pancreatitis through acinar cell stress. We sequenced exon 11 of PNLIPRP2 in a cohort of 256 subjects with chronic pancreatitis (152 alcoholic and 104 non-alcoholic) and 200 controls of Hungarian origin. We observed no significant difference in the distribution of the truncation variant between patients and controls. We analyzed mRNA expression in human pancreatic cDNA samples and found the variant allele markedly reduced. We conclude that the p.W358X truncation variant of PNLIPRP2 is expressed poorly and has no significant effect on the risk of chronic pancreatitis. [ABSTRACT FROM AUTHOR]

Objective: To evaluate the association of enterovirus 71 (EV71) susceptibility and clinical severity with polymorphisms in EV71 receptors, including human scavenger receptor class B member 2 (SCARB2), P-selectin glycoprotein ligand-1 (PSGL-1), and annexin II (ANXA2). Methods: We enrolled laboratory-confirmed EV71 cases and healthy age- and gender-matched controls in Taiwan from 2000 to 2012. We detected genetic polymorphisms in SCARB2, PSGL-1, and ANXA2 and correlated the results with EV71 susceptibility and severity. Results: We collected 599 EV71 cases and 98 controls. Among EV71 patients, the male to female ratio was 1.61, and the mean age was 2.99±2.47 years. For clinical severity, 117 (19.6%) had severe central nervous system involvement with or without cardiopulmonary failure. For outcomes, 46 (7.7%) had sequelae, and 14 (2.3%) died. SCARB2 polymorphisms (rs6824953 and rs11097262) were associated with susceptibility to EV71 infection (OR 1.60, 95% CI 1.07–2.39; and OR 1.64, 95% CI 1.09–2.47, respectively). PSGL-1 polymorphisms (rs7137098 and rs8179137) were significantly associated with severe EV71 infection (OR 1.46, 95% CI 1.1–1.96; and OR 1.47, 95% CI 1.07–2.03, respectively). Conclusions: SCARB2 polymorphisms (rs6824953 and rs11097262) might be associated with EV71 susceptibility. PSGL-1 polymorphisms (rs7137098 and rs8179137) were associated with severe EV71 infection. [ABSTRACT FROM AUTHOR]

The cotton bollworm, Helicoverpa armigera (Hübner) is one of the most serious insect pest species to evolve resistance against many insecticides from different chemical classes. This species has evolved resistance to the pyrethroid insecticides across its native range and is becoming a truly global pest after establishing in South America and having been recently recorded in North America. A chimeric cytochrome P450 gene, CYP337B3, has been identified as a resistance mechanism for resistance to fenvalerate and cypermethrin. Here we show that this resistance mechanism is common around the world with at least eight different alleles. It is present in South America and has probably introgressed into its closely related native sibling species, Helicoverpa zea. The different alleles of CYP337B3 are likely to have arisen independently in different geographic locations from selection on existing diversity. The alleles found in Brazil are those most commonly found in Asia, suggesting a potential origin for the incursion of H. armigera into the Americas. [ABSTRACT FROM AUTHOR]

Background: Pseudogenes are non-functional sequences in the genome with homologous sequences that are functional (i.e. genes). They are abundant in eukaryotes where they have been extensively investigated, while in prokaryotes they are significantly scarcer and less well studied. Here we conduct a comprehensive analysis of the evolution of orthologs of Mycobacterium leprae pseudogenes in prokaryotes. The leprosy pathogen M. leprae is of particular interest since it contains an unusually large number of pseudogenes, comprising approximately 40% of its entire genome. The analysis is conducted in both broad and narrow phylogenetic ranges. Results: We have developed an informatics-based approach to characterize the evolution of pseudogenes. This approach combines tools from phylogenomics, genomics, and transcriptomics. The results we obtain are used to assess the contributions of two mechanisms for pseudogene formation: failed horizontal gene transfer events and disruption of native genes. Conclusions: We conclude that, although it was reported that in most bacteria the former is most likely responsible for the majority of pseudogenization events, in mycobacteria, and in particular in M. leprae with its exceptionally high pseudogene numbers, the latter predominates. We believe that our study sheds new light on the evolution of pseudogenes in bacteria, by utilizing new methodologies that are applied to the unusually abundant M. leprae pseudogenes and their orthologs. [ABSTRACT FROM AUTHOR]

During lytic replication of Kaposi’s sarcoma-associated herpesvirus (KSHV), a nuclear viral long noncoding RNA known as PAN RNA becomes the most abundant polyadenylated transcript in the cell. Knockout or knockdown of KSHV PAN RNA results in loss of late lytic viral gene expression and, consequently, reduction of progeny virion release from the cell. Here, we demonstrate that knockdown of PAN RNA from the related Rhesus macaque rhadinovirus (RRV) phenocopies that of KSHV PAN RNA. These two PAN RNA homologs, although lacking significant nucleotide sequence conservation, can functionally substitute for each other to rescue phenotypes associated with the absence of PAN RNA expression. Because PAN RNA is exclusively nuclear, previous studies suggested that it directly interacts with host and viral chromatin to modulate gene expression. We studied KSHV and RRV PAN RNA homologs using capture hybridization analysis of RNA targets (CHART) and observed their associations with host chromatin, but the loci differ between PAN RNA homologs. Accordingly, we find that KSHV PAN RNA is undetectable in chromatin following cell fractionation. Thus, modulation of gene expression at specific chromatin loci appears not to be the primary, nor the pertinent function of this viral long noncoding RNA. PAN RNA represents a cautionary tale for the investigation of RNA association with chromatin whereby cross-linking of DNA spatially adjacent to an abundant nuclear RNA gives the appearance of specific interactions. Similarly, PAN RNA expression does not affect viral transcription factor complex expression or activity, which is required for generation of the late lytic viral mRNAs. Rather, we provide evidence for an alternative model of PAN RNA function whereby knockdown of KSHV or RRV PAN RNA results in compromised nuclear mRNA export thereby reducing the cytoplasmic levels of viral mRNAs available for production of late lytic viral protein. [ABSTRACT FROM AUTHOR]

Nuclear genes of euglenids and marine diplonemids harbor atypical, nonconventional introns which are not observed in the genomes of other eukaryotes. Nonconventional introns do not have the conserved borders characteristic for spliceosomal introns or the sequence complementary to U1 snRNA at the 5' end. They form a stable secondary structure bringing together both exon/intron junctions, nevertheless, this conformation does not resemble the form of self-splicing or tRNA introns. In the genes studied so far, frequent nonconventional introns insertions at new positions have been observed, whereas conventional introns have been either found at the conserved positions, or simply lost. In this work, we examined the order of intron removal from Euglena gracilis transcripts of the tubA and gapC genes, which contain two types of introns: nonconventional and spliceosomal. The relative order of intron excision was compared for pairs of introns belonging to different types. Furthermore, intermediate products of splicing were analyzed using the PacBio Next Generation Sequencing system. The analysis led to the main conclusion that nonconventional introns are removed in a rapid way but later than spliceosomal introns. Moreover, the observed accumulation of transcripts with conventional introns removed and nonconventional present may suggest the existence of a time gap between the two types of splicing. [ABSTRACT FROM AUTHOR]

The modern cattle was domesticated from aurochs, sharing its physiological traits into two subspecies Bos taurus and Bos indicus. MicroRNAs (miRNAs) are a class of non-coding short RNAs of ~22nt which have a key role in the regulation of many cellular and physiological processes in the animal. The current study was aimed to predict and annotate the potential mutations in indicine miRNAs throughout the genome using de novo and homology-based in silico approaches. Genome-wide mapping was performed in available indicine assembly by the homology-based approach and 768 miRNAs were recovered out of 808 reported taurine miRNAs belonging to 521 unique mature miRNA families. While 42 precursors were dropped due to lack of secondary miRNA structure, increasing stringency or decreasing similarity between the two genomes’ miRNA. Increasing tendency of miRNAs incidence was observed on chr5, chr7, chr8, chr12 and chr21 with 19 polycistronic miRNA within 1-kilobase distance throughout the indicine genome. Notably, 12 miRNAs showed copy number variation. Eighteen miRNAs showed a mutation in their mature sequences in which eight were found in their seed region. Whilst in de novo based approach, 12 novel potential miRNAs on Y chromosome in indicine cattle along with a new miRNA (bind-miR-1264) on chrX were found. The final data set is annotated and explains the impending target genes that are responsible for enhanced immunity, heat tolerance and disease tolerance regulation in indicine. The study conforms to better understanding and perceptive approach towards indicine genome. [ABSTRACT FROM AUTHOR]

Mutations in interleukin-10 receptor (IL-10R) genes are one cause of very early-onset inflammatory bowel disease with perianal lesions, which can be cured by hematopoietic stem cell transplantation. Using a functional test, which assesses responsiveness of peripheral monocytes to IL-10, we identified three unrelated Portuguese patients carrying two novel IL-10RB mutations. In the three patients, sequencing of genomic DNA identified the same large deletion of exon 3 which precluded protein expression. This mutation was homozygous in two patients born from consanguineous families and heterozygous in the third patient born from unrelated parents. Microsatellite analysis of the IL10RB genomic region revealed a common haplotype in the three Portuguese families pointing to a founder deletion inherited from a common ancestor 400 years ago. In the third patient, surface expression of IL-10R was normal but signaling in response to IL-10 was impaired. Complementary DNA sequencing and next-generation sequencing of IL10RB locus with custom-made probes revealed a ≈ 6 Kb duplication encompassing the exon 6 which leads to a frameshift mutation and a loss of the TYK2-interacting Box 2 motif. Altogether, we describe two novel copy number variations in IL10RB, one with founder effect and one preserving cell surface expression but abolishing signaling. [ABSTRACT FROM AUTHOR]

The interaction between Plasmodium vivax Duffy binding protein (PvDBP) and Duffy antigen receptor for chemokines (DARC) has been described as critical for the invasion of human reticulocytes, although increasing reports of P. vivax infections in Duffy-negative individuals questions its unique role. To investigate the genetic diversity of the two main protein ligands for reticulocyte invasion, PvDBP and P. vivax Erythrocyte Binding Protein (PvEBP), we analyzed 458 isolates collected in Cambodia and Madagascar from individuals genotyped as Duffy-positive. First, we observed a high proportion of isolates with multiple copies PvEBP from Madagascar (56%) where Duffy negative and positive individuals coexist compared to Cambodia (19%) where Duffy-negative population is virtually absent. Whether the gene amplification observed is responsible for alternate invasion pathways remains to be tested. Second, we found that the PvEBP gene was less diverse than PvDBP gene (12 vs. 33 alleles) but provided evidence for an excess of nonsynonymous mutations with the complete absence of synonymous mutations. This finding reveals that PvEBP is under strong diversifying selection, and confirms the importance of this protein ligand in the invasion process of the human reticulocytes and as a target of acquired immunity. These observations highlight how genomic changes in parasite ligands improve the fitness of P. vivax isolates in the face of immune pressure and receptor polymorphisms. [ABSTRACT FROM AUTHOR]

Elucidating the complete network of protein-protein interactions, or interactome, is a fundamental goal of the post-genomic era, yet existing interactome maps are far from complete. To increase the throughput and resolution of interactome mapping, methods for protein-protein interaction discovery by co-migration have been introduced. However, accurate identification of interacting protein pairs within the resulting large-scale proteomic datasets is challenging. Consequently, most computational pipelines for co-migration data analysis incorporate external genomic datasets to distinguish interacting from non-interacting protein pairs. The effect of this procedure on interactome mapping is poorly understood. Here, we conduct a rigorous analysis of genomic data integration for interactome recovery across a large number of co-migration datasets, spanning diverse experimental and computational methods. We find that genomic data integration leads to an increase in the functional coherence of the resulting interactome maps, but this comes at the expense of a decrease in power to discover novel interactions. Importantly, putative novel interactions predicted by genomic data integration are no more likely to later be experimentally discovered than those predicted from co-migration data alone. Our results reveal a widespread and unappreciated limitation in a methodology that has been widely used to map the interactome of humans and model organisms. [ABSTRACT FROM AUTHOR]

Endothelial nitric oxide synthase (eNOS) gene polymorphisms are associated with reduced eNOS activity and nitric oxide (NO) production leading to an increase in blood pressure (BP). Regular exercise is the main strategy to minimize the deleterious effects of polymorphisms. However, due to the differences that physical exercise can be performed, some controversial results are found. Therefore it seems reasonable to evaluate the training status (TS). Thus, this study aimed to investigate the association of eNOS gene haplotypes and different levels of TS on nitrite concentrations (NO2-) and BP values in older adult. 424 elderly performed the following assessments: General Functional Fitness Index (GFFI) to estimate TS, systolic and diastolic blood pressure (SBP and DBP), blood collection for analysis of NO2- and g.-786T>C, intron 4b/a (VNTR) and 894G>T polymorphisms. Multivariate logistic regression showed that NO2- was influenced by GFFI and 4b/4a Intron 4. Regarding BP, GFFI influenced SBP and DBP, and just intron 4 was associated with variations in DBP. It can be observed that GFFI affected the NO2-, SBP and DBP independently of haplotypes. Therefore, maintenance of good level of TS can overcome the negative influence of genetics factors (intron 4) by increasing NO2- concentration and decreasing BP values. [ABSTRACT FROM AUTHOR]

The IGFN1 (Immunoglobulin-Like And Fibronectin Type III Domain Containing 1) gene has a role in skeletal muscle function and is also involved in metastatic breast cancer, and the isoforms with three N-terminal globular domains are sufficient for its function in skeletal muscle. Two novel splicing isoforms of IGFN1 have been identified in renal cell carcinoma (RCC), one with 5’exon extension and an isoform with a novel exon. The role of G-quadruplex, a non-B DNA, was explored for the splicing alteration of IGFN1 in RCC. G-quadruplexes are the secondary structures acquired by stacking of G-quartets by Hoogsteen hydrogen bonding in DNA and RNA. IGFN1 has intronic potential G-quadruplex forming sequence (PQS) folding into G-quadruplex and is studied for its involvement in aberrant splicing. A PQS in the intron 15 of IGFN1 gene was observed in our in silico analysis by QGRS mapper and non BdB web servers. We observed PQS folds into stable G-quadruplex structure in gel shift assay and circular dichroism (CD) spectroscopy in the presence of G-quadruplex stabilizing agents Pyridostatin (PDS) and KCl, respectively. G-quadruplex formation site with single base resolution was mapped by Sanger sequencing of the plasmid constructs harbouring the cloned PQS and its mutant. This stable G-quadruplex inhibits reverse transcriptase and taq polymerase in reverse transcriptase & PCR stop assays. PDS changes the different splicing isoforms of IGFN1 in UOK146 cell line, displaying involvement of intronic G-quadruplex in IGFN1 splicing. These results lead us to propose that a stable G-quadruplex structure is formed in IGFN1 intron and a reason behind IGFN1 aberrant splicing which could be targeted for therapeutic intervention. [ABSTRACT FROM AUTHOR]

In plant breeding the use of molecular markers has resulted in tremendous improvement of the speed with which new crop varieties are introduced into the market. Single Nucleotide Polymorphism (SNP) genotyping is routinely used for association studies, Linkage Disequilibrium (LD) and Quantitative Trait Locus (QTL) mapping studies, marker-assisted backcrosses and validation of large numbers of novel SNPs. Here we present the KeyGene SNPSelect technology, a scalable and flexible multiplexed, targeted sequence-based, genotyping solution. The multiplex composition of SNPSelect assays can be easily changed between experiments by adding or removing loci, demonstrating their content flexibility. To demonstrate this versatility, we first designed a 1,056-plex maize assay and genotyped a total of 374 samples originating from an F2 and a Recombinant Inbred Line (RIL) population and a maize germplasm collection. Next, subsets of the most informative SNP loci were assembled in 384-plex and 768-plex assays for further genotyping. Indeed, selection of the most informative SNPs allows cost-efficient yet highly informative genotyping in a custom-made fashion, with average call rates between 88.1% (1,056-plex assay) and 99.4% (384-plex assay), and average reproducibility rates between duplicate samples ranging from 98.2% (1056-plex assay) to 99.9% (384-plex assay). The SNPSelect workflow can be completed from a DNA sample to a genotype dataset in less than three days. We propose SNPSelect as an attractive and competitive genotyping solution to meet the targeted genotyping needs in fields such as plant breeding. [ABSTRACT FROM AUTHOR]

Almost 60 years ago, Severo Ochoa was awarded the Nobel Prize in Physiology or Medicine for his discovery of the enzymatic synthesis of RNA by polynucleotide phosphorylase (PNPase). Although this discovery provided an important tool for deciphering the genetic code, subsequent work revealed that the predominant function of PNPase in bacteria and eukaryotes is catalyzing the reverse reaction, i.e., the release of ribonucleotides from RNA. PNPase has a crucial role in RNA metabolism in bacteria and eukaryotes mainly through its roles in processing and degrading RNAs, but additional functions in RNA metabolism have recently been reported for this enzyme. Here, we discuss these established and noncanonical functions for PNPase and the possibility that the major impact of PNPase on cell physiology is through its unorthodox roles. [ABSTRACT FROM AUTHOR]

Long noncoding RNAs constitute a major fraction of the eukaryotic transcriptome, and together with proteins, they intricately fine-tune various growth regulatory signals to control cellular homeostasis. Here, we describe the functional characterisation of a novel pair of long intergenic noncoding RNAs (lincRNAs) comprised of complementary, fully overlapping sense and antisense transcripts Genomic Instability Inducing RNA (Ginir) and antisense RNA of Ginir (Giniras), respectively, from mouse cells. This transcript pair is expressed in a spatiotemporal manner during embryonic development. The individual levels of the sense and antisense transcripts are finely balanced during embryonic growth and in adult tissues. Functional studies of the individual transcripts performed using overexpression and knock-down strategies in mouse cells has led to the discovery that Ginir RNA is a regulator of cellular proliferation and can act as an oncogene having a preeminent role in malignant transformation. Mechanistically, we demonstrate that the oncogenic function of Ginir is mediated by its interaction with centrosomal protein 112 (Cep112). Additionally, we establish here a specific interaction between Cep112 with breast cancer type 1 susceptibility protein (Brca1), another centrosome-associated protein. Next, we prove that the mutual interaction between Cep112 with Brca1 is significant for mitotic regulation and maintenance of genomic stability. Furthermore, we demonstrate that the Cep112 protein interaction with Brca1 protein is impaired when an elevated level of Ginir RNA is present in the cells, resulting in severe deregulation and abnormality in mitosis, leading to malignant transformation. Inhibiting the Ginir RNA function in transformed cells attenuates transformation and restores genomic stability. Together, these findings unravel, to our knowledge, a hitherto-unknown mechanism of oncogenesis mediated by a long noncoding RNA and establishes a unique role of Cep112–Brca1 interaction being modulated by Ginir RNA in maintaining mitotic fidelity. [ABSTRACT FROM AUTHOR]

The Prenatal Assessment of Genome and Exomes (PAGE) project is a UK-wide study aiming to gain a better understanding of genetic variants causing developmental problems during pregnancy. A further aim of the study is to provide an evidence-base for the introduction of prenatal whole genome and exome sequencing (PWGES) into prenatal diagnostics provided by the NHS, which is expected in 2018. This paper presents the findings of a qualitative interview study undertaken with 20 health professionals and researchers involved in the PAGE project, and explores their implications for understandings of ‘good practice’ in the uses of prenatal genomics clinically. A number of critical issues are identified that will need to be addressed in the development of a model of good ethical practice for prenatal genomics: consent, management of expectations, return of results, and professional duties in the context of PWGES. The analysis presented identifies and illustrates a great deal of complexity and qualitative richness in these issues as they arise in the day-to-day work of genomics professionals. Inclusive, critical discussion of these findings, together with the findings from other empirical studies, normative analysis and scientific discoveries resulting from PAGE, will be required to inform the development of appropriate guidelines of good ethical practice that address the needs and concerns to be encountered in daily clinical practice. [ABSTRACT FROM AUTHOR]

Sexually deceptive orchids of the genus Ophrys may rapidly evolve by adaptation to pollinators. However, understanding of the genetic basis of potential changes and patterns of relationships is hampered by a lack of genomic information. We report the complete plastid genome sequences of Ophrys iricolor and O. sphegodes, representing the two most species-rich lineages of the genus Ophrys. Both plastomes are circular DNA molecules (146754 bp for O. sphegodes and 150177 bp for O. iricolor) with the typical quadripartite structure of plastid genomes and within the average size of photosynthetic orchids. 213 Simple Sequence Repeats (SSRs) (31.5% polymorphic between O. iricolor and O. sphegodes) were identified, with homopolymers and dipolymers as the most common repeat types. SSRs were mainly located in intergenic regions but SSRs located in coding regions were also found, mainly in ycf1 and rpoC2 genes. The Ophrys plastome is predicted to encode 107 distinct genes, 17 of which are completely duplicated in the Inverted Repeat regions. 83 and 87 putative RNA editing sites were detected in 25 plastid genes of the two Ophrys species, all occurring in the first or second codon position. Comparing the rate of nonsynonymous (dN) and synonymous (dS) substitutions, 24 genes (including rbcL and ycf1) display signature consistent with positive selection. When compared with other members of the orchid family, the Ophrys plastome has a complete set of 11 functional ndh plastid genes, with the exception of O. sphegodes that has a truncated ndhF gene. Comparative analysis showed a large co-linearity with other related Orchidinae. However, in contrast to O. iricolor and other Orchidinae, O. sphegodes has a shift of the junction between the Inverted Repeat and Small Single Copy regions associated with the loss of the partial duplicated gene ycf1 and the truncation of the ndhF gene. Data on relative genomic coverage and validation by PCR indicate the presence, with a different ratio, of the two plastome types (i.e. with and without ndhF deletion) in both Ophrys species, with a predominance of the deleted type in O. sphegodes. A search for this deleted plastid region in O. sphegodes nuclear genome shows that the deleted region is inserted in a retrotransposon nuclear sequence. The present study provides useful genomic tools for studying conservation and patterns of relationships of this rapidly radiating orchid genus. [ABSTRACT FROM AUTHOR]

Herpes simplex virus 1 (HSV-1) latency establishment is tightly controlled by promyelocytic leukemia (PML) nuclear bodies (NBs) (or ND10), although their exact contribution is still elusive. A hallmark of HSV-1 latency is the interaction between latent viral genomes and PML NBs, leading to the formation of viral DNA-containing PML NBs (vDCP NBs), and the complete silencing of HSV-1. Using a replication-defective HSV-1-infected human primary fibroblast model reproducing the formation of vDCP NBs, combined with an immuno-FISH approach developed to detect latent/quiescent HSV-1, we show that vDCP NBs contain both histone H3.3 and its chaperone complexes, i.e., DAXX/ATRX and HIRA complex (HIRA, UBN1, CABIN1, and ASF1a). HIRA also co-localizes with vDCP NBs present in trigeminal ganglia (TG) neurons from HSV-1-infected wild type mice. ChIP and Re-ChIP show that vDCP NBs-associated latent/quiescent viral genomes are chromatinized almost exclusively with H3.3 modified on its lysine (K) 9 by trimethylation, consistent with an interaction of the H3.3 chaperones with multiple viral loci and with the transcriptional silencing of HSV-1. Only simultaneous inactivation of both H3.3 chaperone complexes has a significant impact on the deposition of H3.3 on viral genomes, suggesting a compensation mechanism. In contrast, the sole depletion of PML significantly impacts the chromatinization of the latent/quiescent viral genomes with H3.3 without any overall replacement with H3.1. vDCP NBs-associated HSV-1 genomes are not definitively silenced since the destabilization of vDCP NBs by ICP0, which is essential for HSV-1 reactivation in vivo, allows the recovery of a transcriptional lytic program and the replication of viral genomes. Consequently, the present study demonstrates a specific chromatin regulation of vDCP NBs-associated latent/quiescent HSV-1 through an H3.3-dependent HSV-1 chromatinization involving the two H3.3 chaperones DAXX/ATRX and HIRA complexes. Additionally, the study reveals that PML NBs are major actors in latent/quiescent HSV-1 H3.3 chromatinization through a PML NB/histone H3.3/H3.3 chaperone axis. [ABSTRACT FROM AUTHOR]

Background: The present study was conducted to discover genetic imbalances such as DNA copy number variations (CNVs) associated with gastric cancer (GC) and to examine their association with different genes involved in the process of gastric carcinogenesis in Saudi population. Methods: Formalin-fixed paraffin-embedded (FFPE) tissues samples from 33 gastric cancer patients and 15 normal gastric samples were collected. Early and late stages GC samples were genotyped and CNVs were assessed by using Illumina HumanOmni1-Quad v.1.0 BeadChip. Results: Copy number gains were more frequent than losses throughout all GC samples compared to normal tissue samples. The mean number of the altered chromosome per case was 64 for gains and 40 for losses, and the median aberration length was 679115bp for gains and 375889bp for losses. We identified 7 high copy gain, 52 gains, 14 losses, 32 homozygous losses, and 10 copy neutral LOHs (loss of heterozygosities). Copy number gains were frequently detected at 1p36.32, 1q12, 1q22, 2p11.1, 4q23-q25, 5p12-p11, 6p21.33, 9q12-q21.11, 12q11-q12, 14q32.33, 16p13.3, 17p13.1, 17q25.3, 19q13.32, and losses at 1p36.23, 1p36.32, 1p32.1, 1q44, 3q25.2, 6p22.1, 6p21.33, 8p11.22, 10q22.1, 12p11.22, 14q32.12 and 16q24.2. We also identified 2 monosomy at chromosome 14 and 22, 52 partially trisomy and 22 whole chromosome 4 neutral loss of heterozygosities at 13q14.2-q21.33, 5p15.2-p15.1, 5q11.2-q13.2, 5q33.1-q34 and 3p14.2-q13.12. Furthermore, 11 gains and 2 losses at 1p36.32 were detected for 11 different GC samples and this region has not been reported before in other populations. Statistical analysis confirms significant association of H. pylori infection with T4 stage of GC as compare to control and other stages. Conclusions: We found that high frequency of copy number gains and losses at 1p36.23, 1p32.1, 1p36.32, 3q25.2, 6p21.33 and 16q24.2 may be common events in gastric cancer. While novel CNVs at 1p36.32 harbouring PRDM16, TP73 and TP73-AS1 genes showed 11 gains and 2 losses for 11 different GC cases and this region is not reported yet in Database of Genomic Variants may be specific to Saudi population. [ABSTRACT FROM AUTHOR]

Several modern genomic technologies, such as DNA-Methylation arrays, measure spatially registered probes that number in the hundreds of thousands across multiple chromosomes. The measured probes are by themselves less interesting scientifically; instead scientists seek to discover biologically interpretable genomic regions comprised of contiguous groups of probes which may act as biomarkers of disease or serve as a dimension-reducing pre-processing step for downstream analyses. In this paper, we introduce an unsupervised feature learning technique which maps technological units (probes) to biological units (genomic regions) that are common across all subjects. We use ideas from fusion penalties and convex clustering to introduce a method for Spatial Convex Clustering, or SpaCC. Our method is specifically tailored to detecting multi-subject regions of methylation, but we also test our approach on the well-studied problem of detecting segments of copy number variation. We formulate our method as a convex optimization problem, develop a massively parallelizable algorithm to find its solution, and introduce automated approaches for handling missing values and determining tuning parameters. Through simulation studies based on real methylation and copy number variation data, we show that SpaCC exhibits significant performance gains relative to existing methods. Finally, we illustrate SpaCC’s advantages as a pre-processing technique that reduces large-scale genomics data into a smaller number of genomic regions through several cancer epigenetics case studies on subtype discovery, network estimation, and epigenetic-wide association. [ABSTRACT FROM AUTHOR]

Genomic imprinting is the phenomena that leads to silencing of one copy of a gene inherited from a specific parent. Mutations in imprinted regions have been involved in diseases showing parent of origin effects. Identifying genes with evidence of parent of origin expression patterns in family studies allows the detection of more subtle imprinting. Here, we use allele specific expression in lymphoblastoid cell lines from 306 Hutterites related in a single pedigree to provide formal evidence for parent of origin effects. We take advantage of phased genotype data to assign parent of origin to RNA-seq reads in individuals with gene expression data. Our approach identified known imprinted genes, two putative novel imprinted genes, PXDC1 and PWAR6, and 14 genes with asymmetrical parent of origin gene expression. We used gene expression in peripheral blood leukocytes (PBL) to validate our findings, and then confirmed imprinting control regions (ICRs) using DNA methylation levels in the PBLs. [ABSTRACT FROM AUTHOR]

Similar to many diurnal animals, Drosophila melanogaster exhibits a mid-day siesta that is more robust as ambient temperature rises, an adaptive response aimed at minimizing exposure to heat. Mid-day siesta levels are partly regulated by the thermosensitive splicing of a small intron (termed dmpi8) found in the 3’ untranslated region (UTR) of the circadian clock gene period (per). Using the well-studied D. melanogaster latitudinal cline along the eastern coast of Australia, we show that flies from temperate populations sleep less during the day compared to those from tropical regions. We identified combinations of four single nucleotide polymorphisms (SNPs) in the 3’ UTR of per that yield several different haplotypes. The two most abundant of these haplotypes exhibit a reciprocal tropical-temperate distribution in relative frequency. Intriguingly, transgenic flies with the major tropical isoform manifest increased daytime sleep and reduced dmpi8 splicing compared to those carrying the temperate variant. Our results strongly suggest that for a major portion of D. melanogaster in Australia, thermal adaptation of daily sleep behavior included spatially varying selection on ancestrally derived polymorphisms in the per 3’ UTR that differentially control dmpi8 splicing efficiency. Prior work showed that African flies from high altitudes manifest reduced mid-day siesta levels, indicative of parallel latitudinal and altitudinal adaptation across continents. However, geographical variation in per 3’ UTR haplotypes was not observed for African flies, providing a compelling case for inter-continental variation in factors targeted by natural selection in attaining a parallel adaptation. We propose that the ability to calibrate mid-day siesta levels to better match local temperature ranges is a key adaptation contributing to the successful colonization of D. melanogaster beyond its ancestral range in the lowlands of Sub-Saharan Africa. [ABSTRACT FROM AUTHOR]

E3 ubiquitin ligases of the ubiquitin proteasome system (UPS) mediate recognition of substrates and later transfer the ubiquitin (Ub). They are the most expanded components of the system. The Really Interesting New Gene (RING) domain contains 40–60 residues that are highly represented among E3 ubiquitin ligases. The Arabidopsis thaliana E3 ubiquitin ligases with a RING finger primarily contain RING-HC or RING-H2 type domains or less frequently RING-v, RING-C2, RING-D, RING-S/T and RING-G type domains. Our previous work on three E3 ubiquitin ligase families with a RING-H2 type domain, ATL, BTL, and CTL, suggested that a phylogenetic distribution based on the RING domain allowed for the creation a catalog of known domains or unknown conserved motifs. This work provided a useful and comprehensive view of particular families of RING E3 ubiquitin ligases. We updated the annotation of A. thaliana RING proteins and surveyed RING proteins from 30 species across eukaryotes. Based on domain architecture profile of the A. thaliana proteins, we catalogued 4711 RING finger proteins into 107 groups, including 66 previously described gene families or single genes and 36 novel families or undescribed genes. Forty-four groups were specific to a plant lineage while 41 groups consisted of proteins found in all eukaryotic species. Our present study updates the current classification of plant RING finger proteins and reiterates the importance of these proteins in plant growth and adaptation. [ABSTRACT FROM AUTHOR]

Introduction: Pseudogenes are paralogues of functional genes historically viewed as defunct due to either the lack of regulatory elements or the presence of frameshift mutations. Recent evidence, however, suggests that pseudogenes may regulate gene expression, although the functional role of pseudogenes remains largely unknown. We previously reported that MYLKP1, the pseudogene of MYLK that encodes myosin light chain kinase (MLCK), is highly expressed in lung and colon cancer cell lines and tissues but not in normal lung or colon. The MYLKP1 promoter is minimally active in normal bronchial epithelial cells but highly active in lung adenocarcinoma cells. In this study, we further validate MYLKP1 as an oncogene via elucidation of the functional role of MYLKP1 genetic variants in colon cancer risk. Methods: Proliferation and migration assays were performed in MYLKP1-transfected colon and lung cancer cell lines (H441, A549) and commercially-available normal lung and colon cells. Fourteen MYLKP1 SNPs (MAFs >0.01) residing within the 4 kb MYLKP1 promoter region, the core 1.4 kb of MYLKP1 gene, and a 4 kb enhancer region were selected and genotyped in a colorectal cancer cohort. MYLKP1 SNP influences on activity of MYLKP1 promoter (2kb) was assessed by dual luciferase reporter assay. Results: Cancer cell lines, H441 and A549, exhibited increased MYLKP1 expression, increased MYLKP1 luciferase promoter activity, increased proliferation and migration. Genotyping studies identified two MYLKP1 SNPs (rs12490683; rs12497343) that significantly increase risk of colon cancer in African Americans compared to African American controls. Rs12490683 and rs12497343 further increase MYLKP1 promoter activity compared to the wild type MYLKP1 promoter. Conclusion: MYLKP1 is a cancer-promoting pseudogene whose genetic variants differentially enhance cancer risk in African American populations. [ABSTRACT FROM AUTHOR]

Messenger RNA polyadenylation is a universal aspect of gene expression in eukaryotes. In well-established model organisms, this process is mediated by a conserved complex of 15–20 subunits. To better understand this process in apicomplexans, a group of unicellular parasites that causes serious disease in humans and livestock, a computational and high throughput sequencing study of the polyadenylation complex and poly(A) sites in several species was conducted. BLAST-based searches for orthologs of the human polyadenylation complex yielded clear matches to only two—poly(A) polymerase and CPSF73—of the 19 proteins used as queries in this analysis. As the human subunits that recognize the AAUAAA polyadenylation signal (PAS) were not immediately obvious, a computational analysis of sequences adjacent to experimentally-determined apicomplexan poly(A) sites was conducted. The results of this study showed that there exists in apicomplexans an A-rich region that corresponds in position to the AAUAAA PAS. The set of experimentally-determined sites in one species, Sarcocystis neurona, was further analyzed to evaluate the extent and significance of alternative poly(A) site choice in this organism. The results showed that almost 80% of S. neurona genes possess more than one poly(A) site, and that more than 780 sites showed differential usage in the two developmental stages–extracellular merozoites and intracellular schizonts–studied. These sites affected more than 450 genes, and included a disproportionate number of genes that encode membrane transporters and ribosomal proteins. Taken together, these results reveal that apicomplexan species seem to possess a poly(A) signal analogous to AAUAAA even though genes that may encode obvious counterparts of the AAUAAA-recognizing proteins are absent in these organisms. They also indicate that, as is the case in other eukaryotes, alternative polyadenylation is a widespread phenomenon in S. neurona that has the potential to impact growth and development. [ABSTRACT FROM AUTHOR]

Recursive splicing (RS) is an evolutionarily conserved process of removing long introns via multiple steps of splicing. It was first discovered in Drosophila and recently proven to occur also in humans. The detailed mechanism of recursive splicing is not well understood, in particular, whether it is kinetically coupled with transcription. To investigate the dynamic process that underlies recursive splicing, we systematically characterized 342 RS sites in three human cell types using published time-series data that monitored synchronized Pol II elongation and nascent RNA production with 4-thiouridine labeling. We found that half of the RS events occurred post-transcriptionally with long delays. For at least 18–47% RS introns, we detected RS junction reads only after detecting canonical splicing junction reads, supporting the notion that these introns were removed by both recursive splicing and canonical splicing. Furthermore, the choice of which splicing mechanism was used showed cell type specificity. Our results suggest that recursive splicing supplements, rather than replaces, canonical splicing for removing long introns. [ABSTRACT FROM AUTHOR]

Recursive splicing, a process by which a single intron is removed from pre-mRNA transcripts in multiple distinct segments, has been observed in a small subset of Drosophila melanogaster introns. However, detection of recursive splicing requires observation of splicing intermediates that are inherently unstable, making it difficult to study. Here we developed new computational approaches to identify recursively spliced introns and applied them, in combination with existing methods, to nascent RNA sequencing data from Drosophila S2 cells. These approaches identified hundreds of novel sites of recursive splicing, expanding the catalog of recursively spliced fly introns by 4-fold. A subset of recursive sites were validated by RT-PCR and sequencing. Recursive sites occur in most very long (> 40 kb) fly introns, including many genes involved in morphogenesis and development, and tend to occur near the midpoints of introns. Suggesting a possible function for recursive splicing, we observe that fly introns with recursive sites are spliced more accurately than comparably sized non-recursive introns. [ABSTRACT FROM AUTHOR]

Schizophrenia and the affective disorders, here comprising bipolar disorder and major depressive disorder, are psychiatric illnesses that lead to significant morbidity and mortality worldwide. Whilst understanding of their pathobiology remains limited, large case-control studies have recently identified single nucleotide polymorphisms (SNPs) associated with these disorders. However, discerning the functional effects of these SNPs has been difficult as the associated causal genes are unknown. Here we evaluated whether schizophrenia and affective disorder associated-SNPs are correlated with gene expression within human brain tissue. Specifically, to identify expression quantitative trait loci (eQTLs), we leveraged disorder-associated SNPs identified from 11 genome-wide association studies with gene expression levels in post-mortem, neurologically-normal tissue from two independent human brain tissue expression datasets (UK Brain Expression Consortium (UKBEC) and Genotype-Tissue Expression (GTEx)). Utilizing stringent multi-region meta-analyses, we identified 2,224 cis-eQTLs associated with expression of 40 genes, including 11 non-coding RNAs. One cis-eQTL, rs16969968, results in a functionally disruptive missense mutation in CHRNA5, a schizophrenia-implicated gene. Importantly, comparing across tissues, we find that blood eQTLs capture < 10% of brain cis-eQTLs. Contrastingly, > 30% of brain-associated eQTLs are significant in tibial nerve. This study identifies putatively causal genes whose expression in region-specific tissue may contribute to the risk of schizophrenia and affective disorders. [ABSTRACT FROM AUTHOR]

Sophora japonica L. (, Leguminosae) is an important traditional Chinese herb with a long history of cultivation. Its flower buds and fruits contain abundant flavonoids, and therefore, the plants are cultivated for the industrial extraction of rutin. Here, we determined the complete nucleotide sequence of the mitochondrial genome of S. japonica ‘JinhuaiJ2’, the most widely planted variety in Guangxi region of China. The total length of the mtDNA sequence is 484,916 bp, with a GC content of 45.4%. Sophora japonica mtDNA harbors 32 known protein-coding genes, 17 tRNA genes, and three rRNA genes with 17 cis-spliced and five trans-spliced introns disrupting eight protein-coding genes. The gene coding and intron regions, and intergenic spacers account for 7.5%, 5.8% and 86.7% of the genome, respectively. The gene profile of S. japonica mitogenome differs from that of the other Faboideae species by only one or two gene gains or losses. Four of the 17 cis-spliced introns showed distinct length variations in the Faboideae, which could be attributed to the homologous recombination of the short repeats measuring a few bases located precisely at the edges of the putative deletions. This reflects the importance of small repeats in the sequence evolution in Faboideae mitogenomes. Repeated sequences of S. japonica mitogenome are mainly composed of small repeats, with only 20 medium-sized repeats, and one large repeat, adding up to 4% of its mitogenome length. Among the 25 pseudogene fragments detected in the intergenic spacer regions, the two largest ones and their corresponding functional gene copies located in two different sets of medium-sized repeats, point to their origins from homologous recombinations. As we further observed the recombined reads associated with the longest repeats of 2,160 bp with the PacBio long read data set of just 15 × in depth, repeat mediated homologous recombinations may play important role in the mitogenomic evolution of S. japonica. Our study provides insightful knowledge to the genetic background of this important herb species and the mitogenomic evolution in the Faboideae species. [ABSTRACT FROM AUTHOR]

We present ggsashimi, a command-line tool for the visualization of splicing events across multiple samples. Given a specified genomic region, ggsashimi creates sashimi plots for individual RNA-seq experiments as well as aggregated plots for groups of experiments, a feature unique to this software. Compared to the existing versions of programs generating sashimi plots, it uses popular bioinformatics file formats, it is annotation-independent, and allows the visualization of splicing events even for large genomic regions by scaling down the genomic segments between splice sites. ggsashimi is freely available at . It is implemented in python, and internally generates R code for plotting. [ABSTRACT FROM AUTHOR]

Background: Pathophysiological mechanisms are still incompletely understood for leprosy, an urgent public health issue in Brazil. Complement receptor 1 (CR1) binds complement fragments C3b/C4b deposited on mycobacteria, mediating its entrance in macrophages. We investigated CR1 polymorphisms, gene expression and soluble CR1 levels in a case-control study with Brazilian leprosy patients, aiming to understand the role of this receptor in differential susceptibility to the disease. Methodology: Nine polymorphisms were haplotyped by multiplex PCR-SSP in 213 leprosy patients (47% multibacillary) and 297 controls. mRNA levels were measured by qPCR and sCR1 by ELISA, in up to 80 samples. Principal findings: Individuals with the most common recombinant haplotype harboring rs3849266*T in intron 21 and rs3737002*T in exon 26 (encoding p.1408Met of the York Yka+ antigen), presented twice higher susceptibility to leprosy (OR = 2.43, p = 0.017). Paucibacillary patients with these variants presented lower sCR1 levels, thus reducing the anti-inflammatory response (p = 0.040 and p = 0.046, respectively). Furthermore, the most ancient haplotype increased susceptibility to the multibacillary clinical form (OR = 3.04, p = 0.01) and presented the intronic rs12034383*G allele, which was associated with higher gene expression (p = 0.043), probably increasing internalization of the parasite. Furthermore, there was an inverse correlation between the the levels of sCR1 and and mannose-binding lectin (initiator molecule of the lectin pathway of complement, recognized by CR1) (R = -0.52, p = 0.007). Conclusions: The results lead us to suggest a regulatory role for CR1 polymorphisms on mRNA and sCR1 levels, with haplotype-specific effects increasing susceptibility to leprosy, probably by enhancing parasite phagocytosis and inflammation. [ABSTRACT FROM AUTHOR]