An SNP map of human chromosome 22

A biallelic multiple nucleotide length polymorphism explains functional causality at 5p15.33 prostate cancer risk locus

To date, single-nucleotide polymorphisms (SNPs) have been the most intensively investigated class of polymorphisms in genome wide associations studies (GWAS), however, other classes such as insertion-deletion or multiple nucleotide length polymorphism (MNLPs) may also confer disease risk. Multiple reports have shown that the 5p15.33 prostate cancer risk region is a particularly strong expression quantitative trait locus (eQTL) for Iroquois Homeobox 4 (IRX4) transcripts. Here, we demonstrate using epigenome and genome editing that a biallelic (21 and 47 base pairs (bp)) MNLP is the causal variant regulating IRX4 transcript levels. In LNCaP prostate cancer cells (homozygous for the 21 bp short allele), a single copy knock-in of the 47 bp long allele potently alters the chromatin state, enabling de novo functional binding of the androgen receptor (AR) associated with increased chromatin accessibility, Histone 3 lysine 27 acetylation (H3K27ac), and ~3-fold upregulation of IRX4 expression. We further show that an MNLP is amongst the strongest candidate susceptibility variants at two additional prostate cancer risk loci. We estimated that at least 5% of prostate cancer risk loci could be explained by functional non-SNP causal variants, which may have broader implications for other cancers GWAS. More generally, our results underscore the importance of investigating other classes of inherited variation as causal mediators of human traits.

A novel and cost-efficient allele-specific PCR method for multiple SNP genotyping in a single run

Cost-effective methods for DNA genotyping were needed because single nucleotide polymorphisms (SNPs) were essential biomarkers associated with many diseases. Allele-specific PCR (AS-PCR) has the advantages of mature instruments and high sensitivity. But conventional AS-PCR needs to multiply the number of reactions or primers for multiple targets, which complicates the operation and increases the cost. Herein, we proposed a novel AS-PCR method for multiple SNP genotyping in a single run. Wild-type allele-specific primer (WT primer) was designed for each target gene. The sample and WT primers only needed to undergo multiplexed AS-PCR once simultaneously. After AS-PCR, the concentration of remaining primers varied among the samples of each genotype combination, due to the different matching performance between template and WT primers. The remaining primers then triggered multiplexed molecular beacon-rolling circle amplification, and the molecular beacons labelled with different fluorescent dyes corresponded to different targets. The fluorescence ratios of the sample to the positive control were used as the genotyping indexes. This method was able to detect samples with concentrations as low as 10 fM. We successfully applied the method to the multiple genotyping of 23 hair root samples for ADH1B and ALDH2 genes, obtaining completely consistent results with sequencing. The reagent cost was 0.6 dollar for one sample, showing a good cost performance. This proposed approach had a great application prospect in simultaneously rapid and accurate genotyping of multi-SNPs, and provided a new method for personalized health management.

Tristetraprolin Gene Single-Nucleotide Polymorphisms and mRNA Level in Patients With Rheumatoid Arthritis

To observe and evaluate the correlation between single-nucleotide polymorphisms (SNPs) and messenger RNA (mRNA) level related to tristetraprolin (TTP) in Chinese rheumatoid arthritis (RA). TapMan SNP was used for genotyping analysis in 580 RA patients and 554 healthy people. Association between TTP gene polymorphisms (rs251864 and rs3746083) and RA was obtained. Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) technology was applied for the detection of TTP mRNA level in peripheral blood mononuclear cells (PBMCs) in 36 RA patients and 37 healthy people. We observed that the allele T of TTP rs3746083 increased RA susceptibility (p = 0.019). A significant difference was found under the dominant model of rs3746083 (p = 0.037). Further analysis showed the allele distribution of rs3746083 was nominally correlated with RF phenotype of RA patients (p = 0.045). Nevertheless, the association between TTP rs251864 and the incidence of RA was no statistically significant (p > 0.05). The TTP expression level in PBMCs of RA patients was significantly reduced (p < 0.001). In conclusion, the results of this experiment support that TTP may be involved in the pathogenesis of RA.

Electrochemical Impedimetric Study of Non-Watson-Crick Base Pairs of DNA

Electrochemical impedance spectroscopy (EIS) was used to detect non-Watson-Crick base pairs of DNA. Thiol-modified DNA as a probe and mercaptohexanol (MCH) were co-immobilized to form a DNA/MCH mixed self-assembled monolayer on a gold electrode surface and then hybridized with complementary DNAs. The DNA layers were measured by the EIS method and interpreted by equivalent circuits. Every terminal base mismatch of the DNA duplex brought about an increase in the charge-transfer resistance (R_ct), unlike the case with a fully matched DNA duplex. The value of R_ct was highly sensitive to the number of base mismatches for both unpaired and overhang DNA at the terminal. For internal base mismatches, however, no significant increase in R_ct was observed. These experimental results proved that the charge transfer of redox molecules to the electrode surface is largely hindered by an end fraying motion due to base unpairing and dangling overhang. EIS was able to detect these steric properties of DNA strands. Furthermore, an electrode modified with G-quadruplex (G4) DNA demonstrated the influences of bulkiness and loop structure on the accessibility of the redox probe to the electrode.

INDEL detection, the 'Achilles heel' of precise genome editing: a survey of methods for accurate profiling of gene editing induced indels

Advances in genome editing technologies have enabled manipulation of genomes at the single base level. These technologies are based on programmable nucleases (PNs) that include meganucleases, zinc-finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs) and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/CRISPR-associated 9 (Cas9) nucleases and have given researchers the ability to delete, insert or replace genomic DNA in cells, tissues and whole organisms. The great flexibility in re-designing the genomic target specificity of PNs has vastly expanded the scope of gene editing applications in life science, and shows great promise for development of the next generation gene therapies. PN technologies share the principle of inducing a DNA double-strand break (DSB) at a user-specified site in the genome, followed by cellular repair of the induced DSB. PN-elicited DSBs are mainly repaired by the non-homologous end joining (NHEJ) and the microhomology-mediated end joining (MMEJ) pathways, which can elicit a variety of small insertion or deletion (indel) mutations. If indels are elicited in a protein coding sequence and shift the reading frame, targeted gene knock out (KO) can readily be achieved using either of the available PNs. Despite the ease by which gene inactivation in principle can be achieved, in practice, successful KO is not only determined by the efficiency of NHEJ and MMEJ repair; it also depends on the design and properties of the PN utilized, delivery format chosen, the preferred indel repair outcomes at the targeted site, the chromatin state of the target site and the relative activities of the repair pathways in the edited cells. These variables preclude accurate prediction of the nature and frequency of PN induced indels. A key step of any gene KO experiment therefore becomes the detection, characterization and quantification of the indel(s) induced at the targeted genomic site in cells, tissues or whole organisms. In this survey, we briefly review naturally occurring indels and their detection. Next, we review the methods that have been developed for detection of PN-induced indels. We briefly outline the experimental steps and describe the pros and cons of the various methods to help users decide a suitable method for their editing application. We highlight recent advances that enable accurate and sensitive quantification of indel events in cells regardless of their genome complexity, turning a complex pool of different indel events into informative indel profiles. Finally, we review what has been learned about PN-elicited indel formation through the use of the new methods and how this insight is helping to further advance the genome editing field.

