Cell-free DNA testing of an extended range of chromosomal anomalies: clinical experience with 6,388 consecutive cases.

PURPOSE: Cell-free DNA (cfDNA) testing for fetal aneuploidies was broadly implemented for common trisomies and sex-chromosome anomalies (SCAs). However, such an approach identifies only 75 to 85% of clinically relevant aneuploidies. METHODS: We present a consecutive series of 6,388 cases, thus uncovering a broader array of aneuploidies, including the rare autosomal trisomies (RATs) and the maternally inherited deletion and duplication copy-number variations (CNVs), with complete and stratified follow-up by amniocentesis. Combined measurements of z-scores and the fetal fraction, in conjunction with fetal cfDNA enrichment, were used to stratify the likelihood of true and false results. RESULTS: We obtained an incremental diagnostic yield of 50%; RATs and CNVs were found to be significant causes of fetal pathology. Scrutinizing z-scores and the fetal fraction made it possible to distinguish the sources of false-negative results; predict the likelihood of false-positive results for major trisomies and SCAs; classify maternal mosaic SCAs and CNVs, preventing false-positive results; and robustly identify maternally inherited CNVs and detect recurrent genomic disorders as a standardized function of the fetal fraction. CONCLUSION: With the clinical pertinence of this broader detection scheme confirmed, we offer recommendations for its implementation.Genet Med 19 2, 169-175.

Duality of the murine CD8 compartment.

CD8αß plays crucial roles in the thymic selection, differentiation, and activation of some, but not all, CD8(+) T cells, whereas CD8αα does not. To investigate these roles, we produced mice that expressed transgene P14 T-cell receptor ß (TCRß) chain and CD8ß or did not (WT and KO mice, respectively). The primary CD8(+) T-cell response to acute lymphocytic choriomeningitis virus (LCMV) infection was predominantly D(b)/GP33 specific and CD8 independent in KO mice and was mostly CD8 dependent in WT mice. Cytotoxic T lymphocytes (CTL) from KO mice failed to mobilize intracellular Ca(2+) and to kill via perforin/granzyme. Their strong Fas/FasL-mediated cytotoxicity and IFN-γ response were signaled via a Ca(2+)-independent, PI3K-dependent pathway. This was also true for 15-20% of CD8-independent CTL found in WT mice. Conversely, the perforin/granzyme-mediated killing and IFN-γ response of CD8-dependent CTL were signaled via a Ca(2+), p56(lck), and nuclear factor of activated T cells-dependent pathway. Deep sequencing of millions of TCRα chain transcripts revealed that the TCR repertoires of preimmune CD8(+) T cells were highly diverse, but those of LCMV D(b)/GP33-specific CTL, especially from KO mice, were narrow. The immune repertoires exhibited biased use of Vα segments that encoded different complementary-determining region 1α (CDR1α) and CDR2α sequences. We suggest that TCR from WT CD8-independent T cells may engage MHC-peptide complexes in a manner unfavorable for efficient CD8 engagement and Ca(2+) signaling but permissive for Ca(2+)-independent, PI3K-dependent signaling. This duality of the CD8 compartment may provide organisms with broader protective immunity.

Highly diverse TCRα chain repertoire of pre-immune CD8⁺ T cells reveals new insights in gene recombination.

Although the T-cell receptor αδ (TCRαδ) locus harbours large libraries of variable (TRAV) and junctional (TRAJ) gene segments, according to previous studies the TCRα chain repertoire is of limited diversity due to restrictions imposed by sequential coordinate TRAV-TRAJ recombinations. By sequencing tens of millions of TCRα chain transcripts from naive mouse CD8(+) T cells, we observed a hugely diverse repertoire, comprising nearly all possible TRAV-TRAJ combinations. Our findings are not compatible with sequential coordinate gene recombination, but rather with a model in which contraction and DNA looping in the TCRαδ locus provide equal access to TRAV and TRAJ gene segments, similarly to that demonstrated for IgH gene recombination. Generation of the observed highly diverse TCRα chain repertoire necessitates deletion of failed attempts by thymic-positive selection and is essential for the formation of highly diverse TCRαß repertoires, capable of providing good protective immunity.

