High-throughput variation detection and genotyping using microarrays

Large-scale identification of human bone remains via SNP microarray analysis with reference SNP database

Single nucleotide polymorphisms (SNPs) are valuable markers complementary to conventional forensic short tandem repeat (STR) markers in genetic typing, with potential advantages in challenging forensic casework. With the advent of high-throughput technologies, such as microarrays and massively parallel sequencing, the use of SNP typing has now expanded to large-scale forensic applications. Herein, a forensic case is presented to demonstrate the usefulness of SNP typing in identifying large-scale human bone remains with reference database construction. A total of 402 bone remains were recovered from an island in the Jeju Province of Korea where a massive disaster occurred in 1948. The first phase of the identification process was accomplished via conventional DNA typing methods including autosomal and Y-chromosomal STR typing, and mitochondrial DNA sequencing, which resulted in the identification of 74 of 402 remains. The second phase of the identification involved the remaining 327 unidentified remains using SNP typing as a supplementary tool based on Affymetrix resequencing array. The SNP markers of 782 family members were also analyzed and a reference database was constructed for comparison. An additional 51 bone remains were identified in the second phase. SNP data obtained from the supplementary genotyping yielded additional genetic information as well as contributed to kinship testing to determine the second degrees of relationship. In addition SNPs are useful in discriminating ambiguous relationship when only STR data are available. A software program developed for SNP typing system enabled efficient kinship analysis for large-scale forensic identification. The results and the casework are described and discussed.

Investigation of Yersinia pestis and Yersinia pseudotuberculosis strains from Georgia and neighboring countries in the Caucasus by high-density SNP microarray

Yersinia pestis, the causative agent of plague, is a recently evolved clone of the enteropathogenic bacterium Yersinia pseudotuberculosis. Y. pestis has been extensively studied for decades; however, there are insufficient data about the intra-species diversity of this microorganism in certain parts of the world, including the Caucasus region. Using a high-density single-nucleotide polymorphism (SNP) microarray, we genotyped a total of 46 Y. pestis isolates from two plague foci in Georgia and neighboring Caucasus countries together with 12 Y. pseudotuberculosis isolates from Georgia. The genotyping microarray captured a total of 13,525 SNP positions across the Y. pestis and Y. pseudotuberculosis genomes and plasmids with high-throughput capability and superior reproducibility. From this analysis, we confirmed the presence of two independent and relatively distant phylogenetic groups of Y. pestis in the Caucasus region. The signature SNP patterns identified from this study will allow assay development for plague surveillance and pseudotuberculosis diagnostics.

Understanding the Bioinformatics Challenges of Integrating Genomics into Healthcare

Genomic data is paving the way towards personalized healthcare. By unveiling genetic disease-contributing factors, genomic data can aid in the detection, diagnosis, and treatment of a wide range of complex diseases. Integrating genomic data into healthcare is riddled with a wide range of challenges spanning social, ethical, legal, educational, economic, and technical aspects. Bioinformatics is a core integration aspect presenting an overwhelming number of unaddressed challenges. In this paper we tackle the fundamental bioinformatics integration concerns including: genomic data generation, storage, representation, and utilization in conjunction with clinical data. We divide the bioinformatics challenges into a series of seven intertwined integration aspects spanning the areas of informatics, knowledge management, and communication. For each aspect, we provide a detailed discussion of the current research directions, outstanding challenges, and possible resolutions. This paper seeks to help narrow the gap between the genomic applications, which are being predominantly utilized in research settings, and the clinical adoption of these applications.

Surface-Mediated DNA Hybridization: Effects of DNA Conformation, Surface Chemistry, and Electrostatics

Single-molecule Förster Resonance Energy Transfer (FRET) was used to study the dynamic association of mobile donor-labeled ssDNA oligonucleotides ("target") with covalently immobilized complementary acceptor-labeled ssDNA oligonucleotides ("probe"). While probe-target association events were resolved for all experiments, such FRET events were far more likely to occur in systems with complementarity and on hydrophobic, as compared to hydrophilic, surfaces. The distribution of donor-acceptor association-time intervals did not exhibit simple first-order kinetics, and when decomposed into a superposition of first-order processes, only a small fraction of events corresponded to a long-lived state that was presumed to represent true DNA hybridization, while the majority of association events were transient, representing nonspecific associations or partial hybridization. The structure of the DNA target and probe affected both the stability of the hybridized state, as well as the likelihood that an association between the two led to hybridization. In particular, the likelihood of hybridization decreased for longer target strands and for targets with stem-loop secondary structure. The presence of oligonucleotide secondary structure reduced the stability of hybridization, while greater complementarity increased stability of the hybridized state. Interestingly, increased ionic strength (i.e., greater electrostatic screening) increased the probability of hybridization but did not influence the lifetime of the hybridized state. Combined, these observations provide a nuanced view of surface-mediated DNA hybridization, where various factors independently influence the probability and stability of hybridization.

PEMapper and PECaller provide a simplified approach to whole-genome sequencing

The analysis of human whole-genome sequencing data presents significant computational challenges. The sheer size of datasets places an enormous burden on computational, disk array, and network resources. Here, we present an integrated computational package, PEMapper/PECaller, that was designed specifically to minimize the burden on networks and disk arrays, create output files that are minimal in size, and run in a highly computationally efficient way, with the single goal of enabling whole-genome sequencing at scale. In addition to improved computational efficiency, we implement a statistical framework that allows for a base by base error model, allowing this package to perform as well or better than the widely used Genome Analysis Toolkit (GATK) in all key measures of performance on human whole-genome sequences.

Kinship Testing Based on SNPs Using Microarray System

BACKGROUND

Kinship testing using biallelic SNP markers has been demonstrated to be a promising approach as a supplement to standard STR typing, and several systems, such as pyrosequencing and microarray, have been introduced and utilized in real forensic cases. The Affymetrix microarray containing 169 autosomal SNPs developed for forensic application was applied to our practical case for kinship analysis that had remained inconclusive due to partial STR profiles of degraded DNA and possibility of inbreeding within the population.

CASE REPORT

169 autosomal SNPs were typed on array with severely degraded DNA of two bone samples, and the kinship compared to genotypes in a reference database of their putative family members.

RESULTS

Two bone samples remained unidentified through traditional STR typing with partial profiles of 10 or 14 of 16 alleles. Because these samples originated from a geographically isolated population, a cautious approach was required when analyzing and declaring true paternity only based on PI values. In a supplementary SNP typing, 106 and 78 SNPs were obtained, and the match candidates were found in each case with improved PI values than using only STRs and with no discrepant SNPs in comparison.

CONCLUSION

Our case showed that the utility of multiple SNPs on array is expected in practical forensic caseworks with an establishment of reference database.

Advantages of Array-Based Technologies for Pre-Emptive Pharmacogenomics Testing

