Genotyping for human whole-genome scans: past, present, and future

Development of a 17-Plex of Penta- and Tetra-Nucleotide Microsatellites for DNA Profiling and Paternity Testing in Horses

Tetranucleotide and pentanucleotide short tandem repeat (hereafter termed tetraSTR and pentaSTR) polymorphisms have properties that make them desirable for DNA profiling and paternity testing. However, certain species, such as the horse, have far fewer tetraSTRs than other species and for this reason dinucleotide STRs (diSTRs) have become the standard for DNA profiling in horses, despite being less desirable for technical reasons. During our testing of a series of candidate genes as potentially underlying a heritable condition characterized by megaesophagus in the Friesian horse breed, we found that good tetraSTRs do exist in horses but, as expected, at a much lower frequency than in other species, e.g., dogs and humans. Using a series of efficient methods developed in our laboratory for the production of multiplexed tetraSTRs in other species, we identified a set of tetra- and pentaSTRs that we developed into a 17-plex panel for the horse, plus a sex-identifying marker near the amelogenin gene. These markers were tested in 128 horses representing 16 breeds as well as crossbred horses, and we found that these markers have useful genetic variability. Average observed heterozygosities (Ho) ranged from 0.53 to 0.89 for the individual markers (0.66 average Ho for all markers), and 0.62-0.82 for expected heterozygosity (He) within breeds (0.72 average He for all markers). The probability of identity (PI) within breeds for which 10 or more samples were available was at least 1.1 x 10⁻¹¹, and the PI among siblings (PIsib) was 1.5 x 10⁻⁵. Stutter was ≤ 11% (average stutter for all markers combined was 6.9%) compared to the more than 30% typically seen with diSTRs. We predict that it will be possible to develop accurate allelic ladders for this multiplex panel that will make cross-laboratory comparisons easier and will also improve DNA profiling accuracy. Although we were only able to exclude candidate genes for Friesian horse megaesophagus with no unexcluded genes that are possibly causative at this point in time, the study helped us to refine the methods used to develop better tetraSTR multiplexed panels for species such as the horse that have a low frequency of tetraSTRs.

A 13-plex of tetra- and penta-STRs to identify zebrafish

The zebrafish species Danio rerio has become one of the major vertebrate model organisms used in biomedical research. However, there are aspects of the model that need to be improved. One of these is the ability to identify individual fish and fish lines by DNA profiling. Although many dinucleotide short tandem repeat (diSTR) markers are available for this and similar purposes, they have certain disadvantages such as an excessive polymerase slippage ("stutter") that causes difficulties in automated genotyping and cross-laboratory comparisons. Here we report on the development of a 13-plex of tetranucleotide and pentanucleotide STRs (tetraSTRs and pentaSTRs, respectively) that have low stutter. The system uses an inexpensive universal primer labelling system, which can easily be converted to a direct labeling system if desired. This 13-plex was examined in three zebrafish lines (NHGRI-1, kca33Tg, and kca66Tg, originally obtained from ZIRC). The average observed heterozygosity (Ho) and expected heterozygosity (He) in these highly inbred lines were 0.291 and 0.359, respectively, which is very similar to what has been found with diSTRs. The probability of identity (PI) for all fish tested was 2.1 × 10^-5 and the PI for siblings (PIsib) was 6.4 × 10^-3, as calculated by the Genalex package. Ninety percent of the fish tested were correctly identified with their respective strains. It is also demonstrated that this panel can be used to confirm doubled-haploid cell lines. This multiplex should find multiple uses for improving the accuracy and reproducibility of studies using the zebrafish model.

MultiplexSSR: A pipeline for developing multiplex SSR-PCR assays from resequencing data

Next-generation sequencing has greatly promoted the investigation of single nucleotide polymorphisms, while studies of simple sequence repeats are sharply decreasing. However, simple sequence repeats still present some advantages in conservation genetics. In this study, an end-to-end pipeline referred to as MultiplexSSR was established to develop multiplex PCR assays in batches with highly polymorphic simple sequence repeats for capillary platforms from resequencing data. The distribution of single sequence repeats in the genome, the error profiles of genotypes and allelotypes, and the increase in the allele length range depending on the number of individuals were investigated. A total of 98% of single sequence repeats presented lengths of less than 100 bp. The error rate of the genotyping and allelotyping of dimeric patterns was ten times higher than those for other patterns. The error rate of allelotyping was less than that of genotyping. The allele length range reached approximate saturation with 10 individuals. This pipeline uses allele numbers to select highly polymorphic loci, masks loci with variation, and applies in silico PCR to improve primer specificity. The application of the developed multiplex SSR-PCR assays validated the pipeline's robustness, showing higher polymorphism and stability for the developed simple sequence repeats and a lower cost for genotyping and providing low-depth resequencing data from less than a dozen individuals for the development of markers. This pipeline fills the gap between next-generation sequencing and multiplex SSR-PCR.

Identification and characterization of the highly polymorphic locus D14S739 in the Han Chinese population

AIM

To systemically select and evaluate short tandem repeats (STRs) on the chromosome 14 and obtain new STR loci as expanded genotyping markers for forensic application.

METHODS

STRs on the chromosome 14 were filtered from Tandem Repeats Database and further selected based on their positions on the chromosome, repeat patterns of the core sequences, sequence homology of the flanking regions, and suitability of flanking regions in primer design. The STR locus with the highest heterozygosity and polymorphism information content (PIC) was selected for further analysis of genetic polymorphism, forensic parameters, and the core sequence.

RESULTS

Among 26 STR loci selected as candidates, D14S739 had the highest heterozygosity (0.8691) and PIC (0.8432), and showed no deviation from the Hardy-Weinberg equilibrium. 14 alleles were observed, ranging in size from 21 to 34 tetranucleotide units in the core region of (GATA)9-18 (GACA)7-12 GACG (GACA)2 GATA. Paternity testing showed no mutations.

CONCLUSION

D14S739 is a highly informative STR locus and could be a suitable genetic marker for forensic applications in the Han Chinese population.

Abundant contribution of short tandem repeats to gene expression variation in humans

The contribution of repetitive elements to quantitative human traits is largely unknown. Here we report a genome-wide survey of the contribution of short tandem repeats (STRs), which constitute one of the most polymorphic and abundant repeat classes, to gene expression in humans. Our survey identified 2,060 significant expression STRs (eSTRs). These eSTRs were replicable in orthogonal populations and expression assays. We used variance partitioning to disentangle the contribution of eSTRs from that of linked SNPs and indels and found that eSTRs contribute 10-15% of the cis heritability mediated by all common variants. Further functional genomic analyses showed that eSTRs are enriched in conserved regions, colocalize with regulatory elements and may modulate certain histone modifications. By analyzing known genome-wide association study (GWAS) signals and searching for new associations in 1,685 whole genomes from deeply phenotyped individuals, we found that eSTRs are enriched in various clinically relevant conditions. These results highlight the contribution of STRs to the genetic architecture of quantitative human traits.

