Search

Understanding HIV-1 persistence despite antiretroviral therapy (ART) is of paramount importance. Both single-genome sequencing (SGS) and integration site analysis (ISA) provide useful information regarding the structure of persistent HIV DNA populations; however, until recently, there was no way to link integration sites to their cognate proviral sequences. Here, we used multiple-displacement amplification (MDA) of cellular DNA diluted to a proviral endpoint to obtain full-length proviral sequences and their corresponding sites of integration. We applied this method to lymph node and peripheral blood mononuclear cells from 5 ART-treated donors to determine whether groups of identical subgenomic sequences in the 2 compartments are the result of clonal expansion of infected cells or a viral genetic bottleneck. We found that identical proviral sequences can result from both cellular expansion and viral genetic bottlenecks occurring prior to ART initiation and following ART failure. We identified an expanded T cell clone carrying an intact provirus that matched a variant previously detected by viral outgrowth assays and expanded clones with wild-type and drug-resistant defective proviruses. We also found 2 clones from 1 donor that carried identical proviruses except for nonoverlapping deletions, from which we could infer the sequence of the intact parental virus. Thus, MDA-SGS can be used for "viral reconstruction" to better understand intrapatient HIV-1 evolution and to determine the clonality and structure of proviruses within expanded clones, including those with drug-resistant mutations. Importantly, we demonstrate that identical sequences observed by standard SGS are not always sufficient to establish proviral clonality.

We assessed the results of genome sequencing for early-onset dementia. Participants were selected from a memory disorders clinic. Genome sequencing was performed along with C9orf72 repeat expansion testing. All returned sequencing results were Sanger-validated. Prior clinical diagnoses included Alzheimer's disease, frontotemporal dementia, and unspecified dementia. The mean age of onset was 54 (41-76). Fifty percent of patients had a strong family history, 37.5% had some, and 12.5% had no known family history. Nine of 32 patients (28%) had a variant defined as pathogenic or likely pathogenic (P/LP) by American College of Medical Genetics and Genomics standards, including variants in APP, C9orf72, CSF1R, and MAPT Nine patients (including three with P/LP variants) harbored established risk alleles with moderate penetrance (odds ratios of ∼2-5) in ABCA7, AKAP9, GBA, PLD3, SORL1, and TREM2 All six patients harboring these moderate penetrance variants but not P/LP variants also had one or two APOE ε4 alleles. One patient had two APOE ε4 alleles with no other established contributors. In total, 16 patients (50%) harbored one or more genetic variants likely to explain symptoms. We identified variants of uncertain significance (VUSs) in ABI3, ADAM10, ARSA, GRID2IP, MME, NOTCH3, PLCD1, PSEN1, TM2D3, TNK1, TTC3, and VPS13C, also often along with other variants. In summary, genome sequencing for early-onset dementia frequently identified multiple established or possible contributory alleles. These observations add support for an oligogenic model for early-onset dementia.

The detection of thrombin by using CdS nanocrystals (CdS NCs), gold nanoparticles (AuNPs) and luminol is investigated in this work. Thrombin is detected by three methods. One is called the quenching method. It is based on the quenching effect of AuNPs on the yellow fluorescence of CdS NCs (with excitation/emission wavelengths of 355/550 nm) when placed adjacent to CdS NCs. The second method (called amplification method) is based on an amplification mechanism in which the plasmonics on the AuNPs enhance the emission of CdS NCs through distance related Förster resonance energy transfer (FRET). The third method is ratiometric and based on the emission by two luminophores, viz. CdS NCs and luminol. In this method, by increasing the concentration of thrombin, the intensity of CdS NCs decreases, while that of luminol increases. The results showed that ratiometric method was most sensitive (with an LOD of 500 fg.mL-1), followed by the amplification method (6.5 pg.mL-1) and the quenching method (92 pg.mL-1). Hence, the latter is less useful. Graphical abstract Schematic representation of three different methods (quenching, amplification and ratiometric) were applied for detection of thrombin via aptasensor. The CdS nanocrystals, streptavidin (Str) coated AuNPs and also Str-luminol coated AuNPs were used for the construction steps of the electrochemiluminescence (ECL)-based biosensor.

BackgroundIn standard high throughput sequencing analysis, genetic variants are not assigned to a homologous chromosome of origin. This process, called haplotype phasing, can reveal information important for understanding the relationship between genetic variants and biological phenotypes. For example, in genes that carry multiple heterozygous missense variants, phasing resolves whether one or both gene copies are altered. Here, we present a novel approach to phasing variants that takes advantage of unique properties of paired tumor:normal sequencing data from cancer studies.ResultsVAF phasing uses changes in variant allele frequency (VAF) between tumor and normal samples in regions of somatic chromosomal gain or loss to phase germline variants. We apply VAF phasing to 6180 samples from the Cancer Genome Atlas (TCGA) and demonstrate that our method is highly concordant with other standard phasing methods, and can phase an average of 33% more variants than other read-backed phasing methods. Using variant annotation tools designed to score gene haplotypes, we find a suggestive association between carrying multiple missense variants in a single copy of a cancer predisposition gene and earlier age of cancer diagnosis.ConclusionsVAF phasing exploits unique properties of tumor genomes to increase the number of germline variants that can be phased over standard read-backed methods in paired tumor:normal samples. Our phase-informed association testing results call attention to the need to develop more tools for assessing the joint effect of multiple genetic variants.

Pre-mRNA splicing is regulated through multiple trans-acting splicing factors. These regulators interact with the pre-mRNA at intronic and exonic positions. Given that most exons are protein coding, the evolution of exons must be modulated by a combination of selective coding and splicing pressures. It has previously been demonstrated that selective splicing pressures are more easily deconvoluted when phylogenetic comparisons are made for exons of identical size, suggesting that exon size-filtered sequence alignments may improve identification of nucleotides evolved to mediate efficient exon ligation. To test this hypothesis, an exon size database was created, filtering 76 vertebrate sequence alignments based on exon size conservation. In addition to other genomic parameters, such as splice-site strength, gene position, or flanking intron length, this database permits the identification of exons that are size- and/or sequence-conserved. Highly size-conserved exons are always sequence-conserved. However, sequence conservation does not necessitate exon size conservation. Our analysis identified evolutionarily young exons and demonstrated that length conservation is a strong predictor of alternative splicing. A published data set of approximately 5000 exonic SNPs associated with disease was analyzed to test the hypothesis that exon size-filtered sequence comparisons increase detection of splice-altering nucleotides. Improved splice predictions could be achieved when mutations occur at the third codon position, especially when a mutation decreases exon inclusion efficiency. The results demonstrate that coding pressures dominate nucleotide composition at invariable codon positions and that exon size-filtered sequence alignments permit identification of splice-altering nucleotides at wobble positions.

