Mostly natural sequencing-by-synthesis for scRNA-seq using Ultima sequencing.

Here we introduce a mostly natural sequencing-by-synthesis (mnSBS) method for single-cell RNA sequencing (scRNA-seq), adapted to the Ultima genomics platform, and systematically benchmark it against current scRNA-seq technology. mnSBS uses mostly natural, unmodified nucleotides and only a low fraction of fluorescently labeled nucleotides, which allows for high polymerase processivity and lower costs. We demonstrate successful application in four scRNA-seq case studies of different technical and biological types, including 5' and 3' scRNA-seq, human peripheral blood mononuclear cells from a single individual and in multiplex, as well as Perturb-Seq. Benchmarking shows that results from mnSBS-based scRNA-seq are very similar to those using Illumina sequencing, with minor differences in results related to the position of reads relative to annotated gene boundaries, owing to single-end reads of Ultima being closer to gene ends than reads from Illumina. The method is thus compatible with state-of-the-art scRNA-seq libraries independent of the sequencing technology. We expect mnSBS to be of particular utility for cost-effective large-scale scRNA-seq projects.

Publisher Correction: Sex-dependent dominance maintains migration supergene in rainbow trout.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Sex-dependent dominance maintains migration supergene in rainbow trout.

Males and females often differ in their fitness optima for shared traits that have a shared genetic basis, leading to sexual conflict. Morphologically differentiated sex chromosomes can resolve this conflict and protect sexually antagonistic variation, but they accumulate deleterious mutations. However, how sexual conflict is resolved in species that lack differentiated sex chromosomes is largely unknown. Here we present a chromosome-anchored genome assembly for rainbow trout (Oncorhynchus mykiss) and characterize a 55-Mb double-inversion supergene that mediates sex-specific migratory tendency through sex-dependent dominance reversal, an alternative mechanism for resolving sexual conflict. The double inversion contains key photosensory, circadian rhythm, adiposity and sex-related genes and displays a latitudinal frequency cline, indicating environmentally dependent selection. Our results show sex-dependent dominance reversal across a large autosomal supergene, a mechanism for sexual conflict resolution capable of protecting sexually antagonistic variation while avoiding the homozygous lethality and deleterious mutations associated with typical heteromorphic sex chromosomes.

The maize W22 genome provides a foundation for functional genomics and transposon biology.

The maize W22 inbred has served as a platform for maize genetics since the mid twentieth century. To streamline maize genome analyses, we have sequenced and de novo assembled a W22 reference genome using short-read sequencing technologies. We show that significant structural heterogeneity exists in comparison to the B73 reference genome at multiple scales, from transposon composition and copy number variation to single-nucleotide polymorphisms. The generation of this reference genome enables accurate placement of thousands of Mutator (Mu) and Dissociation (Ds) transposable element insertions for reverse and forward genetics studies. Annotation of the genome has been achieved using RNA-seq analysis, differential nuclease sensitivity profiling and bisulfite sequencing to map open reading frames, open chromatin sites and DNA methylation profiles, respectively. Collectively, the resources developed here integrate W22 as a community reference genome for functional genomics and provide a foundation for the maize pan-genome.

Deciphering genetic factors that determine melon fruit-quality traits using RNA-Seq-based high-resolution QTL and eQTL mapping.

Combined quantitative trait loci (QTL) and expression-QTL (eQTL) mapping analysis was performed to identify genetic factors affecting melon (Cucumis melo) fruit quality, by linking genotypic, metabolic and transcriptomic data from a melon recombinant inbred line (RIL) population. RNA sequencing (RNA-Seq) of fruit from 96 RILs yielded a highly saturated collection of > 58 000 single-nucleotide polymorphisms, identifying 6636 recombination events that separated the genome into 3663 genomic bins. Bin-based QTL analysis of 79 RILs and 129 fruit-quality traits affecting taste, aroma and color resulted in the mapping of 241 QTL. Thiol acyltransferase (CmThAT1) gene was identified within the QTL interval of its product, S-methyl-thioacetate, a key component of melon fruit aroma. Metabolic activity of CmThAT1-encoded protein was validated in bacteria and in vitro. QTL analysis of flesh color intensity identified a candidate white-flesh gene (CmPPR1), one of two major loci determining fruit flesh color in melon. CmPPR1 encodes a member of the pentatricopeptide protein family, involved in processing of RNA in plastids, where carotenoid and chlorophyll pigments accumulate. Network analysis of > 12 000 eQTL mapped for > 8000 differentially expressed fruit genes supported the role of CmPPR1 in determining the expression level of plastid targeted genes. We highlight the potential of RNA-Seq-based QTL analysis of small to moderate size, advanced RIL populations for precise marker-assisted breeding and gene discovery. We provide the following resources: a RIL population genotyped with a unique set of SNP markers, confined genomic segments that harbor QTL governing 129 traits and a saturated set of melon eQTLs.

