RNA Loci Research Articles

P309 Clonal outbreak of Trichophyton tonsurans causing tinea capitis among wrestlers in Beijing, China

Poster session 2, September 22, 2022, 12:30 PM - 1:30 PMObjectives Trichophyton tonsurans, an anthropophilic dermatophyte, mostly causes tinea capitis and tinea corporis and is often associated with outbreaks among athletes involved in combat sports, such as wrestlers and judo athletes. Here, we report an outbreak of tinea capitis caused by Trichophyton tonsurans among five juvenile athletes aged ten to fourteen in a wrestling team in Beijing, China.MethodsA total of 5 wrestlers aged 10-14 years presented with alopecia, erythema, scales, and pustule on the scalp (Fig. 1a). Scrapings from the lesions were performed by direct microscopic examination using 10% potassium hydroxide (KOH) with calcofluor white and fungal culture. The fungal pathogens were identified by morphology and sequencing of the internal transcribed spacer (ITS) regions. Multilocus genotyping analysis was performed by sequencing of the five gene loci including the ITS and non-transcribed spacer (NTS) of the ribosomal RNA (rRNA) locus, alkalineprotease-1 (ALP1), metalloprotease-5 (MEP5), carboxypeptidases Y (CarbY), which show intraspecies diversity and can be applied to epidemiological investigations and determining the route of infection transmission. Antifungal susceptibility of terbinafine (TBF), itraconazole (ITC), fluconazole (FLC), ketoconazole (KTC), and amphotericin B (AMB) against the causative fungal isolates was determined by broth microdilution method according to the Clinical and Laboratory Standards Institute (CLSI) M38-A3 document.ResultsThe five patients were diagnosed with tinea capitis because of the scalp lesions, spores, and hyphae observed by direct microscopic examination (Fig. 1b), and positive fungal culture. The causative isolates were all identified as T. tonsurans by macroscopic and microscopic characteristics (Fig. 2) and sequencing of ITS regions. The sequences of ITS, NTS, ALP1, MEP5, and CarbY were identical among the five isolates, revealing a single genotype. The abovementioned sequences have been deposited to Genbank. The minimum inhibitory concentrations (MICs) of TBF, ITC, FLC, KTC, and AMB against five T. tonsurans isolates were 0.015 μg/ml, 0.015 μg/ml, 8 μg/ml, 0.06 μg/ml, and 0.5 μg/ml, respectively.ConclusionWe report an outbreak of tinea capitis caused by T. tonsurans among wrestlers in Beijing, China. Multilocus genotyping analysis revealed that all isolates consisted of a single genotype, suggesting the outbreak may be caused by a single strain of T. tonsurans. The isolates were all susceptible to the common antifungal drugs treating tinea capitis.

Severe Outbreak of Dry Core Rot in Apple Fruits cv. Fuji Caused by Kalmusia variispora During Pre-harvest in Maule Region, Chile.