Ultrasensitive Multiplex Detection of Single Nucleotide Polymorphisms Based on Short-Chain Hybridization Combined with Online Preconcentration of Capillary Electrophoresis

Reliable multiple single nucleotide polymorphisms (SNPs) detection at low abundance is of great significance for disease diagnosis and biomedical research. Herein, we have developed a novel and simple method for multiple SNPs detection combining solid-phase capture by specific hybridization with online preconcentration of capillary gel electrophoresis-laser-induced fluorescence (CGE-LIF). The method presents an excellent performance due to its favorable traits: the solid-phase short-chain hybridization ensures the high specificity of SNP detection; the effective separation ability of CGE can easily achieve multiplex detection; the simple online preconcentration significantly improves the detection sensitivity of fluorescent probe by nearly 100-fold. For a single SNP target, the assay achieves a limit of detection as low as 0.01-0.02% for three different NRAS mutations in the same codon. For multiple SNP targets, as low as 0.05% abundance can be easily realized. Our method is simple, efficient, ultrasensitive, and universal for multiple SNPs detection without complex enzymatic or chemical ligation reaction, which shows great potential in early clinical diagnosis.

Indel detection from Whole Genome Sequencing data and association with lipid metabolism in pigs

The selection in commercial swine breeds for meat-production efficiency has been increasing among the past decades, reducing the intramuscular fat content, which has changed the sensorial and technological properties of pork. Through processes of natural adaptation and selective breeding, the accumulation of mutations has driven the genetic divergence between pig breeds. The most common and well-studied mutations are single-nucleotide polymorphisms (SNPs). However, insertions and deletions (indels) usually represents a fifth part of the detected mutations and should also be considered for animal breeding. In the present study, three different programs (Dindel, SAMtools mpileup, and GATK) were used to detect indels from Whole Genome Sequencing data of Iberian boars and Landrace sows. A total of 1,928,746 indels were found in common with the three programs. The VEP tool predicted that 1,289 indels may have a high impact on protein sequence and function. Ten indels inside genes related with lipid metabolism were genotyped in pigs from three different backcrosses with Iberian origin, obtaining different allelic frequencies on each backcross. Genome-Wide Association Studies performed in the Longissimus dorsi muscle found an association between an indel located in the C1q and TNF related 12 (C1QTNF12) gene and the amount of eicosadienoic acid (C20:2(n-6)).

CXCL12 and CXCR4 polymorphisms and expressions in peripheral blood from patients of hepatocellular carcinoma

AIM

To determine if CXCL12 (rs1801157) and CXCR4 (rs2228014) polymorphisms are associated with hepatocellular carcinoma (HCC) susceptibility, and detect their expressions in peripheral blood.

METHODS

206 HCC patients, 252 chronic hepatitis B patients, 221 liver cirrhosis patients and 275 healthy volunteers were recruited. Genes CXCL12 and CXCR4 were amplified and genotyped. Their expression in peripheral blood were detected.

RESULTS

CXCL12 rs1801157 and CXCR4 rs2228014 polymorphisms were associated with increased susceptibility of HCC, and genotypes GA/AA and CT/TT may be risk factors of HCC (all p < 0.05). Expressions of CXCL12 and CXCR4 in peripheral blood from HCC patients increased significantly (p < 0.05).

CONCLUSION

CXCL12 and CXCR4 polymorphisms may be risk factors for HCC, and they may be potential HCC markers.

Spatial Analysis of Nanofluidic-Embedded Biosensors for Wash-Free Single-Nucleotide Difference Discrimination

In this work, we demonstrate that the analysis of spatially resolved nanofluidic-embedded biosensors permits the fast and direct discrimination of single-nucleotide difference (SND) within oligonucleotide sequences in a single step interaction. We design a sensor with a linear dimension much larger than the channel depth in order to ensure that the reaction over the whole sensor is limited by the convection rate. Thus, the targets are fully collected, inducing a nonuniform spatial hybridization profile. We also use the nanoscale height of the channel, which enables us to minimize the amount of labeled molecules flowing over the sensor and hence to reduce the fluorescence background, to carry out real-time hybridization detection by fluorescence microscopy. Taken together, these design rules allow us to show that the spatial hybridization profile depends on the duplex affinity, and we speculate that the on and off-rate constants can be inferred during target injection, which is not possible in local analysis where the dissociation step through rinsing must be conducted. We finally manage to discriminate a GT mismatch on a microRNA sequence by optimizing the interaction temperature and the probe design after a few minutes of interaction in a single step protocol. This work may be applied to any biosensing transduction scheme with spatial resolution, e.g., surface plasmon resonance imaging, integrated into nanofluidic channels for applications where high oligonucleotide sequence selectivity and short analysis times are required.

Correlation between survivin polymorphism and acute leukemia of children

The correlation between the variations in the polymorphic sites of survivin, rs9904341C/G and rs8073069C/G, and the pathogenesis of acute leukemia, as well as the guiding significance in clinical practice were investigated. We enrolled a total of 182 children with acute leukemia and 200 healthy children as the subjects. In accordance with the case-control method, the polymerase chain reaction was carried out for genetic typing of the two polymorphic sites, rs9904341C/G and rs8073069C/G. In the case group and the healthy group, the frequencies of C and G alleles in rs9904341C/G of survivin were 59.3 and 41.7%, and 46.7 and 50.3%, respectively, and the pairwise comparison showed statistically significant differences (P=0.008). Additionally, the frequencies of genotypes, C/C, C/G and G/G, were 38.5 and 41.7%; 19.8 and 26.5%; 16.5 and 27.0% in the case group and the healthy group, respectively, and the differences in comparisons showed statistical significance (P=0.033). The genotype frequency of C/C in the case group was 38.5%, significantly higher than that in the healthy group (26.5%). Compared with C/C, the risk coefficient of leukemia in patients with genotypes of C/G or G/G was significantly decreased. In the case group and the healthy group, the frequencies of C and G alleles in rs8073069C/G of survivin were 30.5 and 69.5%; 27.7 and 72.3%, respectively, and the pairwise comparison showed no statistically significant differences (P=0.404). Additionally, the frequencies of genotypes, C/C, C/G and G/G, were 11 and 39.0%; 50.0 and 9.0%; 37.5 and 53.5% in the case group and the healthy group, respectively, and the differences in comparisons showed no statistical significance (P=0.62). Compared with the genotype of C/C, we found that the risk of leukemia was not affected in patients with genotypes of C/G and G/G. In conclusion, the SNP of rs9904341C/G in survivin may be correlated with the risk of acute leukemia, and compared with C/C genotype, patients with C/G or G/G may have a decreased risk of acute leukemia. In survivin, rs8073069C/G may have no correlation with the risk of acute leukemia.