Wild emmer genome architecture and diversity elucidate wheat evolution and domestication.

Wheat (Triticum spp.) is one of the founder crops that likely drove the Neolithic transition to sedentary agrarian societies in the Fertile Crescent more than 10,000 years ago. Identifying genetic modifications underlying wheat's domestication requires knowledge about the genome of its allo-tetraploid progenitor, wild emmer (T. turgidum ssp. dicoccoides). We report a 10.1-gigabase assembly of the 14 chromosomes of wild tetraploid wheat, as well as analyses of gene content, genome architecture, and genetic diversity. With this fully assembled polyploid wheat genome, we identified the causal mutations in Brittle Rachis 1 (TtBtr1) genes controlling shattering, a key domestication trait. A study of genomic diversity among wild and domesticated accessions revealed genomic regions bearing the signature of selection under domestication. This reference assembly will serve as a resource for accelerating the genome-assisted improvement of modern wheat varieties.

Draft Assembly of Elite Inbred Line PH207 Provides Insights into Genomic and Transcriptome Diversity in Maize.

Intense artificial selection over the last 100 years has produced elite maize (Zea mays) inbred lines that combine to produce high-yielding hybrids. To further our understanding of how genome and transcriptome variation contribute to the production of high-yielding hybrids, we generated a draft genome assembly of the inbred line PH207 to complement and compare with the existing B73 reference sequence. B73 is a founder of the Stiff Stalk germplasm pool, while PH207 is a founder of Iodent germplasm, both of which have contributed substantially to the production of temperate commercial maize and are combined to make heterotic hybrids. Comparison of these two assemblies revealed over 2500 genes present in only one of the two genotypes and 136 gene families that have undergone extensive expansion or contraction. Transcriptome profiling revealed extensive expression variation, with as many as 10,564 differentially expressed transcripts and 7128 transcripts expressed in only one of the two genotypes in a single tissue. Genotype-specific genes were more likely to have tissue/condition-specific expression and lower transcript abundance. The availability of a high-quality genome assembly for the elite maize inbred PH207 expands our knowledge of the breadth of natural genome and transcriptome variation in elite maize inbred lines across heterotic pools.

High-resolution genetic mapping of maize pan-genome sequence anchors.

In addition to single-nucleotide polymorphisms, structural variation is abundant in many plant genomes. The structural variation across a species can be represented by a 'pan-genome', which is essential to fully understand the genetic control of phenotypes. However, the pan-genome's complexity hinders its accurate assembly via sequence alignment. Here we demonstrate an approach to facilitate pan-genome construction in maize. By performing 18 trillion association tests we map 26 million tags generated by reduced representation sequencing of 14,129 maize inbred lines. Using machine-learning models we select 4.4 million accurately mapped tags as sequence anchors, 1.1 million of which are presence/absence variations. Structural variations exhibit enriched association with phenotypic traits, indicating that it is a significant source of adaptive variation in maize. The ability to efficiently map ultrahigh-density pan-genome sequence anchors enables fine characterization of structural variation and will advance both genetic research and breeding in many crops.

Efficiency and specificity in microRNA biogenesis.

Primary microRNA cleavage by the Drosha-Dgcr8 'Microprocessor' complex is critical for microRNA biogenesis. Yet, the Microprocessor may also cleave other nuclear RNAs in a nonspecific manner. We studied Microprocessor function using mathematical modeling and experiments in mouse and human tissues. We found that the autoregulatory feedback on Microprocessor expression is instrumental for balancing the efficiency and specificity of its activity by effectively tuning Microprocessor levels to those of its pri-miRNA substrate.

Robust selection of sensory organ precursors by the Notch-Delta pathway.

The patterning of multicellular organisms is robust to environmental, genetic, or stochastic fluctuations. Mathematical modeling is instrumental in identifying mechanisms supporting this robustness. The principle of lateral inhibition, whereby a differentiating cell inhibits its neighbors from adopting the same fate, is frequently used for selecting a single cell out of a cluster of equipotent cells. For example, Sensory Organ Precursors (SOP) in the fruit-fly Drosophila implement lateral inhibition by activating the Notch-Delta pathway. We discuss parameters affecting the rate of errors in this process, and the mechanism (inhibitory cis interaction between Notch and Delta) predicted to reduce this error.

