Circulating transcripts in maternal blood reflect a molecular signature of early-onset preeclampsia.

Circulating RNA (C-RNA) is continually released into the bloodstream from tissues throughout the body, offering an opportunity to noninvasively monitor all aspects of pregnancy health from conception to birth. We asked whether C-RNA analysis could robustly detect aberrations in patients diagnosed with preeclampsia (PE), a prevalent and potentially fatal pregnancy complication. As an initial examination, we sequenced the circulating transcriptome from 40 pregnancies at the time of severe, early-onset PE diagnosis and 73 gestational age-matched controls. Differential expression analysis identified 30 transcripts with gene ontology annotations and tissue expression patterns consistent with the placental dysfunction, impaired fetal development, and maternal immune and cardiovascular system dysregulation characteristic of PE. Furthermore, machine learning identified combinations of 49 C-RNA transcripts that classified an independent cohort of patients (early-onset PE, n = 12; control, n = 12) with 85 to 89% accuracy. C-RNA may thus hold promise for improving the diagnosis and identification of at-risk pregnancies.

Targeted or whole genome sequencing of formalin fixed tissue samples: potential applications in cancer genomics.

Current genomic studies are limited by the poor availability of fresh-frozen tissue samples. Although formalin-fixed diagnostic samples are in abundance, they are seldom used in current genomic studies because of the concern of formalin-fixation artifacts. Better characterization of these artifacts will allow the use of archived clinical specimens in translational and clinical research studies. To provide a systematic analysis of formalin-fixation artifacts on Illumina sequencing, we generated 26 DNA sequencing data sets from 13 pairs of matched formalin-fixed paraffin-embedded (FFPE) and fresh-frozen (FF) tissue samples. The results indicate high rate of concordant calls between matched FF/FFPE pairs at reference and variant positions in three commonly used sequencing approaches (whole genome, whole exome, and targeted exon sequencing). Global mismatch rates and C · G > T · A substitutions were comparable between matched FF/FFPE samples, and discordant rates were low (<0.26%) in all samples. Finally, low-pass whole genome sequencing produces similar pattern of copy number alterations between FF/FFPE pairs. The results from our studies suggest the potential use of diagnostic FFPE samples for cancer genomic studies to characterize and catalog variations in cancer genomes.

Whole-genome haplotyping by dilution, amplification, and sequencing.

Standard whole-genome genotyping technologies are unable to determine haplotypes. Here we describe a method for rapid and cost-effective long-range haplotyping. Genomic DNA is diluted and distributed into multiple aliquots such that each aliquot receives a fraction of a haploid copy. The DNA template in each aliquot is amplified by multiple displacement amplification, converted into barcoded sequencing libraries using Nextera technology, and sequenced in multiplexed pools. To assess the performance of our method, we combined two male genomic DNA samples at equal ratios, resulting in a sample with diploid X chromosomes with known haplotypes. Pools of the multiplexed sequencing libraries were subjected to targeted pull-down of a 1-Mb contiguous region of the X-chromosome Duchenne muscular dystrophy gene. We were able to phase the Duchenne muscular dystrophy region into two contiguous haplotype blocks with a mean length of 494 kb. The haplotypes showed 99% agreement with the consensus base calls made by sequencing the individual DNAs. We subsequently used the strategy to haplotype two human genomes. Standard genomic sequencing to identify all heterozygous SNPs in the sample was combined with dilution-amplification-based sequencing data to resolve the phase of identified heterozygous SNPs. Using this procedure, we were able to phase >95% of the heterozygous SNPs from the diploid sequence data. The N50 for a Yoruba male DNA was 702 kb whereas the N50 for a European female DNA was 358 kb. Therefore, the strategy described here is suitable for haplotyping of a set of targeted regions as well as of the entire genome.

A highly scalable peptide-based assay system for proteomics.

We report a scalable and cost-effective technology for generating and screening high-complexity customizable peptide sets. The peptides are made as peptide-cDNA fusions by in vitro transcription/translation from pools of DNA templates generated by microarray-based synthesis. This approach enables large custom sets of peptides to be designed in silico, manufactured cost-effectively in parallel, and assayed efficiently in a multiplexed fashion. The utility of our peptide-cDNA fusion pools was demonstrated in two activity-based assays designed to discover protease and kinase substrates. In the protease assay, cleaved peptide substrates were separated from uncleaved and identified by digital sequencing of their cognate cDNAs. We screened the 3,011 amino acid HCV proteome for susceptibility to cleavage by the HCV NS3/4A protease and identified all 3 known trans cleavage sites with high specificity. In the kinase assay, peptide substrates phosphorylated by tyrosine kinases were captured and identified by sequencing of their cDNAs. We screened a pool of 3,243 peptides against Abl kinase and showed that phosphorylation events detected were specific and consistent with the known substrate preferences of Abl kinase. Our approach is scalable and adaptable to other protein-based assays.

The DExD/H box ATPase Dhh1 functions in translational repression, mRNA decay, and processing body dynamics.

Translation, storage, and degradation of messenger ribonucleic acids (mRNAs) are key steps in the posttranscriptional control of gene expression, but how mRNAs transit between these processes remains poorly understood. In this paper, we functionally characterized the DExD/H box adenosine triphosphatase (ATPase) Dhh1, a critical regulator of the cytoplasmic fate of mRNAs. Using mRNA tethering experiments in yeast, we showed that Dhh1 was sufficient to move an mRNA from an active state to translational repression. In actively dividing cells, translational repression was followed by mRNA decay; however, deleting components of the 5'-3' decay pathway uncoupled these processes. Whereas Dhh1's ATPase activity was not required to induce translational inhibition and mRNA decay when directly tethered to an mRNA, ATP hydrolysis regulated processing body dynamics and the release of Dhh1 from these RNA-protein granules. Our results place Dhh1 at the interface of translation and decay controlling whether an mRNA is translated, stored, or decayed.