Association of interferon gamma gene polymorphism and susceptibility to hepatitis C virus infection in Egyptian patients: A multicenter, family-based study

Background and Aim

Polymorphisms in some genes may influence the persistence of hepatitis C virus (HCV) infection, clinical outcome, HCV replication, and liver damage. This study was conducted to investigate the role of the interferon gamma (IFN‐γ) gene at (+874 T/A, −764 G/C, −179 C/A) single‐nucleotide polymorphisms (SNPs) and its receptor (IFN‐γR2) at (rs 2786067 A/C) SNP in the susceptibility of Egyptian families to HCV infection with high‐resolution techniques.

Methods

In total, 517 Egyptian families, with 2246 subjects, were recruited to this study from the Upper and Lower Egypt governorates and were classified into three groups: 1034 patients with chronic hepatitis C virus, 108 subjects with spontaneous virus clearance (SVC), and 1104 subjects as a healthy control group. All subjects were genotyped for (+874 T/A, rs2430561, −764 G/C, rs2069707, −179 C/A, rs2069709, and rs 27860067, A/C) SNPs of the IFN‐γ gene using the allelic discrimination real‐time polymerase chain reaction technique and were confirmed using sequence‐based typing.

Results

The carriage of T allele of (+874) IFN‐γ is a risky allele and was significantly higher in chronic hepatitis C more than other two groups (odds ratio [OR]: 2.6646, P < 0.0002). On the other hand, the C allele of (−764, rs2069707) is a protective allele and was higher in SVC than the other two groups (OR: 0.2709, P < 0.0001). However, both (−179 C/A, rs 2069709) and (rs 27860067, A/C) SNPs are not polymorphic enough to be studied in the Egyptian population.

Conclusions

HCV infection is associated with the T allele of (+874 rs2430561), while SVC of HCV is associated with the C allele of (−764, rs2069707) of the IFN‐γ gene.

Genome annotation for clinical genomic diagnostics: strengths and weaknesses

The Human Genome Project and advances in DNA sequencing technologies have revolutionized the identification of genetic disorders through the use of clinical exome sequencing. However, in a considerable number of patients, the genetic basis remains unclear. As clinicians begin to consider whole-genome sequencing, an understanding of the processes and tools involved and the factors to consider in the annotation of the structure and function of genomic elements that might influence variant identification is crucial. Here, we discuss and illustrate the strengths and weaknesses of approaches for the annotation and classification of important elements of protein-coding genes, other genomic elements such as pseudogenes and the non-coding genome, comparative-genomic approaches for inferring gene function, and new technologies for aiding genome annotation, as a practical guide for clinicians when considering pathogenic sequence variation. Complete and accurate annotation of structure and function of genome features has the potential to reduce both false-negative (from missing annotation) and false-positive (from incorrect annotation) errors in causal variant identification in exome and genome sequences. Re-analysis of unsolved cases will be necessary as newer technology improves genome annotation, potentially improving the rate of diagnosis.

Interaction research on the antiviral molecule dufulin targeting on southern rice black streaked dwarf virus p9-1 nonstructural protein

Southern rice black streaked dwarf virus (SRBSDV) causes severe harm to rice production. Unfortunately, studies on effective antiviral drugs against SRBSDV and interaction mechanism of antiviral molecule targeting on SRBSDV have not been reported. This study found dufulin (DFL), an ideal anti-SRBSDV molecule, and investigated the interactions of DFL targeting on the nonstructural protein P9-1. The biological sequence information and bonding characterization of DFL to four kinds of P9-1 protein were described with fluorescence titration (FT) and microscale thermophoresis (MST) assays. The sequence analysis indicated that P9-1 had highly-conserved C- and N-terminal amino acid residues and a hypervariable region that differed from 131 aa to 160 aa. Consequently, wild-type (WT-His-P9-1), 23 C-terminal residues truncated (TR-ΔC23-His-P9-1), 6 N-terminal residues truncated (TR-ΔN6-His-P9-1), and Ser138 site-directed (MU-138-His-P9-1) mutant proteins were expressed. The FT and MST assay results indicated that DFL bounded to WT-His-P9-1 with micromole affinity and the 23 C-terminal amino acids were the potential targeting site. This system, which combines a complete sequence analysis, mutant protein expression, and binding action evaluating system, could further advance the understanding of the interaction abilities between antiviral drugs and their targets.

Bioinformatics tools for discovery and functional analysis of single nucleotide polymorphisms

With the high speed DNA sequencing of genome, databases of genome data continue to grow, and the understanding of genetic variation between individuals grows as well. Single nucleotide polymorphisms (SNPs), a main type of genetic variation, are increasingly important resource for understanding the structure and function of the human genome and become a valuable resource for investigating the genetic basis of disease. During the past years, in addition to experimental approaches to characterize specific variants, intense bioinformatics techniques were applied to understand effects of these genetic changes. In the genetics studies, one intends to understand the molecular basis of disease, and computational methods are becoming increasingly important for SNPs selection, prediction and understanding the downstream effects of genetic variation. The review provides systematic information on the available resources and methods for SNPs discovery and analysis. We also report some new results on DNA sequence-based prediction of SNPs in human cytochrome P450, which serves as an example of computational methods to predict and discovery SNPs. Additionally, annotation and prediction of functional SNPs, as well as a comprehensive list of existing tools and online recourses, are reviewed and described.

Brain investigation and brain conceptualization

The brain of a patient with Alzheimer's disease (AD) undergoes changes starting many years before the development of the first clinical symptoms. The recent availability of large prospective datasets makes it possible to create sophisticated brain models of healthy subjects and patients with AD, showing pathophysiological changes occurring over time. However, these models are still inadequate; representations are mainly single-scale and they do not account for the complexity and interdependence of brain changes. Brain changes in AD patients occur at different levels and for different reasons: at the molecular level, changes are due to amyloid deposition; at cellular level, to loss of neuron synapses, and at tissue level, to connectivity disruption. All cause extensive atrophy of the whole brain organ. Initiatives aiming to model the whole human brain have been launched in Europe and the US with the goal of reducing the burden of brain diseases. In this work, we describe a new approach to earlier diagnosis based on a multimodal and multiscale brain concept, built upon existing and well-characterized single modalities.

