Library preparation for next generation sequencing: A review of automation strategies.

Next generation sequencing is in the process of evolving from a technology used for research purposes to one which is applied in clinical diagnostics. Recently introduced high throughput and benchtop instruments offer fully automated sequencing runs at a lower cost per base and faster assay times. In turn, the complex and cumbersome library preparation, starting with isolated nucleic acids and resulting in amplified and barcoded DNA with sequencing adapters, has been identified as a significant bottleneck. Library preparation protocols usually consist of a multistep process and require costly reagents and substantial hands-on-time. Considerable emphasis will need to be placed on standardisation to ensure robustness and reproducibility. This review presents an overview of the current state of automation of library preparation for next generation sequencing. Major challenges associated with library preparation are outlined and different automation strategies are classified according to their functional principle. Pipetting workstations allow high-throughput processing yet offer limited flexibility, whereas microfluidic solutions offer great potential due to miniaturisation and decreased investment costs. For the emerging field of single cell transcriptomics for example, microfluidics enable singularisation of tens of thousands of cells in nanolitre droplets and barcoding of the RNA to assign each nucleic acid sequence to its cell of origin. Finally, two applications, the characterisation of bacterial pathogens and the sequencing within human immunogenetics, are outlined and benefits of automation are discussed.

Review on pneumatic operations in centrifugal microfluidics.

Centrifugal microfluidics allows for miniaturization, automation and parallelization of laboratory workflows. The fact that centrifugal forces are always directed radially outwards has been considered a main drawback for the implementation of complex workflows leading to the requirement of additional actuation forces for pumping, valving and switching. In this work, we review and discuss the combination of centrifugal with pneumatic forces which enables transport of even complex liquids in any direction on centrifugal systems, provides actuation for valving and switching, offers alternatives for mixing and enables accurate and precise metering and aliquoting. In addition, pneumatics can be employed for timing to carry out any of the above listed unit operations in a sequential and cascaded manner. Firstly, different methods to generate pneumatic pressures are discussed. Then, unit operations and applications that employ pneumatics are reviewed. Finally, a tutorial section discusses two examples to provide insight into the design process. The first tutorial explains a comparatively simple implementation of a pneumatic siphon valve and provides a workflow to derive optimum design parameters. The second tutorial discusses cascaded pneumatic operations consisting of temperature change rate actuated valving and subsequent pneumatic pumping. In conclusion, combining pneumatic actuation with centrifugal microfluidics allows for the design of robust fluidic networks with simple fluidic structures that are implemented in a monolithic fashion. No coatings are required and the overall demands on manufacturing are comparatively low. We see the combination of centrifugal forces with pneumatic actuation as a key enabling technology to facilitate compact and robust automation of biochemical analysis.

Production of recombinant human tektin 1, 2, and 4 and in vitro assembly of human tektin 1.

Proteins predicted to be composed of large stretches of coiled-coil structure have often proven difficult to crystallize for structural determination. We have successfully applied EPR spectroscopic techniques to the study of the structure and assembly of full-length human vimentin assembled into native 11 nm filaments, in physiologic solution, circumventing the limitations of crystallizing shorter peptide sequences. Tektins are a small family of highly alpha helical filamentous proteins found in the doublet microtubules of cilia and related structures. Tektins exhibit several similarities to intermediate filaments (IFs): moderate molecular weight, highly alpha helical, hypothesized to be coiled-coil, and homo- and heteromeric assembly into long smooth filaments. In this report, we show the application of IF research methodologies to the study of tektin structure and assembly. To begin in vitro studies, expression constructs for human tektins 1, 2, and 4 were synthesized. Recombinant tektins were produced in E. coli and purified by chromatography. Preparations of tektin 1 successfully formed filaments. The recombinant human tektin 1 was used to produce antibodies which recognized an antigen in mouse testes, most likely present in sperm flagella. Finally, we report the creation of seven mutants to analyze predictions of coiled-coil structure in the rod 1A domain of tektin 1. Although this region is predicted to be coiled-coil, our EPR analysis does not reflect the parallel, in register, coiled-coil structure as demonstrated in vimentin and kinesin. These results document that tektin can be successfully expressed and assembled in vitro, and that SDSL EPR techniques can be used for structural analysis.