Advances and applications of molecular cloning in clinical microbiology

Molecular cloning is based on isolation of a DNA sequence of interest to obtain multiple copies of it in vitro. Application of this technique has become an increasingly important tool in clinical microbiology due to its simplicity, cost effectiveness, rapidity, and reliability. This review entails the recent advances in molecular cloning and its application in the clinical microbiology in the context of polymicrobial infections, recombinant antigens, recombinant vaccines, diagnostic probes, antimicrobial peptides, and recombinant cytokines. Culture-based methods in polymicrobial infection have many limitation, which has been overcome by cloning techniques and provide gold standard technique. Recombinant antigens produced by cloning technique are now being used for screening of HIV, HCV, HBV, CMV, Treponema pallidum, and other clinical infectious agents. Recombinant vaccines for hepatitis B, cholera, influenza A, and other diseases also use recombinant antigens which have replaced the use of live vaccines and thus reduce the risk for adverse effects. Gene probes developed by gene cloning have many applications including in early diagnosis of hereditary diseases, forensic investigations, and routine diagnosis. Industrial application of this technology produces new antibiotics in the form of antimicrobial peptides and recombinant cytokines that can be used as therapeutic agents.

Genome-wide linkage scan identifies two novel genetic loci for coronary artery disease: in GeneQuest families

Coronary artery disease (CAD) is the leading cause of death worldwide. Recent genome-wide association studies (GWAS) identified >50 common variants associated with CAD or its complication myocardial infarction (MI), but collectively they account for <20% of heritability, generating a phenomena of “missing heritability”. Rare variants with large effects may account for a large portion of missing heritability. Genome-wide linkage studies of large families and follow-up fine mapping and deep sequencing are particularly effective in identifying rare variants with large effects. Here we show results from a genome-wide linkage scan for CAD in multiplex GeneQuest families with early onset CAD and MI. Whole genome genotyping was carried out with 408 markers that span the human genome by every 10 cM and linkage analyses were performed using the affected relative pair analysis implemented in GENEHUNTER. Affected only nonparametric linkage (NPL) analysis identified two novel CAD loci with highly significant evidence of linkage on chromosome 3p25.1 (peak NPL  = 5.49) and 3q29 (NPL  = 6.84). We also identified four loci with suggestive linkage on 9q22.33, 9q34.11, 17p12, and 21q22.3 (NPL  = 3.18–4.07). These results identify novel loci for CAD and provide a framework for fine mapping and deep sequencing to identify new susceptibility genes and novel variants associated with risk of CAD.

Value of Mendelian laws of segregation in families: data quality control, imputation, and beyond

When analyzing family data, we dream of perfectly informative data, even whole-genome sequences (WGSs) for all family members. Reality intervenes, and we find that next-generation sequencing (NGS) data have errors and are often too expensive or impossible to collect on everyone. The Genetic Analysis Workshop 18 working groups on quality control and dropping WGSs through families using a genome-wide association framework focused on finding, correcting, and using errors within the available sequence and family data, developing methods to infer and analyze missing sequence data among relatives, and testing for linkage and association with simulated blood pressure. We found that single-nucleotide polymorphisms, NGS data, and imputed data are generally concordant but that errors are particularly likely at rare variants, for homozygous genotypes, within regions with repeated sequences or structural variants, and within sequence data imputed from unrelated individuals. Admixture complicated identification of cryptic relatedness, but information from Mendelian transmission improved error detection and provided an estimate of the de novo mutation rate. Computationally, fast rule-based imputation was accurate but could not cover as many loci or subjects as more computationally demanding probability-based methods. Incorporating population-level data into pedigree-based imputation methods improved results. Observed data outperformed imputed data in association testing, but imputed data were also useful. We discuss the strengths and weaknesses of existing methods and suggest possible future directions, such as improving communication between data collectors and data analysts, establishing thresholds for and improving imputation quality, and incorporating error into imputation and analytical models.

Microsatellite data show recent demographic expansions in sedentary but not in nomadic human populations in Africa and Eurasia

The transition from hunting and gathering to plant and animal domestication was one of the most important cultural and technological revolutions in human history. According to archeologists and paleoanthropologists, this transition triggered major demographic expansions. However, few genetic studies have found traces of Neolithic expansions in the current repartition of genetic polymorphism, pointing rather toward Paleolithic expansions. Here, we used microsatellite autosomal data to investigate the past demographic history of 87 African and Eurasian human populations with contrasted lifestyles (nomadic hunter-gatherers, semi-nomadic herders and sedentary farmers). Likely due to the combination of a higher mutation rate and the possibility to analyze several loci as independent replicates of the coalescent process, the analysis of microsatellite data allowed us to infer more recent expansions than previous genetic studies, potentially resulting from the Neolithic transition. Despite the variability in their location and environment, we found consistent expansions for all sedentary farmers, while we inferred constant population sizes for all hunter-gatherers and most herders that could result from constraints linked to a nomadic or semi-nomadic lifestyle and/or competition for land between herders and farmers. As an exception, we inferred expansions for Central Asian herders. This might be linked with the arid environment of this area that may have been more favorable to nomadic herders than to sedentary farmers. Alternatively, current Central Asian herders may descent from populations who have first experienced a transition from hunter-gathering to sedentary agropastoralism, and then a second transition to nomadic herding.

Next-generation phenotyping: requirements and strategies for enhancing our understanding of genotype-phenotype relationships and its relevance to crop improvement

More accurate and precise phenotyping strategies are necessary to empower high-resolution linkage mapping and genome-wide association studies and for training genomic selection models in plant improvement. Within this framework, the objective of modern phenotyping is to increase the accuracy, precision and throughput of phenotypic estimation at all levels of biological organization while reducing costs and minimizing labor through automation, remote sensing, improved data integration and experimental design. Much like the efforts to optimize genotyping during the 1980s and 1990s, designing effective phenotyping initiatives today requires multi-faceted collaborations between biologists, computer scientists, statisticians and engineers. Robust phenotyping systems are needed to characterize the full suite of genetic factors that contribute to quantitative phenotypic variation across cells, organs and tissues, developmental stages, years, environments, species and research programs. Next-generation phenotyping generates significantly more data than previously and requires novel data management, access and storage systems, increased use of ontologies to facilitate data integration, and new statistical tools for enhancing experimental design and extracting biologically meaningful signal from environmental and experimental noise. To ensure relevance, the implementation of efficient and informative phenotyping experiments also requires familiarity with diverse germplasm resources, population structures, and target populations of environments. Today, phenotyping is quickly emerging as the major operational bottleneck limiting the power of genetic analysis and genomic prediction. The challenge for the next generation of quantitative geneticists and plant breeders is not only to understand the genetic basis of complex trait variation, but also to use that knowledge to efficiently synthesize twenty-first century crop varieties.