As recognised by the National Institutes of Health (NIH) Precision Medicine Initiative (PMI), microarray technology currently provides a rapid, inexpensive means of identifying large numbers of known genomic variants or gene transcripts in experimental and clinical settings. However new generation sequencing techniques are now being introduced in many clinical genetic contexts, particularly where novel mutations are involved. While these methods can be valuable for screening a restricted set of genes for known or novel mutations, implementation of whole genome sequencing in clinical practice continues to present challenges. Even very accurate high-throughput methods with small error rates can generate large numbers of false negative or false positive errors due to the high numbers of simultaneous readings. Additional validation is likely to be required for safe use of any such methods in clinical settings. Custom-designed arrays can offer advantages for screening for common, known mutations and, in this context, may currently be better suited for accredited, quality-controlled clinical genetic screening services, as illustrated by their successful application in several large-scale pre-emptive pharmacogenomics programs now underway. Excessive, inappropriate use of next-generation sequencing may waste scarce research funds and other resources. Microarrays presently remain the technology of choice in applications that require fast, cost-effective genome-wide screening of variants of known importance, particularly for large sample sizes. This commentary considers some of the applications where microarrays continue to offer advantages over next-generation sequencing technologies.

Mango (Mangifera indica L.) germplasm diversity based on single nucleotide polymorphisms derived from the transcriptome

BACKGROUND

Germplasm collections are an important source for plant breeding, especially in fruit trees which have a long duration of juvenile period. Thus, efforts have been made to study the diversity of fruit tree collections. Even though mango is an economically important crop, most of the studies on diversity in mango collections have been conducted with a small number of genetic markers.

RESULTS

We describe a de novo transcriptome assembly from mango cultivar 'Keitt'. Variation discovery was performed using Illumina resequencing of 'Keitt' and 'Tommy Atkins' cultivars identified 332,016 single-nucleotide polymorphisms (SNPs) and 1903 simple-sequence repeats (SSRs). Most of the SSRs (70.1%) were of trinucleotide with the preponderance of motif (GGA/AAG)n and only 23.5% were di-nucleotide SSRs with the mostly of (AT/AT)n motif. Further investigation of the diversity in the Israeli mango collection was performed based on a subset of 293 SNPs. Those markers have divided the Israeli mango collection into two major groups: one group included mostly mango accessions from Southeast Asia (Malaysia, Thailand, Indonesia) and India and the other with mainly of Floridian and Israeli mango cultivars. The latter group was more polymorphic (FS=-0.1 on the average) and was more of an admixture than the former group. A slight population differentiation was detected (FST=0.03), suggesting that if the mango accessions of the western world apparently was originated from Southeast Asia, as has been previously suggested, the duration of cultivation was not long enough to develop a distinct genetic background.

CONCLUSIONS

Whole-transcriptome reconstruction was used to significantly broaden the mango's genetic variation resources, i.e., SNPs and SSRs. The set of SNP markers described in this study is novel. A subset of SNPs was sampled to explore the Israeli mango collection and most of them were polymorphic in many mango accessions. Therefore, we believe that these SNPs will be valuable as they recapitulate and strengthen the history of mango diversity.

Single-Nucleotide Polymorphisms Reveal Spatial Diversity Among Clones of Yersinia pestis During Plague Outbreaks in Colorado and the Western United States

BACKGROUND

In western North America, plague epizootics caused by Yersinia pestis appear to sweep across landscapes, primarily infecting and killing rodents, especially ground squirrels and prairie dogs. During these epizootics, the risk of Y. pestis transmission to humans is highest. While empirical models that include climatic conditions and densities of rodent hosts and fleas can predict when epizootics are triggered, bacterial transmission patterns across landscapes, and the scale at which Y. pestis is maintained in nature during inter-epizootic periods, are poorly defined. Elucidating the spatial extent of Y. pestis clones during epizootics can determine whether bacteria are propagated across landscapes or arise independently from local inter-epizootic maintenance reservoirs.

MATERIAL AND METHODS

We used DNA microarray technology to identify single-nucleotide polymorphisms (SNPs) in 34 Y. pestis isolates collected in the western United States from 1980 to 2006, 21 of which were collected during plague epizootics in Colorado. Phylogenetic comparisons were used to elucidate the hypothesized spread of Y. pestis between the mountainous Front Range and the eastern plains of northern Colorado during epizootics. Isolates collected from across the western United States were included for regional comparisons.

RESULTS

By identifying SNPs that mark individual clones, our results strongly suggest that Y. pestis is maintained locally and that widespread epizootic activity is caused by multiple clones arising independently at small geographic scales. This is in contrast to propagation of individual clones being transported widely across landscapes. Regionally, our data are consistent with the notion that Y. pestis diversifies at relatively local scales following long-range translocation events. We recommend that surveillance and prediction by public health and wildlife management professionals focus more on models of local or regional weather patterns and ecological factors that may increase risk of widespread epizootics, rather than predicting or attempting to explain epizootics on the basis of movement of host species that may transport plague.

Development of a high-throughput resequencing array for the detection of pathogenic mutations in osteogenesis imperfecta

Objective

Osteogenesis imperfecta (OI) is a rare inherited skeletal disease, characterized by bone fragility and low bone density. The mutations in this disorder have been widely reported to be on various exonal hotspots of the candidate genes, including COL1A1, COL1A2, CRTAP, LEPRE1, and FKBP10, thus creating a great demand for precise genetic tests. However, large genome sizes make the process daunting and the analyses, inefficient and expensive. Therefore, we aimed at developing a fast, accurate, efficient, and cheaper sequencing platform for OI diagnosis; and to this end, use of an advanced array-based technique was proposed.

Method

A CustomSeq Affymetrix Resequencing Array was established for high-throughput sequencing of five genes simultaneously. Genomic DNA extraction from 13 OI patients and 85 normal controls and amplification using long-range PCR (LR-PCR) were followed by DNA fragmentation and chip hybridization, according to standard Affymetrix protocols. Hybridization signals were determined using GeneChip Sequence Analysis Software (GSEQ). To examine the feasibility, the outcome from new resequencing approach was validated by conventional capillary sequencing method.

Result

Overall call rates using resequencing array was 96–98% and the agreement between microarray and capillary sequencing was 99.99%. 11 out of 13 OI patients with pathogenic mutations were successfully detected by the chip analysis without adjustment, and one mutation could also be identified using manual visual inspection.

Conclusion

A high-throughput resequencing array was developed that detects the disease-associated mutations in OI, providing a potential tool to facilitate large-scale genetic screening for OI patients. Through this method, a novel mutation was also found.

Resequencing microarray technology for genotyping human papillomavirus in cervical smears

There are more than 40 human papillomaviruses (HPVs) belonging to the alpha genus that cause sexually transmitted infections; these infections are among the most frequent and can lead to condylomas and anogenital intra-epithelial neoplasia. At least 18 of these viruses are causative agents of anogenital carcinomas. We evaluated the performance of a resequencing microarray for the detection and genotyping of alpha HPV of clinical significance using cloned HPV DNA. To reduce the number of HPV genotypes tiled on microarray, we used reconstructed ancestral sequences (RASs) as they are more closely related to the various genotypes than the current genotypes are among themselves. The performance of this approach was tested by genotyping with a set of 40 cervical smears already genotyped using the commercial PapilloCheck kit. The results of the two tests were concordant for 70% (28/40) of the samples and compatible for 30% (12/40). Our findings indicate that RASs were able to detect and identify one or several HPV in clinical samples. Associating RASs with homonym sequences improved the genotyping of HPV present in cases of multiple infection. In conclusion, we demonstrate the diagnostic potential of resequencing technology for genotyping of HPV, and illustrate its value both for epidemiological studies and for monitoring the distribution of HPV in the post-vaccination era.

