Genome Research current issue

Sense-antisense miRNA pairs constitute an elaborate reciprocal regulatory circuit [RESEARCH]

2020-05-27T06:30:21-07:00

Antisense transcription of protein-coding genes has been increasingly recognized as an important regulatory mechanism of gene expression. However, less is known about the extent and importance of antisense transcription of noncoding genes. Here, we investigate the breadth and dynamics of antisense transcription of miRNAs, a class of important noncoding RNAs. Because the antisense transcript of a miRNA is likely to form a hairpin suitable as the substrate of ADARs, which convert adenosine to inosine in double-stranded RNAs, we used A-to-I RNA editing as ultrasensitive readout for antisense transcription of the miRNAs. Through examining the unstranded targeted RNA-seq libraries covering all miRNA loci in 25 types of human tissues, we identified 7275 editing events located in 81% of the antisense strand of the miRNA loci, thus uncovering the previously unknown prevalent antisense transcription of the miRNAs. We found that antisense transcripts are tightly regulated, and a substantial fraction of miRNAs and their antisense transcripts are coexpressed. Sense miRNAs have been shown to down-regulate the coexpressed antisense transcripts, whereas the act of antisense transcription, rather than the transcripts themselves, regulates the expression of sense miRNAs. RNA editing tends to decrease the miRNA accessibility of the antisense transcripts, therefore protecting them from being degraded by the sense-mature miRNAs. Altogether, our study reveals the landscape of antisense transcription and editing of miRNAs, as well as a previously unknown reciprocal regulatory circuit of sense–antisense miRNA pairs.

Dynamic effects of interacting genes underlying rice flowering-time phenotypic plasticity and global adaptation [RESEARCH]

Guo, T., Mu, Q., Wang, J., Vanous, A. E., Onogi, A., Iwata, H., Li, X., Yu, J. — 2020-05-27T06:30:21-07:00

The phenotypic variation of living organisms is shaped by genetics, environment, and their interaction. Understanding phenotypic plasticity under natural conditions is hindered by the apparently complex environment and the interacting genes and pathways. Herein, we report findings from the dissection of rice flowering-time plasticity in a genetic mapping population grown in natural long-day field environments. Genetic loci harboring four genes originally discovered for their photoperiodic effects (Hd1, Hd2, Hd5, and Hd6) were found to differentially respond to temperature at the early growth stage to jointly determine flowering time. The effects of these plasticity genes were revealed with multiple reaction norms along the temperature gradient. By coupling genomic selection and the environmental index, accurate performance predictions were obtained. Next, we examined the allelic variation in the four flowering-time genes across the diverse accessions from the 3000 Rice Genomes Project and constructed haplotypes at both individual-gene and multigene levels. The geographic distribution of haplotypes revealed their preferential adaptation to different temperature zones. Regions with lower temperatures were dominated by haplotypes sensitive to temperature changes, whereas the equatorial region had a majority of haplotypes that are less responsive to temperature. By integrating knowledge from genomics, gene cloning and functional characterization, and environment quantification, we propose a conceptual model with multiple levels of reaction norms to help bridge the gaps among individual gene discovery, field-level phenotypic plasticity, and genomic diversity and adaptation.

Polymorphic centromere locations in the pathogenic yeast Candida parapsilosis [RESEARCH]

2020-05-27T06:30:21-07:00

Centromeres pose an evolutionary paradox: strongly conserved in function but rapidly changing in sequence and structure. However, in the absence of damage, centromere locations are usually conserved within a species. We report here that isolates of the pathogenic yeast species Candida parapsilosis show within-species polymorphism for the location of centromeres on two of its eight chromosomes. Its old centromeres have an inverted-repeat (IR) structure, whereas its new centromeres have no obvious structural features but are located within 30 kb of the old site. Centromeres can therefore move naturally from one chromosomal site to another, apparently spontaneously and in the absence of any significant changes in DNA sequence. Our observations are consistent with a model in which all centromeres are genetically determined, such as by the presence of short or long IRs or by the ability to form cruciforms. We also find that centromeres have been hotspots for genomic rearrangements in the C. parapsilosis clade.