Linkage analysis of a large African family segregating stuttering suggests polygenic inheritance

We describe a pedigree of 71 individuals from the Republic of Cameroon in which at least 33 individuals have a clinical diagnosis of persistent stuttering. The high concentration of stuttering individuals suggests that the pedigree either contains a single highly penetrant gene variant or that assortative mating led to multiple stuttering-associated variants being transmitted in different parts of the pedigree. No single locus displayed significant linkage to stuttering in initial genome-wide scans with microsatellite and SNP markers. By dividing the pedigree into five subpedigrees, we found evidence for linkage to previously reported loci on 3q and 15q, and to novel loci on 2p, 3p, 14q, and a different region of 15q. Using the two-locus mode of Superlink, we showed that combining the recessive locus on 2p and a single-locus additive representation of the 15q loci is sufficient to achieve a two-locus score over 6 on the entire pedigree. For this 2p + 15q analysis, we show LOD scores ranging from 4.69 to 6.57, and the scores are sensitive to which marker is chosen for 15q. Our findings provide strong evidence for linkage at several loci.

DNAJC6 is responsible for juvenile parkinsonism with phenotypic variability

Familial parkinson's disease is both clinically and genetically heterogeneous. By mapping the disease locus with a lod score of 5.13 to a < 3.5 Mbp region at 1p31.3 in a consanguineous family and subsequent exome sequencing analysis, we identified homozygous truncating mutation p.Q734X in DNAJC6. Four members of the family were afflicted with juvenile parkinsonism that presented with mental retardation, pyramidal signs and epilepsy, as well as varying degrees of a progressive neurological disease. Recently a splicing mutation in the same gene was reported in two brothers with juvenile parkinsonism that was not L-Dopa responsive and not accompanied by pyramidal signs or mental retardation. Also, an 80-kb deletion that included DNAJC6 sequences was identified in a boy reported as having obesity, epilepsy and mental retardation but not any signs of parkinsonism. The phenotype of our study family resembles both of those families, which among themselves do not share any clinical features. Our findings further establish DNAJC6 as a juvenile parkinsonism gene, and expand the spectrums of the parkinsonism phenotype and DNAJC6 mutation. DNAJC6 encodes the neuronal co-chaperone auxilin. We found that its transcript is highly significantly more abundant in brain as compared to the non-neural tissues assayed.

Human neutral genetic variation and forensic STR data

The forensic genetics field is generating extensive population data on polymorphism of short tandem repeats (STR) markers in globally distributed samples. In this study we explored and quantified the informative power of these datasets to address issues related to human evolution and diversity, by using two online resources: an allele frequency dataset representing 141 populations summing up to almost 26 thousand individuals; a genotype dataset consisting of 42 populations and more than 11 thousand individuals. We show that the genetic relationships between populations based on forensic STRs are best explained by geography, as observed when analysing other worldwide datasets generated specifically to study human diversity. However, the global level of genetic differentiation between populations (as measured by a fixation index) is about half the value estimated with those other datasets, which contain a much higher number of markers but much less individuals. We suggest that the main factor explaining this difference is an ascertainment bias in forensics data resulting from the choice of markers for individual identification. We show that this choice results in average low variance of heterozygosity across world regions, and hence in low differentiation among populations. Thus, the forensic genetic markers currently produced for the purpose of individual assignment and identification allow the detection of the patterns of neutral genetic structure that characterize the human population but they do underestimate the levels of this genetic structure compared to the datasets of STRs (or other kinds of markers) generated specifically to study the diversity of human populations.

A direct characterization of human mutation based on microsatellites

Mutations are the raw material of evolution, but have been difficult to study directly. We report the largest study of new mutations to date: 2,058 germline changes discovered by analyzing 85,289 Icelanders at 2,477 microsatellites. The paternal-to-maternal mutation rate ratio is 3.3, and the rate in fathers doubles from age 20 to 58 whereas there is no association with age in mothers. Longer microsatellite alleles are more mutagenic and tend to decrease in length, whereas the opposite is seen for shorter alleles. We use these empirical observations to build a model that we apply to individuals for whom we have both genome sequence and microsatellite data, allowing us to estimate key parameters of evolution without calibration to the fossil record. We infer that the sequence mutation rate is 1.4–2.3×10⁻⁸ per base pair per generation (90% credible interval), and that human-chimpanzee speciation occurred 3.7–6.6 million years ago.

lobSTR: A short tandem repeat profiler for personal genomes

Short tandem repeats (STRs) have a wide range of applications, including medical genetics, forensics, and genetic genealogy. High-throughput sequencing (HTS) has the potential to profile hundreds of thousands of STR loci. However, mainstream bioinformatics pipelines are inadequate for the task. These pipelines treat STR mapping as gapped alignment, which results in cumbersome processing times and a biased sampling of STR alleles. Here, we present lobSTR, a novel method for profiling STRs in personal genomes. lobSTR harnesses concepts from signal processing and statistical learning to avoid gapped alignment and to address the specific noise patterns in STR calling. The speed and reliability of lobSTR exceed the performance of current mainstream algorithms for STR profiling. We validated lobSTR's accuracy by measuring its consistency in calling STRs from whole-genome sequencing of two biological replicates from the same individual, by tracing Mendelian inheritance patterns in STR alleles in whole-genome sequencing of a HapMap trio, and by comparing lobSTR results to traditional molecular techniques. Encouraged by the speed and accuracy of lobSTR, we used the algorithm to conduct a comprehensive survey of STR variations in a deeply sequenced personal genome. We traced the mutation dynamics of close to 100,000 STR loci and observed more than 50,000 STR variations in a single genome. lobSTR's implementation is an end-to-end solution. The package accepts raw sequencing reads and provides the user with the genotyping results. It is written in C/C++, includes multi-threading capabilities, and is compatible with the BAM format.

Microsatellite markers for linkage and association studies

During the 1990s and the first several years of this century, microsatellites or short tandem repeats were the workhorse genetic markers for hypothesis-independent studies in human genetics, facilitating genome-wide linkage studies and allelic imbalance studies. However, the rise of higher throughput and cost-effective single-nucleotide polymorphism (SNP) platforms led to the era of the SNP for genome scans. Nevertheless, it is important to note that microsatellites remain highly informative and useful measures of genomic variation for linkage and association studies. Their continued advantage in complementing SNPs lies in their greater allelic diversity than biallelic SNPs as well as in their population history, in which single-step expansion or contraction of the tandem repeat on the background of ancestral SNP haplotypes can break up common haplotypes, leading to greater haplotype diversity within the linkage disequilibrium block of interest. In fact, microsatellites have starred in association studies leading to widely replicated discoveries of type 2 diabetes (TCF7L2) and prostate cancer genes (the 8q21 region). At the end of the day, it will be important to catalog all variation, including SNPs, microsatellites, copy number variations, and polymorphic inversions in human genetic studies. This article describes the utilities of microsatellites and experimental approaches in their use.