De novo finished 2.8 Mbp Staphylococcus aureus genome assembly from 100 bp short and long range paired-end reads.

MOTIVATION: Paired-end sequencing allows circumventing the shortness of the reads produced by second generation sequencers and is essential for de novo assembly of genomes. However, obtaining a finished genome from short reads is still an open challenge. We present an algorithm that exploits the pairing information issued from inserts of potentially any length. The method determines paths through an overlaps graph by using a constrained search tree. We also present a method that automatically determines suited overlaps cutoffs according to the contextual coverage, reducing thus the need for manual parameterization. Finally, we introduce an interactive mode that allows querying an assembly at targeted regions. RESULTS: We assess our methods by assembling two Staphylococcus aureus strains that were sequenced on the Illumina platform. Using 100 bp paired-end reads and minimal manual curation, we produce a finished genome sequence for the previously undescribed isolate SGH-10-168. AVAILABILITY AND IMPLEMENTATION: The presented algorithms are implemented in the standalone Edena software, freely available under the General Public License (GPLv3) at www.genomic.ch/edena.php.

Ultra-deep sequencing of foraminiferal microbarcodes unveils hidden richness of early monothalamous lineages in deep-sea sediments.

Deep-sea floors represent one of the largest and most complex ecosystems on Earth but remain essentially unexplored. The vastness and remoteness of this ecosystem make deep-sea sampling difficult, hampering traditional taxonomic observations and diversity assessment. This problem is particularly true in the case of the deep-sea meiofauna, which largely comprises small-sized, fragile, and difficult-to-identify metazoans and protists. Here, we introduce an ultra-deep sequencing-based metagenetic approach to examine the richness of benthic foraminifera, a principal component of deep-sea meiofauna. We used Illumina sequencing technology to assess foraminiferal richness in 31 unsieved deep-sea sediment samples from five distinct oceanic regions. We sequenced an extremely short fragment (36 bases) of the small subunit ribosomal DNA hypervariable region 37f, which has been shown to accurately distinguish foraminiferal species. In total, we obtained 495,978 unique sequences that were grouped into 1,643 operational taxonomic units, of which about half (841) could be reliably assigned to foraminifera. The vast majority of the operational taxonomic units (nearly 90%) were either assigned to early (ancient) lineages of soft-walled, single-chambered (monothalamous) foraminifera or remained undetermined and yet possibly belong to unknown early lineages. Contrasting with the classical view of multichambered taxa dominating foraminiferal assemblages, our work reflects an unexpected diversity of monothalamous lineages that are as yet unknown using conventional micropaleontological observations. Although we can only speculate about their morphology, the immense richness of deep-sea phylotypes revealed by this study suggests that ultra-deep sequencing can improve understanding of deep-sea benthic diversity considered until now as unknowable based on a traditional taxonomic approach.

Commercially available outbred mice for genome-wide association studies.

Genome-wide association studies using commercially available outbred mice can detect genes involved in phenotypes of biomedical interest. Useful populations need high-frequency alleles to ensure high power to detect quantitative trait loci (QTLs), low linkage disequilibrium between markers to obtain accurate mapping resolution, and an absence of population structure to prevent false positive associations. We surveyed 66 colonies for inbreeding, genetic diversity, and linkage disequilibrium, and we demonstrate that some have haplotype blocks of less than 100 Kb, enabling gene-level mapping resolution. The same alleles contribute to variation in different colonies, so that when mapping progress stalls in one, another can be used in its stead. Colonies are genetically diverse: 45% of the total genetic variation is attributable to differences between colonies. However, quantitative differences in allele frequencies, rather than the existence of private alleles, are responsible for these population differences. The colonies derive from a limited pool of ancestral haplotypes resembling those found in inbred strains: over 95% of sequence variants segregating in outbred populations are found in inbred strains. Consequently it is possible to impute the sequence of any mouse from a dense SNP map combined with inbred strain sequence data, which opens up the possibility of cataloguing and testing all variants for association, a situation that has so far eluded studies in completely outbred populations. We demonstrate the colonies' potential by identifying a deletion in the promoter of H2-Ea as the molecular change that strongly contributes to setting the ratio of CD4+ and CD8+ lymphocytes.

