High-Throughput Massively Parallel Sequencing for Fetal Aneuploidy Detection from Maternal Plasma

Background

Circulating cell-free (ccf) fetal DNA comprises 3–20% of all the cell-free DNA present in maternal plasma. Numerous research and clinical studies have described the analysis of ccf DNA using next generation sequencing for the detection of fetal aneuploidies with high sensitivity and specificity. We sought to extend the utility of this approach by assessing semi-automated library preparation, higher sample multiplexing during sequencing, and improved bioinformatic tools to enable a higher throughput, more efficient assay while maintaining or improving clinical performance.

Methods

Whole blood (10mL) was collected from pregnant female donors and plasma separated using centrifugation. Ccf DNA was extracted using column-based methods. Libraries were prepared using an optimized semi-automated library preparation method and sequenced on an Illumina HiSeq2000 sequencer in a 12-plex format. Z-scores were calculated for affected chromosomes using a robust method after normalization and genomic segment filtering. Classification was based upon a standard normal transformed cutoff value of z = 3 for chromosome 21 and z = 3.95 for chromosomes 18 and 13.

Results

Two parallel assay development studies using a total of more than 1900 ccf DNA samples were performed to evaluate the technical feasibility of automating library preparation and increasing the sample multiplexing level. These processes were subsequently combined and a study of 1587 samples was completed to verify the stability of the process-optimized assay. Finally, an unblinded clinical evaluation of 1269 euploid and aneuploid samples utilizing this high-throughput assay coupled to improved bioinformatic procedures was performed. We were able to correctly detect all aneuploid cases with extremely low false positive rates of 0.09%, <0.01%, and 0.08% for trisomies 21, 18, and 13, respectively.

Conclusions

These data suggest that the developed laboratory methods in concert with improved bioinformatic approaches enable higher sample throughput while maintaining high classification accuracy.     

Yeast ribosomal protein L7 and its homologue Rlp7 are simultaneously present at distinct sites on pre-60S ribosomal particles

Ribosome biogenesis requires >300 assembly factors in Saccharomyces cerevisiae. Ribosome assembly factors Imp3, Mrt4, Rlp7 and Rlp24 have sequence similarity to ribosomal proteins S9, P0, L7 and L24, suggesting that these pre-ribosomal factors could be placeholders that prevent premature assembly of the corresponding ribosomal proteins to nascent ribosomes. However, we found L7 to be a highly specific component of Rlp7-associated complexes, revealing that the two proteins can bind simultaneously to pre-ribosomal particles. Cross-linking and cDNA analysis experiments showed that Rlp7 binds to the ITS2 region of 27S pre-rRNAs, at two sites, in helix III and in a region adjacent to the pre-rRNA processing sites C₁ and E. However, L7 binds to mature 25S and 5S rRNAs and cross-linked predominantly to helix ES7^Lb within 25S rRNA. Thus, despite their predicted structural similarity, our data show that Rlp7 and L7 clearly bind at different positions on the same pre-60S particles. Our results also suggest that Rlp7 facilitates the formation of the hairpin structure of ITS2 during 60S ribosomal subunit maturation.     

Folding and stability of outer membrane protein A (OmpA) from Escherichia coli in an amphipathic polymer, amphipol A8-35

Amphipols are a class of amphipathic polymers designed to maintain membrane proteins in aqueous solutions in the absence of detergents. Denatured β-barrel membrane proteins, like outer membrane proteins OmpA from Escherichia coli and FomA from Fusobacterium nucleatum, can be folded by dilution of the denaturant urea in the presence of amphipol A8-35. Here, the folding kinetics and stability of OmpA in A8-35 have been investigated. Folding is well described by two parallel first-order processes, whose half-times, ~5 and ~70 min, respectively, are independent of A8-35 concentration. The faster process contributed ~55–64 % to OmpA folding. Folding into A8-35 was faster than into dioleoylphosphatidylcholine bilayers and complete at ratios as low as ~0.17 g/g A8-35/OmpA, corresponding to ~1–2 A8-35 molecules per OmpA. Activation energies were determined from the temperature dependence of folding kinetics, monitored both by electrophoresis, which reports on the formation of stable OmpA tertiary structure, and by fluorescence spectroscopy, which reflects changes in the environment of tryptophan side chains. The two methods yielded consistent estimates, namely ~5–9 kJ/mol for the fast process and ~29–37 kJ/mol for the slow one, which is lower than is observed for OmpA folding into dioleoylphosphatidylcholine bilayers. Folding and unfolding titrations with urea demonstrated that OmpA folding into A8-35 is reversible and that amphipol-refolded OmpA is thermodynamically stable at room temperature. Comparison of activation energies for folding and unfolding in A8-35 versus detergent indicates that stabilization of A8-35-trapped OmpA against denaturation by urea is a kinetic, not a thermodynamic phenomenon.     

Regulation of PSII function in Cyanothece sp. ATCC 51142 during a light–dark cycle

Photosynthesis Research - Cyanobacteria, as well as green algae and higher plants, have highly conserved photosynthetic machinery. Cyanothece sp. ATCC 51142 is a unicellular, aerobic, diazotrophic...     

Sequential protein association with nascent 60S ribosomal particles