Genome scan of clot lysis time and its association with thrombosis in a protein C-deficient kindred

BACKGROUND

 Previously, we found increased clot-lysis time (CLT), as measured with a plasma-based assay, to increase the risk of venous thrombosis in two population-based case-control studies. The genes influencing CLT are as yet unknown.

PATIENTS/METHODS

 We tested CLT as risk factor for venous thrombosis in Kindred Vermont II (n = 346), a pedigree suffering from a high thrombosis risk, partially attributable to a type I protein C deficiency. Furthermore, we tested for quantitative trait loci (QTLs) for CLT, using variance component linkage analysis.

RESULTS

 Protein C-deficient family members had shorter CLTs than non-deficient members (median CLT 67 min vs. 75 min). One standard deviation increase in CLT increased the risk of venous thrombosis 2.4-fold in non-deficient family members. Protein C deficiency without elevated CLT increased the risk 6.9-fold. Combining both risk factors yielded a 27.8-fold increased risk. The heritability of CLT was 42-52%. We found suggestive evidence of linkage on chromosome 11 (62 cM), partly explained by the prothrombin 20210A mutation, and on chromosome 13 (52 cM). Thrombin-activatable fibrinolysis inhibitor genotypes did not explain the variation in CLT.

CONCLUSION

Hypofibrinolysis appears to increase thrombosis risk in this family, especially in combination with protein C deficiency. Protein C deficiency is associated with short CLT. CLT is partly genetically regulated. Suggestive QTLs were found on chromosomes 11 and 13.

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.

Early-onset progressive myoclonic epilepsy with dystonia mapping to 16pter-p13.3

The authors present three patients from a consanguineous family afflicted with novel recessive myoclonic epilepsy characterized by very early onset and a steadily progressive course. The onset is in early infancy, and death occurs in the first decade. In addition to various types of myoclonic seizures, episodic phenomena such as dystonias, postictal enduring hemipareses, autonomic involvements, and periods of obtundation and lethargy were also observed. Developmental and neurological retardation, coupled with systemic infections, leads to a full deterioration. The authors designated the disease progressive myoclonic epilepsy with dystonia (PMED). A genome scan for the family and subsequent fine mapping localized the gene responsible for the disease to the most telomeric 6.73 mega base pairs at the p-terminus of chromosome 16, with a maximum multipoint logarithm-of-odds score of 7.83 and a maximum two-point score of 4.25. A candidate gene was analyzed for mutations in patients, but no mutation was found.

A frameshift mutation of ERLIN2 in recessive intellectual disability, motor dysfunction and multiple joint contractures

We present a family afflicted with a novel autosomal recessive disease characterized by progressive intellectual disability, motor dysfunction and multiple joint contractures. No pathology was found by cranial imaging, electromyography and muscle biopsy, but electron microscopy in leukocytes revealed large vacuoles containing flocculent material. We mapped the disease gene by SNP genome scan and linkage analysis to an ∼0.80 cM and 1 Mb region at 8p11.23 with a multipoint logarithm of odds (LOD) score of 12. By candidate gene approach, we identified a homozygous two-nucleotide insertion in ERLIN2, predicted to lead to the truncation of the protein by about 20%. The gene encodes endoplasmic reticulum (ER) lipid raft-associated protein 2 that mediates the ER-associated degradation of activated inositol 1,4,5-trisphosphate receptors and other substrates.

The phenotype caused by PYCR1 mutations corresponds to geroderma osteodysplasticum rather than autosomal recessive cutis laxa type 2

Geroderma osteodysplasticum is a rare autosomal recessive disorder characterized by wrinkled skin on the dorsum of the hands and feet, osteopenia, prognathism, and an elongated and lax face. The mutated gene was identified as GORAB (SCYL1BP1). As well, the PYCR1 gene also was shown to be mutated in a similar disease, designated cutis laxa, autosomal recessive, type IIB (ARCL2B) or cutis laxa with progeroid features. We describe here the clinical findings in four affected individuals in a family with geroderma osteodysplasticum with mental retardation and a homozygous mutation in PYCR1. Although the disease resulting from recessive mutations in that gene has been recently designated ARCL2B, some clinical features, such as prognathism, elongated and lax face, osteopenia and limitation of skin wrinkling to the dorsum of hands and feet, in the patients reported here as well as in others reported with PYCR1 mutations, are generally more common in geroderma osteodysplasticum resulting from recessive GORAB mutations. While the patients with GORAB mutations have severe osteopenia, the patients with PYCR1 mutations have severe mental retardation. In conclusion, the phenotype caused by PYCR1 mutations corresponds to geroderma osteodysplasticum rather than ARCL2B.

A novel recessive GUCY2D mutation causing cone-rod dystrophy and not Leber's congenital amaurosis

Cone-rod dystrophies are inherited retinal dystrophies that are characterized by progressive degeneration of cones and rods, causing an early decrease in central visual acuity and colour vision defects, followed by loss of peripheral vision in adolescence or early adult life. Both genetic and clinical heterogeneity are well known. In a family with autosomal recessive cone-rod dystrophy, genetic analyses comprising genome scan with microsatellite markers, fine mapping and candidate gene approach resulted in the identification of a homozygous missense GUCY2D mutation. This is the first GUCY2D mutation associated with autosomal recessive cone-rod dystrophy rather than Leber's congenital amaurosis (LCA), a severe disease leading to childhood blindness. This study hence establishes GUCY2D, which is a common cause for both recessive LCA and dominant cone-rod dystrophy, as a good candidate for autosomal recessive cone-rod dystrophy.

Two-dimensional, sex-specific autosomal linkage scan of the number of sodium pump sites

OBJECTIVES

The sodium pump consists of the membrane-bound enzyme sodium/potassium-ATPase, which exchanges internal sodium ions for external potassium ions. Obesity, hypertension, and diabetes associate with the activity of the sodium pump, motivating gene discovery for sodium pump number.

METHODS

Variance components linkage analysis was applied to the number of red blood cell sodium pump sites measured by ouabain-binding assays on 1375 members of 46 Utah pedigrees. Both one-dimensional (1D) and two-dimensional (2D) autosome-wide linkage analyses of pump number were performed on the combined sample as well as separately on the male and female subsets.

RESULTS

Two significant 1D linkages were identified: on chromosome 1p13 in the combined sample [1D logarithm of odds (LOD) score = 3.76] and on chromosome 17p21 in the female subset (1D LOD score = 3.24). In addition, two significant 2D linkages were identified in the female subset: on chromosome 10q22 interacting with chromosome 18q11 (2D LOD score = 7.18) and on chromosome 13q21 interacting with chromosome 4q31 (2D LOD score = 6.05). Single-nucleotide polymorphism rs17376826 in neuropeptide Y receptor Y2, an obesity-associated gene and a candidate in the chromosome 4q31 linkage region, is associated with pump number (P = 0.046 in the combined sample and P = 0.042 in the female subset).