BackgroundResistance to immune checkpoint inhibitors (ICIs) has been linked to local immunosuppression independent of major ICI targets (e.g., PD-1). Clinical experience with response prediction based on PD-L1 expression suggests that other factors influence sensitivity to ICIs in non-small cell lung cancer (NSCLC) patients.MethodsTumor specimens from 120 NSCLC patients from 10 institutions were evaluated for PD-L1 expression by immunohistochemistry, and global proliferative profile by targeted RNA-seq.ResultsCell proliferation, derived from the mean expression of 10 proliferation-associated genes (namely BUB1, CCNB2, CDK1, CDKN3, FOXM1, KIAA0101, MAD2L1, MELK, MKI67, and TOP2A), was identified as a marker of response to ICIs in NSCLC. Poorly, moderately, and highly proliferative tumors were somewhat equally represented in NSCLC, with tumors with the highest PD-L1 expression being more frequently moderately proliferative as compared to lesser levels of PD-L1 expression. Proliferation status had an impact on survival in patients with both PD-L1 positive and negative tumors. There was a significant survival advantage for moderately proliferative tumors compared to their combined highly/poorly counterparts (p = 0.021). Moderately proliferative PD-L1 positive tumors had a median survival of 14.6 months that was almost twice that of PD-L1 negative highly/poorly proliferative at 7.6 months (p = 0.028). Median survival in moderately proliferative PD-L1 negative tumors at 12.6 months was comparable to that of highly/poorly proliferative PD-L1 positive tumors at 11.5 months, but in both instances less than that of moderately proliferative PD-L1 positive tumors. Similar to survival, proliferation status has impact on disease control (DC) in patients with both PD-L1 positive and negative tumors. Patients with moderately versus those with poorly or highly proliferative tumors have a superior DC rate when combined with any classification schema used to score PD-L1 as a positive result (i.e., TPS ≥ 50% or ≥ 1%), and best displayed by a DC rate for moderately proliferative tumors of no less than 40% for any classification of PD-L1 as a negative result. While there is an over representation of moderately proliferative tumors as PD-L1 expression increases this does not account for the improved survival or higher disease control rates seen in PD-L1 negative tumors.ConclusionsCell proliferation is potentially a new biomarker of response to ICIs in NSCLC and is applicable to PD-L1 negative tumors.

BackgroundHomozygous Familial Hypercholesterolemia (HoFH) is an inherited recessive condition associated with extremely high levels of low-density lipoprotein (LDL) cholesterol in affected individuals. It is usually caused by homozygous or compound heterozygous functional mutations in the LDL receptor (LDLR). A number of mutations causing FH have been reported in literature and such genetic heterogeneity presents great challenges for disease diagnosis.ObjectiveWe aim to determine the likely genetic defects responsible for three cases of pediatric HoFH in two kindreds.MethodsWe applied whole exome sequencing (WES) on the two probands to determine the likely functional variants among candidate FH genes. We additionally applied 10x Genomics (10xG) Linked-Reads whole genome sequencing (WGS) on one of the kindreds to identify potentially deleterious structural variants (SVs) underlying HoFH. A PCR-based screening assay was also established to detect the LDLR structural variant in a cohort of 641 patients with elevated LDL.ResultsIn the Caucasian kindred, the FH homozygosity can be attributed to two compound heterozygous LDLR damaging variants, an exon 12 p.G592E missense mutation and a novel 3kb exon 1 deletion. By analyzing the 10xG phased data, we ascertained that this deletion allele was most likely to have originated from a Russian ancestor. In the Mexican kindred, the strikingly elevated LDL cholesterol level can be attributed to a homozygous frameshift LDLR variant p.E113fs.ConclusionsWhile the application of WES can provide a cost-effective way of identifying the genetic causes of FH, it often lacks sensitivity for detecting structural variants. Our finding of the LDLR exon 1 deletion highlights the broader utility of Linked-Read WGS in detecting SVs in the clinical setting, especially when HoFH patients remain undiagnosed after WES.

BACKGROUND:The 3-dimensional (3D) conformation of chromatin inside the nucleus is integral to a variety of nuclear processes including transcriptional regulation, DNA replication, and DNA damage repair. Aberrations in 3D chromatin conformation have been implicated in developmental abnormalities and cancer. Despite the importance of 3D chromatin conformation to cellular function and human health, little is known about how 3D chromatin conformation varies in the human population, or whether DNA sequence variation between individuals influences 3D chromatin conformation. RESULTS:To address these questions, we perform Hi-C on lymphoblastoid cell lines from 20 individuals. We identify thousands of regions across the genome where 3D chromatin conformation varies between individuals and find that this variation is often accompanied by variation in gene expression, histone modifications, and transcription factor binding. Moreover, we find that DNA sequence variation influences several features of 3D chromatin conformation including loop strength, contact insulation, contact directionality, and density of local cis contacts. We map hundreds of quantitative trait loci associated with 3D chromatin features and find evidence that some of these same variants are associated at modest levels with other molecular phenotypes as well as complex disease risk. CONCLUSION:Our results demonstrate that common DNA sequence variants can influence 3D chromatin conformation, pointing to a more pervasive role for 3D chromatin conformation in human phenotypic variation than previously recognized.

BackgroundRhodosporidium toruloides is a promising host for the production of bioproducts from lignocellulosic biomass. A key prerequisite for efficient pathway engineering is the availability of robust genetic tools and resources. However, there is a lack of characterized promoters to drive expression of heterologous genes for strain engineering in R. toruloides.ResultsThis data describes a set of native R. toruloides promoters, characterized over time in four different media commonly used for cultivation of this yeast. The promoter sequences were selected using transcriptional analysis and several of them were found to drive expression bidirectionally. Promoter expression strength was determined by measurement of EGFP and mRuby2 reporters by flow cytometry. A total of 20 constitutive promoters (12 monodirectional and 8 bidirectional) were found, and are expected to be of potential value for genetic engineering of R. toruloides.ConclusionsA set of robust and constitutive promoters to facilitate genetic engineering of R. toruloides is presented here, ranging from a promoter previously used for this purpose (P7, glyceraldehyde 3-phosphate dehydrogenase, GAPDH) to stronger monodirectional (e.g., P15, mitochondrial adenine nucleotide translocator, ANT) and bidirectional (e.g., P9 and P9R, histones H3 and H4, respectively) promoters. We also identified promoters that may be useful for specific applications such as late-stage expression (e.g., P3, voltage-dependent anion channel protein 2, VDAC2). This set of characterized promoters significantly expands the range of engineering tools available for this yeast and can be applied in future metabolic engineering studies.

BackgroundInexpensive high-throughput DNA sequencing has democratized access to genetic information for most organisms so that research utilizing a genome or transcriptome of an organism is not limited to model systems. However, the quality of the assemblies of sampled genomes can vary greatly which hampers utility for comparisons and meaningful interpretation. The uncertainty of the completeness of a given genome sequence can limit feasibility of asserting patterns of high rates of gene loss reported in many lineages.ResultsWe propose a computational framework and sequence resource for assessing completeness of fungal genomes called FGMP (Fungal Genome Mapping Project). Our approach is based on evolutionary conserved sets of proteins and DNA elements and is applicable to various types of genomic data. We present a comparison of FGMP and state-of-the-art methods for genome completeness assessment utilizing 246 genome assemblies of fungi. We discuss genome assembly improvements/degradations in 57 cases where assemblies have been updated, as recorded by NCBI assembly archive.ConclusionFGMP is an accurate tool for quantifying level of completion from fungal genomic data. It is particularly useful for non-model organisms without reference genomes and can be used directly on unassembled reads, which can help reducing genome sequencing costs.

The CRISPR/Cas system is a highly specific genome editing tool capable of distinguishing alleles differing by even a single base pair. Target sites might carry genetic variations that are not distinguishable by sgRNA designing tools based on one reference genome. AlleleAnalyzer is an open-source software that incorporates single-nucleotide variants and short insertions and deletions to design sgRNAs for precisely editing 1 or multiple haplotypes of a sequenced genome, currently supporting 11 Cas proteins. It also leverages patterns of shared genetic variation to optimize sgRNA design for different human populations. AlleleAnalyzer is available at https://github.com/keoughkath/AlleleAnalyzer .