Intragenic repeat expansion in the cell wall protein gene HPF1 controls yeast chronological aging [RESEARCH]

2020-05-27T06:30:21-07:00

Aging varies among individuals due to both genetics and environment, but the underlying molecular mechanisms remain largely unknown. Using a highly recombined Saccharomyces cerevisiae population, we found 30 distinct quantitative trait loci (QTLs) that control chronological life span (CLS) in calorie-rich and calorie-restricted environments and under rapamycin exposure. Calorie restriction and rapamycin extended life span in virtually all genotypes but through different genetic variants. We tracked the two major QTLs to the cell wall glycoprotein genes FLO11 and HPF1. We found that massive expansion of intragenic tandem repeats within the N-terminal domain of HPF1 was sufficient to cause pronounced life span shortening. Life span impairment by HPF1 was buffered by rapamycin but not by calorie restriction. The HPF1 repeat expansion shifted yeast cells from a sedentary to a buoyant state, thereby increasing their exposure to surrounding oxygen. The higher oxygenation altered methionine, lipid, and purine metabolism, and inhibited quiescence, which explains the life span shortening. We conclude that fast-evolving intragenic repeat expansions can fundamentally change the relationship between cells and their environment with profound effects on cellular lifestyle and longevity.

Global fitness landscapes of the Shine-Dalgarno sequence [RESEARCH]

2020-05-27T06:30:21-07:00

Shine-Dalgarno sequences (SD) in prokaryotic mRNA facilitate protein translation by pairing with rRNA in ribosomes. Although conventionally defined as AG-rich motifs, recent genomic surveys reveal great sequence diversity, questioning how SD functions. Here, we determined the molecular fitness (i.e., translation efficiency) of 4⁹ synthetic 9-nt SD genotypes in three distinct mRNA contexts in Escherichia coli. We uncovered generic principles governing the SD fitness landscapes: (1) Guanine contents, rather than canonical SD motifs, best predict the fitness of both synthetic and endogenous SD; (2) the genotype-fitness correlation of SD promotes its evolvability by steadily supplying beneficial mutations across fitness landscapes; and (3) the frequency and magnitude of deleterious mutations increase with background fitness, and adjacent nucleotides in SD show stronger epistasis. Epistasis results from disruption of the continuous base pairing between SD and rRNA. This "chain-breaking" epistasis creates sinkholes in SD fitness landscapes and may profoundly impact the evolution and function of prokaryotic translation initiation and other RNA-mediated processes. Collectively, our work yields functional insights into the SD sequence variation in prokaryotic genomes, identifies a simple design principle to guide bioengineering and bioinformatic analysis of SD, and illuminates the fundamentals of fitness landscapes and molecular evolution.

Determining the impact of uncharacterized inversions in the human genome by droplet digital PCR [METHOD]

2020-05-27T06:30:21-07:00

Despite the interest in characterizing genomic variation, the presence of large repeats at the breakpoints hinders the analysis of many structural variants. This is especially problematic for inversions, since there is typically no gain or loss of DNA. Here, we tested novel linkage-based droplet digital PCR (ddPCR) assays to study 20 inversions ranging from 3.1 to 742 kb flanked by inverted repeats (IRs) up to 134 kb long. Of those, we validated 13 inversions predicted by different genome-wide techniques. In addition, we obtained new experimental human population information across 95 African, European, and East Asian individuals for 16 inversions, including four already validated variants without high-throughput genotyping methods. Through comparison with previous data, independent replicates and both inversion breakpoints, we demonstrate that the technique is highly accurate and reproducible. Most studied inversions are widespread across continents, and their frequency is negatively correlated with genetic length. Moreover, all except two show clear signs of being recurrent, and we could better define the factors affecting recurrence levels and estimate the inversion rate across the genome. Finally, the generated genotypes have allowed us to check inversion functional effects, validating gene expression differences reported before for two inversions and finding new candidate associations. Therefore, the developed methodology makes it possible to screen these and other complex genomic variants quickly in a large number of samples for the first time, highlighting the importance of direct genotyping to assess their potential consequences and clinical implications.