CONCLUSION

Pump number is influenced by multiple genes, possibly including neuropeptide Y receptor Y2.

Sequence determinants of human microsatellite variability

BACKGROUND

Microsatellite loci are frequently used in genomic studies of DNA sequence repeats and in population studies of genetic variability. To investigate the effect of sequence properties of microsatellites on their level of variability we have analyzed genotypes at 627 microsatellite loci in 1,048 worldwide individuals from the HGDP-CEPH cell line panel together with the DNA sequences of these microsatellites in the human RefSeq database.

RESULTS

Calibrating PCR fragment lengths in individual genotypes by using the RefSeq sequence enabled us to infer repeat number in the HGDP-CEPH dataset and to calculate the mean number of repeats (as opposed to the mean PCR fragment length), under the assumption that differences in PCR fragment length reflect differences in the numbers of repeats in the embedded repeat sequences. We find the mean and maximum numbers of repeats across individuals to be positively correlated with heterozygosity. The size and composition of the repeat unit of a microsatellite are also important factors in predicting heterozygosity, with tetra-nucleotide repeat units high in G/C content leading to higher heterozygosity. Finally, we find that microsatellites containing more separate sets of repeated motifs generally have higher heterozygosity.

CONCLUSIONS

These results suggest that sequence properties of microsatellites have a significant impact in determining the features of human microsatellite variability.

A novel locus for adolescent idiopathic scoliosis on chromosome 12p

Adolescent idiopathic scoliosis (AIS) is a common disorder with strong evidence for genetic predisposition. Quantitative trait loci (QTLs) for AIS susceptibility have been identified on chromosomes. We performed a genome-wide genetic linkage scan in seven multiplex families using 400 marker loci with a mean spacing of 8.6 cM. We used Genehunter Plus to generate linkage statistics, expressed as homogeneity (HLOD) scores, under dominant and recessive genetic models. We found a significant linkage signal on chromosome 12p, whose support interval extends from near 12 pter, spanning approximately 10 million bases or 31 cM. Fine mapping within the region using 20 additional markers reveals maximum HLOD = 3.7 at 5 cM under a dominant inheritance model, and a split peak maximum HLOD = 3.2 at 8 and 18 cM under a recessive inheritance model. The linkage support interval contains 95 known genes. We found evidence suggestive of linkage on chromosomes 1, 6, 7, 8, and 14. This study is the first to find evidence of an AIS susceptibility locus on chromosome 12. Detection of AIS susceptibility QTLs on multiple chromosomes in this and other studies demonstrate that the condition is genetically heterogeneous.

Autosomal recessive spastic paraplegia (SPG45) with mental retardation maps to 10q24.3-q25.1

Hereditary spastic paraplegias (HSPs) are characterized by progressive spasticity in the lower limbs. They are clinically heterogeneous, and pure forms as well as complicated forms with other accompanying clinical findings are known. HSPs are also genetically heterogeneous. We performed clinical and genetic studies in a consanguineous family with five affected members. A genome scan using 405 microsatellite markers for eight members of the family identified candidate gene loci, and subsequent fine mapping in 16 members identified the gene locus responsible for the HSP. The clinical manifestations were very early onset spastic paraplegia (SPG) accompanied by mental retardation and ocular signs. The gene locus was identified as the interval 102.05-106.64 Mbp on chromosome 10. Gene MRPL43 was analyzed in the patients. No mutation but high levels of mRNA were detected. We have mapped a novel autosomal recessive complicated form of HSP (SPG45) to a 4.6-Mbp region at 10q24.3-q25.1 with multipoint logarithm of odds scores >4.5.

The impact of data quality on the identification of complex disease genes: experience from the Family Blood Pressure Program

The application of genome-wide linkage scans to uncover susceptibility loci for complex diseases offers great promise for the risk assessment, treatment, and understanding of these diseases. However, for most published studies, linkage signals are typically modest and vary considerably from one study to another. The multicenter Family Blood Pressure Program has analyzed genome-wide linkage scans of over 12 000 individuals. Based on this experience, we developed a protocol for large linkage studies that reduces two sources of data error: pedigree structure and marker genotyping errors. We then used the linkage signals, before and after data cleaning, to illustrate the impact of missing and erroneous data. A comprehensive error-checking protocol is an important part of complex disease linkage studies and enhances gene mapping. The lack of significant and reproducible linkage findings across studies is, in part, due to data quality.

A genomewide comparison of population structure at STRPs and nearby SNPs in humans

Patterns of population structure provide insights into evolutionary processes and help identify groups of individuals for genotype-phenotype association studies. With increasing availability of polymorphic molecular markers across genomes, the examination of population structure using large numbers of unlinked loci has become a common practice in evolutionary biology and human genetics. The two classes of molecular variation most widely used for this purpose, short tandem repeat polymorphisms (STRPs) and single-nucleotide polymorphisms (SNPs), differ in mutational properties expected to affect population structure. To measure the relative ability of these loci to describe population structure, we compared diversity at neighboring STRPs and SNPs from 720 genomic regions in the four populations that comprise the Human HapMap. Comparing loci from the same genomic regions allowed us to focus on the contribution of mutational differences (rather than variation in genealogical history) to disparities in population structure between STRPs and SNPs. Relative to average values for SNPs from the same regions, STRPs had lower F(st), but higher G(st)' and I(n) values. STRP-SNP correlations in population structure across genomic regions were statistically significant but weak in magnitude. Separate analyses by repeat type showed that these correlations were driven primarily by tetranucleotide and trinucleotide STRPs; measures of population structure at dinucleotides and SNPs were not significantly correlated. Pairwise comparisons among populations revealed effects of divergence time on differences in population structure between STRPs and SNPs. Collectively, these results confirm that individual STRPs can provide more information about population structure than individual SNPs, but suggest that the difference in structure at STRPs and SNPs depends on local genealogical history. Our study motivates theoretical comparisons of population structure at loci with different mutational properties.

Microsatellites are molecular clocks that support accurate inferences about history