Dynamic profiling of mRNA turnover reveals gene-specific and system-wide regulation of mRNA decay.

RNA levels are determined by the rates of both transcription and decay, and a mechanistic understanding of the complex networks regulating gene expression requires methods that allow dynamic measurements of transcription and decay in living cells with minimal perturbation. Here, we describe a metabolic pulse-chase labeling protocol using 4-thiouracil combined with large-scale RNA sequencing to determine decay rates of all mRNAs in Saccharomyces cerevisiae. Profiling in various growth and stress conditions reveals that mRNA turnover is highly regulated both for specific groups of transcripts and at the system-wide level. For example, acute glucose starvation induces global mRNA stabilization but increases the degradation of all 132 detected ribosomal protein mRNAs. This effect is transient and can be mimicked by inhibiting the target-of-rapamycin kinase. Half-lives of mRNAs critical for galactose (GAL) metabolism are also highly sensitive to changes in carbon source. The fast reduction of GAL transcripts in glucose requires their dramatically enhanced turnover, highlighting the importance of mRNA decay in the control of gene expression. The approach described here provides a general platform for the global analysis of mRNA turnover and transcription and can be applied to dissect gene expression programs in a wide range of organisms and conditions.

Interactions between muscle fibers and segment boundaries in zebrafish.

The most obvious segmental structures in the vertebrate embryo are somites: transient structures that give rise to vertebrae and much of the musculature. In zebrafish, most somitic cells give rise to long muscle fibers that are anchored to intersegmental boundaries. Therefore, this boundary is analogous to the mammalian tendon in that it transduces muscle-generated force to the skeletal system. We have investigated interactions between somite boundaries and muscle fibers. We define three stages of segment boundary formation. The first stage is the formation of the initial epithelial somite boundary. The second "transition" stage involves both the elongation of initially round muscle precursor cells and somite boundary maturation. The third stage is myotome boundary formation, where the boundary becomes rich in extracellular matrix and all muscle precursor cells have elongated to form long muscle fibers. It is known that formation of the initial epithelial somite boundary requires Notch signaling; vertebrate Notch pathway mutants show severe defects in somitogenesis. However, many zebrafish Notch pathway mutants are homozygous viable suggesting that segmentation of their larval and adult body plans at least partially recovers. We show that epithelial somite boundary formation and slow-twitch muscle morphogenesis are initially disrupted in after eight (aei) mutant embryos (which lack function of the Notch ligand, DeltaD); however, myotome boundaries form later ("recover") in a Hedgehog-dependent fashion. Inhibition of Hedgehog-induced slow muscle induction in aei/deltaD and deadly seven (des)/notch1a mutant embryos suggests that slow muscle is necessary for myotome boundary recovery in the absence of initial epithelial somite boundary formation. Because we have previously demonstrated that slow muscle migration triggers fast muscle cell elongation in zebrafish, we hypothesize that migrating slow muscle facilitates myotome boundary formation in aei/deltaD mutant embryos by patterning coordinated fast muscle cell elongation. In addition, we utilized genetic mosaic analysis to show that somite boundaries also function to limit the extent to which fast muscle cells can elongate. Combined, our results indicate that multiple interactions between somite boundaries and muscle fibers mediate zebrafish segmentation.

Yeast poly(A)-binding protein Pab1 shuttles between the nucleus and the cytoplasm and functions in mRNA export.

Pab1 is the major poly(A)-binding protein in yeast. It is a multifunctional protein that mediates many cellular functions associated with the 3'-poly(A)-tail of messenger RNAs. Here, we characterize Pab1 as an export cargo of the protein export factor Xpo1/Crm1. Pab1 is a major Xpo1/Crm1-interacting protein in yeast extracts and binds directly to Xpo1/Crm1 in a RanGTP-dependent manner. Pab1 shuttles rapidly between the nucleus and the cytoplasm and partially accumulates in the nucleus when the function of Xpo1/Crm1 is inhibited. However, Pab1 can also be exported by an alternative pathway, which is dependent on the MEX67-mRNA export pathway. Import of Pab1 is mediated by the import receptor Kap108/Sxm1 through a nuclear localization signal in its fourth RNA-binding domain. Interestingly, inhibition of Pab1's nuclear import causes a kinetic delay in the export of mRNA. Furthermore, the inviability of a pab1 deletion strain is suppressed by a mutation in the 5'-3' exoribonuclease RRP6, a component of the nuclear exosome. Therefore, nuclear Pab1 may be required for efficient mRNA export and may function in the quality control of mRNA in the nucleus.

Identification, function and structure of the mycobacterial sulfotransferase that initiates sulfolipid-1 biosynthesis.

Sulfolipid-1 (SL-1) is an abundant sulfated glycolipid and potential virulence factor found in Mycobacterium tuberculosis. SL-1 consists of a trehalose-2-sulfate (T2S) disaccharide elaborated with four lipids. We identified and characterized a conserved mycobacterial sulfotransferase, Stf0, which generates the T2S moiety of SL-1. Biochemical studies demonstrated that the enzyme requires unmodified trehalose as substrate and is sensitive to small structural perturbations of the disaccharide. Disruption of stf0 in Mycobacterium smegmatis and M. tuberculosis resulted in the loss of T2S and SL-1 formation, respectively. The structure of Stf0 at a resolution of 2.6 A reveals the molecular basis of trehalose recognition and a unique dimer configuration that encloses the substrate into a bipartite active site. These data provide strong evidence that Stf0 carries out the first committed step in the biosynthesis of SL-1 and establish a system for probing the role of SL-1 in M. tuberculosis infection.