Combination of an unbiased amplification method and a resequencing microarray for detecting and genotyping equine arteritis virus

This study shows that an unbiased amplification method applied to equine arteritis virus RNA significantly improves the sensitivity of the real-time reverse transcription-quantitative PCR (RT-qPCR) recommended by the World Organization for Animal Health. Twelve viral RNAs amplified using this method were hybridized on a high-density resequencing microarray for effective viral characterization.

Mitochondriome and cholangiocellular carcinoma

Cholangiocellular carcinoma (CCA) of the liver was the target of more interest, recently, due mainly to its increased incidence and possible association to new environmental factors. Somatic mitochondrial DNA (mtDNA) mutations have been found in several cancers. Some of these malignancies contain changes of mtDNA, which are not or, very rarely, found in the mtDNA databases. In terms of evolutionary genetics and oncology, these data are extremely interesting and may be considered a sign of poor fitness, which may conduct in some way to different cellular processes, including carcinogenesis. MitoChip analysis is a strong tool for investigations in experimental oncology and was carried out on three CCA cell lines (HuCCT1, Huh-28 and OZ) with different outcome in human and a Papova-immortalized normal hepatocyte cell line (THLE-3). Real time quantitative PCR, western blot analysis, transmission electron microscopy, confocal laser microscopy, and metabolic assays including L-Lactate and NAD⁺/NADH assays were meticulously used to identify mtDNA copy number, oxidative phosphorylation (OXPHOS) content, ultrastructural morphology, mitochondrial membrane potential (ΔΨ_m), and differential composition of metabolites, respectively. Among 102 mtDNA changes observed in the CCA cell lines, 28 were non-synonymous coding region alterations resulting in an amino acid change. Thirty-eight were synonymous and 30 involved ribosomal RNA (rRNA) and transfer RNA (tRNA) regions. We found three new heteroplasmic mutations in two CCA cell lines (HuCCT1 and Huh-28). Interestingly, mtDNA copy number was decreased in all three CCA cell lines, while complexes I and III were decreased with depolarization of mitochondria. L-Lactate and NAD+/NADH assays were increased in all three CCA cell lines. MtDNA alterations seem to be a common event in CCA. This is the first study using MitoChip analysis with comprehensive metabolic studies in CCA cell lines potentially creating a platform for future studies on the interactions between normal and neoplastic cells.

Mechanisms of surface-mediated DNA hybridization

Single-molecule total internal reflection fluorescence microscopy was employed in conjunction with resonance energy transfer (RET) to observe the dynamic behavior of donor-labeled ssDNA at the interface between aqueous solution and a solid surface decorated with complementary acceptor-labeled ssDNA. At least 100,000 molecular trajectories were determined for both complementary strands and negative control ssDNA. RET was used to identify trajectory segments corresponding to the hybridized state. The vast majority of molecules from solution adsorbed nonspecifically to the surface, where a brief two-dimensional search was performed with a 7% chance of hybridization. Successful hybridization events occurred with a characteristic search time of ∼0.1 s, and unsuccessful searches resulted in desorption from the surface, ultimately repeating the adsorption and search process. Hybridization was reversible, and two distinct modes of melting (i.e., dehybridization) were observed, corresponding to long-lived (∼15 s) and short-lived (∼1.4 s) hybridized time intervals. A strand that melted back onto the surface could rehybridize after a brief search or desorb from the interface. These mechanistic observations provide guidance for technologies that involve DNA interactions in the near-surface region, suggesting a need to design surfaces that both enhance the complex multidimensional search process and stabilize the hybridized state.

Characterization of sal(A), a novel gene responsible for lincosamide and streptogramin A resistance in Staphylococcus sciuri

Natural resistance to lincosamides and streptogramins A (LSA), which is a species characteristic of Bacillus subtilis and Enterococcus faecalis, has never been documented in the Staphylococcus genus. We investigate here the molecular basis of the LSA phenotype exhibited by seven reference strains of Staphylococcus sciuri, including the type strains of the three described subspecies. By whole-genome sequencing of strain ATCC 29059, we identified a candidate gene that encodes an ATP-binding cassette protein similar to the Lsa and VmlR resistance determinants. Isolation and reverse transcription-quantitative PCR (qRT-PCR) expression studies confirmed that Sal(A) can confer a moderate resistance to lincosamides (8 times the MIC of lincomycin) and a high-level resistance to streptogramins A (64 times the MIC of pristinamycin II). The chromosomal location of sal(A) between two housekeeping genes of the staphylococcal core genome supports the gene's ancient origins and thus innate resistance to these antimicrobials within S. sciuri subspecies.

Molecular epidemiology and clinical spectrum of hereditary spastic paraplegia in the Japanese population based on comprehensive mutational analyses

Hereditary spastic paraplegia (HSP) is one of the most genetically heterogeneous neurodegenerative disorders characterized by progressive spasticity and pyramidal weakness of lower limbs. Because >30 causative genes have been identified, screening of multiple genes is required for establishing molecular diagnosis of individual patients with HSP. To elucidate molecular epidemiology of HSP in the Japanese population, we have conducted mutational analyses of 16 causative genes of HSP (L1CAM, PLP1, ATL1, SPAST, CYP7B1, NIPA1, SPG7, KIAA0196, KIF5A, HSPD1, BSCL2, SPG11, SPG20, SPG21, REEP1 and ZFYVE27) using resequencing microarrays, array-based comparative genomic hybridization and Sanger sequencing. The mutational analysis of 129 Japanese patients revealed 49 mutations in 46 patients, 32 of which were novel. Molecular diagnosis was accomplished for 67.3% (33/49) of autosomal dominant HSP patients. Even among sporadic HSP patients, mutations were identified in 11.1% (7/63) of them. The present study elucidated the molecular epidemiology of HSP in the Japanese population and further broadened the mutational and clinical spectra of HSP.

Identification of novel sequence variations in microRNAs in chronic lymphocytic leukemia

We have analyzed the miRNA sequence variations in patients with CLL and the effect of these variations on their secondary structure and expression.

MicroRNA (miRNA) expression is deregulated in many tumors including chronic lymphocytic leukemia (CLL). Although the particular mechanism(s) responsible for their aberrant expression is not well characterized, the presence of mutations and single-nucleotide polymorphisms (SNPs) in miRNA genes, possibly affecting their secondary structure and expression, has been described. In CLL; however, the impact and frequency of such variations have yet to be elucidated. Using a custom resequencing microarray, we screened sequence variations in 109 cancer-related pre-miRNAs in 98 CLL patients. Additionally, the primary regions of miR-29b-2/29c and miR-16-1 were analyzed by Sanger sequencing in another cohort of 213 and 193 CLL patients, respectively. Altogether, we describe six novel miR-sequence variations and the presence of SNPs (n = 27), most of which changed the miR-secondary structure. Moreover, some of the identified SNPs have a significantly different frequency in CLL when compared with a control population. Additionally, we identified a novel variation in miR-16-1 that had not been described previously in CLL patients. We show that this variation affects the expression of mature miR-16-1. We also show that the expression of another miRNA with pathogenetic relevance for CLL, namely miR-29b-2, is influenced by the presence of a polymorphic insertion, which is more frequent in CLL than in a control population. Altogether, these data suggest that sequence variations may occur during CLL development and/or progression.