Association between IL-10 gene promoter polymorphism and hepatitis B viral infection in an Egyptian population

Cytokines play critical roles in the pathogenesis of hepatitis B virus infection (HBV). This work was designed to study the effect of IL-10 gene polymorphisms (-1082G/A and -819C/T) on susceptibility of Egyptians to HBV. Genotyping was performed using single-stranded polymorphism-polymerase chain reaction in 118 Egyptian hepatitis B patients and 119 healthy controls, and IL-10 serum levels were measured using ELISA. The frequency of IL-10 -1082G/G was significantly higher in HBV patients than in healthy controls, and G/A and A/A were not significantly different between groups. The distribution of IL-10 -819 genotypes was not significantly different between the HBV and healthy control groups. Although AT was significantly different between controls and patients, the distribution of the other haplotypes was not. IL-10 levels were significantly lower among hepatitis B patients. Our data stress the importance of IL-10 gene polymorphism in HBV infection. Depending on our preliminary work, IL-10 -1082G/G may act as a host genetic factor in the susceptibility to HBV infection in Egyptians.

Sequencing of the core MHC region of black grouse (Tetrao tetrix) and comparative genomics of the galliform MHC

BACKGROUND

The MHC, which is regarded as the most polymorphic region in the genomes of jawed vertebrates, plays a central role in the immune system by encoding various proteins involved in the immune response. The chicken MHC-B genomic region has a highly streamlined gene content compared to mammalian MHCs. Its core region includes genes encoding Class I and Class IIB molecules but is only ~92Kb in length. Sequences of other galliform MHCs show varying degrees of similarity as that of chicken. The black grouse (Tetrao tetrix) is a wild galliform bird species which is an important model in conservation genetics and ecology. We sequenced the black grouse core MHC-B region and combined this with available data from related species (chicken, turkey, gold pheasant and quail) to perform a comparative genomics study of the galliform MHC. This kind of analysis has previously been severely hampered by the lack of genomic information on avian MHC regions, and the galliformes is still the only bird lineage where such a comparison is possible.

RESULTS

In this study, we present the complete genomic sequence of the MHC-B locus of black grouse, which is 88,390 bp long and contains 19 genes. It shows the same simplicity as, and almost perfect synteny with, the corresponding genomic region of chicken. We also use 454-transcriptome sequencing to verify expression in 17 of the black grouse MHC-B genes. Multiple sequence inversions of the TAPBP gene and TAP1-TAP2 gene block identify the recombination breakpoints near the BF and BLB genes. Some of the genes in the galliform MHC-B region also seem to have been affected by selective forces, as inferred from deviating phylogenetic signals and elevated rates of non-synonymous nucleotide substitutions.

CONCLUSIONS

We conclude that there is large synteny between the MHC-B region of the black grouse and that of other galliform birds, but that some duplications and rearrangements have occurred within this lineage. The MHC-B sequence reported here will provide a valuable resource for future studies on the evolution of the avian MHC genes and on links between immunogenetics and ecology of black grouse.

Exploring the implications of INDELs in neuropsychiatric genetics: challenges and perspectives

The decade passed after publishing the Human Genome first draft faced an enormous growth at the understanding of the genomic variation among different subjects, populations, and groups of patients. Single nucleotide polymorphisms (SNPs) and insertion or deletions (INDELs) have been increasingly recognized as a major type of genetic variations, with potential impact in protein activities and gene expression changes observed in complex genetic traits, like neuropsychiatric diseases. INDELs represent the second most common class of variations after SNPs, but there is still an important gap between the number of INDELs reported and the actual knowledge about the functional implications of such variations. There are approximately 10 million SNPs already reported, and the human populations are expected to collectively harbor at least 1.6-2.5 million INDELs. One of the major challenges is to find better platforms to screen for INDELs in a high throughput manner. The discordance in between the data from different studies might be explained by the diverse approaches employed to sequence the genomes with variable platforms. Short INDEL variations increased the scope of genetic markers in human genetic diseases, and various studies showed that common microdeletions and smaller INDELs might be highly associated with neuropsychiatric diseases such as schizophrenia, autism, mental retardation, and Alzheimer disease. The rapidly increasing amount of resequencing, genotyping, and personal genome data generated by large-scale genetic human projects require the development of integrated bioinformatics tools able to efficiently manage and analyze these genetic data. Our group is currently dealing with different approaches that might optimize sequencing and bioinformatics analyses of short INDELs to broaden our research capabilities of identifying those intriguing genetic variations. Hopefully, INDELs might become a new trend in association studies in neuropsychiatric genetics since so far the level of significant and positive associations with the standard SNPs reported presents limited predictive application.

Strategies for genotyping

The identification of genomic loci linked to or associated with human disease has been greatly facilitated by the evolution of genotyping strategies and techniques. The success of these strategies continues to be based upon clear clinical assessment, accurate sample handling, and careful data management, but also increasingly upon experimental design. Technological advances in the field of genotyping have permitted increasingly complex and large population studies to be performed. An understanding of publicly available genetic variation databases, including an awareness of the limitations of these data, and an appreciation of the strategic approaches that should be used to exploit this information will provide tremendous insight for researchers are aiming to utilize this accessible technology. As genome-wide association studies (GWAS) and Next Generation (NextGen) sequencing become the mainstays of genetic analyses, it is important that their technical strengths and limitations, as well as their impact on study design, be understood before use in a linkage or genetic association study.

Small insertions and deletions (INDELs) in human genomes

In this review, we focus on progress that has been made with detecting small insertions and deletions (INDELs) in human genomes. Over the past decade, several million small INDELs have been discovered in human populations and personal genomes. The amount of genetic variation that is caused by these small INDELs is substantial. The number of INDELs in human genomes is second only to the number of single nucleotide polymorphisms (SNPs), and, in terms of base pairs of variation, INDELs cause similar levels of variation as SNPs. Many of these INDELs map to functionally important sites within human genes, and thus, are likely to influence human traits and diseases. Therefore, small INDEL variation will play a prominent role in personalized medicine.

Mining SNPs from EST sequences using filters and ensemble classifiers