A biological solution to a fundamental distributed computing problem.

Computational and biological systems are often distributed so that processors (cells) jointly solve a task, without any of them receiving all inputs or observing all outputs. Maximal independent set (MIS) selection is a fundamental distributed computing procedure that seeks to elect a set of local leaders in a network. A variant of this problem is solved during the development of the fly's nervous system, when sensory organ precursor (SOP) cells are chosen. By studying SOP selection, we derived a fast algorithm for MIS selection that combines two attractive features. First, processors do not need to know their degree; second, it has an optimal message complexity while only using one-bit messages. Our findings suggest that simple and efficient algorithms can be developed on the basis of biologically derived insights.

miRNA-based mechanism for the commitment of multipotent progenitors to a single cellular fate.

When stem cells and multipotent progenitors differentiate, they undergo fate restriction, enabling a single fate and blocking differentiation along alternative routes. We herein present a mechanism whereby such unequivocal commitment is achieved, based on microRNA (miRNA)-dependent repression of an alternative cell fate. We show that the commitment of monocyte RAW264.7 progenitors to active macrophage differentiation involves rapid up-regulation of miR-155 expression, which leads to the suppression of the alternative pathway, namely RANK ligand-induced osteoclastogenesis, by repressing the expression of MITF, a transcription factor essential for osteoclast differentiation. A temporal asymmetry, whereby miR-155 expression precedes and overrides the activation of the osteoclast transcriptional program, provides the means for coherent macrophage differentiation, even in the presence of osteoclastogenic signals. Based on these findings, we propose that miRNA may provide a general mechanism for the unequivocal commitment underlying stem cell differentiation.

Error minimization in lateral inhibition circuits.

The pattern of the sensory bristles in the fruit fly Drosophila is remarkably reproducible. Each bristle arises from a sensory organ precursor (SOP) cell that is selected, through a lateral inhibition process, from a cluster of proneural cells. Although this process is well characterized, the mechanism ensuring its robustness remains obscure. Using probabilistic modeling, we defined the sources of error in SOP selection and examined how they depend on the underlying molecular circuit. We found that rapid inhibition of the neural differentiation of nonselected cells, coupled with high cell-to-cell variability in the timing of selection, is crucial for accurate SOP selection. Cell-autonomous interactions (cis interactions) between the Notch receptor and its ligands Delta or Serrate facilitate accurate SOP selection by shortening the effective delay between the time when the inhibitory signal is initiated in one cell and the time when it acts on neighboring cells, suggesting that selection relies on competition between cis and trans interactions of Notch with its ligands. The cis interaction model predicts that the increase in ectopic SOP selections observed with reduced Notch abundance can be compensated for by reducing the abundance of the Notch ligands Delta and Serrate. We validated this prediction experimentally by quantifying the frequency of ectopic bristles in flies carrying heterozygous null mutations of Notch, Delta, or Serrate or combinations of these alleles. We propose that susceptibility to errors distinguishes seemingly equivalent designs of developmental circuits regulating pattern formation.

MicroRNA expression patterns and function in endodermal differentiation of human embryonic stem cells.

BACKGROUND/AIMS: microRNAs (miRNAs) are small noncoding RNAs that regulate cognate mRNAs post-transcriptionally. Human embryonic stem cells (hESC), which exhibit the characteristics of pluripotency and self-renewal, may serve as a model to study the role of miRNAs in early human development. We aimed to determine whether endodermally-differentiated hESC demonstrate a unique miRNA expression pattern, and whether overexpression of endoderm-specific miRNA may affect hESC differentiation. METHODS: miRNA expression was profiled in undifferentiated and NaButyrate-induced differentiated hESC of two lines, using microarray and quantitative RT-PCR. Then, the effect of lentiviral-based overexpression of liver-specific miR-122 on hESC differentiation was analyzed, using genomewide gene microarrays. RESULTS: The miRNA profiling revealed expression of three novel miRNAs in undifferentiated and differentiated hESC. Upon NaButyrate induction, two of the most upregulated miRNAs common to both cell lines were miR-24 and miR-10a, whose target genes have been shown to inhibit endodermal differentiation. Furthermore, induction of several liver-enriched miRNAs, including miR-122 and miR-192, was observed in parallel to induction of endodermal gene expression. Stable overexpression of miR-122 in hESC was unable to direct spontaneous differentiation towards a clear endodermal fate, but rather, delayed general differentiation of these cells. CONCLUSIONS: Our results demonstrate that expression of specific miRNAs correlates with that of specific genes upon differentiation, and highlight the potential role of miRNAs in endodermal differentiation of hESC.