Microsatellite length mutations are often modeled using the generalized stepwise mutation process, which is a type of random walk. If this model is sufficiently accurate, one can estimate the coalescence time between alleles of a locus after a mathematical transformation of the allele lengths. When large-scale microsatellite genotyping first became possible, there was substantial interest in using this approach to make inferences about time and demography, but that interest has waned because it has not been possible to empirically validate the clock by comparing it with data in which the mutation process is well understood. We analyzed data from 783 microsatellite loci in human populations and 292 loci in chimpanzee populations, and compared them with up to one gigabase of aligned sequence data, where the molecular clock based upon nucleotide substitutions is believed to be reliable. We empirically demonstrate a remarkable linearity (r(2) > 0.95) between the microsatellite average square distance statistic and sequence divergence. We demonstrate that microsatellites are accurate molecular clocks for coalescent times of at least 2 million years (My). We apply this insight to confirm that the African populations San, Biaka Pygmy, and Mbuti Pygmy have the deepest coalescent times among populations in the Human Genome Diversity Project. Furthermore, we show that microsatellites support unbiased estimates of population differentiation (F(ST)) that are less subject to ascertainment bias than single nucleotide polymorphism (SNP) F(ST). These results raise the prospect of using microsatellite data sets to determine parameters of population history. When genotyped along with SNPs, microsatellite data can also be used to correct for SNP ascertainment bias.

Genetic Effect on Blood Pressure Is Modulated by Age

Genome-wide linkage analysis was performed for systolic and diastolic blood pressures in the Hypertension Genetic Epidemiology Network. We investigated the role of gene–age interactions using a recently developed variance components method that incorporates age variation in genetic effects. Substantially improved linkage evidence, in terms of both the number of linkage peaks and their significance levels, was observed. Twenty-six linkage peaks were identified with maximum logarithm of odds scores ranging between 3.0 and 4.6, 15 of which were cross-validated by the literature. The chromosomal region 1p36 that showed the highest logarithm of odds score in our study was found to be supported by evidence from 3 studies. The new method also led to vastly improved validation across ethnic groups. Ten of the 15 supported linkage peaks were cross-validated between 2 different ethnic groups, and 2 peaks on chromosomal region 1q31 and 16p11 were validated in 3 ethnic groups. In conclusion, this investigation demonstrates that genetic effects on blood pressure vary by age. The improved genetic linkage results presented here should help to identify the specific genetic variants that explain the observed results.

Air pollutants, oxidative stress and human health

Air pollutants have, and continue to be, major contributing factors to chronic diseases and mortality, subsequently impacting public health. Chronic diseases include: chronic obstructive pulmonary diseases (COPD), cardiovascular diseases (CVD), asthma, and cancer. Byproducts of oxidative stress found in air pollutants are common initiators or promoters of the damage produced in such chronic diseases. Such air pollutants include: ozone, sulfur oxides, carbon monoxide, nitrogen oxides, and particulate matter. Interaction between oxidative stress byproducts and certain genes within our population may modulate the expression of specific chronic diseases. In this brief review we attempt to provide some insight into what we currently know about the health problems associated with various air pollutants and their relationship in promoting chronic diseases through changes in oxidative stress and modulation of gene expression. Such insight eventually may direct the means for effective public health prevention and treatment of diseases associated with air pollution and treatment of diseases associated with air pollution.

Predisposition of genetic polymorphism with the risk of urolithiasis

Urolithiasis is a relevant clinical problem with a subsequent burden for health system. The aim of this review is to provide recent progress made using genetic polymorphisms to define pathophysiology, to identify persons at risk for kidney stone disease and to predict treatment response. Population case-control studies are useful both as an alternative and an adjunct as compared to family studies. These involve either whole genome scanning or candidate gene approaches. While whole genome scanning is likely to be widely used in future, at present, candidate gene studies are more feasible. When performing candidate gene case-control studies factors such as study design, methods for recruitment of case and controls, selection of candidate genes, functional significance of polymorphisms chosen for study and statistical analysis require close attention to ensure that only genuine associations are detected. Some of the significant genes that play role in stone formation include calcitonin receptor gene (CTR), vitamin D receptor (VDR), Urokinase, Interleukin, (IL-1β, IL-Ra), E-Cadherin, Androgen & oestrogen receptor gene, vascular endothelial growth factor (VEGF) and Arginine p21. In our case-control study we studied CTR, VDR, Urokinase, IL-1β(-511 and +3954), IL-Ra from north India and predict that VDR, IL-β (-511) and IL-1Ra gene may be used as a possible genetic marker for earlier detection in patients who are at risk for calcium oxalate stone disease. Further, linkage disequilibrium and haplotype structure of a certain candidate gene is important for association analysis. When a certain polymorphic allele has been found to be associated with disease, it is further explained on basis of LD and haplotype structure by one or more other alleles. Once it is determined which haplotype carries the risk allele, by means of molecular biological functional analyses, the variants on that haplotype allele truly causing the effect can be determined.

A first stage genome-wide screen for regions shared identical-by-descent in Hutterite families with multiple sclerosis

The complexity of multiple sclerosis (MS) genetics has made the search for novel genes using traditional sharing methods problematic. In order to minimize the genetic heterogeneity present in the MS population we have screened the Canadian MS population for individuals belonging to the Hutterite Brethren. Seven Hutterites with clinically definite MS were ascertained and are related to a common founder by eight generations. Six of the 7 affected individuals and 21 of their unaffected family members (total = 27) were genotyped for 807 markers. Haplotypes were then inspected for sharing among the six MS patients. There were three haplotypes shared among all six MS patients. The haplotypes were located at 2q34-35, 4q31-32, and 17p13. An additional 15 haplotypes were shared among five of the six Hutterites MS patients. The HLA Class II region was one of the highlighted regions; however, the shared MHC haplotype bore the DRB1*04 allele and not the MS-associated DRB1*15 allele providing further evidence of the complexity of the MHC. Additional genotyping to refine the haplotypes followed by screening for potential variants may lead to the identification of a novel MS susceptibility gene(s) in this unique population.

Homozygous WNT10b mutation and complex inheritance in Split-Hand/Foot Malformation

Split-Hand/Foot Malformation (SHFM) is a complex limb malformation affecting the central rays of the autopod. We studied a large consanguineous kindred afflicted with autosomal recessive SHFM. Twelve affected members had central feet reductions with or without hand involvement while the remaining one had the mildest phenotype and atypical SHFM. We identified by homozygosity mapping a novel SHFM locus at 12q13.11-q13 with a maximum multipoint lod score of 5.47 and by subsequent candidate gene approach a homozygous missense WNT10b mutation (p.R332W) in all affected individuals but the atypical case plus in an asymptomatic female. We propose that either a second locus contributes to the manifestation of SHFM phenotype or a suppressor locus prevented trait manifestation in the non-penetrant female. We also investigated linkage to the five known SHFM loci. Four of the loci were excluded, while in TP63 [tumor protein p63 (SHFM4)], the only known gene responsible for SHFM, we detected in most affected subjects a rare insertion variant (rs34201045) at the alternate promoter used for transcription of the N-terminal-truncated p63 isotype. This is the first reported WNT10b mutation on the pathogenesis of limb development and recessive mutation in SHFM.

Linkage disequilibrium between STRPs and SNPs across the human genome