Molecular profiling of a 6-hydroxydopamine model of Parkinson's disease.

Convection enhanced delivery of 6-hydroxydopamine (6-OHDA) to the rat striatum results in a model of Parkinson's disease. An important feature of this unilateral model is the progressive loss of dopaminergic (DA) neurons over the course of several weeks. To improve the understanding of this model, gene expression changes in the substantia nigra, which contains the DA neuron cell bodies, and the striatum, which contains the DA neuron synaptic terminals, were examined using DNA microarrays. Samples were collected and behavior was analyzed from vehicle and toxin treated animals at 3 days, 1 week, 2 weeks and 4 weeks following 6-OHDA treatment. Tissue DA content was determined and samples from animals which exhibited a substantial depletion of striatal DA were included in the subsequent gene expression analysis. The results of the gene expression analysis indicated that 6-OHDA elicits a vigorous inflammatory response, comprised of several distinct pathways, in the striatum at the earliest time point tested. In contrast, relatively few gene expression changes were observed in the SN at the 3-day time point. In both tissues examined there was evidence for a vigorous inflammatory response at the 1- and 2-week time points, which was substantially diminished by the 4-week time point. Inflammation plays a prominent role in the 6-OHDA model of Parkinson's disease.

Phylogenetic analysis of the tenascin gene family: evidence of origin early in the chordate lineage.

BACKGROUND: Tenascins are a family of glycoproteins found primarily in the extracellular matrix of embryos where they help to regulate cell proliferation, adhesion and migration. In order to learn more about their origins and relationships to each other, as well as to clarify the nomenclature used to describe them, the tenascin genes of the urochordate Ciona intestinalis, the pufferfish Tetraodon nigroviridis and Takifugu rubripes and the frog Xenopus tropicalis were identified and their gene organization and predicted protein products compared with the previously characterized tenascins of amniotes. RESULTS: A single tenascin gene was identified in the genome of C. intestinalis that encodes a polypeptide with domain features common to all vertebrate tenascins. Both pufferfish genomes encode five tenascin genes: two tenascin-C paralogs, a tenascin-R with domain organization identical to mammalian and avian tenascin-R, a small tenascin-X with previously undescribed GK repeats, and a tenascin-W. Four tenascin genes corresponding to tenascin-C, tenascin-R, tenascin-X and tenascin-W were also identified in the X. tropicalis genome. Multiple sequence alignment reveals that differences in the size of tenascin-W from various vertebrate classes can be explained by duplications of specific fibronectin type III domains. The duplicated domains are encoded on single exons and contain putative integrin-binding motifs. A phylogenetic tree based on the predicted amino acid sequences of the fibrinogen-related domains demonstrates that tenascin-C and tenascin-R are the most closely related vertebrate tenascins, with the most conserved repeat and domain organization. Taking all lines of evidence together, the data show that the tenascins referred to as tenascin-Y and tenascin-N are actually members of the tenascin-X and tenascin-W gene families, respectively. CONCLUSION: The presence of a tenascin gene in urochordates but not other invertebrate phyla suggests that tenascins may be specific to chordates. Later genomic duplication events led to the appearance of four family members in vertebrates: tenascin-C, tenascin-R, tenascin-W and tenascin-X.

Cloning and characterization of a Moraxella bovis cytotoxin gene.