Association study of mitochondrial genetic polymorphisms in asthmatic children

It has been suggested that mitochondrial dysfunction plays a role in the pathogenesis of asthma. To test whether mitochondrial variants influence the risk of asthma, we analyzed 16,158 mtSNPs in a sample of 372 asthmatic children and 395 healthy children using the DNA pooling technique and genome wide association analysis. Stratified analysis by sex was performed to explain the differences observed between sexes in the etiology of asthma. Different variants were detected to be significant in the sample of girls and boys with the smallest adjusted p values being 1.4 × 10(-09) (mt5295) and 3.6 × 10(-12) (mt16158), respectively. Most of the significant locations found in boys are within the CYB gene and the non-coding region. For girls, most of the significant mtSNPs lie within NADH-dehydrogenase-subunits. The variants reported here have not previously been described in connection with asthma. Although further studies in other cohorts are needed to confirm these findings our study highlights the importance of the mitochondria among the factors that contribute to the risk of asthma.

Cell-material interactions revealed via material techniques of surface patterning

Cell-material interactions constitute a key fundamental topic in biomaterials study. Various cell cues and matrix cues as well as soluble factors regulate cell behaviors on materials. These factors are coupled with each other as usual, and thus it is very difficult to unambiguously elucidate the role of each regulator. The recently developed material techniques of surface patterning afford unique ways to reveal the underlying science. This paper reviews the pertinent material techniques to fabricate patterns of microscale and nanoscale resolutions, and corresponding cell studies. Some issues are emphasized, such as cell localization on patterned surfaces of chemical contrast, and effects of cell shape, cell size, cell-cell contact, and seeding density on differentiation of stem cells. Material cues to regulate cell adhesion, cell differentiation and other cell events are further summed up. Effects of some physical properties, such as surface topography and matrix stiffness, on cell behaviors are also discussed; nanoscaled features of substrate surfaces to regulate cell fate are summarized as well. The pertinent work sheds new insight into the cell-material interactions, and is stimulating for biomaterial design in regenerative medicine, tissue engineering, and high-throughput detection, diagnosis, and drug screening.

Genomic profiling of mantle cell lymphoma

Genomic profiling of mantle cell lymphoma (MCL) cells has enabled a better understanding of the complex mechanisms underlying the pathogenesis of disease. Besides the t(11;14)(q13;q32) leading to cyclin D1 overexpression, MCL exhibits a characteristic pattern of DNA copy number aberrations that differs from those detected in other B-cell lymphomas. These genomic changes disrupt selected oncogenes and suppressor genes that are required for lymphoma development and progression, many of which are components of cell cycle, DNA damage response and repair, apoptosis, and cell-signaling pathways. Additionally, some of them may represent effective therapeutic targets. A number of genomic and molecular abnormalities have been correlated with the clinical outcome of patients with MCL and are considered prognostic factors. However, only a few genomic markers have been shown to predict the response to current or novel targeted therapies. One representative example is the high-level amplification of the BCL2 gene, which predicts a good response to pro-apoptotic BH3 mimetic drugs. In summary, genomic analyses have contributed to the substantial advances made in the comprehension of the pathogenesis of MCL, providing a solid basis for the identification of optimal therapeutic targets and for the design of new molecular therapies aiming to cure this fatal disease.

A model of binding on DNA microarrays: understanding the combined effect of probe synthesis failure, cross-hybridization, DNA fragmentation and other experimental details of affymetrix arrays

BACKGROUND

DNA microarrays are used both for research and for diagnostics. In research, Affymetrix arrays are commonly used for genome wide association studies, resequencing, and for gene expression analysis. These arrays provide large amounts of data. This data is analyzed using statistical methods that quite often discard a large portion of the information. Most of the information that is lost comes from probes that systematically fail across chips and from batch effects. The aim of this study was to develop a comprehensive model for hybridization that predicts probe intensities for Affymetrix arrays and that could provide a basis for improved microarray analysis and probe development. The first part of the model calculates probe binding affinities to all the possible targets in the hybridization solution using the Langmuir isotherm. In the second part of the model we integrate details that are specific to each experiment and contribute to the differences between hybridization in solution and on the microarray. These details include fragmentation, wash stringency, temperature, salt concentration, and scanner settings. Furthermore, the model fits probe synthesis efficiency and target concentration parameters directly to the data. All the parameters used in the model have a well-established physical origin.

RESULTS

For the 302 chips that were analyzed the mean correlation between expected and observed probe intensities was 0.701 with a range of 0.88 to 0.55. All available chips were included in the analysis regardless of the data quality. Our results show that batch effects arise from differences in probe synthesis, scanner settings, wash strength, and target fragmentation. We also show that probe synthesis efficiencies for different nucleotides are not uniform.

CONCLUSIONS

To date this is the most complete model for binding on microarrays. This is the first model that includes both probe synthesis efficiency and hybridization kinetics/cross-hybridization. These two factors are sequence dependent and have a large impact on probe intensity. The results presented here provide novel insight into the effect of probe synthesis errors on Affymetrix microarrays; furthermore, the algorithms developed in this work provide useful tools for the analysis of cross-hybridization, probe synthesis efficiency, fragmentation, wash stringency, temperature, and salt concentration on microarray intensities.

Molecular diagnostic and genetic characterization of highly pathogenic viruses: application during Crimean-Congo haemorrhagic fever virus outbreaks in Eastern Europe and the Middle East

Several haemorrhagic fevers are caused by highly pathogenic viruses that must be handled in Biosafety level 4 (BSL–4) containment. These zoonotic infections have an important impact on public health and the development of a rapid and differential diagnosis in case of outbreak in risk areas represents a critical priority. We have demonstrated the potential of a DNA resequencing microarray (PathogenID v2.0) for this purpose. The microarray was first validated in vitro using supernatants of cells infected with prototype strains from five different families of BSL-4 viruses (e.g. families Arenaviridae, Bunyaviridae, Filoviridae, Flaviviridae and Paramyxoviridae). RNA was amplified based on isothermal amplification by Phi29 polymerase before hybridization. We were able to detect and characterize Nipah virus and Crimean–Congo haemorrhagic fever virus (CCHFV) in the brains of experimentally infected animals. CCHFV was finally used as a paradigm for epidemics because of recent outbreaks in Turkey, Kosovo and Iran. Viral variants present in human sera were characterized by BLASTN analysis. Sensitivity was estimated to be 10⁵–10⁶ PFU/mL of hybridized cDNA. Detection specificity was limited to viral sequences having ∼13–14% of global divergence with the tiled sequence, or stretches of ∼20 identical nucleotides. These results highlight the benefits of using the PathogenID v2.0 resequencing microarray to characterize geographical variants in the follow-up of haemorrhagic fever epidemics; to manage patients and protect communities; and in cases of bioterrorism.