Abundant single nucleotide polymorphisms (SNPs) provide the most complete information for genome-wide association studies. However, due to the bottleneck of manual discovery of putative SNPs and the inaccessibility of the original sequencing reads, it is essential to develop a more efficient and accurate computational method for automated SNP detection. We propose a novel computational method to rapidly find true SNPs in public-available EST (expressed sequence tag) databases; this method is implemented as SNPDigger. EST sequences are clustered and aligned. SNP candidates are then obtained according to a measure of redundant frequency. Several new informative biological features, such as the structural neighbor profiles and the physical position of the SNP, were extracted from EST sequences, and the effectiveness of these features was demonstrated. An ensemble classifier, which employs a carefully selected feature set, was included for the imbalanced training data. The sensitivity and specificity of our method both exceeded 80% for human genetic data in the cross validation. Our method enables detection of SNPs from the user's own EST dataset and can be used on species for which there is no genome data. Our tests showed that this method can effectively guide SNP discovery in ESTs and will be useful to avoid and save the cost of biological analyses.

The Schistosoma japonicum genome reveals features of host-parasite interplay

Schistosoma japonicum is a parasitic flatworm that causes human schistosomiasis, which is a significant cause of morbidity in China and the Philippines. Here we present a draft genomic sequence for the worm. The genome provides a global insight into the molecular architecture and host interaction of this complex metazoan pathogen, revealing that it can exploit host nutrients, neuroendocrine hormones and signalling pathways for growth, development and maturation. Having a complex nervous system and a well-developed sensory system, S. japonicum can accept stimulation of the corresponding ligands as a physiological response to different environments, such as fresh water or the tissues of its intermediate and mammalian hosts. Numerous proteases, including cercarial elastase, are implicated in mammalian skin penetration and haemoglobin degradation. The genomic information will serve as a valuable platform to facilitate development of new interventions for schistosomiasis control.

Sequences, annotation and single nucleotide polymorphism of the major histocompatibility complex in the domestic cat

Two sequences of major histocompatibility complex (MHC) regions in the domestic cat, 2.976 and 0.362 Mbps, which were separated by an ancient chromosome break (55-80 MYA) and followed by a chromosomal inversion were annotated in detail. Gene annotation of this MHC was completed and identified 183 possible coding regions, 147 human homologues, possible functional genes and 36 pseudo/unidentified genes) by GENSCAN and BLASTN, BLASTP RepeatMasker programs. The first region spans 2.976 Mbp sequence, which encodes six classical class II antigens (three DRA and three DRB antigens) lacking the functional DP, DQ regions, nine antigen processing molecules (DOA/DOB, DMA/DMB, TAPASIN, and LMP2/LMP7,TAP1/TAP2), 52 class III genes, nineteen class I genes/gene fragments (FLAI-A to FLAI-S). Three class I genes (FLAI-H, I-K, I-E) may encode functional classical class I antigens based on deduced amino acid sequence and promoter structure. The second region spans 0.362 Mbp sequence encoding no class I genes and 18 cross-species conserved genes, excluding class I, II and their functionally related/associated genes, namely framework genes, including three olfactory receptor genes. One previously identified feline endogenous retrovirus, a baboon retrovirus derived sequence (ECE1) and two new endogenous retrovirus sequences, similar to brown bat endogenous retrovirus (FERVmlu1, FERVmlu2) were found within a 140 Kbp interval in the middle of class I region. MHC SNPs were examined based on comparisons of this BAC sequence and MHC homozygous 1.9x WGS sequences and found that 11,654 SNPs in 2.84 Mbp (0.00411 SNP per bp), which is 2.4 times higher rate than average heterozygous region in the WGS (0.0017 SNP per bp genome), and slightly higher than the SNP rate observed in human MHC (0.00337 SNP per bp).

Construction of a dense SNP map for bovine chromosome 6 to assist the assembly of the bovine genome sequence

A linkage map was constructed for bovine chromosome 6 (BTA6), using 399 single nucleotide polymorphisms (SNPs) detected primarily from PCR-resequencing. The efficiency of SNP detection was highly dependent on the source of sequence information chosen for primer design (BAC-end sequences, introns or promoters). The SNPs were used to build a linkage map comprising 104 cM on BTA6. The SNP order in the linkage map corresponded very well with radiation hybrid (RH) maps available for BTA6 as well as with expected positions in the human comparative map, but diverged significantly from the current assembly of the bovine genome (Btau_3.1). When performing linkage analysis with the marker order suggested from the Btau_3.1 we observed an expansion of the genetic map from 104 cM to 137 cM, strongly suggesting a reordering of scaffolds in the current version of the bovine genome assembly. The extent of LD on BTA6 was evaluated by calculating the average r(2) for SNP pairs separated by given distances. The decline of LD was rapid with distance, such that r(2) was 0.1 at 100 kb. Our results indicate that linkage mapping will be a valuable source of information for correcting errors in the current bovine assembly. These errors were sufficiently frequent to be of concern for the accuracy of mapping QTL with panels of SNPs whose positions are based on the current assembly.

Challenges and standards in integrating surveys of structural variation

There has been an explosion of data describing newly recognized structural variants in the human genome. In the flurry of reporting, there has been no standard approach to collecting the data, assessing its quality or describing identified features. This risks becoming a rampant problem, in particular with respect to surveys of copy number variation and their application to disease studies. Here, we consider the challenges in characterizing and documenting genomic structural variants. From this, we derive recommendations for standards to be adopted, with the aim of ensuring the accurate presentation of this form of genetic variation to facilitate ongoing research.

Identification and characterization of single nucleotide polymorphisms in 6 growth-correlated genes in porcine by denaturing high performance liquid chromatography

The growth-correlated genes that are part of the neuroendocrine growth axis play crucial roles in the regulation of growth and development of pig. The identification of genetic polymorphisms in these genes will enable the scientist to evaluate the biological relevance of such polymorphisms and to gain a better understanding of quantitative traits like growth. In the present study, seven pairs of primers were designed to obtain unknown sequences of growth-correlated genes, and other 25 pairs of primers were designed to identify single nucleotide polymorphisms (SNP) using the denaturing high-performance liquid chromatography (DHPLC) technology in four pig breeds (Duroc, Landrace, Lantang and Wuzhishan), significantly differing in growth and development characteristics. A total of 101 polymorphisms were discovered in 10,707 base pairs (bp) from six genes of the ghrelin (GHRL), leptin (LEP), insulin-like growth factor II (IGF-II), insulin-like growth factor binding protein 2 (IGFBP-2), insulin-like growth factor binding protein 3 (IGFBP-3), and somatostatin (SS). The observed average distances between the SNP in the 5'UTR, coding regions, introns and 3'UTR were 134, 521, 81 and 92 bp, respectively. Four SNPs were found in the coding regions of IGF-II, IGFBP-2 and LEP, respectively. Two synonymous mutations were obtained in IGF-II and LEP genes respectively, and two non-synonymous were found in IGFBP-2 and LEP genes, respectively. Seven other mutations were also observed. Thirty-two PCR-RFLP markers were found among 101 polymorphisms of the six genes. The SNP discovered in this study would provide suitable markers for association studies of candidate genes with growth related traits in pig.

