Revisiting rustrela virus - new cases of encephalitis and a solution to the capsid enigma

Rustrela virus (RusV, species Rubivirus strelense) is a recently discovered relative of rubella virus (RuV) that has been detected in cases of encephalitis across a wide spectrum of mammals, including placental and marsupial animals. Here we diagnosed two additional cases of fatal RusV-associated meningoencephalitis in a South American coati (Nasua nasua) and a Eurasian otter (Lutra lutra) that were detected in a zoological garden with history of prior RusV infections. Both animals showed abnormal movement or unusual behaviour and their brains tested positive for RusV using specific RT-qPCR and RNA in situ hybridization. As previous sequencing of RusV proved to be very challenging, we employed a sophisticated target-specific capture enrichment with specifically designed RNA baits to generate complete RusV genome sequences from both detected encephalitic animals and apparently healthy wild yellow-necked field mice (Apodemus flavicollis). Furthermore, the technique was used to revise three previously published RusV genomes from two encephalitic animals and a wild yellow-necked field mouse. Virus-to-host sequence ratio and thereby sequence coverage improved markedly using the enrichment method as compared to standard procedures. When comparing the newly generated RusV sequences to the previously published RusV genomes, we identified a previously undetected stretch of 309 nucleotides predicted to represent the intergenic region and the sequence encoding the N-terminus of the capsid protein. This indicated that the original RusV sequence was likely incomplete due to misassembly of the genome at a region with an exceptionally high G+C content of >80 mol%, which could not be resolved even by enormous sequencing efforts with standard methods. The updated capsid protein amino acid sequence now resembles those of RuV and ruhugu virus in size and harbours a predicted RNA binding domain that was not encoded in the original RusV genome version. The new sequence data indicate that RusV has the largest overall genome (9,631 nucleotides), intergenic region (290 nucleotides) and capsid protein-encoding sequence (331 codons) within the genus Rubivirus.

Structural Characterization of the Highly Restricted Down Syndrome Critical Region on 21q22.13: New KCNJ6 and DSCR4 Transcript Isoforms

Down syndrome (DS) is caused by trisomy of chromosome 21 and it is the most common genetic cause of intellectual disability (ID) in humans. Subjects with DS show a typical phenotype marked by facial dysmorphisms and ID. Partial trisomy 21 (PT21) is a rare genotype characterized by the duplication of a delimited chromosome 21 (Hsa21) portion and it may or may not be associated with DS diagnosis. The highly restricted Down syndrome critical region (HR-DSCR) is a region of Hsa21 present in three copies in all individuals with PT21 and a diagnosis of DS. This region, located on distal 21q22.13, is 34 kbp long and does not include characterized genes. The HR-DSCR is annotated as an intergenic region between KCNJ6-201 transcript encoding for potassium inwardly rectifying channel subfamily J member 6 and DSCR4-201 transcript encoding Down syndrome critical region 4. Two transcripts recently identified by massive RNA-sequencing (RNA-Seq) and automatically annotated on Ensembl database reveal that the HR-DSCR seems to be partially crossed by KCNJ6-202 and DSCR4-202 isoforms. KCNJ6-202 shares the coding sequence with KCNJ6-201 which is involved in many physiological processes, including heart rate in cardiac cells and circuit activity in neuronal cells. DSCR4-202 transcript has the first two exons in common with DSCR4-201, the only experimentally verified gene uniquely present in Hominidae. In this study, we performed in silico and in vitro analyses of the HR-DSCR. Bioinformatic data, obtained using Sequence Read Archive (SRA) and SRA-BLAST software, were confirmed by Reverse Transcription-Polymerase Chain Reaction (RT-PCR) and Sanger sequencing on a panel of human tissues. Our data demonstrate that the HR-DSCR cannot be defined as an intergenic region. Further studies are needed to investigate the functional role of the new transcripts, likely involved in DS phenotypes.

Lupus susceptibility region containing CTLA4 rs17268364 functionally reduces CTLA4 expression by binding EWSR1 and correlates IFN-α signature