Effects of molecular size and surface hydrophobicity on oligonucleotide interfacial dynamics

Single-molecule total internal reflection fluorescence microscopy was used to observe the dynamic behavior of polycytosine single-stranded DNA (ssDNA) (1-50 nucleotides long) at the interface between aqueous solution and hydrophilic (oligoethylene glycol-modified fused silica, OEG) and hydrophobic (octadecyltriethoxysilane-modified fused silica, OTES) solid surfaces. High throughput molecular tracking was used to determine >75,000 molecular trajectories for each molecular length, which were then used to calculate surface residence time and squared displacement (i.e., "step-size") distributions. On hydrophilic OEG surfaces, the surface residence time increased systematically with ssDNA chain length, as expected due to increasing molecule-surface interactions. Interestingly, the residence time decreased with increasing ssDNA length on the hydrophobic OTES surface, particularly for longer chains. Similarly, the interfacial mobility of polynucleotides slowed with increasing chain length on OEG, but became faster on OTES. On OTES surfaces, the rates associated with desorption and surface diffusion exhibited the distinctive anomalous temperature dependence that is characteristic of hydrophobic interactions for short-chain species but not for longer chains. These combined observations suggest that long oligonucleotides adopt conformations minimizing hydrophobic interactions, e.g., by internal sequestration of hydrophobic nucleobases.

Evaluating mitochondrial DNA variation in autism spectrum disorders

Despite the increasing speculation that oxidative stress and abnormal energy metabolism may play a role in Autism Spectrum Disorders (ASD), and the observation that patients with mitochondrial defects have symptoms consistent with ASD, there are no comprehensive published studies examining the role of mitochondrial variation in autism. Therefore, we have sought to comprehensively examine the role of mitochondrial DNA (mtDNA) variation with regard to ASD risk, employing a multi-phase approach. In phase 1 of our experiment, we examined 132 mtDNA single-nucleotide polymorphisms (SNPs) genotyped as part of our genome-wide association studies of ASD. In phase 2 we genotyped the major European mitochondrial haplogroup-defining variants within an expanded set of autism probands and controls. Finally in phase 3, we resequenced the entire mtDNA in a subset of our Caucasian samples (∼400 proband-father pairs). In each phase we tested whether mitochondrial variation showed evidence of association to ASD. Despite a thorough interrogation of mtDNA variation, we found no evidence to suggest a major role for mtDNA variation in ASD susceptibility. Accordingly, while there may be attractive biological hints suggesting the role of mitochondria in ASD our data indicate that mtDNA variation is not a major contributing factor to the development of ASD.

High-density oligonucleotide-based resequencing assay for mutations causing syndromic and non-syndromic forms of thoracic aortic aneurysms and dissections

Thoracic aortic aneurysm and dissection is associated with increasing mortality rate that may occur as part of a syndrome or as an isolated familial condition. Several genes have been implicated in causing TAAD, though an appropriate genetic test for their parallel testing is not yet available. Herein, we describe the novel 117-kb "MFSTAAD chip" that may help to understand the genetic basis of TAAD. A custom duplicate resequencing assay was developed to cover eight genes previously described in TAAD; FBN1, TGFBR1&2, COL3A1, MYH11, ACTA2, SLC2A10 and NOTCH1. GSEQ and SeqC software were used for data analysis. The analytical sensitivity of the assay was validated by the recognition of 182 known mutations (153 point mutations, 21 deletions, 7 insertions and 1 duplication) and a cohort of 28 patients were selected to determine the mutation yield, whereby 18 of them were previously negative for mutations in the genes FBN1 and TGFBR2. The assay had significantly higher sensitivity for point mutations (100%) and the largest deletion of 16 bp was detectable through a decline in the hybridization strength. The overall analytical sensitivity was 85%. Mutation testing of 28 unrelated TAAD patients revealed 4 known and 6 possibly pathogenic mutations with a mutation yield of 32%. The MFSTAAD chip is an alternative tool to next-generation sequencing that allows parallel analysis of several genes on a single platform. Refinements in the probe design and data analysis software will increase the analytical sensitivity of insertions and deletions making this assay even more applicable for clinical testing.

Mitochondrial DNA mutations in pancreatic cancer

BACKGROUND

Somatic mutations of mitochondrial DNA (mtDNA) are increasingly being recognized in many human cancers, but automated sequencing of 16.5 kb of DNA poses an onerous task. We have recently described an oligonucleotide microarray (MitoChip) for rapid and accurate sequencing of the entire mitochondrial genome (Zhou et al., J Mol Diagnostics, 8: 9_14, 2006), greatly facilitating the analysis of mtDNA mutations in cancer. In this report, we perform a comprehensive cataloging of somatic mutations in the mitochondrial genome of human pancreatic cancers using our novel array-based approach.

MATERIALS AND METHODS

MitoChip analysis was performed on DNA isolated from 15 histologically confirmed resection specimens of pancreatic ductal adenocarcinomas. In all cases, matched nonneoplastic pancreatic tissue was obtained as germline control for mtDNA sequence. DNA was extracted from snap-frozen cryostat-embedded specimens and hybridized to the sequencing microarray after appropriate polymerase chain reaction amplification and labeling steps. The vast majority of somatic mutational analyses of mtDNA in human cancers utilize lymphocyte DNA as germline control, without excluding the potential for organ-specific polymorphisms. Therefore, we also examined a series of 15 paired samples of DNA obtained from nonneoplastic pancreata and corresponding EBV-immortalized lymphoblastoid cell lines to determine whether lymphocyte DNA provides an accurate surrogate for the mtDNA sequence of pancreatic tissue.

RESULTS

We sequenced 497,070 base pairs of mtDNA in the 15 matched samples of pancreatic cancer and nonneoplastic pancreatic tissue, and 467,269 base pairs (94.0%) were assigned by the automated genotyping software. All 15 pancreatic cancers demonstrated at least one somatic mtDNA mutation compared to the control germline DNA with a range of 1-14 alterations. Of the 71 somatic mutations observed in our series, 18 were nonsynonymous coding region alterations (i.e., resulting in an amino acid change), 22 were synonymous coding region alterations, and 31 involved noncoding mtDNA segments (including ribosomal and transfer RNAs). Overall, somatic mutations in the coding region most commonly involved the ND4, COI, and CYTB genes; of note, an A-G transition at nucleotide position 841 in the 12sRNA was observed in three independent samples. In the paired analysis of nonneoplastic pancreata and lymphoblastoid cell line DNA, 14 nucleotide discrepancies were observed out of 226,876 nucleotide sequences (a concordance rate of 99.99%), with 9 samples demonstrating a perfect match across all bases assigned.

CONCLUSIONS

Our findings confirm that somatic mtDNA mutations are common in pancreatic cancers, and therefore, have the potential to be a clinically useful biomarker for early detection. Further, our studies confirm that lymphocyte DNA is an excellent, albeit not perfect, surrogate for nonneoplastic pancreatic tissues in terms of being utilized as a germline control. Finally, our report confirms the utility of a high-throughput array-based platform for mtDNA mutational analyses of human cancers.

Application of resequencing microarrays in microbial detection and characterization