Human ADAR3 is a catalytically inactive member of the Adenosine Deaminase Acting on RNA (ADAR) protein family, whose active members catalyze A-to-I RNA editing in metazoans. Until now, the reasons for the catalytic incapability of ADAR3 has not been defined and its biological function rarely explored. Yet, its exclusive expression in the brain and involvement in learning and memory suggest a central role in the nervous system. Here we describe the engineering of a catalytically active ADAR3 enzyme using a combination of computational design and functional screening. Five mutations (A389V, V485I, E527Q, Q549R and Q733D) engender RNA deaminase in human ADAR3. By way of its catalytic activity, the ADAR3 pentamutant was used to identify potential binding targets for wild type ADAR3 in a human glioblastoma cell line. Novel ADAR3 binding sites discovered in this manner include the 3'-UTRs of the mRNAs encoding early growth response 1 (EGR1) and dual specificity phosphatase 1 (DUSP1); both known to be activity-dependent immediate early genes that respond to stimuli in the brain. Further studies reveal that the wild type ADAR3 protein can regulate transcript levels for DUSP1 and EGR1, suggesting a novel role ADAR3 may play in brain function.

RNA-chromatin interactions represent an important aspect of the transcriptional regulation of genes and transposable elements. However, analyses of chromatin-associated RNAs (caRNAs) are often limited to one caRNA at a time. Here, we describe the iMARGI (in situ mapping of RNA-genome interactome) technique, which is used to discover caRNAs and reveal their respective genomic interaction loci. iMARGI starts with in situ crosslinking and genome fragmentation, followed by converting each proximal RNA-DNA pair into an RNA-linker-DNA chimeric sequence. These chimeric sequences are subsequently converted into a sequencing library suitable for paired-end sequencing. A standardized bioinformatic software package, iMARGI-Docker, is provided to decode the paired-end sequencing data into caRNA-DNA interactions. Compared to its predecessor MARGI (mapping RNA-genome interactions), the number of input cells for iMARGI is 3-5 million (a 100-fold reduction), experimental time is reduced, and clear checkpoints have been established. It takes a few hours a day and a total of 8 d to complete the construction of an iMARGI sequencing library and 1 d to carry out data processing with iMARGI-Docker.

Gliomatosis peritonei is a rare pathologic finding that is associated with ovarian teratomas and malignant mixed germ cell tumors. The occurrence of gliomatosis as a mature glial implant can impart an improved prognosis to patients with immature ovarian teratoma, making prompt and accurate diagnosis important. We describe a case of recurrent immature teratoma in a 10-yr-old female patient, in which comparative analysis of the RNA sequencing gene expression data from the patient's tumor was used effectively to aid in the diagnosis of gliomatosis peritonei.

Feline infectious peritonitis (FIP) is a highly fatal disease caused by a virulent feline coronavirus in domestic and wild cats. We have previously reported the synthesis of potent coronavirus 3C-like protease (3CLpro) inhibitors and the efficacy of a protease inhibitor, GC376, in client-owned cats with FIP. In this study, we studied the effect of the amino acid changes in 3CLpro of feline coronavirus from a feline patient who received antiviral treatment for prolonged duration. We generated recombinant 3CLpro containing the identified amino acid changes (N25S, A252S or K260 N) and determined their susceptibility to protease inhibitors in the fluorescence resonance energy transfer assay. The assay showed that N25S in 3CLpro confers a small change (up to 1.68-fold increase in the 50% inhibitory concentration) in susceptibility to GC376, but other amino acid changes do not affect susceptibility. Modelling of 3CLpro carrying the amino acid changes was conducted to probe the structural basis for these findings. The results of this study may explain the observed absence of clinical resistance to the long-term antiviral treatment in the patients.

Eye regression occurs across cave-dwelling populations of many species and is often coupled with a decrease or loss in eye function. Teleost fishes are among the few vertebrates to undergo widespread colonization of caves and often exhibit eye regression with blindness. Cave populations of the poeciliid fish Poecilia mexicana (cave molly) exhibit reduced-albeit functional-eyes, offering the opportunity to investigate partial eye regression. We sequenced eye transcriptomes of cave and surface populations of P. mexicana to identify differentially expressed genes that potentially underlie eye regression in cave mollies. We identified 28 significantly differentially expressed genes, 20 of which were directly related to light sensitivity, eye structure and visual signaling. Twenty-six of these genes were downregulated in cave compared to surface populations. Functional enrichment analysis revealed eye-related gene ontologies that were under-represented in cave mollies. In addition, a set of co-expressed genes related to vision and circadian rhythm was correlated with habitat type (cave versus surface). Our study suggests that differential gene expression plays a key role in the beginning evolutionary stages of eye regression in P. mexicana, shedding further light on regressive evolution in cavefish.

Fine mapping and validation of genes causing β cell failure from susceptibility loci identified in type 2 diabetes genome-wide association studies (GWAS) poses a significant challenge. The VPS13C-C2CD4A-C2CD4B locus on chromosome 15 confers diabetes susceptibility in every ethnic group studied to date. However, the causative gene is unknown. FoxO1 is involved in the pathogenesis of β cell dysfunction, but its link to human diabetes GWAS has not been explored. Here we generated a genome-wide map of FoxO1 superenhancers in chemically identified β cells using 2-photon live-cell imaging to monitor FoxO1 localization. When parsed against human superenhancers and GWAS-derived diabetes susceptibility alleles, this map revealed a conserved superenhancer in C2CD4A, a gene encoding a β cell/stomach-enriched nuclear protein of unknown function. Genetic ablation of C2cd4a in β cells of mice phenocopied the metabolic abnormalities of human carriers of C2CD4A-linked polymorphisms, resulting in impaired insulin secretion during glucose tolerance tests as well as hyperglycemic clamps. C2CD4A regulates glycolytic genes, and notably represses key β cell "disallowed" genes, such as lactate dehydrogenase A We propose that C2CD4A is a transcriptional coregulator of the glycolytic pathway whose dysfunction accounts for the diabetes susceptibility associated with the chromosome 15 GWAS locus.