MEDEA: analysis of transcription factor binding motifs in accessible chromatin [METHOD]

Mariani, L., Weinand, K., Gisselbrecht, S. S., Bulyk, M. L. — 2020-05-27T06:30:21-07:00

Deciphering the interplay between chromatin accessibility and transcription factor (TF) binding is fundamental to understanding transcriptional regulation, control of cellular states, and the establishment of new phenotypes. Recent genome-wide chromatin accessibility profiling studies have provided catalogs of putative open regions, where TFs can recognize their motifs and regulate gene expression programs. Here, we present motif enrichment in differential elements of accessibility (MEDEA), a computational tool that analyzes high-throughput chromatin accessibility genomic data to identify cell-type-specific accessible regions and lineage-specific motifs associated with TF binding therein. To benchmark MEDEA, we used a panel of reference cell lines profiled by ENCODE and curated by the ENCODE Project Consortium for the ENCODE-DREAM Challenge. By comparing results with RNA-seq data, ChIP-seq peaks, and DNase-seq footprints, we show that MEDEA improves the detection of motifs associated with known lineage specifiers. We then applied MEDEA to 610 ENCODE DNase-seq data sets, where it revealed significant motifs even when absolute enrichment was low and where it identified novel regulators, such as NRF1 in kidney development. Finally, we show that MEDEA performs well on both bulk and single-cell ATAC-seq data. MEDEA is publicly available as part of our Glossary-GENRE suite for motif enrichment analysis.

Exploring dimension-reduced embeddings with Sleepwalk [METHOD]

Ovchinnikova, S., Anders, S. — 2020-05-27T06:30:21-07:00

Dimension-reduction methods, such as t-SNE or UMAP, are widely used when exploring high-dimensional data describing many entities, for example, RNA-seq data for many single cells. However, dimension reduction is commonly prone to introducing artifacts, and we hence need means to see where a dimension-reduced embedding is a faithful representation of the local neighborhood and where it is not. We present Sleepwalk, a simple but powerful tool that allows the user to interactively explore an embedding, using color to depict original or any other distances from all points to the cell under the mouse cursor. We show how this approach not only highlights distortions but also reveals otherwise hidden characteristics of the data, and how Sleepwalk's comparative modes help integrate multisample data and understand differences between embedding and preprocessing methods. Sleepwalk is a versatile and intuitive tool that unlocks the full power of dimension reduction and will be of value not only in single-cell RNA-seq but also in any other area with matrix-shaped big data.

Genome-wide CRISPR screening reveals genes essential for cell viability and resistance to abiotic and biotic stresses in Bombyx mori [METHOD]

2020-05-27T06:30:21-07:00

High-throughput genetic screens are powerful methods to interrogate gene function on a genome-wide scale and identify genes responsible to certain stresses. Here, we developed a piggyBac strategy to deliver pooled sgRNA libraries stably into cell lines. We used this strategy to conduct a screen based on genome-wide clustered regularly interspaced short palindromic repeat technology (CRISPR)-Cas9 in Bombyx mori cells. We first constructed a single guide RNA (sgRNA) library containing 94,000 sgRNAs, which targeted 16,571 protein-coding genes. We then generated knockout collections in BmE cells using the piggyBac transposon. We identified 1006 genes that are essential for cell viability under normal growth conditions. Of the identified genes, 82.4% (829 genes) were homologous to essential genes in seven animal species. We also identified 838 genes whose loss facilitated cell growth. Next, we performed context-specific positive screens for resistance to biotic or nonbiotic stresses using temperature and baculovirus separately, which identified several key genes and pathways from each screen. Collectively, our results provide a novel and versatile platform for functional annotations of B. mori genomes and deciphering key genes responsible for various conditions. This study also shows the effectiveness, practicality, and convenience of genome-wide CRISPR screens in nonmodel organisms.