Patterns of linkage disequilibrium (LD) reveal the action of evolutionary processes and provide crucial information for association mapping of disease genes. Although recent studies have described the landscape of LD among single nucleotide polymorphisms (SNPs) from across the human genome, associations involving other classes of molecular variation remain poorly understood. In addition to recombination and population history, mutation rate and process are expected to shape LD. To test this idea, we measured associations between short-tandem-repeat polymorphisms (STRPs), which can mutate rapidly and recurrently, and SNPs in 721 regions across the human genome. We directly compared STRP-SNP LD with SNP-SNP LD from the same genomic regions in the human HapMap populations. The intensity of STRP-SNP LD, measured by the average of D', was reduced, consistent with the action of recurrent mutation. Nevertheless, a higher fraction of STRP-SNP pairs than SNP-SNP pairs showed significant LD, on both short (up to 50 kb) and long (cM) scales. These results reveal the substantial effects of mutational processes on LD at STRPs and provide important measures of the potential of STRPs for association mapping of disease genes.

Molecular pathology in anatomic pathology practice: a review of basic principles

Molecular testing in pathology emerged shortly after polymerase chain reaction became a standard molecular biology assay. Testing efforts began in the clinical laboratories primarily with assays for genetically inherited diseases and assays for clonality in hematologic malignancies. Today, the field has evolved into "molecular diagnostics," which encompasses testing in almost every area of anatomic pathology. Molecular testing is now even making its way definitively into both surgical pathology and cytopathology, although molecular anatomic pathology is still young with few standard tissue-based molecular assays. As more clinically valuable information is gained from molecular pathology testing of tissues, unique challenges are also becoming apparent at the intersection between tissue diagnosis and DNA diagnosis. This review focuses on basic molecular pathology concepts, with particular emphasis on the challenge of tissue-based testing in anatomic pathology.

Ignoring intermarker linkage disequilibrium induces false-positive evidence of linkage for consanguineous pedigrees when genotype data is missing for any pedigree member

Missing genotype data can increase false-positive evidence for linkage when either parametric or nonparametric analysis is carried out ignoring intermarker linkage disequilibrium (LD). Previously it was demonstrated by Huang et al. [1] that no bias occurs in this situation for affected sib-pairs with unrelated parents when either both parents are genotyped or genotype data is available for two additional unaffected siblings when parental genotypes are missing. However, this is not the case for autosomal recessive consanguineous pedigrees, where missing genotype data for any pedigree member within a consanguinity loop can increase false-positive evidence of linkage. False-positive evidence for linkage is further increased when cryptic consanguinity is present. The amount of false-positive evidence for linkage, and which family members aid in its reduction, is highly dependent on which family members are genotyped. When parental genotype data is available, the false-positive evidence for linkage is usually not as strong as when parental genotype data is unavailable. For a pedigree with an affected proband whose first-cousin parents have been genotyped, further reduction in the false-positive evidence of linkage can be obtained by including genotype data from additional affected siblings of the proband or genotype data from the proband's sibling-grandparents. For the situation, when parental genotypes are unavailable, false-positive evidence for linkage can be reduced by including genotype data from either unaffected siblings of the proband or the proband's married-in-grandparents in the analysis.

Genome-wide scan of bipolar disorder and investigation of population stratification effects on linkage: support for susceptibility loci at 4q21, 7q36, 9p21, 12q24, 14q24, and 16p13

Bipolar disorder (BPD) is a complex genetic disorder with cycling symptoms of depression and mania. Despite the extreme complexity of this psychiatric disorder, attempts to localize genes which confer vulnerability to the disorder have had some success. Chromosomal regions including 4p16, 12q24, 18p11, 18q22, and 21q21 have been repeatedly linked to BPD in different populations. Here we present the results of a whole genome scan for linkage to BPD in an Irish population. Our most significant result was at 14q24 which yielded a non-parametric LOD (NPL) score of 3.27 at the D14S588 marker with a nominal P-value of 0.0006 under a narrow (bipolar type I only) model of affection. We previously reported linkage to 14q22-24 in a subset of the families tested in this analysis. We also obtained suggestive evidence for linkage at 4q21, 9p21, 12q24, and 16p13, chromosomal regions that have all been previously linked to BPD. Additionally, we report on a novel approach to linkage analysis, STRUCTURE-Guided Linkage Analysis (SGLA), which is designed to reduce genetic heterogeneity and increase the power to detect linkage. Application of this technique resulted in more highly significant evidence for linkage of BPD to three regions including 16p13, a locus that has been repeatedly linked to numerous psychiatric disorders.

Multistage Sampling for Genetic Studies

In the past, to study Mendelian diseases, segregating families have been carefully ascertained for segregation analysis, followed by collecting extended multiplex families for linkage analysis. This would then be followed by association studies, using independent case-control samples and/or additional family data. Recently, for complex diseases, the initial sampling has been for a genome-wide linkage analysis, often using independent sib-pairs or nuclear families, to identify candidate regions for follow-up with association studies, again using case-control samples and/or additional family data. We now have the ability to conduct genome-wide association studies using 100,000-500,000 diallelic genetic markers. For such studies we focus especially on efficient two-stage association sampling designs, which can retain nearly optimal statistical power at about half the genotyping cost. Similarly, beginning an association study by genotyping pooled samples may also be a viable option if the cost of accurately pooling DNA samples outweighs genotyping costs. Finally, we note that the sampling of family data for linkage analysis is not a practice that should be automatically discontinued.

A novel method combining linkage disequilibrium information and imputed functional knowledge for tagSNP selection

Analyses of high-density SNPs in genetic studies have the potential problems of prohibitive genotyping costs and inflated false discovery rates. Current methods select subsets of representative SNPs (tagSNPs) using information either on potential biologic functionality of the SNPs or on the underlying linkage disequilibrium (LD) structure, but not both. Combining the two types of information may lead to more effective tagSNP selection. The proposed method combines both functional and LD information using a weighted factor analysis (WFA) model. The WFA was applied to the dense SNP collection from 129 genes sequenced by the SeattleSNPs Program for Genomic Application. TagSNPs selected by WFA were compared with those selected by an LD-based method. WFA allowed prioritization of SNPs that would otherwise share equivalent ranking due to underlying LD structure alone. Furthermore, WFA consistently included SNPs not selected by function or by LD alone. A literature review of a subset of genes revealed that SNPs selected by WFA were more likely represented in published reports.

A genome-wide scan of loci linked to serum adiponectin in two populations of African descent

OBJECTIVE

The objectives were to identify quantitative trait loci linked to serum adiponectin concentration and to estimate heritability in two populations of African descent.

RESEARCH METHODS AND PROCEDURES

We conducted a genome scan for serum adiponectin concentration in two populations of African descent. Genome-wide microsatelitte markers were typed in an African-American population consisting of 203 families from the Chicago area and in a Nigerian Yoruba population consisting of 146 families. Linkage analysis was performed to identify loci. Variance component model was used to estimate heritability.