The black morel (Morchella importuna Kuo, O'Donnell and Volk) was once an uncultivable wild mushroom, until the development of exogenous nutrient bag (ENB), making its agricultural production quite feasible and stable. To date, how the nutritional acquisition of the morel mycelium is fulfilled to trigger its fruiting remains unknown. To investigate the mechanisms involved in ENB decomposition, the genome of a cultivable morel strain (M. importuna SCYDJ1-A1) was sequenced and the genes coding for the decay apparatus were identified. Expression of the encoded carbohydrate-active enzymes (CAZymes) was then analyzed by metatranscriptomics and metaproteomics in combination with biochemical assays. The results show that a diverse set of hydrolytic and redox CAZymes secreted by the morel mycelium is the main force driving the substrate decomposition. Plant polysaccharides such as starch and cellulose present in ENB substrate (wheat grains plus rice husks) were rapidly degraded, whereas triglycerides were accumulated initially and consumed later. ENB decomposition led to a rapid increase in the organic carbon content in the surface soil of the mushroom bed, which was thereafter consumed during morel fruiting. In contrast to the high carbon consumption, no significant acquisition of nitrogen was observed. Our findings contribute to an increasingly detailed portrait of molecular features triggering morel fruiting.

Electrochemical sensors are major players in the race for improved molecular diagnostics due to their convenience, temporal resolution, manufacturing scalability, and their ability to support real-time measurements. This is evident in the ever-increasing number of health-related electrochemical sensing platforms, ranging from single-measurement point-of-care devices to wearable devices supporting immediate and continuous monitoring. In support of the need for such systems to rapidly process large data volumes, we describe here an open-source, easily customizable, multiplatform compatible program for the real-time control, processing, and visualization of electrochemical data. The software's architecture is modular and fully documented, allowing the easy customization of the code to support the processing of voltammetric (e.g., square-wave and cyclic) and chronoamperometric data. The program, which we have called Software for the Analysis and Continuous Monitoring of Electrochemical Systems (SACMES), also includes a graphical interface allowing the user to easily change analysis parameters (e.g., signal/noise processing, baseline correction) in real-time. To demonstrate the versatility of SACMES we use it here to analyze the real-time data output by (1) the electrochemical, aptamer-based measurement of a specific small-molecule target, (2) a monoclonal antibody-detecting DNA-scaffold sensor, and (3) the determination of the folding thermodynamics of an electrode-attached, redox-reporter-modified protein.

Most current methods for detecting natural selection from DNA sequence data are limited in that they are either based on summary statistics or a composite likelihood, and as a consequence, do not make full use of the information available in DNA sequence data. We here present a new importance sampling approach for approximating the full likelihood function for the selection coefficient. Our method CLUES treats the ancestral recombination graph (ARG) as a latent variable that is integrated out using previously published Markov Chain Monte Carlo (MCMC) methods. The method can be used for detecting selection, estimating selection coefficients, testing models of changes in the strength of selection, estimating the time of the start of a selective sweep, and for inferring the allele frequency trajectory of a selected or neutral allele. We perform extensive simulations to evaluate the method and show that it uniformly improves power to detect selection compared to current popular methods such as nSL and SDS, and can provide reliable inferences of allele frequency trajectories under many conditions. We also explore the potential of our method to detect extremely recent changes in the strength of selection. We use the method to infer the past allele frequency trajectory for a lactase persistence SNP (MCM6) in Europeans. We also infer the trajectory of a SNP (EDAR) in Han Chinese, finding evidence that this allele's age is much older than previously claimed. We also study a set of 11 pigmentation-associated variants. Several genes show evidence of strong selection particularly within the last 5,000 years, including ASIP, KITLG, and TYR. However, selection on OCA2/HERC2 seems to be much older and, in contrast to previous claims, we find no evidence of selection on TYRP1.

Non-homologous end-joining (NHEJ) plays an important role in double-strand break (DSB) repair of DNA. Recent studies have shown that the error patterns of NHEJ are strongly biased by sequence context, but these studies were based on relatively few templates. To investigate this more thoroughly, we systematically profiled ∼1.16 million independent mutational events resulting from CRISPR/Cas9-mediated cleavage and NHEJ-mediated DSB repair of 6872 synthetic target sequences, introduced into a human cell line via lentiviral infection. We find that: (i) insertions are dominated by 1 bp events templated by sequence immediately upstream of the cleavage site, (ii) deletions are predominantly associated with microhomology and (iii) targets exhibit variable but reproducible diversity with respect to the number and relative frequency of the mutational outcomes to which they give rise. From these data, we trained a model that uses local sequence context to predict the distribution of mutational outcomes. Exploiting the bias of NHEJ outcomes towards microhomology mediated events, we demonstrate the programming of deletion patterns by introducing microhomology to specific locations in the vicinity of the DSB site. We anticipate that our results will inform investigations of DSB repair mechanisms as well as the design of CRISPR/Cas9 experiments for diverse applications including genome-wide screens, gene therapy, lineage tracing and molecular recording.

The 17th International HLA and Immunogenetics Workshop (IHIW) conducted a project entitled "The Study of Haplotypes in Families by NGS HLA". We investigated the HLA haplotypes of 1017 subjects in 263 nuclear families sourced from five US clinical immunogenetics laboratories, primarily as part of the evaluation of related donor candidates for hematopoietic stem cell and solid organ transplantation. The parents in these families belonged to five broad groups - African (72 parents), Asian (115), European (210), Hispanic (118) and "Other" (11). High-resolution HLA genotypes were generated for each subject using next-generation sequencing (NGS) HLA typing systems. We identified the HLA haplotypes in each family using HaplObserve, software that builds haplotypes in families by reviewing HLA allele segregation from parents to children. We calculated haplotype frequencies within each broad group, by treating the parents in each family as unrelated individuals. We also calculated standard measures of global linkage disequilibrium (LD) and conditional asymmetric LD for each ethnic group, and used untruncated and two-field allele names to investigate LD patterns. Finally we demonstrated the utility of consensus DNA sequences in identifying novel variants, confirming them using HLA allele segregation at the DNA sequence level.

Bacteriophages are abundant within the human gastrointestinal tract, yet their interactions with gut bacteria remain poorly understood, particularly with respect to CRISPR-Cas immunity. Here, we show that the type I-C CRISPR-Cas system in the prevalent gut Actinobacterium Eggerthella lenta is transcribed and sufficient for specific targeting of foreign and chromosomal DNA. Comparative analyses of E. lenta CRISPR-Cas systems across (meta)genomes revealed 2 distinct clades according to cas sequence similarity and spacer content. We assembled a human virome database (HuVirDB), encompassing 1,831 samples enriched for viral DNA, to identify protospacers. This revealed matches for a majority of spacers, a marked increase over other databases, and uncovered "hyper-targeted" phage sequences containing multiple protospacers targeted by several E. lenta strains. Finally, we determined the positional mismatch tolerance of observed spacer-protospacer pairs. This work emphasizes the utility of merging computational and experimental approaches for determining the function and targets of CRISPR-Cas systems.

Notholaenids are an unusual group of ferns that have adapted to, and diversified within, the deserts of Mexico and the southwestern United States. With approximately 40 species, this group is noted for being desiccation-tolerant and having "farina"-powdery exudates of lipophilic flavonoid aglycones-that occur on both the gametophytic and sporophytic phases of their life cycle. The most recent circumscription of notholaenids based on plastid markers surprisingly suggests that several morphological characters, including the expression of farina, are homoplasious. In a striking case of convergence, Notholaena standleyi appears to be distantly related to core Notholaena, with several taxa not before associated with Notholaena nested between them. Such conflicts can be due to morphological homoplasy resulting from adaptive convergence or, alternatively, the plastid phylogeny itself might be misleading, diverging from the true species tree due to incomplete lineage sorting, hybridization, or other factors. In this study, we present a species phylogeny for notholaenid ferns, using four low-copy nuclear loci and concatenated data from three plastid loci. A total of 61 individuals (49 notholaenids and 12 outgroup taxa) were sampled, including 31 out of 37 recognized notholaenid species. The homeologous/allelic nuclear sequences were retrieved using PacBio sequencing and the PURC bioinformatics pipeline. Each dataset was first analyzed individually using maximum likelihood and Bayesian inference, and the species phylogeny was inferred using *BEAST. Although we observed several incongruences between the nuclear and plastid phylogenies, our principal results are broadly congruent with previous inferences based on plastid data. By mapping the presence of farina and their biochemical constitutions on our consensus phylogenetic tree, we confirmed that the characters are indeed homoplastic and have complex evolutionary histories. Hybridization among recognized species of the notholaenid clade appears to be relatively rare compared to that observed in other well-studied fern genera.

Pacemaker cardiomyocytes that create the sinoatrial node are essential for the initiation and maintenance of proper heart rhythm. However, illuminating developmental cues that direct their differentiation has remained particularly challenging due to the unclear cellular origins of these specialized cardiomyocytes. By discovering the origins of pacemaker cardiomyocytes, we reveal an evolutionarily conserved Wnt signaling mechanism that coordinates gene regulatory changes directing mesoderm cell fate decisions, which lead to the differentiation of pacemaker cardiomyocytes. We show that in zebrafish, pacemaker cardiomyocytes derive from a subset of Nkx2.5+ mesoderm that responds to canonical Wnt5b signaling to initiate the cardiac pacemaker program, including activation of pacemaker cell differentiation transcription factors Isl1 and Tbx18 and silencing of Nkx2.5. Moreover, applying these developmental findings to human pluripotent stem cells (hPSCs) notably results in the creation of hPSC-pacemaker cardiomyocytes, which successfully pace three-dimensional bioprinted hPSC-cardiomyocytes, thus providing potential strategies for biological cardiac pacemaker therapy.

Localized arrays of proteins cooperatively assemble onto chromosomes to control DNA activity in many contexts. Binding cooperativity is often mediated by specific protein-protein interactions, but cooperativity through DNA structure is becoming increasingly recognized as an additional mechanism. During the site-specific DNA recombination reaction that excises phage λ from the chromosome, the bacterial DNA architectural protein Fis recruits multiple λ-encoded Xis proteins to the attR recombination site. Here, we report X-ray crystal structures of DNA complexes containing Fis + Xis, which show little, if any, contacts between the two proteins. Comparisons with structures of DNA complexes containing only Fis or Xis, together with mutant protein and DNA binding studies, support a mechanism for cooperative protein binding solely by DNA allostery. Fis binding both molds the minor groove to potentiate insertion of the Xis β-hairpin wing motif and bends the DNA to facilitate Xis-DNA contacts within the major groove. The Fis-structured minor groove shape that is optimized for Xis binding requires a precisely positioned pyrimidine-purine base-pair step, whose location has been shown to modulate minor groove widths in Fis-bound complexes to different DNA targets.

Clinical persistence of Chlamydia trachomatis (Ct) sexually transmitted infections (STIs) is a major public health concern. In vitro persistence is known to develop through interferon gamma (IFN-γ) induction of indoleamine 2,3-dioxygenase (IDO), which catabolizes tryptophan, an essential amino acid for Ct replication. The organism can recover from persistence by synthesizing tryptophan from indole, a substrate for the enzyme tryptophan synthase. The majority of Ct strains, except for reference strain B/TW-5/OT, contain an operon comprised of α and β subunits that encode TrpA and TrpB, respectively, and form a functional αββα tetramer. However, trpA mutations in ocular Ct strains, which are responsible for the blinding eye disease known as trachoma, abrogate tryptophan synthesis from indole. We examined serial urogenital samples from a woman who had recurrent Ct infections over 4 years despite antibiotic treatment. The Ct isolates from each infection episode were genome sequenced and analyzed for phenotypic, structural, and functional characteristics. All isolates contained identical mutations in trpA and developed aberrant bodies within intracellular inclusions, visualized by transmission electron microscopy, even when supplemented with indole following IFN-γ treatment. Each isolate displayed an altered αββα structure, could not synthesize tryptophan from indole, and had significantly lower trpBA expression but higher intracellular tryptophan levels compared with those of reference Ct strain F/IC-Cal3. Our data indicate that emergent mutations in the tryptophan operon, which were previously thought to be restricted only to ocular Ct strains, likely resulted in in vivo persistence in the described patient and represents a novel host-pathogen adaptive strategy for survival.IMPORTANCEChlamydia trachomatis (Ct) is the most common sexually transmitted bacterium with more than 131 million cases occurring annually worldwide. Ct infections are often asymptomatic, persisting for many years despite treatment. In vitro recovery from persistence occurs when indole is utilized by the organism's tryptophan synthase to synthesize tryptophan, an essential amino acid for replication. Ocular but not urogenital Ct strains contain mutations in the synthase that abrogate tryptophan synthesis. Here, we discovered that the genomes of serial isolates from a woman with recurrent, treated Ct STIs over many years were identical with a novel synthase mutation. This likely allowed long-term in vivo persistence where active infection resumed only when tryptophan became available. Our findings indicate an emerging adaptive host-pathogen evolutionary strategy for survival in the urogenital tract that will prompt the field to further explore chlamydial persistence, evaluate the genetics of mutant Ct strains and fitness within the host, and their implications for disease pathogenesis.

The faster a disease can be diagnosed, the sooner effective treatment can be initiated, motivating a drive to replace standard laboratory techniques with point-of-care technologies that return answers in minutes rather than hours. Thus motivated, we describe the development of an E-DNA scaffold sensor for the rapid and convenient measurement of antibodies diagnostic of syphilis. To achieve this (and in contrast to previous sensors of this class, which relied on single, linear epitopes for detection), we utilized a near full-length antigen as the sensor's recognition element, allowing us to simultaneously display multiple epitopes. The resultant sensor is able to detect antibodies against Treponema pallidum pallidum, the causative agent of syphilis, at clinically relevant concentrations in samples in less than 10 min. Preliminary results obtained using sero-positive and sero-negative human samples suggest the clinical sensitivity and specificity of the approach compare well to current gold-standard tests, while being simple and rapid enough to deploy at the point of care.

Mammalian sperm RNA is increasingly recognized as an additional source of paternal hereditary information beyond DNA. Environmental inputs, including an unhealthy diet, mental stresses and toxin exposure, can reshape the sperm RNA signature and induce offspring phenotypes that relate to paternal environmental stressors. Our understanding of the categories of sperm RNAs (such as tRNA-derived small RNAs, microRNAs, ribosomal RNA-derived small RNAs and long non-coding RNAs) and associated RNA modifications is expanding and has begun to reveal the functional diversity and information capacity of these molecules. However, the coding mechanism endowed by sperm RNA structures and by RNA interactions with DNA and other epigenetic factors remains unknown. How sperm RNA-encoded information is decoded in early embryos to control offspring phenotypes also remains unclear. Complete deciphering of the 'sperm RNA code' with regard to metabolic control could move the field towards translational applications and precision medicine, and this may lead to prevention of intergenerational transmission of obesity and type 2 diabetes mellitus susceptibility.

MicroRNAs (miRNAs) are important posttranscriptional regulators of metabolism and energy homeostasis. Dysregulation of certain miRNAs in the liver has been shown to contribute to the pathogenesis of Type 2 diabetes (T2D), in part by impairing hepatic insulin sensitivity. By small RNA-sequencing analysis, we identified seven hepatic miRNAs (including miR-29b) that are consistently aberrantly expressed across five different rodent models of metabolic dysfunction that share the feature of insulin resistance (IR). We also showed that hepatic miR-29b exhibits persistent dysregulation during disease progression in a rat model of diabetes, UCD-T2DM. Furthermore, we observed that hepatic levels of miR-29 family members are attenuated by interventions known to improve IR in rodent and rhesus macaque models. To examine the function of the miR-29 family in modulating insulin sensitivity, we used locked nucleic acid (LNA) technology and demonstrated that acute in vivo suppression of the miR-29 family in adult mice leads to significant reduction of fasting blood glucose (in both chow-fed lean and high-fat diet-fed obese mice) and improvement in insulin sensitivity (in chow-fed lean mice). We carried out whole transcriptome studies and uncovered candidate mechanisms, including regulation of DNA methyltransferase 3a (Dnmt3a) and the hormone-encoding gene Energy homeostasis associated (Enho). In sum, we showed that IR/T2D is linked to dysregulation of hepatic miR-29b across numerous models and that acute suppression of the miR-29 family in adult mice leads to improved glycemic control. Future studies should investigate the therapeutic utility of miR-29 suppression in different metabolic disease states.Enho; insulin resistance; liver; microRNA-29 (miR-29); UCD-T2DM.

Allohexaploid (2n = 6x = 42) intermediate wheatgrass (Thinopyrum intermedium), abbreviated IWG, is an outcrossing perennial grass belonging to the tertiary gene pool of wheat. Perenniality would be valuable option for grain production, but attempts to introgress this complex trait from wheat-Thinopyrum hybrids have not been commercially successful. Efforts to breed IWG itself as a dual-purpose forage and grain crop have demonstrated useful progress and applications, but grain yields are significantly less than wheat. Therefore, genetic and physical maps have been developed to accelerate domestication of IWG. Herein, these maps were used to identify quantitative trait loci (QTLs) and candidate genes associated with IWG grain production traits in a family of 266 full-sib progenies derived from two heterozygous parents, M26 and M35. Transgressive segregation was observed for 17 traits related to seed size, shattering, threshing, inflorescence capacity, fertility, stem size, and flowering time. A total of 111 QTLs were detected in 36 different regions using 3826 genotype-by-sequence markers in 21 linkage groups. The most prominent QTL had a LOD score of 15 with synergistic effects of 29% and 22% over the family means for seed retention and percentage of naked seeds, respectively. Many QTLs aligned with one or more IWG gene models corresponding to 42 possible domestication orthogenes including the wheat Q and RHT genes. A cluster of seed-size and fertility QTLs showed possible alignment to a putative Z self-incompatibility gene, which could have detrimental grain-yield effects when genetic variability is low. These findings elucidate pathways and possible hurdles in the domestication of IWG.

DLX transcription factors (TFs) are master regulators of the developing vertebrate brain, driving forebrain GABAergic neuronal differentiation. Ablation of Dlx1&2 alters expression of genes that are critical for forebrain GABAergic development. We integrated epigenomic and transcriptomic analyses, complemented with in situ hybridization (ISH), and in vivo and in vitro studies of regulatory element (RE) function. This revealed the DLX-organized gene regulatory network at genomic, cellular, and spatial levels in mouse embryonic basal ganglia. DLX TFs perform dual activating and repressing functions; the consequences of their binding were determined by the sequence and genomic context of target loci. Our results reveal and, in part, explain the paradox of widespread DLX binding contrasted with a limited subset of target loci that are sensitive at the epigenomic and transcriptomic level to Dlx1&2 ablation. The regulatory properties identified here for DLX TFs suggest general mechanisms by which TFs orchestrate dynamic expression programs underlying neurodevelopment.

The California poppy (Eschscholzia californica) is renowned for its brilliant golden-orange flowers, though white petal variants have been described. By whole-transcriptome sequencing, we have discovered in multiple white petal varieties a single deletion leading to altered splicing and C-terminal truncation of phytoene synthase (PSY), a key enzyme in carotenoid biosynthesis. Our findings underscore the diverse roles of phytoene synthase in shaping horticultural traits, and resolve a longstanding mystery of the regaled golden poppy.

The availability of genome-wide epigenomic datasets enables in-depth studies of epigenetic modifications and their relationships with chromatin structures and gene expression. Various alignment tools have been developed to align nucleotide or protein sequences in order to identify structurally similar regions. However, there are currently no alignment methods specifically designed for comparing multi-track epigenomic signals and detecting common patterns that may explain functional or evolutionary similarities. We propose a new local alignment algorithm, EpiAlign, designed to compare chromatin state sequences learned from multi-track epigenomic signals and to identify locally aligned chromatin regions. EpiAlign is a dynamic programming algorithm that novelly incorporates varying lengths and frequencies of chromatin states. We demonstrate the efficacy of EpiAlign through extensive simulations and studies on the real data from the NIH Roadmap Epigenomics project. EpiAlign is able to extract recurrent chromatin state patterns along a single epigenome, and many of these patterns carry cell-type-specific characteristics. EpiAlign can also detect common chromatin state patterns across multiple epigenomes, and it will serve as a useful tool to group and distinguish epigenomic samples based on genome-wide or local chromatin state patterns.

The proliferation, differentiation and survival of mononuclear phagocytes depend on signals from the receptor for macrophage colony-stimulating factor, CSF1R. The mammalian Csf1r locus contains a highly conserved super-enhancer, the fms-intronic regulatory element (FIRE). Here we show that genomic deletion of FIRE in mice selectively impacts CSF1R expression and tissue macrophage development in specific tissues. Deletion of FIRE ablates macrophage development from murine embryonic stem cells. Csf1rΔFIRE/ΔFIRE mice lack macrophages in the embryo, brain microglia and resident macrophages in the skin, kidney, heart and peritoneum. The homeostasis of other macrophage populations and monocytes is unaffected, but monocytes and their progenitors in bone marrow lack surface CSF1R. Finally, Csf1rΔFIRE/ΔFIRE mice are healthy and fertile without the growth, neurological or developmental abnormalities reported in Csf1r-/- rodents. Csf1rΔFIRE/ΔFIRE mice thus provide a model to explore the homeostatic, physiological and immunological functions of tissue-specific macrophage populations in adult animals.

Mutations in the transcription factor FOXA1 define a unique subset of prostate cancers but the functional consequences of these mutations and whether they confer gain or loss of function is unknown1-9. Here, by annotating the landscape of FOXA1 mutations from 3,086 human prostate cancers, we define two hotspots in the forkhead domain: Wing2 (around 50% of all mutations) and the highly conserved DNA-contact residue R219 (around 5% of all mutations). Wing2 mutations are detected in adenocarcinomas at all stages, whereas R219 mutations are enriched in metastatic tumours with neuroendocrine histology. Interrogation of the biological properties of wild-type FOXA1 and fourteen FOXA1 mutants reveals gain of function in mouse prostate organoid proliferation assays. Twelve of these mutants, as well as wild-type FOXA1, promoted an exaggerated pro-luminal differentiation program, whereas two different R219 mutants blocked luminal differentiation and activated a mesenchymal and neuroendocrine transcriptional program. Assay for transposase-accessible chromatin using sequencing (ATAC-seq) of wild-type FOXA1 and representative Wing2 and R219 mutants revealed marked, mutant-specific changes in open chromatin at thousands of genomic loci and exposed sites of FOXA1 binding and associated increases in gene expression. Of note, ATAC-seq peaks in cells expressing R219 mutants lacked the canonical core FOXA1-binding motifs (GTAAAC/T) but were enriched for a related, non-canonical motif (GTAAAG/A), which was preferentially activated by R219-mutant FOXA1 in reporter assays. Thus, FOXA1 mutations alter its pioneering function and perturb normal luminal epithelial differentiation programs, providing further support for the role of lineage plasticity in cancer progression.

Determining the carbon sources for active microbial populations in the subsurface is a challenging but highly informative component of subsurface microbial ecology. This work developed a method to provide ecological insights into groundwater microbial communities by characterizing community RNA through its radiocarbon and ribosomal RNA (rRNA) signatures. RNA was chosen as the biomolecule of interest because rRNA constitutes the majority of RNA in prokaryotes, represents recently active organisms, and yields detailed taxonomic information. The method was applied to a groundwater filter collected from a shallow alluvial aquifer in Colorado. RNA was extracted, radiometrically dated, and the 16S rRNA was analyzed by RNA-Seq. The RNA had a radiocarbon signature (Δ14C) of -193.4 ± 5.6‰. Comparison of the RNA radiocarbon signature to those of potential carbon pools in the aquifer indicated that at least 51% of the RNA was derived from autotrophy, in close agreement with the RNA-Seq data, which documented the prevalence of autotrophic taxa, such as Thiobacillus and Gallionellaceae. Overall, this hybrid method for RNA analysis provided cultivation-independent information on the in-situ carbon sources of active subsurface microbes and reinforced the importance of autotrophy and the preferential utilization of dissolved over sedimentary organic matter in alluvial aquifers.

Alu retroelements, whose retrotransposition requires prior transcription by RNA polymerase III to generate Alu RNAs, represent the most numerous non-coding RNA (ncRNA) gene family in the human genome. Alu transcription is generally kept to extremely low levels by tight epigenetic silencing, but it has been reported to increase under different types of cell perturbation, such as viral infection and cancer. Alu RNAs, being able to act as gene expression modulators, may be directly involved in the mechanisms determining cellular behavior in such perturbed states. To directly address the regulatory potential of Alu RNAs, we generated IMR90 fibroblasts and HeLa cell lines stably overexpressing two slightly different Alu RNAs, and analyzed genome-wide the expression changes of protein-coding genes through RNA-sequencing. Among the genes that were upregulated or downregulated in response to Alu overexpression in IMR90, but not in HeLa cells, we found a highly significant enrichment of pathways involved in cell cycle progression and mitotic entry. Accordingly, Alu overexpression was found to promote transition from G1 to S phase, as revealed by flow cytometry. Therefore, increased Alu RNA may contribute to sustained cell proliferation, which is an important factor of cancer development and progression.

Transfer RNAs (tRNAs) are ubiquitous across the tree of life. Although tRNA structure is highly conserved, there is still significant variation in sequence features between clades, isotypes and even isodecoders. This variation not only impacts translation, but as shown by a variety of recent studies, nontranslation-associated functions are also sensitive to small changes in tRNA sequence. Despite the rapidly growing number of sequenced genomes, there is a lack of tools for both small- and large-scale comparative genomics analysis of tRNA sequence features. Here, we have integrated over 150 000 tRNAs spanning all domains of life into tRNAviz, a web application for exploring and visualizing tRNA sequence features. tRNAviz implements a framework for determining consensus sequence features and can generate sequence feature distributions by isotypes, clades and anticodons, among other tRNA properties such as score. All visualizations are interactive and exportable. The web server is publicly available at http://trna.ucsc.edu/tRNAviz/.

Sample multiplexing facilitates scRNA-seq by reducing costs and identifying artifacts such as cell doublets. However, universal and scalable sample barcoding strategies have not been described. We therefore developed MULTI-seq: multiplexing using lipid-tagged indices for single-cell and single-nucleus RNA sequencing. MULTI-seq reagents can barcode any cell type or nucleus from any species with an accessible plasma membrane. The method involves minimal sample processing, thereby preserving cell viability and endogenous gene expression patterns. When cells are classified into sample groups using MULTI-seq barcode abundances, data quality is improved through doublet identification and recovery of cells with low RNA content that would otherwise be discarded by standard quality-control workflows. We use MULTI-seq to track the dynamics of T-cell activation, perform a 96-plex perturbation experiment with primary human mammary epithelial cells and multiplex cryopreserved tumors and metastatic sites isolated from a patient-derived xenograft mouse model of triple-negative breast cancer.

Deletions in the 16.6 kb mitochondrial genome have been implicated in numerous disorders that often display muscular and/or neurological symptoms due to the high-energy demands of these tissues. We describe a catalogue of 4489 putative mitochondrial DNA (mtDNA) deletions, including their frequency and relative read rate, using a combinatorial approach of mitochondria-targeted PCR, next-generation sequencing, bioinformatics, post-hoc filtering, annotation, and validation steps. Our bioinformatics pipeline uses MapSplice, an RNA-seq splice junction detection algorithm, to detect and quantify mtDNA deletion breakpoints rather than mRNA splices. Analyses of 93 samples from postmortem brain and blood found (i) the 4977 bp 'common deletion' was neither the most frequent deletion nor the most abundant; (ii) brain contained significantly more deletions than blood; (iii) many high frequency deletions were previously reported in MitoBreak, suggesting they are present at low levels in metabolically active tissues and are not exclusive to individuals with diagnosed mitochondrial pathologies; (iv) many individual deletions (and cumulative metrics) had significant and positive correlations with age and (v) the highest deletion burdens were observed in major depressive disorder brain, at levels greater than Kearns-Sayre Syndrome muscle. Collectively, these data suggest the Splice-Break pipeline can detect and quantify mtDNA deletions at a high level of resolution.

Deoxyribozymes capable of catalyzing sequence-specific RNA cleavage have found broad applications in biotechnology, DNA computing and environmental sensing. Among these, deoxyribozyme 8-17 is the most common small DNA motif capable of catalyzing RNA cleavage. However, the extent to which other DNA molecules with similar catalytic motifs exist remains elusive. Here we report a novel RNA-cleaving deoxyribozyme called 10-12opt that functions with an equally small catalytic motif and an unusually short binding arm. This deoxyribozyme contains a 14-nucleotide catalytic core that preferentially catalyzes RNA cleavage at UN dinucleotide junctions (kobs = 0.9 h-1 for UU cleavage). Surprisingly, the left binding arm contains only three nucleotides and forms two canonical base pairs with the RNA substrate. Mutational analysis reveals that a riboguanosine residue 3-nucleotide downstream of cleavage site must not form canonical base pairing for the optimal catalysis, and this nucleobase likely participates in catalysis with its carbonyl O6 atom. Furthermore, we demonstrate that deoxyribozyme 10-12opt can be utilized to cleave certain microRNA sequences which are not preferentially cleaved by 8-17. Together, these results suggest that this novel RNA-cleaving deoxyribozyme forms a distinct catalytic structure than 8-17 and that sequence space may contain additional examples of DNA molecules that can cleave RNA at site-specific locations.

G protein-coupled receptors (GPCRs) are targets for ∼35% of approved drugs but only ∼15% of the ∼800 human GPCRs are currently such targets. GPCRomics, the use of unbiased, hypothesis-generating methods [e.g., RNA-sequencing (RNA-seq)], with tissues and cell types to identify and quantify GPCR expression, has led to the discovery of previously unrecognized GPCRs that contribute to functional responses and pathophysiology and that may be therapeutic targets. The combination of GPCR expression data with validation studies (e.g., signaling and functional activities) provides opportunities for the discovery of disease-relevant GPCR targets and therapeutics. Here, we review insights from GPCRomic approaches, gaps in knowledge, and future directions by which GPCRomics can advance GPCR biology and the discovery of new GPCR-targeted drugs.

Compared to coding sequences, untranslated regions of the transcriptome are not well conserved, and functional annotation of these sequences is challenging. Global relationships between nucleotide composition of 3' UTR sequences and their sequence conservation have been appreciated since mammalian genomes were first sequenced, but the functional relevance of these patterns remain unknown. We systematically measured the effect on gene expression of the sequences of more than 25,000 RNA-binding protein (RBP) binding sites in primary mouse T cells using a massively parallel reporter assay. GC-rich sequences were destabilizing of reporter mRNAs and come from more rapidly evolving regions of the genome. These sequences were more likely to be folded in vivo and contain a number of structural motifs that reduced accumulation of a heterologous reporter protein. Comparison of full-length 3' UTR sequences across vertebrate phylogeny revealed that strictly conserved 3' UTRs were GC-poor and enriched in genes associated with organismal development. In contrast, rapidly evolving 3' UTRs tended to be GC-rich and derived from genes involved in metabolism and immune responses. Cell-essential genes had lower GC content in their 3' UTRs, suggesting a connection between unstructured mRNA noncoding sequences and optimal protein production. By reducing gene expression, GC-rich RBP-occupied sequences act as a rapidly evolving substrate for gene regulatory interactions.

The piRNA pathway is an adaptive mechanism that maintains genome stability by repression of selfish genomic elements. In the male germline of Drosophila melanogaster repression of Stellate genes by piRNAs generated from Supressor of Stellate (Su(Ste)) locus is required for male fertility, but both Su(Ste) piRNAs and their targets are absent in other Drosophila species. We found that D. melanogaster genome contains multiple X-linked non-coding genomic repeats that have sequence similarity to the protein-coding host gene vasa. In the male germline, these vasa-related AT-chX repeats produce abundant piRNAs that are antisense to vasa; however, vasa mRNA escapes silencing due to imperfect complementarity to AT-chX piRNAs. Unexpectedly, we discovered AT-chX piRNAs target vasa of Drosophila mauritiana in the testes of interspecies hybrids. In the majority of hybrid flies, the testes were strongly reduced in size and germline content. A minority of hybrids maintained wild-type array of premeiotic germ cells in the testes, but in them harmful Stellate genes were derepressed due to the absence of Su(Ste) piRNAs, and meiotic failures were observed. Thus, the piRNA pathway contributes to reproductive isolation between D. melanogaster and closely related species, causing hybrid male sterility via misregulation of two different host protein factors.

ObjectiveGenome-wide association studies (GWASs) for epithelial ovarian cancer (EOC) have focused largely on populations of European ancestry. We aimed to identify common germline variants associated with EOC risk in Asian women.MethodsGenotyping was performed as part of the OncoArray project. Samples with >60% Asian ancestry were included in the analysis. Genotyping was performed on 533,631 SNPs in 3238 Asian subjects diagnosed with invasive or borderline EOC and 4083 unaffected controls. After imputation, genotypes were available for 11,595,112 SNPs to identify associations.ResultsAt chromosome 6p25.2, SNP rs7748275 was associated with risk of serous EOC (odds ratio [OR] = 1.34, P = 8.7 × 10-9) and high-grade serous EOC (HGSOC) (OR = 1.34, P = 4.3 × 10-9). SNP rs6902488 at 6p25.2 (r2 = 0.97 with rs7748275) lies in an active enhancer and is predicted to impact binding of STAT3, P300 and ELF1. We identified additional risk loci with low Bayesian false discovery probability (BFDP) scores, indicating they are likely to be true risk associations (BFDP

Eosinophilic esophagitis (EoE) is a recently described disease in which exposure to specific foods and allergens leads to type 2 inflammation, epithelial barrier dysfunction, and difficulty in swallowing. In the current issue of the JCI, Wen and colleagues investigate tissue T cell heterogeneity in patients with EoE using single-cell RNA sequencing (scRNA-seq). Esophageal epithelium from individuals with EoE convtained a prominent population of Th2 cells not seen in controls. The short-chain fatty acid (SCFA) receptor FFAR3 was found to be highly expressed in EoE Th2 cells. Experiments presented here provide evidence that SCFAs may promote type 2 inflammation in allergic diseases such as EoE and asthma. This study provides an early example of scRNA-seq for identifying relevant cell populations and mechanisms underlying allergic diseases.

The majority of the human genome is made of transposable elements, giving rise to interspaced repeats, including Long INterspersed Element-1s (LINE-1s or L1s). L1s are active human transposable elements involved in genomic diversity and evolution; however, they can also contribute to genomic instability and diseases. L1s require host factors to complete their life cycles, whereas the host has evolved numerous mechanisms to restrict L1-induced mutagenesis. Restriction mechanisms in somatic cells include methylation of the L1 promoter, anti-viral factors and RNA-mediated processes such as small RNAs. microRNAs (miRNAs or miRs) are small non-coding RNAs that post-transcriptionally repress multiple target genes often found in the same cellular pathways. We have recently established that miR-128 functions as a novel restriction factor inhibiting L1 mobilization in somatic cells. We have further demonstrated that miR-128 functions through a dual mechanism; by directly targeting L1 RNA for degradation and indirectly by inhibiting a cellular co-factor which L1 is dependent on to transpose to new genomic locations (TNPO1). Here, we add another piece to the puzzle of the enigmatic L1 lifecycle. We show that miR-128 also inhibits another key cellular factor, hnRNPA1 (heterogeneous nuclear ribonucleoprotein A1), by significantly reducing mRNA and protein levels through direct interaction with the coding sequence (CDS) of hnRNPA1 mRNA. In addition, we demonstrate that repression of hnRNPA1 using hnRNPA1-shRNA significantly decreases de novo L1 retro-transposition and that induced hnRNPA1 expression enhances L1 mobilization. Furthermore, we establish that hnRNPA1 is a functional target of miR-128. Finally, we determine that induced hnRNPA1 expression in miR-128-overexpressing cells can partly rescue the miR-128-induced repression of L1's ability to transpose to different genomic locations. Thus, we have identified an additional mechanism by which miR-128 represses L1 retro-transposition and mediates genomic stability.

Patients with neurofibromatosis type 1 (NF1) are predisposed to develop neurofibromas, but the underlying molecular mechanisms of neurofibromagenesis are not fully understood. We showed dual genetic deletion of Runx1 and Runx3 in Schwann cells (SCs) and SC precursors delayed neurofibromagenesis and prolonged mouse survival. We identified peripheral myelin protein 22 (Pmp22/Gas3) related to neurofibroma initiation. Knockdown of Pmp22 with short hairpin RNAs increased Runx1fl/fl;Runx3fl/fl;Nf1fl/fl;DhhCre tumor-derived sphere numbers and enabled significantly more neurofibroma-like microlesions on transplantation. Conversely, overexpression of Pmp22 in mouse neurofibroma SCs decreased cell proliferation. Mechanistically, RUNX1/3 regulated alternative promoter usage and induced levels of protein expression of Pmp22 to control SC growth. Last, pharmacological inhibition of RUNX/core-binding factor β (CBFB) activity significantly reduced neurofibroma volume in vivo. Thus, we identified a signaling pathway involving RUNX1/3 suppression of Pmp22 in neurofibroma initiation and/or maintenance. Targeting disruption of RUNX/CBFB interaction might provide a novel therapy for patients with neurofibroma.

The immunoglobulin heavy variable (IGHV) and T cell beta variable (TRBV) loci are among the most complex and variable regions in the human genome. Generated through a process of gene duplication/deletion and diversification, these loci can vary extensively between individuals in copy number and contain genes that are highly similar, making their analysis technically challenging. Here, we present a comprehensive study of the functional gene segments in the IGHV and TRBV loci, quantifying their copy number and single-nucleotide variation in a globally diverse sample of 109 (IGHV) and 286 (TRBV) humans from over a 100 populations. We find that the IGHV and TRBV gene families exhibit starkly different patterns of variation. In addition to providing insight into the different evolutionary paths of the IGHV and TRBV loci, our results are also important to the adaptive immune repertoire sequencing community, where the lack of frequencies of common alleles and copy number variants is hampering existing analytical pipelines.