Microarrays are powerful, highly parallel assays that are transforming microbiological diagnostics and research. The adaptation of microarray-based resequencing technology for microbial detection and characterization resulted in the development of a number assays that have unique advantages over other existing technologies. This technological platform seems to be especially useful for sensitive and high-resolution multiplexed diagnostics for clinical syndromes with similar symptoms, screening environmental samples for biothreat agents, as well as genotyping and whole-genome analysis of single pathogens.

Alanyl-tRNA synthetase mutation in a family with dominant distal hereditary motor neuropathy

OBJECTIVE

To identify a new genetic cause of distal hereditary motor neuropathy (dHMN), which is also known as a variant of Charcot-Marie-Tooth disease (CMT), in a Chinese family.

METHODS

We investigated a Chinese family with dHMN clinically, electrophysiologically, and genetically. We screened for the mutations of 28 CMT or related pathogenic genes using an originally designed microarray resequencing DNA chip.

RESULTS

Investigation of the family history revealed an autosomal dominant transmission pattern. The clinical features of the family included mild weakness and wasting of the distal muscles of the lower limb and foot deformity, without clinical sensory involvement. Electrophysiologic studies revealed motor neuropathy. MRI of the lower limbs showed accentuated fatty infiltration of the gastrocnemius and vastus lateralis muscles. All 4 affected family members had a heterozygous missense mutation c.2677G>A (p.D893N) of alanyl-tRNA synthetase (AARS), which was not found in the 4 unaffected members and control subjects.

CONCLUSION

An AARS mutation caused dHMN in a Chinese family. AARS mutations result in not only a CMT phenotype but also a dHMN phenotype.

Sequencing by hybridization of long targets

Sequencing by Hybridization (SBH) reconstructs an n-long target DNA sequence from its biochemically determined l-long subsequences. In the standard approach, the length of a uniformly random sequence that can be unambiguously reconstructed is limited to due to repetitive subsequences causing reconstruction degeneracies. We present a modified sequencing method that overcomes this limitation without the need for different types of biochemical assays and is robust to error.

A high-throughput resequencing microarray for autosomal dominant spastic paraplegia genes

Hereditary spastic paraplegias (HSP) are a heterogeneous group of neurological disorders. Insidiously progressive spastic weakness of the lower extremities is the common criterion in all forms described. Clinically, HSP is differentiated into pure (uncomplicated) and complex (complicated) forms. While pure HSP is predominantly characterized by signs and symptoms of pyramidal tract dysfunction, additional neurological and non-neurological symptoms occur in complicated forms. Autosomal dominant, autosomal recessive, and X-linked modes of inheritance have been described and at least 48 subtypes, termed SPG1-48, have been genetically defined. Although in autosomal dominant HSP families 50-60% of etiologies can be established by genetic testing, genotype predictions based on the phenotype are limited. In order to realize high-throughput genotyping for dominant HSP, we designed a resequencing microarray for six autosomal dominant genes on the Affymetrix CustomSEQ array platform. For validation purposes, 10 previously Sanger sequenced patients with autosomal dominant HSP and 40 positive controls with known mutations in ATL1, SPAST, NIPA1, KIF5A, and BSCL2 (32 base exchanges, eight small indels) were resequenced on this array. DNA samples of 45 additional patients with AD spastic paraplegia were included in the study. With two different sequencing analysis software modules (GSEQ, SeqC), all missense/nonsense mutations in the positive controls were identified while indels had a detection rate of only 50%. In total, 244 common synonymous single-nucleotide polymorphisms (SNPs) annotated in dbSNP (build 132) corresponding to 22 distinct sequence variations were found in the 53 analyzed patients. Among the 22 different sequence variations (SPAST n = 15, ATL1 n = 3, KIF5A n = 2, HSPD1 n = 1, BSCL2 n = 1, NIPA1 n = 0), 12 were rare variants that have not been previously described and whose clinical significance is unknown. In SPAST-negative cases, a genetic diagnosis could be established in 11% by resequencing. Resequencing microarray technology can therefore efficiently be used to study genotypes and mutations in large patient cohorts.

DNA microarray-based mutation discovery and genotyping

DNA microarrays provide an efficient means of identifying single-nucleotide polymorphisms (SNPs) in DNA samples and characterizing their frequencies in individual and mixed samples. We have studied the parameters that determine the sensitivity of DNA probes to SNPs and found that the melting temperature (T (m)) of the probe is the primary determinant of probe sensitivity. An isothermal-melting temperature DNA microarray design, in which the T (m) of all probes is tightly distributed, can be implemented by varying the length of DNA probes within a single DNA microarray. I describe guidelines for designing isothermal-melting temperature DNA microarrays and protocols for labeling and hybridizing DNA samples to DNA microarrays for SNP discovery, genotyping, and quantitative determination of allele frequencies in mixed samples.

Pneumococcal genome sequencing tracks a vaccine escape variant formed through a multi-fragment recombination event

Streptococcus pneumoniae ('pneumococcus') causes an estimated 14.5 million cases of serious disease and 826,000 deaths annually in children under 5 years of age(1). The highly effective introduction of the PCV7 pneumococcal vaccine in 2000 in the United States(2,3) provided an unprecedented opportunity to investigate the response of an important pathogen to widespread, vaccine-induced selective pressure. Here, we use array-based sequencing of 62 isolates from a US national monitoring program to study five independent instances of vaccine escape recombination(4), showing the simultaneous transfer of multiple and often large (up to at least 44 kb) DNA fragments. We show that one such new strain quickly became established, spreading from east to west across the United States. These observations clarify the roles of recombination and selection in the population genomics of pneumococcus and provide proof of principle of the considerable value of combining genomic and epidemiological information in the surveillance and enhanced understanding of infectious diseases.

Evaluating Japanese patients with the Marfan syndrome using high-throughput microarray-based mutational analysis of fibrillin-1 gene

Marfan syndrome (MS) is an inherited connective tissue disorder, and detailed evaluations of multiple organ systems are required for its diagnosis. Genetic testing of the disease-causing fibrillin-1 gene (FBN1) is also important in this diagnostic scheme. The aim of this study was to define the clinical characteristics of Japanese patients with MS and enable the efficient and accurate diagnosis of MS with mutational analysis using a high-throughput microarray-based resequencing system. Fifty-three Japanese probands were recruited, and their clinical characteristics were evaluated using the Ghent criteria. For mutational analysis, an oligonucleotide microarray was designed to interrogate FBN1, and the entire exon and exon-intron boundaries of FBN1 were sequenced. Clinical evaluation revealed more pulmonary phenotypes and fewer skeletal phenotypes in Japanese patients with MS compared to Caucasians. The microarray-based resequencing system detected 35 kinds of mutations, including 23 new mutations. The mutation detection rate for patients who fulfilled the Ghent criteria reached 71%. Of note, splicing mutations accounted for 19% of all mutations, which is more than previously reported. In conclusion, this comprehensive approach successfully detected clinical phenotypes of Japanese patients with MS and demonstrated the usefulness and feasibility of this microarray-based high-throughput resequencing system for mutational analysis of MS.

Mutation screening of multiple genes in Spanish patients with autosomal recessive retinitis pigmentosa by targeted resequencing