Spliced leader trapping reveals widespread alternative splicing patterns in the highly dynamic transcriptome of Trypanosoma brucei.

Trans-splicing of leader sequences onto the 5'ends of mRNAs is a widespread phenomenon in protozoa, nematodes and some chordates. Using parallel sequencing we have developed a method to simultaneously map 5'splice sites and analyze the corresponding gene expression profile, that we term spliced leader trapping (SLT). The method can be applied to any organism with a sequenced genome and trans-splicing of a conserved leader sequence. We analyzed the expression profiles and splicing patterns of bloodstream and insect forms of the parasite Trypanosoma brucei. We detected the 5' splice sites of 85% of the annotated protein-coding genes and, contrary to previous reports, found up to 40% of transcripts to be differentially expressed. Furthermore, we discovered more than 2500 alternative splicing events, many of which appear to be stage-regulated. Based on our findings we hypothesize that alternatively spliced transcripts present a new means of regulating gene expression and could potentially contribute to protein diversity in the parasite. The entire dataset can be accessed online at TriTrypDB or through: http://splicer.unibe.ch/.

Profile of small interfering RNAs from cotton plants infected with the polerovirus Cotton leafroll dwarf virus.

BACKGROUND: In response to infection, viral genomes are processed by Dicer-like (DCL) ribonuclease proteins into viral small RNAs (vsRNAs) of discrete sizes. vsRNAs are then used as guides for silencing the viral genome. The profile of vsRNAs produced during the infection process has been extensively studied for some groups of viruses. However, nothing is known about the vsRNAs produced during infections of members of the economically important family Luteoviridae, a group of phloem-restricted viruses. Here, we report the characterization of a population of vsRNAs from cotton plants infected with Cotton leafroll dwarf virus (CLRDV), a member of the genus Polerovirus, family Luteoviridae. RESULTS: Deep sequencing of small RNAs (sRNAs) from leaves of CLRDV-infected cotton plants revealed that the vsRNAs were 21- to 24-nucleotides (nt) long and that their sequences matched the viral genome, with higher frequencies of matches in the 3- region. There were equivalent amounts of sense and antisense vsRNAs, and the 22-nt class of small RNAs was predominant. During infection, cotton Dcl transcripts appeared to be up-regulated, while Dcl2 appeared to be down-regulated. CONCLUSIONS: This is the first report on the profile of sRNAs in a plant infected with a virus from the family Luteoviridae. Our sequence data strongly suggest that virus-derived double-stranded RNA functions as one of the main precursors of vsRNAs. Judging by the profiled size classes, all cotton DCLs might be working to silence the virus. The possible causes for the unexpectedly high accumulation of 22-nt vsRNAs are discussed. CLRDV is the causal agent of Cotton blue disease, which occurs worldwide. Our results are an important contribution for understanding the molecular mechanisms involved in this and related diseases.

Whole genome sequencing of Saccharomyces cerevisiae: from genotype to phenotype for improved metabolic engineering applications.

BACKGROUND: The need for rapid and efficient microbial cell factory design and construction are possible through the enabling technology, metabolic engineering, which is now being facilitated by systems biology approaches. Metabolic engineering is often complimented by directed evolution, where selective pressure is applied to a partially genetically engineered strain to confer a desirable phenotype. The exact genetic modification or resulting genotype that leads to the improved phenotype is often not identified or understood to enable further metabolic engineering. RESULTS: In this work we performed whole genome high-throughput sequencing and annotation can be used to identify single nucleotide polymorphisms (SNPs) between Saccharomyces cerevisiae strains S288c and CEN.PK113-7D. The yeast strain S288c was the first eukaryote sequenced, serving as the reference genome for the Saccharomyces Genome Database, while CEN.PK113-7D is a preferred laboratory strain for industrial biotechnology research. A total of 13,787 high-quality SNPs were detected between both strains (reference strain: S288c). Considering only metabolic genes (782 of 5,596 annotated genes), a total of 219 metabolism specific SNPs are distributed across 158 metabolic genes, with 85 of the SNPs being nonsynonymous (e.g., encoding amino acid modifications). Amongst metabolic SNPs detected, there was pathway enrichment in the galactose uptake pathway (GAL1, GAL10) and ergosterol biosynthetic pathway (ERG8, ERG9). Physiological characterization confirmed a strong deficiency in galactose uptake and metabolism in S288c compared to CEN.PK113-7D, and similarly, ergosterol content in CEN.PK113-7D was significantly higher in both glucose and galactose supplemented cultivations compared to S288c. Furthermore, DNA microarray profiling of S288c and CEN.PK113-7D in both glucose and galactose batch cultures did not provide a clear hypothesis for major phenotypes observed, suggesting that genotype to phenotype correlations are manifested post-transcriptionally or post-translationally either through protein concentration and/or function. CONCLUSIONS: With an intensifying need for microbial cell factories that produce a wide array of target compounds, whole genome high-throughput sequencing and annotation for SNP detection can aid in better reducing and defining the metabolic landscape. This work demonstrates direct correlations between genotype and phenotype that provides clear and high-probability of success metabolic engineering targets. The genome sequence, annotation, and a SNP viewer of CEN.PK113-7D are deposited at http://www.sysbio.se/cenpk.