PolyScan: an automatic indel and SNP detection approach to the analysis of human resequencing data

Small insertions and deletions (indels) and single nucleotide polymorphisms (SNPs) are common genetic variants that are thought to be associated with a wide variety of human diseases. Owing to the genome's size and complexity, manually characterizing each one of these variations in an individual is not practical. While significant progress has been made in automated single-base mutation discovery from the sequences of diploid PCR products, automated and reliable detection of indels continues to pose difficult challenges. In this paper, we present PolyScan, an algorithm and software implementation designed to provide de novo heterozygous indel detection and improved SNP identification in the context of high-throughput medical resequencing. Tests on a human diploid PCR-based sequence data set, consisting of 90,270 traces from 13 genes, indicate that PolyScan identified approximately 90% of the 151 consensus indel sites and approximately 84% of the 1546 heterozygous indels previously identified by manual inspection. Tests on tumor-derived data show that PolyScan better identifies high-quality, low-level mutations as compared with other mutation detection software. Moreover, SNP identification improves when reprocessing the results of other programs. These results suggest that PolyScan may play a useful role in the post human genome project research era.

An initial map of insertion and deletion (INDEL) variation in the human genome

Although many studies have been conducted to identify single nucleotide polymorphisms (SNPs) in humans, few studies have been conducted to identify alternative forms of natural genetic variation, such as insertion and deletion (INDEL) polymorphisms. In this report, we describe an initial map of human INDEL variation that contains 415,436 unique INDEL polymorphisms. These INDELs were identified with a computational approach using DNA re-sequencing traces that originally were generated for SNP discovery projects. They range from 1 bp to 9989 bp in length and are split almost equally between insertions and deletions, relative to the chimpanzee genome sequence. Five major classes of INDELs were identified, including (1) insertions and deletions of single-base pairs, (2) monomeric base pair expansions, (3) multi-base pair expansions of 2-15 bp repeat units, (4) transposon insertions, and (5) INDELs containing random DNA sequences. Our INDELs are distributed throughout the human genome with an average density of one INDEL per 7.2 kb of DNA. Variation hotspots were identified with up to 48-fold regional increases in INDEL and/or SNP variation compared with the chromosomal averages for the same chromosomes. Over 148,000 INDELs (35.7%) were identified within known genes, and 5542 of these INDELs were located in the promoters and exons of genes, where gene function would be expected to be influenced the greatest. All INDELs in this study have been deposited into dbSNP and have been integrated into maps of human genetic variation that are available to the research community.

Relationship between IFN-γ gene polymorphism and susceptibility to intrauterine HBV infection

AIM: To explore the susceptibility of children to intrauterine HBV infection by studying the relationship between IFN-γ gene polymorphism, including IFN-γ+874A/T single nucleotide polymorphism(SNP) and CA repeat microsatellite polymorphism and intrauterine HBV infection.

METHODS: A TaqMan fluorescence polymerase chain reaction in the IFN-γ+874A/T single nucleotide polymorphism was tested in the intrauterine HBV infection group(group I) and the normal immune children group(group II). Capillary electrophoresis was performed in the above two groups to assay the IFN-γ CA repeat microsatellite polymorphism.

RESULTS: Frequencies of AA, AT and TT genotypes were 67.4%, 19.6% and 13.0% in the intrauterine HBV infection group, and 45.2%, 30.1% and 24.7% in the normal immune children group, respectively. A significant difference was found in the frequency distribution of IFN-γ+874 genotype between the two groups (χ² = 5.102, P = 0.02389). In the intrauterine HBV infection group the AA genotype was more common than in the normal immune group. Frequency of IFN-γ+874A allele was 77.17% in the intrauterine HBV infection group, and 60.27% in the normal immune children group. In the intrauterine HBV infection group the IFN-γ+874A allele was more common than in normal immune group. A significant difference was found in the frequency distribution between the two groups (χ² = 7.238, P = 0.02389, OR = 2.228, 95% CI = 1.244-3.992). (CA₁₂)⁺/(CA₁₂)⁺ of IFN-γ CA microsatellite polymorphism was 11.90% in the intrauterine HBV infection group and 26.47% in the normal immune children group. A significant difference was found in the frequency distribution between the two groups (χ² = 5.64, P = 0.0176). Frequency of IFN-γ CA repeat was 25% in the intrauterine HBV infection group and 43.38% in the normal immune children group. The frequency of IFN-γ CA repeat was less in the intrauterine HBV infection group than in normal immune group. A significant difference was found in the frequency distribution between the two groups (χ² = 7.548, P = 0.0060).

CONCLUSION: There is a relationship between IFN-γ+874A/T SNP and intrauterine HBV infection as well as between IFN-γ CA microsatellite polymorphism and intrauterine HBV infection. IFN-γ gene polymorphism might be important in determining individual’s susceptibility to intrauterine HBV infection.

Common single nucleotide polymorphism of hypoxia-inducible factor-1alpha and its impact on the clinicopathological features of esophageal squamous cell carcinoma

OBJECTIVE

Angiogenesis is one of the most important molecular events in solid tumor development and growth, in which hypoxia-inducible factor (HIF)-1alpha is a key regulator and plays an important role. Studies have shown that a single nucleotide polymorphism (C1772T) in the HIF-1alpha gene exerts a large effect on the phenotype of human head and neck squamous cell carcinoma and renal cell carcinoma. But the impact of the C1772T polymorphism on the clinicopathological features of human esophageal squamous cell carcinoma (ESCC) remains unknown, and thus it is the main focus of our study.

METHODS

The C1772T genotype of 95 ESCC patients and 104 healthy controls were studied by using the polymerase chain reaction and restriction fragment length polymorphism. Mutations were confirmed by direct DNA sequencing. The impact of C1772T on tumor size, invasive depth, lymph node metastasis, distant metastasis, histological grade and TNM stage was also studied.

RESULTS

The genotype frequency observed in the patients and controls was 11.58% versus 10.58%, respectively, for genotype C/T (P > 0.05). Genotype T/T was not found in our study. Larger tumors and a higher rate of lymph node metastasis was found for the C/T group.

CONCLUSIONS

Although there is no significant difference of genotype distribution between ESCC patients and healthy controls, genotype C/T is associated with larger tumor and higher rate of lymph node metastasis.

SNPdetector: a software tool for sensitive and accurate SNP detection