OBJECTIVE: To identify the Moraxella bovis cytotoxin gene. PROCEDURE: Hemolytic and nonhemolytic strains of M. bovis were compared by use of western blotting to identify proteins unique to hemolytic strains. Oligonucleotide primers, designed on the basis of amino acid sequences of 2 tryptic peptides derived from 1 such protein and conserved regions of the C and B genes from members of the repeats in the structural toxin (RTX) family of bacterial toxins, were used to amplify cytotoxin-specific genes from M. bovis genomic DNA. Recombinant proteins were expressed, and antisera against these proteins were produced in rabbits. RESULTS: Several proteins ranging in molecular mass from 55 to 75 kd were unique to the hemolytic strain. An open reading frame encoding a 927-amino acid protein with a predicted molecular mass of 98.8 kd was amplified from M. bovis genomic DNA. The deduced amino acid sequence encoded by this open reading frame was homologous to RTX toxins. Antisera against the recombinant carboxy terminus encoded by this open reading frame neutralized hemolytic and cytolytic activities of native M. bovis cytotoxin. CONCLUSIONS AND CLINICAL RELEVANCE: A gene was identified in M bovis that encodes a protein with sequence homology to other RTX toxins. Results of cytotoxin neutralization assays support the hypothesis that M. bovis cytotoxin is encoded by this gene and belongs in the RTX family of bacterial exoproteins. Identification of this gene and expression of recombinant cytotoxin could facilitate the development of improved vaccines against infectious bovine keratoconjunctivitis.

The sequence and gene characterization of a 400-kb candidate region for IDDM4 on chromosome 11q13.

Type 1 diabetes is a complex disorder with interaction of both genetic and environmental factors. One of the loci, IDDM4, has been mapped to chromosome 11q13, with evidence of association to two markers, D11S1917 and H0570polyA. To identify putative candidate genes for IDDM4, we have constructed a 400-kb clone contig in this region and sequenced the clones. We have also sequenced the orthologous DNA from mouse. Previously, we identified a cDNA for the low-density lipoprotein receptor-related protein 5 gene (LRP5) 3 kb distal to H0570polyA. We have now determined the exon-intron structure of this gene. Detailed sequence analysis has identified a further three genes in this region: the CGI-85 gene (previously identified by W.-C. Lin) and two novel genes, C11orf24 and C11orf23. The C11orf24 gene has no known similarity to other genes, and its function is unknown. C11orf23 has similarity to the SIT4 (sporulation-induced transcript 4)-associated protein (SAP) family of yeast proteins, which are involved in regulation of the cell cycle. The full-length C11orf23 cDNA is the first mammalian orthologue of the yeast SAP family to be identified. Identification of these four genes in a 400-kb region of the IDDM4 region underpins our strategy to identify the IDDM4 locus.

Molecular cloning and pharmacological characterization of the canine B1 and B2 bradykinin receptors.

The dog is a valuable animal model in the study of the physiological role of both the B1 and B2 bradykinin receptors. To more thoroughly characterize the pharmacological properties of the canine kinin receptors we isolated the cDNA sequence encoding the B1 and B2 bradykinin receptor subtypes and overexpressed them in Chinese hamster ovary (CHO) cells. The cDNA sequence of the canine B1 bradykinin receptor encodes a protein comprised of 350 amino acids that is 76% identical to the human B1 bradykinin receptor. The cDNA sequence of the canine B2 bradykinin receptor encodes a protein of 392 amino acids that is 81% identical to the human B2 bradykinin receptor. The amino acid sequence of the canine B1 and B2 receptors are 35% identical. Pharmacological studies of the cloned receptors revealed that the agonist affinity of the dog B1 receptor is similar to the rodent B1 receptors, and differs from the human form in that there is no preference for the presence of the N-terminal Lys residue of [des-Arg10]Lys-bradykinin. Significantly, the B1 receptor antagonist [des-Arg9,Leu8]BK behaves as partial agonist on the cloned dog B1 receptor. The dog B2 receptor exhibits the 'classical' pharmacological properties of this receptor subtype.

Development- and differentiation-dependent reorganization of intermediate filaments in fiber cells.