CRISPR-Cas9-mediated pinpoint microbial genome editing aided by target-mismatched sgRNAs [METHOD]

Lee, H. J., Kim, H. J., Lee, S. J. — 2020-05-27T06:30:21-07:00

Genome editing has been revolutionized by the CRISPR-Cas9 system. CRISPR-Cas9 is composed of single-molecular guide RNA (sgRNA) and a proteinaceous Cas9 nuclease, which recognizes a specific target sequence and a protospacer adjacent motif (PAM) sequence and, subsequently, cleaves the targeted DNA sequence. This CRISPR-Cas9 system has been used as an efficient negative-selection tool to cleave unedited or unchanged target DNAs during site-specific mutagenesis and, consequently, obtain microbial cells with desired mutations. This study aimed to investigate the genome editing efficiency of the CRISPR-Cas9 system for in vivo oligonucleotide-directed mutagenesis in bacteria. This system successfully introduced two- to four-base mutations in galK in Escherichia coli with high editing efficiencies (81%–86%). However, single-point mutations (T504A or C578A) were rarely introduced with very low editing efficiencies (<3%), probably owing to mismatch tolerance. To resolve this issue, we designed one- or two-base mismatches in the sgRNA sequence to recognize target sequences in galK in E. coli. A single-point nucleotide mutation (T504A or C578A in the galK gene) was successfully introduced in 36%–95% of negatively selected E. coli cells using single-base mismatched sgRNAs. Sixteen targets were randomly selected through genome-wide single-base editing experiments using mismatched sgRNAs. Consequently, out of 48 desired single-base mutations, 25 single bases were successfully edited, using mismatched sgRNAs. Finally, applicable design rules for target-mismatched sgRNAs were provided for single-nucleotide editing in microbial genomes.

Single-cell-resolution transcriptome map of human, chimpanzee, bonobo, and macaque brains [RESOURCES]

2020-05-27T06:30:21-07:00

Identification of gene expression traits unique to the human brain sheds light on the molecular mechanisms underlying human evolution. Here, we searched for uniquely human gene expression traits by analyzing 422 brain samples from humans, chimpanzees, bonobos, and macaques representing 33 anatomical regions, as well as 88,047 cell nuclei composing three of these regions. Among 33 regions, cerebral cortex areas, hypothalamus, and cerebellar gray and white matter evolved rapidly in humans. At the cellular level, astrocytes and oligodendrocyte progenitors displayed more differences in the human evolutionary lineage than the neurons. Comparison of the bulk tissue and single-nuclei sequencing revealed that conventional RNA sequencing did not detect up to two-thirds of cell-type-specific evolutionary differences.

Comprehensive analyses of 723 transcriptomes enhance genetic and biological interpretations for complex traits in cattle [RESOURCES]

2020-05-27T06:30:21-07:00

By uniformly analyzing 723 RNA-seq data from 91 tissues and cell types, we built a comprehensive gene atlas and studied tissue specificity of genes in cattle. We demonstrated that tissue-specific genes significantly reflected the tissue-relevant biology, showing distinct promoter methylation and evolution patterns (e.g., brain-specific genes evolve slowest, whereas testis-specific genes evolve fastest). Through integrative analyses of those tissue-specific genes with large-scale genome-wide association studies, we detected relevant tissues/cell types and candidate genes for 45 economically important traits in cattle, including blood/immune system (e.g., CCDC88C) for male fertility, brain (e.g., TRIM46 and RAB6A) for milk production, and multiple growth-related tissues (e.g., FGF6 and CCND2) for body conformation. We validated these findings by using epigenomic data across major somatic tissues and sperm. Collectively, our findings provided novel insights into the genetic and biological mechanisms underlying complex traits in cattle, and our transcriptome atlas can serve as a primary source for biological interpretation, functional validation, studies of adaptive evolution, and genomic improvement in livestock.