Apple (Malus × domestica) is an important fruit crop in Chile, with a cultivated area of 32,313 ha concentrated (63%) in Maule Region (35°25' S). Unusual core rot on 'Fuji' apples was observed at harvest in a commercial orchard in Curicó, Maule Region, with an incidence ranged between 22 to 35% in 2018 and 2019. Previously, in 2017, an incidence of 30% was estimated on 'Fuji' fruits destined to the Asian market. Internal decay symptoms consisted of dry, corky light to dark-brown tissue, within the seed locules initially. In moderate to severe cases, the necrotic lesion progresses deeper into the mesoderm. External symptoms were quite subtle, and typically, the disease goes unnoticed. However, infected fruit ripen earlier. Small pieces (2-3 mm) from the internal lesion margin of symptomatic apples (n = 50) were placed on Potato Dextrose Agar (PDA) (2%) and incubated at 20°C for 10 days. Pure cultures (n = 41) were obtained and transferred to Malt Extract Agar (MEA) (2%). Colonies on MEA produced an even to slight undulating buff margin with white woolly aerial mycelium, and immerse ochreous in the center, changing gray to olivaceous aerial mycelium with age. On the underside, colonies were umber and buff in the center and margin, respectively. After 10 days, numerous densely aggregate dark-brown mature pycnidia were observed. Aseptate conidia were subglobose to cylindrical, straight, and sometimes curved with rounded at both ends, that was initially hyaline to pale olive, thin, smooth wall with mean dimensions of (2.9-) 3.4 (-4.4) x (1.5-) 1.8 (2.2) µm (n=50). Based on morphology, the fungus was identified as Kalmusia variispora (Verkley et al. 2014). The internal transcribed space (ITS), portion of β-tubulin (TUB), and large subunits of the nuclear ribosomal RNA (LSU) loci were used for molecular identification, using primers ITS4/ITS5, Bt2a/Bt2b, and LR0R/LR5 (Ariyawansa et al., 2014). BLAST searches indicated 100% identity with K. variispora (ex-type CBS 121517). The maximum parsimony phylogenetic analysis placed Chilean isolates in the K. variispora clade. Sequences were deposited in GenBank (OL711706 to OL711709, OL739499 to OL739502 and OL711710 to OL711713 for ITS, TUB and LSU, respectively). Pathogenicity tests were conducted using four K. variispora isolates. 'Fuji' apples (n = 20) were surface disinfested (75% ethanol, 30 s) and then wounded and inoculated with conidial suspension (50 L of 106 conidia/mL) deposited in the middle and into the core region using a sterile fine-tipped micropipette. Additionally, 20 one-year dormant rooted cuttings 'Fuji' and 'Cripps Pink' were pruned and immediately inoculated on the pruning wound. An equal number of apples and rooted cuttings treated with sterile water were used as controls. The experiments were repeated once. All inoculated fruits developed lateral lesions (22 to 37 mm) and dry core rot (18 to 36 mm) symptoms identical to those described in the original outbreak, after 20 days at 20°C in a commercial packing box. The inoculated cuttings produced canker lesions of 10 to 21 mm in length, and dieback symptoms were observed after 3 months. No symptoms were observed on the negative controls. Koch's postulates were fulfilled by 100% reisolating K. variispora. Previously, Alternaria spp. have been reported as the primary pathogen associated with moldy core and dry core rot of apples worldwide (McLeod et al., 2014) and in Chile (Elfar et al., 2018). However, Kalmusia spp. have been associated with dry core rot in apples (McLeod et al., 2014) and have been isolated from canker symptoms on apples in Chile (Díaz et al. 2021). To our knowledge, this is the first report of a severe outbreak of K. variispora causing dry core rot in apples in Chile and worldwide.

The developmentally timed decay of an essential microRNA family is seed-sequence dependent.

SUMMARYMicroRNA (miRNA) abundance is tightly controlled by regulation of biogenesis and decay. Here, we show that the mir-35 miRNA family undergoes selective decay at the transition from embryonic to larval development in C. elegans. The seed sequence of the miRNA is necessary and largely sufficient for this regulation. Sequences outside the seed (3′ end) regulate mir-35 abundance in the embryo but are not necessary for sharp decay at the transition to larval development. Enzymatic modifications of the miRNA 3′ end are neither prevalent nor correlated with changes in decay, suggesting that miRNA 3′ end display is not a core feature of this mechanism and further supporting a seed-driven decay model. Our findings demonstrate that seed-sequence-specific decay can selectively and coherently regulate all redundant members of a miRNA seed family, a class of mechanism that has great biological and therapeutic potential for dynamic regulation of a miRNA family’s target repertoire.

Population genetics of Cryptosporidium parvum subtypes in cattle in Poland: the geographical change of strain prevalence and circulation over time

BackgroundCryptosporidium parvum (C. parvum) is a cosmopolitan parasite that infects various livestock animals including cattle. Microsatellite typing tools for identification of C. parvum subtypes are currently employed to better understand the species-specific epidemiology of cattle cryptosporidiosis. The aim of this study was to analyse the population genetics of C. parvum strains infecting cattle and recognise geographical distribution and time-span correlations in subtype prevalence in Poland. In total, 1601 faecal samples were collected from 2014 to 2018 from healthy cattle from dairy, meat and mixed breeds at the age of 1 week to 4 months. The 267 farms visited were randomly selected and represented all Polish provinces. PCR–RFLP based identification of C. parvum at the 18 small subunit ribosomal RNA (SSU rRNA) locus was performed, followed by strain subtyping by GP60-PCR. ResultsThe overall prevalence of C. parvum in Polish cattle was estimated at 6.2% (100/1601). Animals below the age of 1 month were the major host for this parasite. Excluding one breed, that of dairy-meat mixed, there were no significant differences observed between breed and presence of C. parvum infections (95% TPIAll breeds: 1.67–73.53%; POPR = 0.05—0.95). Infected animals were detected in 15 out of 16 Polish provinces, with significant regional prevalence diffrences (Kruskal–Wallis rank sum test, Kruskal–Wallis χ2 = 13.46, p < 0.001). When the population genetics of C. parvum strains were analysed, 11 parasite subtypes from the IIa and IId genetic families were identified. Compared to other parasite strains, IIaA17G1R1 and IIaA17G2R1 appeared at statistically significantly higher frequency (F-test, F = 3.39; p = 0.0003). The prevalence of C. parvum subtypes in cattle was breed-related (Chi-squared test, χ2 = 143.6; p < 0.001).ConclusionsThe analysis of the population genetics of C. parvum subtypes showed that strains from the IIa subtype family predominated in the tested cattle population. However, relations in changes of subtype prevalence and circulation over time were observed. They were associated with the disappearance of some strains and emergence of new variants from the same genetic family in different geographical locations.

Dissecting super-enhancer heterogeneity: time to re-examine cancer subtypes?

The heterogeneity of transcriptional regulations by super-enhancers (SEs) is poorly understood in human cancers. Herein, we summarize a bioinformatics workflow for genome-wide SE profiling and identification of subtype-specific SEs and regulatory networks. Dissecting SE heterogeneity provides new insights into cancer biology and alternative therapeutic strategies for cancer precision medicine.

Multiple Cases of Bacterial Sequence Erroneously Incorporated Into Publicly Available Chloroplast Genomes.

Public sequencing databases are invaluable resources to biological researchers, but assessing data veracity as well as the curation and maintenance of such large collections of data can be challenging. Genomes of eukaryotic organelles, such as chloroplasts and other plastids, are particularly susceptible to assembly errors and misrepresentations in these databases due to their close evolutionary relationships with bacteria, which may co-occur within the same environment, as can be the case when sequencing plants. Here, based on sequence similarities with bacterial genomes, we identified several suspicious chloroplast assemblies present in the National Institutes of Health (NIH) Reference Sequence (RefSeq) collection. Investigations into these chloroplast assemblies reveal examples of erroneous integration of bacterial sequences into chloroplast ribosomal RNA (rRNA) loci, often within the rRNA genes, presumably due to the high similarity between plastid and bacterial rRNAs. The bacterial lineages identified within the examined chloroplasts as the most likely source of contamination are either known associates of plants, or co-occur in the same environmental niches as the examined plants. Modifications to the methods used to process untargeted ‘raw’ shotgun sequencing data from whole genome sequencing efforts, such as the identification and removal of bacterial reads prior to plastome assembly, could eliminate similar errors in the future.

From genotype to phenotype: genetics of mammalian long non-coding RNAs in vivo.

Genome-wide sequencing has led to the discovery of thousands of long non-coding RNA (lncRNA) loci in the human genome, but evidence of functional significance has remained controversial for many lncRNAs. Genetically engineered model organisms are considered the gold standard for linking genotype to phenotype. Recent advances in CRISPR-Cas genome editing have led to a rapid increase in the use of mouse models to more readily survey lncRNAs for functional significance. Here, we review strategies to investigate the physiological relevance of lncRNA loci by highlighting studies that have used genetic mouse models to reveal key in vivo roles for lncRNAs, from fertility to brain development. We illustrate how an investigative approach, starting with whole-gene deletion followed by transcription termination and/or transgene rescue strategies, can provide definitive evidence for the in vivo function of mammalian lncRNAs.

Myeloid Leukemia Dependencies at CTCF-Enriched Long Noncoding RNA Loci

The noncoding genome presents a largely untapped source of biological insights, including tens of thousands of long noncoding RNA (lncRNA) loci. While some produce bona fide lncRNAs, others exert transcript-independent cis-regulatory effects, and the lack of predictive features renders their mechanistic dissection highly challenging. Here, we describe CTCF-enriched lncRNA loci (C-LNC) as a putative new subclass of functional genetic elements exemplified by MYNRL15 - myeloid leukemia noncoding regulatory locus on chromosome 15.Initially identified by an expression-guided CRISPRi screen of hematopoietic stem and progenitor (HSPC) / acute myeloid leukemia (AML) lncRNA signatures (480 genes, 1545 sgRNAs), we found MYNRL15 dependency in myeloid leukemia cells of diverse genetic backgrounds. Interestingly, cis and trans perturbation approaches revealed both the MYNRL15 transcript and its flanking protein-coding genes to be dispensable. High density CRISPR tiling of a 15 kb area centered on MYNRL15 (1613 sgRNAs) instead uncovered two crucial, candidate cis-regulatory DNA elements in the locus, which drive the MYNRL15 perturbation phenotype. To determine the molecular basis of MYNRL15 dependence, we performed transcriptome, chromatin conformation, chromatin accessibility, and CTCF profiling. RNA-sequencing established MYNRL15's involvement in maintaining key cancer dependency pathways (e.g. cell cycle, ribosome, spliceosome). Further, MYNRL15 perturbation associated with the coordinated dysregulation of several chromosome 15 neighbourhoods, and formation of a long-range chromatin interaction between the locus and the base of a distal loop, as detected via next-generation Capture-C. The gained interaction was accompanied by diffuse gains in chromatin accessibility across the distal interaction sites (ATAC-seq) as well as reduced CTCF occupancy at the MYRNL15 locus (CTCF CUT&RUN), altogether indicating the 3D re-organization of chromosome 15 following MYNRL15 perturbation. Integrative analysis of the chromatin conformation and transcriptome data, combined with a small CRISPR-Cas9 knockout screen of protein-coding genes from the gained interaction region (29 genes, 149 sgRNAs), pinpointed two potent cancer dependency genes that are located in the region and downregulated following MYNRL15 perturbation: namely, WDR61 and IMP3. Individual knockout of both genes robustly depleted myeloid leukemia cells, recapitulating the MYNRL15 perturbation phenotype and positioning WDR61 and IMP3 as its regulatory targets. Importantly, in primary cells, MYNRL15 perturbation eradicated AML blasts while sparing 50-60% of CD34 + HSPCs in vitro, and reduced patient-derived AML xenografts up to 10-fold in vivo, indicating a potential therapeutic window.Having implicated MYNRL15 in 3D genome organization and demonstrated its role in myeloid leukemia cells, we explored whether MYNRL15 may belong to a sub-category of biologically relevant lncRNA loci that have thus far been overlooked due to their lack of transcript-specific functions. Remarkably, elevated CTCF density (e.g. number of CTCF binding sites per kb of gene length) distinguishes MYNRL15 and 531 other lncRNA loci in K562 cells, of which 43-54% associate with genetic subgroups and/or survival in AML patient cohorts, and 18.4% are functionally required for leukemia maintenance as determined by CRISPR-Cas9 screening. The latter hit identification rate represents a substantial improvement over typical lncRNA essentiality screens (which range from 2-6%) - illustrating the effectiveness of CTCF density metrics in refining functional lncRNA candidate lists, and underlining the relevance such loci hold for AML and cancer pathophysiology in general. Curated C-LNC catalogs in other cell types will facilitate the search for noncoding oncogenic vulnerabilities in AML and other malignancies. [Display omitted] DisclosuresReinhardt: Celgene Corporation: Consultancy; Novartis: Consultancy; Bluebird Bio: Consultancy; Janssen: Consultancy; CLS Behring: Research Funding; Roche: Research Funding. Klusmann: Bluebird Bio: Consultancy; Novartis: Consultancy; Roche: Consultancy; Jazz Pharmaceuticals: Consultancy.

Regulation of Fetal Hemoglobin Expression By the VHL-HIF1α Oxygen Sensing System

Defining the mechanisms that control the perinatal switch from γ-globin (HBG1 and HBG2) to β-globin (HBB) gene expression in human red blood cells (RBCs) has informed novel approaches to reactivate fetal hemoglobin (HbF, α2γ2) therapeutically for sickle cell disease and β-thalassemia. However, one longstanding unsolved problem is to explain how HbF becomes elevated in conditions such as blood loss, hypoxia and hemolysis. These conditions are associated with accelerated RBC production, also referred to as stress erythropoiesis, driven by activation of hypoxia-inducible factor (HIF) via a canonical O 2 sensing pathway. At high O 2 levels (“normoxia”), O 2-dependent prolyl hydroxylase domain (PHD) enzymes hydroxylate HIFα, thereby targeting it for ubiquitination by the von Hippel-Lindau (VHL) E3 ubiquitin ligase complex, followed by proteasomal degradation. At low O 2 tension (hypoxia), PHD activity is reduced, causing HIF1α to accumulate, dimerize with constitutively expressed HIF1β, and bind hypoxia response elements (HREs) to activate a broad array of genes that facilitate hypoxic adaptation. We identified VHL and HIF1α, as negative and positive regulators of HbF expression, respectively.Disruption of the VHL gene in CD34 + hematopoietic stem and progenitor cells (HSPCs) by transfection with ribonucleoprotein (RNP) consisting of Cas9 and VHL-targeting guide RNA increased HbF expression from 7.5% ±1.2% in control cells to 30.9% ± 4.8% (mean ± SD, P<0.0001) in RBC progeny generated by in vitro differentiation. Similarly, γ-globin mRNA was induced by 5 -fold after disruption of VHL in HUDEP-2 cells, an immortalized erythroid line that expresses mainly adult hemoglobin (HbA, α2/β2). Mass spectrometry and transcriptome analysis of VHL-disrupted CD34 + HSPCs revealed increased HIF1α protein expression with no change in the corresponding mRNA. Inhibition of HIF1A mRNA by RNA interference suppressed γ-globin induction in VHL-/- HUDEP-2 clones and in RBCs generated from VHL RNP-treated CD34 + cells, indicating VHL knockout induced HbF through HIF1α protein accumulation. CUT&RUN analysis of VHL-depleted erythroblasts revealed HIF1α/HIF1β heterodimer occupancy at BGLT3, a long-noncoding RNA locus approximately 2.7kb 3' of HBG1. The HIF1α-bound BGLT3 locus contains two canonical HREs (ACGTG) separated by 13 bp. Disruption of each HRE motif in VHL-/- HUDEP-2 cells by base editing caused additive reductions in BGLT3 HIF1α occupancy and γ-globin expression. Mechanistically, VHL depletion caused the accumulation of HIFα/β heterodimers at BGLT3, recruitment of the transcriptional activators GATA1 and P300, and targeted chromatin accessibility to establish the active enhancer mark H3K27ac. These changes were accompanied by altered chromosome conformation to favor long-range interactions between the γ-globin loci (HBG1 and HBG2) and the locus control region, a powerful upstream enhancer. Treatment of healthy donor CD34 + HSPCs-derived erythroblasts with the clinically approved PHD inhibitor FG-4592 (Roxadustat) caused HbF to increase from 4.76%±1.22% at baseline to 12.86 ± 2.40% (mean ± SD in 3 biological replicates, P<0.01). Treatment of the same erythroblasts with FG-4592 and hydroxyurea, a widely used SCD drug that acts partly by inducing HbF, caused HbF to increase from 4.76% ± 1.22% at baseline to 24.1% ± 5.3% (mean ± SD in 3 biological replicates, P<0.01), indicating an additive effect. Our findings link developmental globin gene regulation with O 2 sensing, provide a mechanism for HbF induction during stress erythropoiesis, and identify a novel therapeutic approach for β-hemoglobinopathies. DisclosuresBlobel: Pfizer: Consultancy; Fulcrum Therapeutics, Inc.: Consultancy. Weiss: Beam Therapeutics: Current holder of stock options in a privately-held company; Forma Therapeutics: Consultancy; Cellarity Inc.: Consultancy; Novartis: Consultancy.

Genome-wide CRISPR interference screen identifies long non-coding RNA loci required for differentiation and pluripotency.

Although many long non-coding RNAs (lncRNAs) exhibit lineage-specific expression, the vast majority remain functionally uncharacterized in the context of development. Here, we report the first described human embryonic stem cell (hESC) lines to repress (CRISPRi) or activate (CRISPRa) transcription during differentiation into all three germ layers, facilitating the modulation of lncRNA expression during early development. We performed an unbiased, genome-wide CRISPRi screen targeting thousands of lncRNA loci expressed during endoderm differentiation. While dozens of lncRNA loci were required for proper differentiation, most differentially expressed lncRNAs were not, supporting the necessity for functional screening instead of relying solely on gene expression analyses. In parallel, we developed a clustering approach to infer mechanisms of action of lncRNA hits based on a variety of genomic features. We subsequently identified and validated FOXD3-AS1 as a functional lncRNA essential for pluripotency and differentiation. Taken together, the cell lines and methodology described herein can be adapted to discover and characterize novel regulators of differentiation into any lineage.

Prevalence of Trypanosome Species in Cattle Near Ruma National Park, Lambwe Valley, Kenya: An Update From the Historical Focus for African Trypanosomosis.

The effective control of diseases in areas shared with wildlife depends on the validity of the epidemiologic parameters that guide interventions. Epidemiologic data on animal trypanosomosis in Lambwe valley are decades old, and the recent suspected outbreaks of the disease in the valley necessitate the urgent bridging of this data gap. This cross-sectional study estimated the prevalence of bovine trypanosomosis, identified risk factors, and investigated the occurrence of species with zoonotic potential in Lambwe valley. The area is ~324 km2, of which 120 km2 is the Ruma National Park. Blood was sampled from the jugular and marginal ear veins of 952 zebu cattle between December 2018 and February 2019 and tested for trypanosomes using the Buffy Coat Technique (BCT) and PCR-High-Resolution Melting (HRM) analysis of the 18S RNA locus. Risk factors for the disease were determined using logistic regression. The overall trypanosome prevalence was 11.0% by BCT [95% confidence interval (CI): 9.0–13.0] and 27.9% by PCR-HRM (95% CI: 25.1–30.8). With PCR-HRM as a reference, four species of trypanosomes were detected at prevalences of 12.7% for T. congolense savannah (95% CI: 10.6–14.8), 7.7% for T. brucei brucei (CI: 6.0–9.4), 8.7% for T. vivax (CI: 6.9–10.5), and 1.3% for T. theileri (CI: 0.6–2.0). About 2.4% of cattle had mixed infections (CI: 1.4–3.41). No human-infective trypanosomes were found. Infections clustered across villages but were not associated with animal age, sex, herd size, and distance from the park. Approximately 85% of infections occurred within 2 km of the park. These findings add to evidence that previous interventions eliminated human trypanosomosis but not bovine trypanosomosis. Risk-tailored intervention within 2 km of Ruma Park, especially in the north and south ends, coupled with stringent screening with molecular tools, could significantly reduce bovine trypanosomosis.

High-quality reference genome and annotation aids understanding of berry development for evergreen blueberry (Vaccinium darrowii)

Vaccinium darrowii Camp (2n = 2x = 24) is a native North American blueberry species and an important source of traits such as low chill requirement in commercial southern highbush blueberry breeding (Vaccinium corymbosum, 2n = 4x = 48). We present a chromosomal-scale genome of V. darrowii generated by the combination of PacBio sequencing and high throughput chromatin conformation capture (Hi–C) scaffolding technologies, yielding a total length of 1.06 Gigabases (Gb). Over 97.8% of the genome sequences are scaffolded into 24 chromosomes representing the two haplotypes. The primary haplotype assembly of V. darrowii contains 34,809 protein-coding genes. Comparison to a V. corymbosum haplotype assembly reveals high collinearity between the two genomes with small intrachromosomal rearrangements in eight chromosome pairs. With small RNA sequencing, the annotation was further expanded to include more than 200,000 small RNA loci and 638 microRNAs expressed in berry tissues. Transcriptome analysis across fruit development stages indicates that genes involved in photosynthesis are downregulated, while genes involved in flavonoid and anthocyanin biosynthesis are significantly increased at the late stage of berry ripening. A high-quality reference genome and accompanying annotation of V. darrowii is a significant new resource for assessing the evergreen blueberry contribution to the breeding of southern highbush blueberries.

Genome-Wide Analysis of the Regulatory Landscape in Different Pre-B-ALL Subtypes

BackgroundAcute pre-B leukemia results from a failure in cell differentiation. In normal development, cells commit to their lineage and differentiated phenotype state by the activation of transcription factors (TFs), or TF modules, in a series of decision points at which a choice is made between alternative lineage fates. Previously, we showed that the ETV6-RUNX1 gene translocation that occurs in 25% of pediatric leukemia patients silences regulatory regions containing RUNX binding sites (enhancers) at target gene loci implicated in lymphoid lineage differentiation (Teppo et al. 2016). However, the genome-wide changes in enhancer activities at binding sites of other TFs functioning in B-cell differentiation remain poorly understood. Moreover, the genome-wide identification of enhancers altered during leukemogenesis is incomplete.Material & MethodsTo measure transcripts generated both at enhancer and gene regions, global run-on sequencing (GRO-seq) technique was used to measure nascent RNAs (Heinäniemi et al. 2016) from human pre-B-ALL patient samples (N=35) and several cell lines representing different subtypes. For validation, additional samples were generated from cell line models with overexpression or silencing of leukemic TFs. To this collection of 151 nascent transcriptomes, additional 92 GRO-seq samples representing different cell lineages (Liivrand et al 2017) were used to compare the leukemic cells to primary cells and other cancers. Enhancer RNA (eRNA) levels were compared at TF binding motif centered enhancer regions. TF ChIP-seq data from leukemic cells and HSC were acquired from the NCBI Gene Expression Omnibus for validation of this approach.ResultsTo gain insight into the derailed regulatory landscape in leukemic cells, we identified eRNA regions both from normal and leukemic samples and performed statistical analysis distinguishing regulatory regions with differential activity between pre-B-ALL-subtypes or compared to normal CD34+ bone marrow precursors. Next, association of the significantly altered enhancers with gene transcripts was utilized to build regulatory maps for differentially expressed gene loci. For this purpose, we developed a transcript identification and quantification pipeline suitable for detecting altered expression at different coding and non-coding RNA loci from nascent transcriptomes (Liivrand et al 2017; unpublished). To globally evaluate altered TF activities across pre-B-ALL subtypes underlying the differences observed, we present here an unbiased approach based on eRNA quantification at TF motif-centered enhancers. The analysis of eRNA levels at binding sites of TFs with established role in leukemia was used as a benchmark: we could confirm low eRNA transcription at RUNX1 binding sites in the ETV6-RUNX1 subtype, elevated eRNA levels at HOXA9 and MEIS1 binding sites in the MLL rearranged subtype, while TCF3, BCL6 and PBX1 bound sites were most active in TCF3-PBX1 samples.ConclusionsOur results provide the first genome-wide transcriptional regulatory landscape of different pre-B-ALL subtypes. We identified novel subtype-specific active enhancers and annotated regulatory non-coding regions that can be utilized to evaluate the functional impact of genetic alterations outside gene regions. DisclosuresNo relevant conflicts of interest to declare.

Sporadic isolation of Tetratrichomonas species from the cattle urogenital tract.

Tritrichomonas foetus is a venereal trichomonad parasite which causes reproductive issues in cattle. No other trichomonads are known to be urogenital pathogens in cattle, but there are several reports of Tetratrichomonas and Pentatrichomonas isolates of unclear origin from the cattle urogenital tract (UGT) in the Americas. This study reports the first case of a non-T. foetus cattle urogenital trichomonad isolate in Europe. Molecular analysis of the internal transcribed spacer (ITS) 1-5.8S ribosomal RNA-ITS 2 and 18S ribosomal RNA loci suggest that the isolate is a Tetratrichomonas species from a lineage containing other previously described bull preputial isolates. We identified close sequence similarity between published urogenital and gastrointestinal Tetratrichomonas spp., and this is reviewed alongside further evidence regarding the gastrointestinal origin of non-T. foetus isolates. Routine screening for T. foetus is based on culture and identification by microscopy, and so considering other trichomonad parasites of the bovine UGT is important to avoid misdiagnosis.

Expansion of rDNA and pericentromere satellite repeats in the genomes of bank voles Myodes glareolus exposed to environmental radionuclides.

Altered copy number of certain highly repetitive regions of the genome, such as satellite DNA within heterochromatin and ribosomal RNA loci (rDNA), is hypothesized to help safeguard the genome against damage derived from external stressors. We quantified copy number of the 18S rDNA and a pericentromeric satellite DNA (Msat‐160) in bank voles (Myodes glareolus) inhabiting the Chernobyl Exclusion Zone (CEZ), an area that is contaminated by radionuclides and where organisms are exposed to elevated levels of ionizing radiation. We found a significant increase in 18S rDNA and Msat‐160 content in the genomes of bank voles from contaminated locations within the CEZ compared with animals from uncontaminated locations. Moreover, 18S rDNA and Msat‐160 copy number were positively correlated in the genomes of bank voles from uncontaminated, but not in the genomes of animals inhabiting contaminated, areas. These results show the capacity for local‐scale geographic variation in genome architecture and are consistent with the genomic safeguard hypothesis. Disruption of cellular processes related to genomic stability appears to be a hallmark effect in bank voles inhabiting areas contaminated by radionuclides.

A Small RNA Is Linking CRISPR-Cas and Zinc Transport.

The function and mode of action of small regulatory RNAs is currently still understudied in archaea. In the halophilic archaeon Haloferax volcanii, a plethora of sRNAs have been identified; however, in-depth functional analysis is missing for most of them. We selected a small RNA (s479) from Haloferax volcanii for detailed characterization. The sRNA gene is encoded between a CRISPR RNA locus and the Cas protein gene cluster, and the s479 deletion strain is viable and was characterized in detail. Transcriptome studies of wild-type Haloferax cells and the deletion mutant revealed upregulation of six genes in the deletion strain, showing that this sRNA has a clearly defined function. Three of the six upregulated genes encode potential zinc transporter proteins (ZnuA1, ZnuB1, and ZnuC1) suggesting the involvement of s479 in the regulation of zinc transport. Upregulation of these genes in the deletion strain was confirmed by northern blot and proteome analyses. Furthermore, electrophoretic mobility shift assays demonstrate a direct interaction of s479 with the target znuC1 mRNA. Proteome comparison of wild-type and deletion strains further expanded the regulon of s479 deeply rooting this sRNA within the metabolism of H. volcanii especially the regulation of transporter abundance. Interestingly, s479 is not only encoded next to CRISPR–cas genes, but the mature s479 contains a crRNA-like 5′ handle, and experiments with Cas protein deletion strains indicate maturation by Cas6 and interaction with Cas proteins. Together, this might suggest that the CRISPR–Cas system is involved in s479 function.

Intergenic RNA mainly derives from nascent transcripts of known genes

BackgroundEukaryotic genomes undergo pervasive transcription, leading to the production of many types of stable and unstable RNAs. Transcription is not restricted to regions with annotated gene features but includes almost any genomic context. Currently, the source and function of most RNAs originating from intergenic regions in the human genome remain unclear.ResultsWe hypothesize that many intergenic RNAs can be ascribed to the presence of as-yet unannotated genes or the “fuzzy” transcription of known genes that extends beyond the annotated boundaries. To elucidate the contributions of these two sources, we assemble a dataset of more than 2.5 billion publicly available RNA-seq reads across 5 human cell lines and multiple cellular compartments to annotate transcriptional units in the human genome. About 80% of transcripts from unannotated intergenic regions can be attributed to the fuzzy transcription of existing genes; the remaining transcripts originate mainly from putative long non-coding RNA loci that are rarely spliced. We validate the transcriptional activity of these intergenic RNAs using independent measurements, including transcriptional start sites, chromatin signatures, and genomic occupancies of RNA polymerase II in various phosphorylation states. We also analyze the nuclear localization and sensitivities of intergenic transcripts to nucleases to illustrate that they tend to be rapidly degraded either on-chromatin by XRN2 or off-chromatin by the exosome.ConclusionsWe provide a curated atlas of intergenic RNAs that distinguishes between alternative processing of well-annotated genes from independent transcriptional units based on the combined analysis of chromatin signatures, nuclear RNA localization, and degradation pathways.

Chromosome-level genome assembly of the hard-shelled mussel Mytilus coruscus, a widely distributed species from the temperate areas of East Asia.

BackgroundThe hard-shelled mussel (Mytilus coruscus) is widely distributed in the temperate seas of East Asia and is an important commercial bivalve in China. Chromosome-level genome information of this species will contribute not only to the development of hard-shelled mussel genetic breeding but also to studies on larval ecology, climate change biology, marine biology, aquaculture, biofouling, and antifouling.FindingsWe applied a combination of Illumina sequencing, Oxford Nanopore Technologies sequencing, and high-throughput chromosome conformation capture technologies to construct a chromosome-level genome of the hard-shelled mussel, with a total length of 1.57 Gb and a median contig length of 1.49 Mb. Approximately 90.9% of the assemblies were anchored to 14 linkage groups. We assayed the genome completeness using BUSCO. In the metazoan dataset, the present assemblies have 89.4% complete, 1.9% incomplete, and 8.7% missing BUSCOs. Gene modeling enabled the annotation of 37,478 protein-coding genes and 26,917 non-coding RNA loci. Phylogenetic analysis showed that M. coruscus is the sister taxon to the clade including Modiolus philippinarum and Bathymodiolus platifrons. Conserved chromosome synteny was observed between hard-shelled mussel and king scallop, suggesting that this is shared ancestrally. Transcriptomic profiling indicated that the pathways of catecholamine biosynthesis and adrenergic signaling in cardiomyocytes might be involved in metamorphosis.ConclusionsThe chromosome-level assembly of the hard-shelled mussel genome will provide novel insights into mussel genome evolution and serve as a fundamental platform for studies regarding the planktonic-sessile transition, genetic diversity, and genomic breeding of this bivalve.

Genomic Space of MGMT in Human Glioma Revisited: Novel Motifs, Regulatory RNAs, NRF1, 2, and CTCF Involvement in Gene Expression.

Background: The molecular regulation of increased MGMT expression in human brain tumors, the associated regulatory elements, and linkages of these to its epigenetic silencing are not understood. Because the heightened expression or non-expression of MGMT plays a pivotal role in glioma therapeutics, we applied bioinformatics and experimental tools to identify the regulatory elements in the MGMT and neighboring EBF3 gene loci. Results: Extensive genome database analyses showed that the MGMT genomic space was rich in and harbored many undescribed RNA regulatory sequences and recognition motifs. We extended the MGMT’s exon-1 promoter to 2019 bp to include five overlapping alternate promoters. Consensus sequences in the revised promoter for (a) the transcriptional factors CTCF, NRF1/NRF2, GAF, (b) the genetic switch MYC/MAX/MAD, and (c) two well-defined p53 response elements in MGMT intron-1, were identified. A putative protein-coding or non-coding RNA sequence was located in the extended 3′ UTR of the MGMT transcript. Eleven non-coding RNA loci coding for miRNAs, antisense RNA, and lncRNAs were identified in the MGMT-EBF3 region and six of these showed validated potential for curtailing the expression of both MGMT and EBF3 genes. ChIP analysis verified the binding site in MGMT promoter for CTCF which regulates the genomic methylation and chromatin looping. CTCF depletion by a pool of specific siRNA and shRNAs led to a significant attenuation of MGMT expression in human GBM cell lines. Computational analysis of the ChIP sequence data in ENCODE showed the presence of NRF1 in the MGMT promoter and this occurred only in MGMT-proficient cell lines. Further, an enforced NRF2 expression markedly augmented the MGMT mRNA and protein levels in glioma cells. Conclusions: We provide the first evidence for several new regulatory components in the MGMT gene locus which predict complex transcriptional and posttranscriptional controls with potential for new therapeutic avenues.

Composition and distribution of fish environmental DNA in an Adirondack watershed.

BackgroundEnvironmental DNA (eDNA) surveys are appealing options for monitoring aquatic biodiversity. While factors affecting eDNA persistence, capture and amplification have been heavily studied, watershed-scale surveys of fish communities and our confidence in such need further exploration.MethodsWe characterized fish eDNA compositions using rapid, low-volume filtering with replicate and control samples scaled for a single Illumina MiSeq flow cell, using the mitochondrial 12S ribosomal RNA locus for taxonomic profiling. Our goals were to determine: (1) spatiotemporal variation in eDNA abundance, (2) the filtrate needed to achieve strong sequencing libraries, (3) the taxonomic resolution of 12S ribosomal sequences in the study environment, (4) the portion of the expected fish community detectable by 12S sequencing, (5) biases in species recovery, (6) correlations between eDNA compositions and catch per unit effort (CPUE) and (7) the extent that eDNA profiles reflect major watershed features. Our bioinformatic approach included (1) estimation of sequencing error from unambiguous mappings and simulation of taxonomic assignment error under various mapping criteria; (2) binning of species based on inferred assignment error rather than by taxonomic rank; and (3) visualization of mismatch distributions to facilitate discovery of distinct haplotypes attributed to the same reference. Our approach was implemented within the St. Regis River, NY, USA, which supports tribal and recreational fisheries and has been a target of restoration activities. We used a large record of St. Regis-specific observations to validate our assignments.ResultsWe found that 300 mL drawn through 25-mm cellulose nitrate filters yielded greater than 5 ng/µL DNA at most sites in summer, which was an approximate threshold for generating strong sequencing libraries in our hands. Using inferred sequence error rates, we binned 12S references for 110 species on a state checklist into 85 single-species bins and seven multispecies bins. Of 48 bins observed by capture survey in the St. Regis, we detected eDNA consistent with 40, with an additional four detections flagged as potential contaminants. Sixteen unobserved species detected by eDNA ranged from plausible to implausible based on distributional data, whereas six observed species had no 12S reference sequence. Summed log-ratio compositions of eDNA-detected taxa correlated with log(CPUE) (Pearson’s R = 0.655, P < 0.001). Shifts in eDNA composition of several taxa and a genotypic shift in channel catfish (Ictalurus punctatus) coincided with the Hogansburg Dam, NY, USA. In summary, a simple filtering apparatus operated by field crews without prior expertise gave useful summaries of eDNA composition with minimal evidence of field contamination. 12S sequencing achieved useful taxonomic resolution despite the short marker length, and data exploration with standard bioinformatic tools clarified taxonomic uncertainty and sources of error.