Retinitis Pigmentosa (RP) is a heterogeneous group of inherited retinal dystrophies characterised ultimately by the loss of photoreceptor cells. RP is the leading cause of visual loss in individuals younger than 60 years, with a prevalence of about 1 in 4000. The molecular genetic diagnosis of autosomal recessive RP (arRP) is challenging due to the large genetic and clinical heterogeneity. Traditional methods for sequencing arRP genes are often laborious and not easily available and a screening technique that enables the rapid detection of the genetic cause would be very helpful in the clinical practice. The goal of this study was to develop and apply microarray-based resequencing technology capable of detecting both known and novel mutations on a single high-throughput platform. Hence, the coding regions and exon/intron boundaries of 16 arRP genes were resequenced using microarrays in 102 Spanish patients with clinical diagnosis of arRP. All the detected variations were confirmed by direct sequencing and potential pathogenicity was assessed by functional predictions and frequency in controls. For validation purposes 4 positive controls for variants consisting of previously identified changes were hybridized on the array. As a result of the screening, we detected 44 variants, of which 15 are very likely pathogenic detected in 14 arRP families (14%). Finally, the design of this array can easily be transformed in an equivalent diagnostic system based on targeted enrichment followed by next generation sequencing.

High-throughput retina-array for screening 93 genes involved in inherited retinal dystrophy

PURPOSE

Retinal dystrophy (RD) is a broad group of hereditary disorders with heterogeneous genotypes and phenotypes. Current available genetic testing for these diseases is complicated, time consuming, and expensive. This study was conducted to develop and apply a microarray-based, high-throughput resequencing system to detect sequence alterations in genes related to inherited RD.

METHODS

A customized 300-kb resequencing chip, Retina-Array, was developed to detect sequence alterations of 267,550 bases of both sense and antisense sequence in 1470 exons spanning 93 genes involved in inherited RD. Retina-Array was evaluated in 19 patient samples with inherited RD provided by the eyeGENE repository and four Centre d'Etudes du Polymorphisme Humaine reference samples through a high-throughput experimental approach that included an automated PCR assay setup and quantification, efficient post-quantification data processing, optimized pooling and fragmentation, and standardized chip processing.

RESULTS

The performance of the chips demonstrated that the average base pair call rate and accuracy were 93.56% and 99.86%, respectively. In total, 304 candidate variations were identified using a series of customized screening filters. Among 174 selected variations, 123 (70.7%) were further confirmed by dideoxy sequencing. Analysis of patient samples using Retina-Array resulted in the identification of 10 known mutations and 12 novel variations with high probability of deleterious effects.

CONCLUSIONS

This study suggests that Retina-Array might be a valuable tool for the detection of disease-causing mutations and disease severity modifiers in a single experiment. Retinal-Array may provide a powerful and feasible approach through which to study genetic heterogeneity in retinal diseases.

Genetic variation and linkage disequilibrium in Bacillus anthracis

We performed whole-genome amplification followed by hybridization of custom-designed resequencing arrays to resequence 303 kb of genomic sequence from a worldwide panel of 39 Bacillus anthracis strains. We used an efficient algorithm contained within a custom software program, UniqueMER, to identify and mask repetitive sequences on the resequencing array to reduce false-positive identification of genetic variation, which can arise from cross-hybridization. We discovered a total of 240 single nucleotide variants (SNVs) and showed that B. anthracis strains have an average of 2.25 differences per 10,000 bases in the region we resequenced. Common SNVs in this region are found to be in complete linkage disequilibrium. These patterns of variation suggest there has been little if any historical recombination among B. anthracis strains since the origin of the pathogen. This pattern of common genetic variation suggests a framework for recognizing new or genetically engineered strains.

Mitochondrial genome sequence analysis: a custom bioinformatics pipeline substantially improves Affymetrix MitoChip v2.0 call rate and accuracy

Background

Mitochondrial genome sequence analysis is critical to the diagnostic evaluation of mitochondrial disease. Existing methodologies differ widely in throughput, complexity, cost efficiency, and sensitivity of heteroplasmy detection. Affymetrix MitoChip v2.0, which uses a sequencing-by-genotyping technology, allows potentially accurate and high-throughput sequencing of the entire human mitochondrial genome to be completed in a cost-effective fashion. However, the relatively low call rate achieved using existing software tools has limited the wide adoption of this platform for either clinical or research applications. Here, we report the design and development of a custom bioinformatics software pipeline that achieves a much improved call rate and accuracy for the Affymetrix MitoChip v2.0 platform. We used this custom pipeline to analyze MitoChip v2.0 data from 24 DNA samples representing a broad range of tissue types (18 whole blood, 3 skeletal muscle, 3 cell lines), mutations (a 5.8 kilobase pair deletion and 6 known heteroplasmic mutations), and haplogroup origins. All results were compared to those obtained by at least one other mitochondrial DNA sequence analysis method, including Sanger sequencing, denaturing HPLC-based heteroduplex analysis, and/or the Illumina Genome Analyzer II next generation sequencing platform.

Results

An average call rate of 99.75% was achieved across all samples with our custom pipeline. Comparison of calls for 15 samples characterized previously by Sanger sequencing revealed a total of 29 discordant calls, which translates to an estimated 0.012% for the base call error rate. We successfully identified 4 known heteroplasmic mutations and 24 other potential heteroplasmic mutations across 20 samples that passed quality control.

Conclusions

Affymetrix MitoChip v2.0 analysis using our optimized MitoChip Filtering Protocol (MFP) bioinformatics pipeline now offers the high sensitivity and accuracy needed for reliable, high-throughput and cost-efficient whole mitochondrial genome sequencing. This approach provides a viable alternative of potential utility for both clinical diagnostic and research applications to traditional Sanger and other emerging sequencing technologies for whole mitochondrial genome analysis.

Large scale mtDNA sequencing reveals sequence and functional conservation as major determinants of homoplasmic mtDNA variant distribution

The mitochondrial DNA (mtDNA) is highly variable, containing large numbers of pathogenic mutations and neutral polymorphisms. The spectrum of homoplasmic mtDNA variation was characterized in 730 subjects and compared with known pathogenic sites. The frequency and distribution of variants in protein coding genes were inversely correlated with conservation at the amino acid level. Analysis of tRNA secondary structures indicated a preference of variants for the loops and some acceptor stem positions. This comprehensive overview of mtDNA variants distinguishes between regions and positions which are likely not critical, mainly conserved regions with pathogenic mutations and essential regions containing no mutations at all.

A new mitochondria-related disease showing myopathy with episodic hyper-creatine kinase-emia

OBJECTIVE

To elucidate the relationship between mitochondrial DNA (mtDNA) alterations and a mitochondrial disease with a distinct combination of characteristic symptoms, namely episodic hyper-creatine kinase (CK)-emia and mild myopathy.

METHODS

We selected 9 patients with mtDNA np8291 alteration from 586 patients suspected to have a mitochondrial disease, and assessed them clinically, pathologically, and genetically. These 9 patients had undiagnosed mitochondrial myopathy with episodic hyper-CK-emia, all showing similar symptoms and progression.

RESULTS