A Key Regulator of the Glycolytic and Gluconeogenic Central Metabolic Pathways in Sinorhizobium meliloti.

The order Rhizobiales contains numerous agriculturally, biotechnologically, and medically important bacteria, including the rhizobia, and the genera Agrobacterium, Brucella, and Methylobacterium, among others. These organisms tend to be metabolically versatile, but there has been relatively little investigation into the regulation of their central carbon metabolic pathways. Here, RNA-sequencing and promoter fusion data are presented to show that the PckR protein is a key regulator of central carbon metabolism in Sinorhizobium meliloti; during growth with gluconeogenic substrates, PckR represses expression of the complete Entner-Doudoroff glycolytic pathway and induces expression of the pckA and fbaB gluconeogenic genes. Electrophoretic mobility shift assays indicate that PckR binds an imperfect palindromic sequence that overlaps the promoter or transcriptional start site in the negatively regulated promoters, or is present in tandem upstream the promoter motifs in the positively regulated promoters. Genetic and in vitro electrophoretic mobility shift assay experiments suggest that elevated concentrations of a PckR effector ligand results in the dissociation of PckR from its target binding site, and evidence is presented that suggests phosphoenolpyruvate may function as the effector. Characterization of missense pckR alleles identified three conserved residues important for increasing the affinity of PckR for its cognate effector molecule. Bioinformatics analyses illustrates that PckR is limited to a narrow phylogenetic range consisting of the Rhizobiaceae, Phyllobacteriaceae, Brucellaceae, and Bartonellaceae families. These data provide novel insights into the regulation of the core carbon metabolic pathways of this pertinent group of α-proteobacteria.

Technical Insights into Highly Sensitive Isolation and Molecular Characterization of Fixed and Live Circulating Tumor Cells for Early Detection of Tumor Invasion.

Circulating Tumor Cells (CTC) and Circulating Tumor Microemboli (CTM) are Circulating Rare Cells (CRC) which herald tumor invasion and are expected to provide an opportunity to improve the management of cancer patients. An unsolved technical issue in the CTC field is how to obtain highly sensitive and unbiased collection of these fragile and heterogeneous cells, in both live and fixed form, for their molecular study when they are extremely rare, particularly at the beginning of the invasion process. We report on a new protocol to enrich from blood live CTC using ISET® (Isolation by SizE of Tumor/Trophoblastic Cells), an open system originally developed for marker-independent isolation of fixed tumor cells. We have assessed the impact of our new enrichment method on live tumor cells antigen expression, cytoskeleton structure, cell viability and ability to expand in culture. We have also explored the ISET® in vitro performance to collect intact fixed and live cancer cells by using spiking analyses with extremely low number of fluorescent cultured cells. We describe results consistently showing the feasibility of isolating fixed and live tumor cells with a Lower Limit of Detection (LLOD) of one cancer cell per 10 mL of blood and a sensitivity at LLOD ranging from 83 to 100%. This very high sensitivity threshold can be maintained when plasma is collected before tumor cells isolation. Finally, we have performed a comparative next generation sequencing (NGS) analysis of tumor cells before and after isolation from blood and culture. We established the feasibility of NGS analysis of single live and fixed tumor cells enriched from blood by our system. This study provides new protocols for detection and characterization of CTC collected from blood at the very early steps of tumor invasion.