Identification of single nucleotide polymorphisms (SNPs) and mutations is important for the discovery of genetic predisposition to complex diseases. PCR resequencing is the method of choice for de novo SNP discovery. However, manual curation of putative SNPs has been a major bottleneck in the application of this method to high-throughput screening. Therefore it is critical to develop a more sensitive and accurate computational method for automated SNP detection. We developed a software tool, SNPdetector, for automated identification of SNPs and mutations in fluorescence-based resequencing reads. SNPdetector was designed to model the process of human visual inspection and has a very low false positive and false negative rate. We demonstrate the superior performance of SNPdetector in SNP and mutation analysis by comparing its results with those derived by human inspection, PolyPhred (a popular SNP detection tool), and independent genotype assays in three large-scale investigations. The first study identified and validated inter- and intra-subspecies variations in 4,650 traces of 25 inbred mouse strains that belong to either the Mus musculus species or the M. spretus species. Unexpected heterozygosity in CAST/Ei strain was observed in two out of 1,167 mouse SNPs. The second study identified 11,241 candidate SNPs in five ENCODE regions of the human genome covering 2.5 Mb of genomic sequence. Approximately 50% of the candidate SNPs were selected for experimental genotyping; the validation rate exceeded 95%. The third study detected ENU-induced mutations (at 0.04% allele frequency) in 64,896 traces of 1,236 zebra fish. Our analysis of three large and diverse test datasets demonstrated that SNPdetector is an effective tool for genome-scale research and for large-sample clinical studies. SNPdetector runs on Unix/Linux platform and is available publicly (http://lpg.nci.nih.gov).

An interactive bovine in silico SNP database (IBISS)

An interactive bovine in silico SNP (IBISS) database has been created through the clustering and aligning of bovine EST and mRNA sequences. Approximately 324,000 EST and mRNA sequences were clustered to produce 29,965 clusters (producing 48,679 consensus sequences) and 48,565 singletons. A SNP screening regime was placed on variations detected in the multiple sequence alignment files to determine which SNPs are more likely to be real rather than sequencing errors. A small subset of predicted SNPs was validated on a diverse set of bovine DNA samples using PCR amplification and sequencing. Fifty percent of the predicted SNPs in the "putative >1" category were polymorphic in the population sampled. The IBISS database represents more than just a SNP database; it is also a genomic database containing uniformly annotated predicted gene mRNA and protein sequences, gene structure, and genomic organization information.

Genetics of common forms of obesity: a brief overview

The obesity epidemic is attributable to dietary and behavioral trends acting on a person's genetic makeup to determine body mass and susceptibility to obesity-related disease. Common forms of obesity have a strong hereditary component, yet genetic pathways that contribute to obesity have not yet been elucidated. Many genetic association studies have been reported, but few have been successfully replicated. New research tools and large studies will lead to an understanding of genes and their interaction to cause obesity, which may help guide successful interventions and treatments.

Annotation: velo-cardio-facial syndrome

BACKGROUND

Velo-cardio-facial syndrome (VCFS), the most frequent known interstitial deletion identified in man, is associated with chromosomal microdeletions in the q11 band of chromosome 22. Individuals with VCFS are reported to have a characteristic behavioural phenotype with high rates of behavioural, psychiatric, neuropsychological and linguistic disorders.

METHODS

A selective literature review was undertaken.

RESULTS

Children and adults with VCFS have high rates of behavioural, psychiatric and communication disorders. While VCFS children have high rates of ADHD, anxiety and affective disorders, adults have high rates of psychotic disorders, particularly schizophrenia. In addition, the presence of a chromosome 22q11 deletion is associated with specific neuropsychological and neuroanatomical abnormalities.

CONCLUSIONS

People with VCFS have a characteristic behavioural phenotype with high rates of behavioural, psychiatric, neuropsychological and communication disorders. Early diagnosis and treatment within a multidisciplinary framework is of paramount importance for VCFS individuals as this will have a major effect in determining the long-term outcome in affected individuals. Longitudinal studies of VCFS children are currently under way to identify precursor symptoms and areas of dysfunction which precede the later development of major psychiatric disorder. Identification of such prodromal features in VCFS may have enormous implications for the clinical management of major psychiatric disorder in VCFS and in the wider population.

Geographic stratification of linkage disequilibrium: a worldwide population study in a region of chromosome 22

Recent studies of haplotype diversity in a number of genomic regions have suggested that long stretches of DNA are preserved in the same chromosome, with little evidence of recombination events. The knowledge of the extent and strength of these haplotypes could become a powerful tool for future genetic analysis of complex traits. Different patterns of linkage disequilibrium (LD) have been found when comparing individuals of African and European descent, but there is scarce knowledge about the worldwide population stratification. Thus, the study of haplotype composition and the pattern of LD from a global perspective are relevant for elucidating their geographical stratification, as it may have implications in the future analysis of complex traits. We have typed 12 single nucleotide polymorphisms in a chromosome 22 region--previously described as having high LD levels in European populations -- in 39 different world populations. Haplotype structure has a clear continental structure with marked heterogeneity within some continents (Africa, America). The pattern of LD among neighbouring markers exhibits a strong clustering of all East Asian populations on the one hand and of Western Eurasian populations (including Europe) on the other, revealing only two major LD patterns, but with some very specific outliers due to specific demographic histories. Moreover, it should be taken into account that African populations are highly heterogeneous. The present results support the existence of a wide (but not total) communality in LD patterns in human populations from different continental regions, despite differences in their demographic histories, as population factors seem to be less relevant compared with genomic forces in shaping the patterns of LD.

Single nucleotide polymorphisms: aging and diseases

Differences of more than 3 million nucleotides can bee seen comparing the genomes of two individuals as a result of single nucleotide polymorphism (SNP). More and more SNPs can be identified and it seems that these alterations are behind of several biological phenomena. Personal differences in these nucleotides result for example in elevated disease susceptibilities, that is, certain nucleotides are more frequent in patients suffering from different diseases comparing to the healthy population. SNPs may cause substantial alterations in the cells, e.g. the enzyme activity of the respective gene changes, but in other cases the effects of the SNPs are not so pronounced. Later results indicate that SNPs can be rendered to individuals living a longer life than the average. Perhaps these results will not directly lead to the lengthening of the maximal life span; however, genes that play an important role in the aging process could be identified. In this respect SNPs are important factors in determining the information level of the cells of individuals which determines the maximal life span (I. Semsei On the nature of aging. Mech. Ageing Dev . 2000; 117: 93-108), in turn SNP is one of the factors that determine the aging process. Since there are certain age-related diseases, the discovery and the description of the SNPs as a function of age and diseases may result in a better understanding of the common roots of aging and those diseases.

Linkage disequilibrium and association studies in higher plants: present status and future prospects

During the last two decades, DNA-based molecular markers have been extensively utilized for a variety of studies in both plant and animal systems. One of the major uses of these markers is the construction of genome-wide molecular maps and the genetic analysis of simple and complex traits. However, these studies are generally based on linkage analysis in mapping populations, thus placing serious limitations in using molecular markers for genetic analysis in a variety of plant systems. Therefore, alternative approaches have been suggested, and one of these approaches makes use of linkage disequilibrium (LD)-based association analysis. Although this approach of association analysis has already been used for studies on genetics of complex traits (including different diseases) in humans, its use in plants has just started. In the present review, we first define and distinguish between LD and association mapping, and then briefly describe various measures of LD and the two methods of its depiction. We then give a list of different factors that affect LD without discussing them, and also discuss the current issues of LD research in plants. Later, we also describe the various uses of LD in plant genomics research and summarize the present status of LD research in different plant genomes. In the end, we discuss briefly the future prospects of LD research in plants, and give a list of softwares that are useful in LD research, which is available as electronic supplementary material (ESM).

Reduced variation on the chicken Z chromosome

Understanding the population genetic factors that shape genome variability is pivotal to the design and interpretation of studies using large-scale polymorphism data. We analyzed patterns of polymorphism and divergence at Z-linked and autosomal loci in the domestic chicken (Gallus gallus) to study the influence of mutation, effective population size, selection, and demography on levels of genetic diversity. A total of 14 autosomal introns (8316 bp) and 13 Z-linked introns (6856 bp) were sequenced in 50 chicken chromosomes from 10 highly divergent breeds. Genetic variation was significantly lower at Z-linked than at autosomal loci, with one segregating site every 39 bp at autosomal loci (theta(W) = 5.8 +/- 0.8 x 10(-3)) and one every 156 bp on the Z chromosome (theta(W) = 1.4 +/- 0.4 x 10(-3)). This difference may in part be due to a low male effective population size arising from skewed reproductive success among males, evident both in the wild ancestor-the red jungle fowl-and in poultry breeding. However, this effect cannot entirely explain the observed three- to fourfold reduction in Z chromosome diversity. Selection, in particular selective sweeps, may therefore have had an impact on reducing variation on the Z chromosome, a hypothesis supported by the observation of heterogeneity in diversity levels among loci on the Z chromosome and the lower recombination rate on Z than on autosomes. Selection on sex-linked genes may be particularly important in organisms with female heterogamety since the heritability of sex-linked sexually antagonistic alleles advantageous to males is improved when fathers pass a Z chromosome to their sons.

Mutant G-protein-coupled receptors as a cause of human diseases

G-protein-coupled receptors (GPCR) are involved in directly and indirectly controlling an extraordinary variety of physiological functions. Their key roles in cellular communication have made them the target for more than 60% of all currently prescribed drugs. Mutations in GPCR can cause acquired and inherited diseases such as retinitis pigmentosa (RP), hypo- and hyperthyroidism, nephrogenic diabetes insipidus, several fertility disorders, and even carcinomas. To date, over 600 inactivating and almost 100 activating mutations in GPCR have been identified which are responsible for more than 30 different human diseases. The number of human disorders is expected to increase given the fact that over 160 GPCR have been targeted in mice. Herein, we summarize the current knowledge relevant to understanding the molecular basis of GPCR function, with primary emphasis on the mechanisms underlying GPCR malfunction responsible for different human diseases.

The allele frequency spectrum in genome-wide human variation data reveals signals of differential demographic history in three large world populations

We have studied a genome-wide set of single-nucleotide polymorphism (SNP) allele frequency measures for African-American, East Asian, and European-American samples. For this analysis we derived a simple, closed mathematical formulation for the spectrum of expected allele frequencies when the sampled populations have experienced nonstationary demographic histories. The direct calculation generates the spectrum orders of magnitude faster than coalescent simulations do and allows us to generate spectra for a large number of alternative histories on a multidimensional parameter grid. Model-fitting experiments using this grid reveal significant population-specific differences among the demographic histories that best describe the observed allele frequency spectra. European and Asian spectra show a bottleneck-shaped history: a reduction of effective population size in the past followed by a recent phase of size recovery. In contrast, the African-American spectrum shows a history of moderate but uninterrupted population expansion. These differences are expected to have profound consequences for the design of medical association studies. The analytical methods developed for this study, i.e., a closed mathematical formulation for the allele frequency spectrum, correcting the ascertainment bias introduced by shallow SNP sampling, and dealing with variable sample sizes provide a general framework for the analysis of public variation data.

Mining single nucleotide polymorphisms from EST data of silkworm, Bombyx mori, inbred strain Dazao

We made use of 81,635 expressed sequence tags (ESTs) derived from 12 different cDNA libraries of the silkworm, Bombyx mori, inbred strain Dazao (P50), to identify high-quality candidate single nucleotide polymorphisms (SNPs). By PHRAP assembling, 12,980 contigs containing 11,537 contigs assembled by more than one read were obtained, and 101 candidate SNPs and 27 single base insertions/deletions were identified from 117 contigs assembled from 1576 high-quality reads base-called with PHRED and screened on the basis of the neighborhood quality standard (NQS). Simultaneously, we also predicted 40 SNPs in coding regions (cSNPs), of which 26 were predicted to lead to amino acid non-synonymous variations and 14 synonymous substitutions. Also, the 1.66:1 ratio of transition/transversion is different from that of other insects. As the first SNP analysis of a Lepidoptera, B. mori, the single nucleotide polymorphic density is estimated to be 1.3 x 10(-3) by sequence diversity. This analysis shows that expressed sequences from multiple libraries may provide an abundant source of comparative reads to mine for cSNPs from the silkworm genome.

Genome coverage and sequence fidelity of phi29 polymerase-based multiple strand displacement whole genome amplification

Major efforts are underway to systematically define the somatic and germline genetic variations causally associated with disease. Genome-wide genetic analysis of actual clinical samples is, however, limited by the paucity of genomic DNA available. Here we have tested the fidelity and genome representation of phi29 polymerase-based genome amplification (phi29MDA) using direct sequencing and high density oligonucleotide arrays probing >10,000 SNP alleles. Genome representation was comprehensive and estimated to be 99.82% complete, although six regions encompassing a maximum of 5.62 Mb failed to amplify. There was no degradation in the accuracy of SNP genotyping and, in direct sequencing experiments sampling 500,000 bp, the estimated error rate (9.5 x 10(-6)) was the same as in paired unamplified samples. The detection of cancer-associated loss of heterozygosity and copy number changes, including homozygous deletion and gene amplification, were similarly robust. These results suggest that phi29MDA yields high fidelity, near-complete genome representation suitable for high resolution genetic analysis.