Patients had mild muscle weakness and episodic hyper-CK-emia triggered by infections or drugs. Five of 9 patients were initially diagnosed with other conditions, such as myasthenia gravis, polymyositis, viral myositis, and drug-induced myopathy, because these conditions were acute or subacute, and 9 patients showed the same 16 mtDNA alterations, which have been reported to be nonpathological polymorphisms. Muscle biopsy revealed ragged-red fibers, highly expressed succinate dehydrogenase staining fibers, and cytochrome c oxidase-deficient fibers. Because their mitochondrial sequence data was almost the same, and 9 patients live in widely separated cities in Japan, the alterations may have arisen from a single source.

INTERPRETATION

These findings suggest that mild myopathy with episodic hyper-CK-emia associated with some of the 16 mtDNA alterations or at least with their mitochondria, could be a novel mitochondrial disease. Therefore, we propose that this disease be named as "mitochondrial myopathy with episodic hyper-CK-emia (MIMECK)." These alterations could work concomitantly and probably modify the impact of medications or other environmental factors. We believe these findings provide an insight into a novel aspect of mitochondrial disease pathogenesis.

Comprehensive mutational analysis of LRRK2 reveals variants supporting association with autosomal dominant Parkinson's disease

Parkinson's disease (PD) is a neurodegenerative disorder characterized by neurodegeneration, most notably of dopaminergic neurons in the substantia nigra. To date, six causative genes have been identified including LRRK2, whose mutations are the most frequent in autosomal dominant PD (Ad-PD). We conducted a comprehensive mutational analysis of LRRK2 in 30 Ad-PD (11 Japanese and 19 Caucasian) families employing a DNA microarray-based resequencing system and direct nucleotide sequence analysis, and identified 23 variants including two known mutations, p.G2019S and p.I1371V, in three Caucasian families and one Caucasian family, respectively, a novel putative pathogenic mutation, p.N1221K, in one Japanese family, and a known nonsynonymous variant, p.G2385R, in two Japanese families. Detailed analysis of the frequency of p.G2385R among 100 Japanese Ad-PD, 73 sporadic PD (sPD) and 238 controls revealed that the frequency of the p.G2385R variant was significantly higher in Ad-PD than in controls (allele frequency, 9.0 vs 2.1%) (χ(2)=16.32, P=5.34 × 10(-5)). The p.G2385R variant, however, did not show complete cosegregation with PD. In addition, the frequency of p.G2385R was also higher in sPD than in controls, although not significant (allele frequency, 3.4 vs 2.1%) (χ(2)=0.76, P=0.38). These observations support the possibility that p.G2385R is associated with an increased risk of PD.

Targeted sequencing using Affymetrix CustomSeq Arrays

This unit provides a basic protocol for oligo hybridization-based sequencing technology and resulting data analysis specific to the Affymetrix GeneChip CustomSeq Resequencing Array platform. All steps and critical aspects related to array design, experimental protocols, data management, and base-calling algorithms are addressed. This unit is particularly appropriate for sequencing targeted regions of the genome of up to 300 kilobases. The basic technology is most suitable for detecting substitution mutations, unless targeted indel probes are added.

Investigation of parameters that affect the success rate of microarray-based allele-specific hybridization assays

Background

The development of microarray-based genetic tests for diseases that are caused by known mutations is becoming increasingly important. The key obstacle to developing functional genotyping assays is that such mutations need to be genotyped regardless of their location in genomic regions. These regions include large variations in G+C content, and structural features like hairpins.

Methods/Findings

We describe a rational, stable method for screening and combining assay conditions for the genetic analysis of 42 Phenylketonuria-associated mutations in the phenylalanine hydroxylase gene. The mutations are located in regions with large variations in G+C content (20–75%). Custom-made microarrays with different lengths of complementary probe sequences and spacers were hybridized with pooled PCR products of 12 exons from each of 38 individual patient DNA samples. The arrays were washed with eight buffers with different stringencies in a custom-made microfluidic system. The data were used to assess which parameters play significant roles in assay development.

Conclusions

Several assay development methods found suitable probes and assay conditions for a functional test for all investigated mutation sites. Probe length, probe spacer length, and assay stringency sufficed as variable parameters in the search for a functional multiplex assay. We discuss the optimal assay development methods for several different scenarios.

Potentiometric Detection of DNA Molecules Hybridization Using Gene Field Effect Transistor and Intercalator

We propose a new concept of a gene field effect transistor (FET) for detection of allele specific oligonucleotide hybridization, which is in principle based on charge density change at the gate insulator. The electrical characteristics of the gene FET were found to shift after specific binding of biomolecules at the surface of the gate insulator. Allele specific oligonucleotide hybridization and reaction between double-stranded DNA and intercalator were successfully detected with gene FETs because they have intrinsic charges in an aqueous solution. Ability to discriminate single nucleotide polymorphism (SNP) was also examined using the gene FET. Our results show that control of hybridization temperature and utilization of intercalator lead to more precise SNP analysis using the gene FET.

Reconstructed ancestral sequences improve pathogen identification using resequencing DNA microarrays

We describe the benefit of using reconstructed ancestral sequences (RAS) on resequencing microarrays for rapid pathogen identification, with Enterobacteriaceae rpoB sequences as a model. Our results demonstrate a sharp improvement of call rate and accuracy when using RASs as compared to extant sequences. This improvement was attributed to the lower sequence divergence of RASs, which also expanded the sequence space covered by the microarray. Extension of this novel microarray design strategy to viruses, antimicrobial resistance elements or toxins is straightforward.

The genetics of blood pressure and hypertension: the role of rare variation

The role of heredity in influencing blood pressure and risk of hypertension is well recognized. However, progress in identifying specific genetic variation that contributes to heritability is very limited. This is in spite of completion of the human genome sequence, the development of extraordinary amounts of information about genome sequence variation and the investigation of blood pressure inheritance in linkage analysis, candidate gene studies and, most recently genome-wide association studies. This paper considers the progress of this research and the obstacles that have been encountered. This work has made clear that the genetic architecture of blood pressure regulation in the population is not likely to be shaped by commonly occurring genetic variation in a discrete set of blood pressure-influencing genes. Rather heritability may be accounted for by rare variation that has its biggest impact within pedigrees rather than on the population at large. Rare variants in a wide range of genes are likely to be the focus of high blood pressure genetics for the next several years and the emerging strategies that can be applied to uncover this genetic variation and the problems that must confronted are considered.

The immunosuppressive pipeline: meeting unmet needs in liver transplantation

Liver transplantation is now recognized as the treatment of choice for end-stage liver failure. Its success can be attributed largely to the generation of selective immunosuppressive agents, which have resulted in a dramatic reduction in the incidence of acute rejection and improvements in the short- and long-term outcomes of patients. However, the unresolved limitation of current immunosuppressive agents is long-term toxicity, which results in increases in the incidence and severity of cardiovascular, neurological, and renal diseases. Our recent understanding of the pathways of cell activation has resulted in the development of a new generation of immunosuppressive agents that may address the challenges facing transplantation today and allow the minimization or substitution of existing agents. Furthermore, advances in our understanding of the mechanisms of tolerance and the identification of biomarker signatures hold the promise that in some patients transplantation may be able to be performed without the need for long-term immunosuppression (tolerance).