PURPOSE: To define the remodeling of lens fiber cell intermediate filaments (IF) that occurs with both development and differentiation. METHODS: Prenatal and postnatal mice were probed for the IF proteins phakosin, filensin, and vimentin, using light microscope immunocytochemical methodology. RESULTS: The pattern of vimentin accumulation in elongating fiber cells changed with development. Early in development vimentin first emerged predominantly as focal accumulations in the basal region of both epithelial and primary fiber cells. A light diffuse cytoplasmic staining was also noted. Later in embryonic development, and through maturity, vimentin in fiber cells was predominantly associated with the plasma membrane with no anterior-posterior polarity. Phakosin and filensin were first detected in the very latest stages of primary fiber elongation and continued to accumulate well after cells had completed elongation. Initially, these proteins accumulated in the anterior half of the fiber cells and were cytoplasmic in distribution. After P13, the pattern of initial distribution in differentiating fiber cells changed to a predominantly plasma membrane localization. Neither beaded filament protein showed focal basal accumulations. In mature lenses, all three proteins ultimately disappeared from the nuclear fiber cells. CONCLUSIONS: Beaded filament protein accumulation lags significantly behind both primary and secondary fiber cell elongation, suggesting a functional role subsequent to elongation. The subcellular distribution of vimentin and the beaded filament proteins showed marked differences within the cell, with differentiation, and with development. The differences in time of initial synthesis and in distribution of these IF proteins may bear on hypotheses about the role of IFs in fiber cell elongation and in structural-functional polarity of the fiber cell.

Expression of the type I diabetes-associated gene LRP5 in macrophages, vitamin A system cells, and the Islets of Langerhans suggests multiple potential roles in diabetes.

LRP5 is a novel member of the low-density lipoprotein receptor family that is genetically associated with Type 1 diabetes. As a start to defining the normal function of LRP5 and to generate testable hypotheses of its potential role in Type 1 diabetes pathogenesis, we carried out an extensive expression analysis of this gene at the mRNA and protein levels in normal human, monkey, and mouse, as well as in non-obese diabetic (NOD) mice at several stages of diabetes development. In all species, expression of LRP5 was found in four functionally important cell types: the distributed mononuclear phagocyte system, the islets of Langerhans, vitamin A-metabolizing cells, and CNS neurons. Given the critical role of macrophages in the onset and progression of islet cell destruction in Type 1 diabetes and the hypothesized role of retinoids as modifiers of diabetes progression, these findings suggest that LRP5 may confer Type 1 diabetes risk by altering the normal functioning of one or more of these regulatory systems. Specifically, given that the LRP5 polymorphisms associated with diabetes are in the promoter region of the gene, alterations in LRP5 expression may be responsible for diabetes susceptibility and therefore may be potential targets for therapeutic intervention. (J Histochem Cytochem 48:1357-1368, 2000)

The NOD Idd9 genetic interval influences the pathogenicity of insulitis and contains molecular variants of Cd30, Tnfr2, and Cd137.

Previous analyses of NOD mice have shown that some genes control the development of both insulitis and diabetes, while other loci influence diabetes without reducing insulitis. Evidence for the existence of a gene only influencing diabetes, Idd9 on mouse chromosome 4, is provided here by the development of a novel congenic mouse strain, NOD.B10 Idd9. NOD.B10 Idd9 mice display profound resistance to diabetes even though nearly all develop insulitis. Subcongenic analysis has demonstrated that alleles of at least three B10 genes, Idd9.1, Idd9.2, and Idd9.3 are required to produce Idd9-mediated diabetes resistance. Candidate genes with amino acid differences between the NOD and B10 strains have been localized to the 5.6 cM Idd9.2 interval (Tnfr2, Cd30) and to the 2.0 cM Idd9.3 interval (Cd137).

A juvenile-onset, progressive cataract locus on chromosome 3q21-q22 is associated with a missense mutation in the beaded filament structural protein-2.

Juvenile-onset cataracts are distinguished from congenital cataracts by the initial clarity of the lens at birth and the gradual development of lens opacity in the second and third decades of life. Genomewide linkage analysis in a multigenerational pedigree, segregating for autosomal dominant juvenile-onset cataracts, identified a locus in chromosome region 3q21.2-q22.3. Because of the proximity of the gene coding for lens beaded filament structural protein-2 (BFSP2) to this locus, we screened for mutations in the coding sequence of BFSP2. We observed a unique C-->T transition, one that was not observed in 200 normal chromosomes. We predicted that this led to a nonconservative R287W substitution in exon 4 that cosegregated with cataracts. This mutation alters an evolutionarily conserved arginine residue in the central rod domain of the intermediate filament. On consideration of the proposed function of BFSP2 in the lens cytoskeleton, it is likely that this alteration is the cause of cataracts in the members of the family we studied. This is the first example of a mutation in a noncrystallin structural gene that leads to a juvenile-onset, progressive cataract.