The Sinorhizobium meliloti glycine betaine biosynthetic genes (betlCBA) are induced by choline and highly expressed in bacteroids.

The symbiotic soil bacterium Sinorhizobium meliloti has the capacity to synthesize the osmoprotectant glycine betaine from choline-O-sulfate and choline. This pathway is encoded by the betICBA locus, which comprises a regulatory gene, betI, and three structural genes, betC (choline sulfatase), betB (betaine aldehyde dehydrogenase), and betA (choline dehydrogenase). Here, we report that betICBA genes constitute a single operon, despite the existence of intergenic regions containing mosaic elements between betI and betC, and betB and betA. The regulation of the bet operon was investigated by using transcriptional lacZ (beta-galactosidase) fusions and has revealed a strong induction by choline at concentrations as low as 25 microM and to a lesser extent by choline-O-sulfate and acetylcholine but not by osmotic stress or oxygen. BetI is a repressor of the bet transcription in the absence of choline, and a nucleotide sequence of dyad symmetry upstream of betI was identified as a putative betI box. Measurements of intracellular pools of choline, well correlated with beta-galactosidase activities, strongly suggested that BetI senses the endogenous choline pool that modulates the intensity of BetI repression. In contrast to Escherichia coli, BetI did not repress choline transport. During symbiosis with Medicago sativa, S. meliloti bet gene expression was observed within the infection threads, in young and in mature nodules. The existence of free choline in nodule cytosol, peribacteroid space, and bacteroids was demonstrated, and the data suggest that bet regulation in planta is mediated by BetI repression, as in free-living cells. Neither Nod nor Fix phenotypes were significantly impaired in a betI::omega mutant, indicating that glycine betaine biosynthesis from choline is not crucial for nodulation and nitrogen fixation.

A new gene in A. rubens: A sea star Ig kappa gene.

The sea star Asterias rubens reacts specifically to the antigen:HRP (horse-radish peroxydase) and produces an antibody anti-HRP. We previously identified a candidate Ig kappa gene corresponding to this manuscript. We show now the gene referred to as: "sea star Ig kappa gene in its specificity".

Whole-Genome Shotgun Sequence of Arthrospira platensis Strain Paraca, a Cultivated and Edible Cyanobacterium.

Here we report the whole-genome shotgun sequence of a Peruvian strain of Arthrospira platensis (Paraca), a cultivated and edible haloalkaliphilic cyanobacterium of great scientific, technical, and economic potential.

Whole-Genome Shotgun Sequence of Pseudomonas viridiflava, a Bacterium Species Pathogenic to Ararabidopsis thaliana.

We report here the first whole-genome shotgun sequence of Pseudomonas viridiflava strain UASWS38, a bacterium species pathogenic to the biological model plant Ararabidopsis thaliana but also usable as a biological control agent and thus of great scientific interest for understanding the genetics of plant-microbe interactions.

Analysis of the salivary microbiome using culture-independent techniques.

BACKGROUND: The salivary microbiota is a potential diagnostic indicator of several diseases. Culture-independent techniques are required to study the salivary microbial community since many of its members have not been cultivated. METHODS: We explored the bacterial community composition in the saliva sample using metagenomic whole genome shotgun (WGS) sequencing, the extraction of 16S rRNA gene fragments from metagenomic sequences (16S-WGS) and high-throughput sequencing of PCR-amplified bacterial 16S rDNA gene (16S-HTS) regions V1 and V3. RESULTS: The hierarchical clustering of data based on the relative abundance of bacterial genera revealed that distances between 16S-HTS datasets for V1 and V3 regions were greater than those obtained for the same V region with different numbers of PCR cycles. Datasets generated by 16S-HTS and 16S-WGS were even more distant. Finally, comparison of WGS and 16S-based datasets revealed the highest dissimilarity.The analysis of the 16S-HTS, WGS and 16S-WGS datasets revealed 206, 56 and 39 bacterial genera, respectively, 124 of which have not been previously identified in salivary microbiomes. A large fraction of DNA extracted from saliva corresponded to human DNA. Based on sequence similarity search against completely sequenced genomes, bacterial and viral sequences represented 0.73% and 0.0036% of the salivary metagenome, respectively. Several sequence reads were identified as parts of the human herpesvirus 7. CONCLUSIONS: Analysis of the salivary metagenome may have implications in diagnostics e.g. in detection of microorganisms and viruses without designing specific tests for each pathogen.