The colorectal microRNAome.

MicroRNAs (miRNAs) are a class of small noncoding RNAs that have important regulatory roles in multicellular organisms. The public miRNA database contains 321 human miRNA sequences, 234 of which have been experimentally verified. To explore the possibility that additional miRNAs are present in the human genome, we have developed an experimental approach called miRNA serial analysis of gene expression (miRAGE) and used it to perform the largest experimental analysis of human miRNAs to date. Sequence analysis of 273,966 small RNA tags from human colorectal cells allowed us to identify 200 known mature miRNAs, 133 novel miRNA candidates, and 112 previously uncharacterized miRNA* forms. To aid in the evaluation of candidate miRNAs, we disrupted the Dicer locus in three human colorectal cancer cell lines and examined known and novel miRNAs in these cells. These studies suggest that the human genome contains many more miRNAs than currently identified and provide an approach for the large-scale experimental cloning of novel human miRNAs in human tissues.

Identification of hundreds of conserved and nonconserved human microRNAs.

MicroRNAs are noncoding RNAs of approximately 22 nucleotides that suppress translation of target genes by binding to their mRNA and thus have a central role in gene regulation in health and disease. To date, 222 human microRNAs have been identified, 86 by random cloning and sequencing, 43 by computational approaches and the rest as putative microRNAs homologous to microRNAs in other species. To prove our hypothesis that the total number of microRNAs may be much larger and that several have emerged only in primates, we developed an integrative approach combining bioinformatic predictions with microarray analysis and sequence-directed cloning. Here we report the use of this approach to clone and sequence 89 new human microRNAs (nearly doubling the current number of sequenced human microRNAs), 53 of which are not conserved beyond primates. These findings suggest that the total number of human microRNAs is at least 800.

Design principle of gene expression used by human stem cells: implication for pluripotency.

Human embryonic stem cells (ESC) are undifferentiated and are endowed with the capacities of self-renewal and pluripotential differentiation. Adult stem cells renew their own tissue, but whether they can transdifferentiate to other tissues is still controversial. To understand the genetic program that underlies the pluripotency of stem cells, we compared the transcription profile of ESC with that of progenitor/stem cells of human hematopoietic and keratinocytic origins, along with their mature cells to be viewed as snapshots along tissue differentiation. ESC gene profiles show higher complexity with significantly more highly expressed genes than adult cells. We hypothesize that ESC use a strategy of expressing genes that represent various differentiation pathways and selection of only a few for continuous expression upon differentiation to a particular target. Such a strategy may be necessary for the pluripotency of ESC. The progenitors of either hematopoietic or keratinocytic cells also follow the same design principle. Using advanced clustering, we show that many of the ESC expressed genes are turned off in the progenitors/stem cells followed by a further down-regulation in adult tissues. Concomitantly, genes specific to the target tissue are up-regulated toward mature cells of skin or blood.

MicroRNA expression detected by oligonucleotide microarrays: system establishment and expression profiling in human tissues.

MicroRNAs (MIRs) are a novel group of conserved short approximately 22 nucleotide-long RNAs with important roles in regulating gene expression. We have established a MIR-specific oligonucleotide microarray system that enables efficient analysis of the expression of the human MIRs identified so far. We show that the 60-mer oligonucleotide probes on the microarrays hybridize with labeled cRNA of MIRs, but not with their precursor hairpin RNAs, derived from amplified, size-fractionated, total RNA of human origin. Signal intensity is related to the location of the MIR sequences within the 60-mer probes, with location at the 5' region giving the highest signals, and at the 3' end, giving the lowest signals. Accordingly, 60-mer probes harboring one MIR copy at the 5' end gave signals of similar intensity to probes containing two or three MIR copies. Mismatch analysis shows that mutations within the MIR sequence significantly reduce or eliminate the signal, suggesting that the observed signals faithfully reflect the abundance of matching MIRs in the labeled cRNA. Expression profiling of 150 MIRs in five human tissues and in HeLa cells revealed a good overall concordance with previously published results, but also with some differences. We present novel data on MIR expression in thymus, testes, and placenta, and have identified MIRs highly enriched in these tissues. Taken together, these results highlight the increased sensitivity of the DNA microarray over other methods for the detection and study of MIRs, and the immense potential in applying such microarrays for the study of MIRs in health and disease.