Autosomal-dominant congenital cataract associated with a deletion mutation in the human beaded filament protein gene BFSP2.

Congenital cataracts are a common major abnormality of the eye that frequently cause blindness in infants. At least one-third of all cases are familial; autosomal-dominant congenital cataract appears to be the most-common familial form in the Western world. Elsewhere, in family ADCC-3, we mapped an autosomal-dominant cataract gene to chromosome 3q21-q22, near the gene that encodes a lens-specific beaded filament protein gene, BFSP2. By sequencing the coding regions of BFSP2, we found that a deletion mutation, DeltaE233, is associated with cataracts in this family. This is the first report of an inherited cataract that is caused by a mutation in a cytoskeletal protein.

Cloning of a novel member of the low-density lipoprotein receptor family.

A gene encoding a novel transmembrane protein was identified by DNA sequence analysis within the insulin-dependent diabetes mellitus (IDDM) locus IDDM4 on chromosome 11q13. Based on its chromosomal position, this gene is a candidate for conferring susceptibility to diabetes. The gene, termed low-density lipoprotein receptor related protein 5 (LRP5), encodes a protein of 1615 amino acids that contains conserved modules which are characteristic of the low-density lipoprotein (LDL) receptor family. These modules include a putative signal peptide for protein export, four epidermal growth factor (EGF) repeats with associated spacer domains, three LDL-receptor (LDLR) repeats, a single transmembrane spanning domain, and a cytoplasmic domain. The encoded protein has a unique organization of EGF and LDLR repeats; therefore, LRP5 likely represents a new category of the LDLR family. Both human and mouse LRP5 cDNAs have been isolated and the encoded mature proteins are 95% identical, indicating a high degree of evolutionary conservation.

Isolation and characterization of LRP6, a novel member of the low density lipoprotein receptor gene family.

A novel member of the low density lipoprotein receptor (LDLR) gene family has been identified and characterized. This gene, termed LDL receptor-related protein 6 (LRP6), encodes a transmembrane protein which has 71% identity and is structurally similar to the protein encoded by LRP5, a proposed candidate gene for type 1 diabetes located on human chromosome 11q13. LRP6 maps to human chromosome 12p11-p13. Mouse Lrp6 encodes a protein that has 98% identity to human LRP6 and maps to chromosome 6. Unlike other members of the LDLR family, LRP6 and LRP5 display a unique pattern of four epidermal growth factor (EGF) and three LDLR repeats in the extracellular domain. The cytoplasmic domain of LRP6 is not similar to other members of the LDLR family, while comparison with LRP5 reveals proline-rich motifs that may mediate protein-protein interactions. Thus, it is likely that LRP6 and LRP5 comprise a new class of the LDLR family.

Primary sequence, secondary structure, gene structure, and assembly properties suggests that the lens-specific cytoskeletal protein filensin represents a novel class of intermediate filament protein.