Background Dysregulation of T cells mediated immune responses is a hallmark in the development of systemic lupus erythematosus (SLE). Recent genome wide association study (GWAS) revealed the genetic contribution of variants located in the cytotoxic T lymphocyte-associated protein-4 (CTLA4)-inducible T cell co-stimulator (ICOS) intergenic region to SLE susceptibility. Our aim is to find a functional variant in this region. Methods The genetic association results in the CTLA4-ICOS region from previous GWAS were adopted to select the potential variant which was further replicated in two independent cohorts (Henan cohort 2053 SLE patients and 1845 healthy controls, Beijing cohort 2303 SLE patients and 19,262 healthy). In order to explore the functional significance in SLE, bioinformatics with validation experiments (including electrophoretic mobility shift assay and luciferase reporter assay) and mRNA expression analysis were also performed. Results A variant located in the CTLA4-ICOS intergenic region, rs17268364, was associated with susceptibility to SLE patients in Chinese populations (risk allele, pmeta = 7.02×10−11, OR 1.19, 95%CI 1.13–1.26). The bioinformatics suggested that rs17268364 might affect the expression of CTLA4, not ICOS. The rs17268364 risk G allele containing sequence reduced the expression of the reporter gene by binding transcriptional repressor Ewing sarcoma breakpoint region 1 (EWSR1). Following genotype-mRNA expression, the analysis also showed the risk allele of rs17268364 was associated with low CTLA4 expression in lupus nephritis (LN) patients. Healthy individuals carrying rs17268364 risk G allele was significantly correlated with higher levels of IFN-α signature including increased lymphocyte antigen 6E (LY6E) (p=0.031), interferon-stimulated gene 15 (ISG15) (p=0.038), interferon regulatory factor 9 (IRF9) (p=0.028), and interferon regulatory factor 5 (IRF5) (p=0.040) mRNA expression. Conclusions The present study confirmed the functional role of rs17268364 in the CTLA4-ICOS intergenic region that increased SLE susceptibility in the Chinese population.

199. Polymorphisms in Key Regulatory Regions of the bla operon Correlate with the Cefazolin Inoculum Effect in Methicillin-Susceptible Staphylococcus aureus (MSSA)

Background The cefazolin inoculum effect (CzIE), defined as Cz minimum inhibitory concentration ≥ 16 µg/ml at high inoculum (HI-MIC), has been associated with poor clinical outcomes in patients with MSSA bacteremia or osteomyelitis. The CzIE is correlated with the presence of the blaZ gene, one of the components of the bla operon encoding the BlaZ β-lactamase (type A, B, C or D). Other portions of the bla operon include blaR and blaI (encoding the antibiotic sensor and transcriptional repressor, respectively) and the intergenic region with operator and promoter sequences (Figure 1). In BlaR, residue 293 mediates signal transduction, and the Z and R dyads in the intergenic region are the DNA-binding sites for BlaI (Figure 2). Previous experiments have shown that the regulatory portions of the bla operon play a key role in the CzIE. Here, we investigated the association between the CzIE and specific variations in the regulatory sequences of the bla operon. Figure 1. Functioning of the bla operon and the production of the staphylococcal β-lactamase BlaZ. Figure 2. Structure and key regions of the intergenic region of the bla operon, incluiding the promoter and the BlaI DNA-binding regions (Z and R dyads). Methods A total of 437 MSSA containing blaZ were evaluated for the CzIE using broth microdilution at high inoculum. Using whole genome sequencing, the sequences of the bla operons were classified into cassettes based on unique changes in predicted amino acid sequences of BlaZ, BlaR and BlaI paired with specific nucleotide alterations in the intergenic region. The bla operon sequence of S. aureus ATCC29213 was used as reference (cassette 0). Results Among 437 MSSA isolates, 46% exhibited the CzIE. We identified 55 unique bla cassettes. The bla cassettes were phylogenetically grouped in 7 clusters (Figure 3) which grouped cassettes with different BlaZ types and variations in the Z dyad, the -35 box, residue 293 of BlaR, and the blaI ribosomal binding site. Each cluster had an association to the CzIE and distinct Cz HI-MICs. The combination of: a BlaZ type A, C or D, an adenine in the position -66 of blaZ (-35 box of blaZ), a cytosine in the position -22 of blaZ(Z dyad), and either an arginine or a serine in position 293 of BlaR was a very strong predictor of the CzIE (Figure 4). Figure 3. Phylogenetical organization of bla operon cassettes into clusters, their association with polymorphisms in key regulatory regions and the CzIE. GM Cz-MIC: Geometric mean of the Cefazolin MIC at high Inoculum. Figure 4. Variations of the bla operon, their association with the CzIE, their GM (Geometric Mean) of the Cefazolin MIC and the MIC distribution of the strains with each specific combination of polymorphisms. Conclusion Specific variations in regulatory portions of the bla operon, which are likely to influence BlaZ expression, are highly associated with the CzIE, supporting the notion that regulation of blaZ is the key factor responsible for the CzIE in MSSA. Disclosures Lorena Diaz, PhD , Nothing to disclose William R. Miller, MD , Entasis Therapeutics (Scientific Research Study Investigator)Merck (Grant/Research Support) William R. Miller, MD , Entasis (Individual(s) Involved: Self): Scientific Research Study Investigator; Merck (Individual(s) Involved: Self): Grant/Research Support Cesar A. Arias, M.D., MSc, Ph.D., FIDSA, Entasis Therapeutics (Grant/Research Support)MeMed Diagnostics (Grant/Research Support)Merk (Grant/Research Support)

Mitochondrial DNA Variation in Peruvian Honey Bee (Apis mellifera L.) Populations Using the tRNAleu-cox2 Intergenic Region