RESULTS

Estimates of heritability adjusted for age, gender, and BMI were 0.45 and 0.70 for the African-American and Nigerian families, respectively. In both populations, adiponectin was significantly negatively correlated with BMI, height, and weight. After adjusting for age, gender, and BMI, we found evidence of genetic linkage to adiponectin on chromosomes 11 [limit of detection (LOD) score = 2.89] and 17 (LOD score = 1.35) in the Nigerian sample. Among the African-Americans, we found genetic linkage on chromosomes 2 (LOD score = 1.82), 4 (LOD score = 2.12), and 11 (LOD score = 2.33). Analysis based on combined data yielded a maximum LOD score of 3.21 on chromosome 11.

DISCUSSION

Consistency of the finding on chromosome 11 suggests that this region is likely to be involved in regulation of adiponectin, either through a primary influence on hormone levels or through pathways influencing body composition. These results suggest that adiponectin could be a potential therapeutic target for obesity.

Identification of a locus modulating serum C-reactive protein levels on chromosome 5p15

OBJECTIVE

Individual propensity to chronic, low-grade inflammation--a determinant of atherosclerosis-is in part under the control of genetic factors. To identify genes involved in this modulation, we performed a 10cM genome screen for linkage with plasma C-reactive protein in 38 extended families including 317 non-diabetic and 177 type 2 diabetic family members (2547 relative pairs).

METHODS AND RESULTS

In a variance component analysis, heritability of CRP values was significant (h(2)=0.39, p<0.0001). This effect was independent of BMI and was present in both diabetic (h(2)=0.42, p=0.003) and non-diabetic (h(2)=0.34, p=0.0015) relatives. The strongest evidence of linkage with CRP was on chromosome 5p15, where the LOD score reached genome-wide significance (LOD=3.41, genome-wide p=0.013). Both diabetic and non-diabetic family members contributed to linkage at this location. Smaller linkage peaks were detected on chromosomes 5q35 (LOD=1.35) and 17p11 (LOD=1.33). When the analysis was restricted to diabetic family members, another peak of moderate intensity (LOD=2.17) was evident at 3p21.

CONCLUSIONS

A major gene influencing CRP levels appears to be located on chromosome 5p15, with an effect that is independent of diabetes. Another gene on 3p21 may control CRP variation but only in the presence of a diabetic or insulin-resistant environment.

Genetic characterization of a new set of recombinant inbred lines (LGXSM) formed from the inter-cross of SM/J and LG/J inbred mouse strains

A new set of LGXSM recombinant inbred (RI) strains is presented. The RI strain panel consists of 18 remaining strains of the original 55 founding strains. Strain characterization is based on 506 polymorphic microsatellites and 4,289 single nucleotide polymorphisms (SNPs) distributed across the genome. Average microsatellite inter-marker distance is 4.80+/-4.84 Mb or 2.91+/-3.21 F(2) cM. SNPs are more densely spaced at 0.57+/-1.27 Mb. Ninety-five percent of all microsatellite inter-marker intervals are separated by less than 15.00 Mb or 8.50 F(2) cM, while 95% of the SNPs are less than 0.95 Mb apart. Strains show expected low levels of nonsyntenic association among loci and complete genomic independence. During inbreeding, the RI strains went through strong natural selection on the agouti locus on Chromosome 2, especially when the epistatically interacting tyrosinase locus on Chromosome 7 carried the wild-type allele. The LG/J and SM/J strains differ in a large number of biomedically important traits, and they and their inter-cross progeny have been used in multiple mapping studies. The LGxSM RI strain panel provides a powerful new resource for mapping the genetic bases of complex traits and should prove to be of great biomedical utility in modeling complex human diseases such as obesity and diabetes.

Evidence for a gene influencing heart rate on chromosome 5p13-14 in a meta-analysis of genome-wide scans from the NHLBI Family Blood Pressure Program

Background

Elevated resting heart rate has been shown in multiple studies to be a strong predictor of cardiovascular disease. Previous family studies have shown a significant heritable component to heart rate with several groups conducting genomic linkage scans to identify quantitative trait loci.

Methods

We performed a genome-wide linkage scan to identify quantitative trait loci influencing resting heart rate among 3,282 Caucasians and 3,989 African-Americans in three independent networks comprising the Family Blood Pressure Program (FBPP) using 368 microsatellite markers. Mean heart rate measurements were used in a regression model including covariates for age, body mass index, pack-years, currently drinking alcohol (yes/no), hypertension status and medication usage to create a standardized residual for each gender/ethnic group within each study network. This residual was used in a nonparametric variance component model to generate a LOD score and a corresponding P value for each ethnic group within each study network. P values from each ethnic group and study network were merged using an adjusted Fisher's combining P values method and the resulting P values were converted to LOD scores. The entire analysis was redone after individuals currently taking beta-blocker medication were removed.

Results

We identified significant evidence of linkage (LOD = 4.62) to chromosome 10 near 142.78 cM in the Caucasian group of HyperGEN. Between race and network groups we identified a LOD score of 1.86 on chromosome 5 (between 39.99 and 45.34 cM) in African-Americans in the GENOA network and the same region produced a LOD score of 1.12 among Caucasians within a different network (HyperGEN). Combining all network and race groups we identified a LOD score of 1.92 (P = 0.0013) on chromosome 5p13-14. We assessed heterogeneity for this locus between networks and ethnic groups and found significant evidence for low heterogeneity (P ≤ 0.05).

Conclusion

We found replication (LOD > 1) between ethnic groups and between study networks with low heterogeneity on chromosome 5p13-14 suggesting that a gene in this region influences resting heart rate.

Advances in statistical human genetics over the last 25 years

The past 25 years has seen an explosion in the number of genetic markers that can be measured on DNA samples at an ever decreasing cost. Although basic statistical methods for analysing such data gathered on samples of either independent individuals or family members, one or two markers at a time, were already well developed before this explosion occurred, there has been a corresponding burst in activity to develop multiple marker models to find disease-causing gene variants, capitalizing on the data that have become available, to increase the power of such methods. This has required the concomitant development of faster algorithms to speed up the computation of various likelihoods. For linkage analysis, to obtain the approximate locations for genes of interest, Mendelian segregation models have been extended to be more realistic and statistical models that do not assume specific modes of inheritance have been extended to allow for the analysis of larger pedigree structures. For association analysis, to obtain more precise locations for genes of interest, the recent completion of the first stage of the HapMap project has spurred the development, still underway, of novel experimental designs and analytical methods to combat the curse of dimensionality and the resulting multiple testing problem. Perhaps the greatest current challenge concerns how best to gather and synthesize the many lines of evidence possible in order to discover the genetic determinants underlying complex diseases.