The ocular lens fiber cell assembles a novel cytoskeletal element, the Beaded Filament, from CP49 and filensin, two proteins expressed only in the differentiated lens fiber cell. We report the primary sequence, secondary structural analysis, gene structure and Yeast Two Hybrid interaction data for human filensin, and develop a consensus model of filensin from the human and previously reported bovine and chicken filensin sequences. This consensus model, combined with gene structure and Yeast Two Hybrid studies establish that filensin is a member of the Intermediate Filament family of proteins. Specifically, filensin exhibits (1) divergence at amino acid sequence motifs otherwise highly conserved among intermediate filament proteins, (2) a loss of 29 amino acids from the central rod domain which is unique among cytoplasmic intermediate filament proteins, (3) an absence of sequence identity with any existing class of intermediate filament protein, (4) a gene structure unique among intermediate filament family, (5) an inability to dimerize with representatives of Type I, II, and III intermediate filament proteins. Thus, at each level of analysis, we find that filensin is similar to the consensus model of intermediate filament proteins, supporting our conclusion that filensin's relatedness to the IF family is not the consequence of convergent evolution. However, filensin also shows unique or extreme distinctions from the consensus intermediate filament protein at each level of analysis, indicating that filensin constitutes a novel class of IF protein. Some of filensin's unique features are incompatible with current models of IF assembly. Analysis of filensin gene structure suggests that the 29 amino acid reduction in the central rod domain was not the result of a single splice site mutation, the mechanism suggested for the transition between nuclear lamins and cytoplasmic intermediate filament proteins.

Effects of kinins on isolated stomachs of control and transgenic knockout B2 receptor mice.

The aim of this study was to investigate the pharmacological profile of the kinin B1 and B2 receptors in isolated stomachs from wild-type control and B2 receptor knockout mice. Isometric contractions evoked by bradykinin (BK) (9 nM) and desArg9BK (28 nM) were shown to be different. The contraction induced by desArg9BK had a longer duration than that evoked by BK and increased during incubation in vitro in stomachs of wild-type controls, while in the transgenic B2 receptor knockout mice, the contractions evoked by desArg9BK and BK were similar and followed the B1 receptor agonist pattern. BK but not the carboxypeptidase-resistant analog, [Phe8psi(CH2-NH)Arg9]BK, was found to be active in the stomach of B2 receptor knockout mice. BK-induced contractions were prevented by mergetpa (a carboxypeptidase M inhibitor) (10 microM) and by a the B receptor antagonist, AcLys[DbetaNal7,Ile8]desArg9BK (R 715) (0.88 microM), while not being influenced by the B2 receptor antagonist HOE 140 (0.38 microM). BK and [Phe8psi(CH2-NH)Arg9]BK were potent contractants of the wild-type mice stomach and their effects were not influenced by mergetpa or by the B receptor antagonist: they were reduced by HOE 140. After incubation in vitro for 3-4 hours, the tissues were treated with HOE 140 (4 microM) and FR-173657 (17 microM) to eliminate B2 receptor function. In these tissues, BK evoked a B1-like contraction which was inhibited by mergetpa (10 microM) and antagonized by R 715 (8 microM). The results indicate that BK acts primarily on B2 receptors. However, after intramural conversion to desArg9BK, activation of B1 receptors of the mice stomach occurs. In the tissues of B2 receptor knockout mice, BK behaves as a pure B1 receptor agonist while in stomachs of control animals, the B2 receptor contribution is overwhelming. After complete blockade of the B2 receptor, BK is able to evoke B1-mediated responses similar to those observed in tissues of B2 receptor knockout mice. It is concluded that the disruption of the B2 receptor gene eliminates the B2 receptor without influencing the B1 receptor system.

Protection of a restriction enzyme from heat inactivation by [alpha]-crystallin.

PURPOSE: To determine whether the chaperone activity of human alpha-crystallin can protect a restriction enzyme from heat inactivation. METHODS: The restriction enzyme Nde I was heated in the presence or absence of purified bovine alpha-crystallin. Following heat treatment, the enzymatic activity of the heat treated samples was assayed by cleavage of plasmid DNA. The extent of digestion was monitored by agarose gel electrophoresis and visualization of DNA fragments by ethidium bromide staining. RESULTS: Heating of Nde I in the absence of alpha-crystallin resulted in inactivation. However, Nde I heated in the presence of alpha-crystallin remained active. Furthermore, an increased amount of alpha-crystallin provided a longer period of thermal protection. CONCLUSIONS: The chaperone activity and thermo-protective effect of alpha-crystallin extend to protection of enzymatic activity, not merely the protection from thermally induced aggregation/denaturation. In addition, inclusion of alpha-crystallin during some enzymatic reactions may be beneficial.