Chromatin immunoprecipitation (ChIP) of plant transcription factors followed by sequencing (ChIP-SEQ) or hybridization to whole genome arrays (ChIP-CHIP).

Chromatin immunoprecipitation (ChIP) is a powerful technique to study interactions between transcription factors (TFs) and DNA in vivo. For genome-wide de novo discovery of TF-binding sites, the DNA that is obtained in ChIP experiments needs to be processed for sequence identification. The sequences can be identified by direct sequencing (ChIP-SEQ) or hybridization to microarrays (ChIP-CHIP). Given the small amounts of DNA that are usually obtained in ChIP experiments, successful and reproducible sample processing is challenging. Here we provide a detailed procedure for ChIP of plant TFs, as well as protocols for sample preparation for ChIP-SEQ and for ChIP-CHIP. Our ChIP procedure is optimized for high signal-to-noise ratio starting with tissue fixation, followed by nuclei isolation, immunoprecipitation, DNA amplification and purification. We also provide a guide for primary data analysis of ChIP-SEQ data. The complete protocol for ChIP-SEQ/ChIP-CHIP sample preparation starting from plant harvest takes approximately 7 d.

Metagenomic study of the oral microbiota by Illumina high-throughput sequencing.

To date, metagenomic studies have relied on the utilization and analysis of reads obtained using 454 pyrosequencing to replace conventional Sanger sequencing. After extensively scanning the 16S ribosomal RNA (rRNA) gene, we identified the V5 hypervariable region as a short region providing reliable identification of bacterial sequences available in public databases such as the Human Oral Microbiome Database. We amplified samples from the oral cavity of three healthy individuals using primers covering an approximately 82-base segment of the V5 loop, and sequenced using the Illumina technology in a single orientation. We identified 135 genera or higher taxonomic ranks from the resulting 1,373,824 sequences. While the abundances of the most common phyla (Firmicutes, Proteobacteria, Actinobacteria, Fusobacteria and TM7) are largely comparable to previous studies, Bacteroidetes were less present. Potential sources for this difference include classification bias in this region of the 16S rRNA gene, human sample variation, sample preparation and primer bias. Using an Illumina sequencing approach, we achieved a much greater depth of coverage than previous oral microbiota studies, allowing us to identify several taxa not yet discovered in these types of samples, and to assess that at least 30,000 additional reads would be required to identify only one additional phylotype. The evolution of high-throughput sequencing technologies, and their subsequent improvements in read length enable the utilization of different platforms for studying communities of complex flora. Access to large amounts of data is already leading to a better representation of sample diversity at a reasonable cost.

De novo bacterial genome sequencing: millions of very short reads assembled on a desktop computer.

Novel high-throughput DNA sequencing technologies allow researchers to characterize a bacterial genome during a single experiment and at a moderate cost. However, the increase in sequencing throughput that is allowed by using such platforms is obtained at the expense of individual sequence read length, which must be assembled into longer contigs to be exploitable. This study focuses on the Illumina sequencing platform that produces millions of very short sequences that are 35 bases in length. We propose a de novo assembler software that is dedicated to process such data. Based on a classical overlap graph representation and on the detection of potentially spurious reads, our software generates a set of accurate contigs of several kilobases that cover most of the bacterial genome. The assembly results were validated by comparing data sets that were obtained experimentally for Staphylococcus aureus strain MW2 and Helicobacter acinonychis strain Sheeba with that of their published genomes acquired by conventional sequencing of 1.5- to 3.0-kb fragments. We also provide indications that the broad coverage achieved by high-throughput sequencing might allow for the detection of clonal polymorphisms in the set of DNA molecules being sequenced.