Ribosome biogenesis in eukaryotes depends on the coordinated action of ribosomal and nonribosomal proteins that guide the assembly of preribosomal particles. These intermediate particles follow a maturation pathway in which important changes in their protein composition occur. The mechanisms involved in the coordinated assembly of the ribosomal particles are poorly understood. We show here that the association of preribosomal factors with pre-60S complexes depends on the presence of earlier factors, a phenomenon essential for ribosome biogenesis. The analysis of the composition of purified preribosomal complexes blocked in maturation at specific steps allowed us to propose a model of sequential protein association with, and dissociation from, early pre-60S complexes for several preribosomal factors such as Mak11, Ssf1, Rlp24, Nog1, and Nog2. The presence of either Ssf1 or Nog2 in complexes that contain the 27SB pre-rRNA defines novel, distinct pre-60S particles that contain the same pre-rRNA intermediates and that differ only by the presence or absence of specific proteins. Physical and functional interactions between Rlp24 and Nog1 revealed that the assembly steps are, at least in part, mediated by direct protein-protein interactions.     

Interaction of discoidin domain receptor 1 with collagen type 1

Discoidin domain receptor 1 (DDR1) is a widely expressed tyrosine kinase receptor which binds to and gets activated by collagens including collagen type 1. Little is understood about the interaction of DDR1 with collagen and its possible functional implications. Here, we elucidate the binding pattern of the DDR1 extracellular domain (ECD) to collagen type 1 and its impact on collagen fibrillogenesis. Our in vitro assays utilized DDR1-Fc fusion proteins, which contain only the ECD of DDR1. Using surface plasmon resonance, we confirmed that further oligomerization of DDR1-Fc (by means of anti-Fc antibody) greatly enhances its binding to immobilized collagen type 1. Single-molecule imaging by means of atomic force microscopy revealed that DDR1 oligomers bound at overlapping or adjacent collagen molecules and were nearly absent on isolated collagen molecules. Interaction of DDR1 oligomers with collagen was found to modulate collagen fibrillogenesis both in vitro and in cell-based assays. Collagen fibers formed in the presence of DDR1 had a larger average diameter, were more cross-linked and lacked the native banded structure. The presence of DDR1 ECD resulted in "locking" of collagen molecules in an incomplete fibrillar state both in vitro and on surfaces of cells overexpressing DDR1. Our results signify an important functional role of the DDR1 ECD, which occurs naturally in kinase-dead isoforms of DDR1 and as a shedded soluble protein. The modulation of collagen fibrillogenesis by the DDR1 ECD elucidates a novel mechanism of collagen regulation by DDR1.     

The p21-activated protein kinase inhibitor Skb15 and its budding yeast homologue are 60S ribosome assembly factors

Ribosome biogenesis is driven by a large number of preribosomal factors that associate with and dissociate from the preribosomal particles along the maturation pathway. We have previously shown that budding yeast Mak11, whose homologues in other eukaryotes were described as modulating a p21-activated protein kinase function, accumulates in Rlp24-associated pre-60S complexes when their maturation is impeded in Saccharomyces cerevisiae. The functional inactivation of WD40 repeat protein Mak11 interfered with the 60S rRNA maturation, led to a cell cycle delay in G(1), and blocked green fluorescent protein-tagged Rpl25 in the nucleoli of yeast cells, indicating an early role of Mak11 in ribosome assembly. Surprisingly, Mak11 inactivation also led to a dramatic destabilization of Rlp24. The suppression of the thermosensitive phenotype of a mak11 mutant by RLP24 overexpression and a direct in vitro interaction between Rlp24 and Mak11 suggest that Mak11 acts as an Rlp24 cofactor during early steps of 60S ribosomal subunit assembly. Moreover, we found that Skb15, the Mak11 homologue in Schizosaccharomyces pombe, also associated with preribosomes and affected 60S biogenesis in fission yeast. It is thus likely that the previously observed phenotypes for MAK11 homologues in other eukaryotes are secondary to the main function of these proteins in ribosome formation.     

Whole genome bisulfite sequencing of cell-free DNA and its cellular contributors uncovers placenta hypomethylated domains

BackgroundCirculating cell-free fetal DNA has enabled non-invasive prenatal fetal aneuploidy testing without direct discrimination of the maternal and fetal DNA. Testing may be improved by specifically enriching the sample material for fetal DNA. DNA methylation may allow for such a separation of DNA; however, this depends on knowledge of the methylomes of circulating cell-free DNA and its cellular contributors.ResultsWe perform whole genome bisulfite sequencing on a set of unmatched samples including circulating cell-free DNA from non-pregnant and pregnant female donors and genomic DNA from maternal buffy coat and placenta samples. We find CpG cytosines within longer fragments are more likely to be methylated. Comparison of the methylomes of placenta and non-pregnant circulating cell-free DNA reveal many of the 51,259 identified differentially methylated regions are located in domains exhibiting consistent placenta hypomethylation across millions of consecutive bases. We find these placenta hypomethylated domains are consistently located within regions exhibiting low CpG and gene density. Differentially methylated regions identified when comparing placenta to non-pregnant circulating cell-free DNA are recapitulated in pregnant circulating cell-free DNA, confirming the ability to detect differential methylation in circulating cell-free DNA mixtures.ConclusionsWe generate methylome maps for four sample types at single-base resolution, identify a link between DNA methylation and fragment length in circulating cell-free DNA, identify differentially methylated regions between sample groups, and uncover the presence of megabase-size placenta hypomethylated domains.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0645-x) contains supplementary material, which is available to authorized users.     

inSilicoDb: an R/Bioconductor package for accessing human Affymetrix expert-curated datasets from GEO

Microarray technology has become an integral part of biomedical research and increasing amounts of datasets become available through public repositories. However, re-use of these datasets is severely hindered by unstructured, missing or incorrect biological samples information; as well as the wide variety of preprocessing methods in use. The inSilicoDb R/Bioconductor package is a command-line front-end to the InSilico DB, a web-based database currently containing 86 104 expert-curated human Affymetrix expression profiles compiled from 1937 GEO repository series. The use of this package builds on the Bioconductor project's focus on reproducibility by enabling a clear workflow in which not only analysis, but also the retrieval of verified data is supported.