Mitochondrial DNA variations of Peruvian honey bee populations were surveyed by using the tRNAleu-cox2 intergenic region. Only two studies have characterized these populations, indicating the presence of Africanized honey bee colonies in different regions of Peru and varied levels of Africanization, but the current status of its genetic diversity is unknown. A total of 512 honey bee colonies were sampled from three regions to characterize them. Our results revealed the presence of European and African haplotypes: the African haplotypes identified belong to sub-lineage AI (13) and sub-lineage AIII (03), and the European haplotypes to lineages C (06) and M (02). Of 24 haplotypes identified, 15 new sequences are reported here (11 sub-lineage AI, 2 sub-lineage AIII, and 2 lineage M). Peruvian honey bee populations presented a higher proportion from African than European haplotypes. High proportions of African haplotype were reported for Piura and Junín, unlike Lima, which showed more European haplotypes from lineage C. Few colonies belonging to lineage M would represent accidental purchase or traces of the introduction into Peru in the 19th century.

Novel EDGE encoding method enhances ability to identify genetic interactions

Assumptions are made about the genetic model of single nucleotide polymorphisms (SNPs) when choosing a traditional genetic encoding: additive, dominant, and recessive. Furthermore, SNPs across the genome are unlikely to demonstrate identical genetic models. However, running SNP-SNP interaction analyses with every combination of encodings raises the multiple testing burden. Here, we present a novel and flexible encoding for genetic interactions, the elastic data-driven genetic encoding (EDGE), in which SNPs are assigned a heterozygous value based on the genetic model they demonstrate in a dataset prior to interaction testing. We assessed the power of EDGE to detect genetic interactions using 29 combinations of simulated genetic models and found it outperformed the traditional encoding methods across 10%, 30%, and 50% minor allele frequencies (MAFs). Further, EDGE maintained a low false-positive rate, while additive and dominant encodings demonstrated inflation. We evaluated EDGE and the traditional encodings with genetic data from the Electronic Medical Records and Genomics (eMERGE) Network for five phenotypes: age-related macular degeneration (AMD), age-related cataract, glaucoma, type 2 diabetes (T2D), and resistant hypertension. A multi-encoding genome-wide association study (GWAS) for each phenotype was performed using the traditional encodings, and the top results of the multi-encoding GWAS were considered for SNP-SNP interaction using the traditional encodings and EDGE. EDGE identified a novel SNP-SNP interaction for age-related cataract that no other method identified: rs7787286 (MAF: 0.041; intergenic region of chromosome 7)–rs4695885 (MAF: 0.34; intergenic region of chromosome 4) with a Bonferroni LRT p of 0.018. A SNP-SNP interaction was found in data from the UK Biobank within 25 kb of these SNPs using the recessive encoding: rs60374751 (MAF: 0.030) and rs6843594 (MAF: 0.34) (Bonferroni LRT p: 0.026). We recommend using EDGE to flexibly detect interactions between SNPs exhibiting diverse action.

A Translation-Aborting Small Open Reading Frame in the Intergenic Region Promotes Translation of a Mg2+ Transporter in Salmonella Typhimurium

ABSTRACT Bacterial mRNAs often harbor upstream open reading frames (uORFs) in the 5′ untranslated regions (UTRs). Translation of the uORF usually affects downstream gene expression at the levels of transcription and/or translation initiation. Unlike other uORFs mostly located in the 5′ UTR, we discovered an 8-amino-acid ORF, designated mgtQ, in the intergenic region between the mgtC virulence gene and the mgtB Mg2+ transporter gene in the Salmonella mgtCBRU operon. Translation of mgtQ promotes downstream mgtB Mg2+ transporter expression at the level of translation by releasing the ribosome-binding sequence of the mgtB gene that is sequestered in a translation-inhibitory stem-loop structure. Interestingly, mgtQ Asp2 and Glu5 codons that induce ribosome destabilization are required for mgtQ-mediated mgtB translation. Moreover, the mgtQ Asp and Glu codons-mediated mgtB translation is counteracted by the ribosomal subunit L31 that stabilizes ribosome. Substitution of the Asp2 and Glu5 codons in mgtQ decreases MgtB Mg2+ transporter production and thus attenuates Salmonella virulence in mice, likely by limiting Mg2+ acquisition during infection. IMPORTANCE Translation initiation regions in mRNAs that include the ribosome-binding site (RBS) and the start codon are often sequestered within a secondary structure. Therefore, to initiate protein synthesis, the mRNA secondary structure must be unfolded to allow the RBS to be accessible to the ribosome. Such unfolding can be achieved by various mechanisms that include translation of a small upstream open reading frame (uORF). In the intracellular pathogen Salmonella enterica serovar Typhimurium, translation of the Mg2+ transporter mgtB gene is enhanced by an 8-amino-acid upstream ORF, namely, mgtQ, that harbors Asp and Glu codons, which are likely to destabilize ribosome during translation. Translation of the mgtQ ORF promotes the formation of a stem-loop mRNA structure sequestering anti-RBS and thus releases the mgtB RBS. Because mgtQ-mediated MgtB Mg2+ transporter production is required for Salmonella virulence, this pathogen seems to control the virulence determinant production exquisitely via this uORF during infection.