Linkage disequilibrium between STRPs and SNPs across the human genome

Characterization and visualization of tandem repeats at genome scale

Tandem repeat (TR) variation is associated with gene expression changes and numerous rare monogenic diseases. Although long-read sequencing provides accurate full-length sequences and methylation of TRs, there is still a need for computational methods to profile TRs across the genome. Here we introduce the Tandem Repeat Genotyping Tool (TRGT) and an accompanying TR database. TRGT determines the consensus sequences and methylation levels of specified TRs from PacBio HiFi sequencing data. It also reports reads that support each repeat allele. These reads can be subsequently visualized with a companion TR visualization tool. Assessing 937,122 TRs, TRGT showed a Mendelian concordance of 98.38%, allowing a single repeat unit difference. In six samples with known repeat expansions, TRGT detected all expansions while also identifying methylation signals and mosaicism and providing finer repeat length resolution than existing methods. Additionally, we released a database with allele sequences and methylation levels for 937,122 TRs across 100 genomes.

YHP: Y-chromosome Haplogroup Predictor for predicting male lineages based on Y-STRs

Human Y chromosome reflects the evolutionary process of males. Male lineage tracing by Y chromosome is of great use in evolutionary, forensic, and anthropological studies. Identifying the male lineage based on the specific distribution of Y haplogroups narrows down the investigation scope, which has been used in forensic scenarios. However, existing software aids in familial searching using Y-STRs (Y-chromosome short tandem repeats) to predict Y-SNP (Y-chromosome single nucleotide polymorphism) haplogroups, they often lack resolution. In this study, we developed YHP (Y Haplogroup Predictor), a novel software offering high-resolution haplogroup inference without requiring extensive Y-SNP sequencing. Leveraging existing datasets (219 haplogroups, 4064 samples in total), YHP predicts haplogroups with 0.923 accuracy under the highest haplogroup resolution, employing a random forest algorithm. YHP, available on Github (https://github.com/cissy123/YHP-Y-Haplogroup-Predictor-), facilitates high-resolution haplogroup prediction, haplotype mismatch analysis, and haplotype similarity comparison. Notably, it demonstrates efficacy in East Asian populations, benefiting from training data from eight distinct East Asian ethnic populations. Moreover, it enables seamless integration of additional training sets, extending its utility to diverse populations.

Microsatellites used in forensics are in regions enriched for trait-associated variants

The 20 short tandem repeat (STR) loci of the combined DNA index system (CODIS) are the basis of the vast majority of forensic genetics in the United States. One argument for permissive rules about the collection of CODIS genotypes is that the CODIS loci are thought to contain little information about ancestry or traits. However, in the past 20 years, a growing field has identified hundreds of thousands of genotype-trait associations. Here, we conduct a survey of the landscape of such associations surrounding the CODIS loci as compared with non-CODIS STRs. Although this study cannot establish or quantify associations between CODIS genotypes and phenotypes, we find that the regions around the CODIS loci are enriched for both known pathogenic variants (> 90th percentile) and for trait-associated SNPs identified in genome-wide association studies (GWAS) (≥ 95th percentile in 10kb and 100kb flanking regions), compared with other random sets of autosomal tetranucleotide-repeat STRs.

Microsatellites used in forensics are located in regions unusually rich in trait-associated variants

The 20 short tandem repeat (STR) markers of the combined DNA index system (CODIS) are the basis of the vast majority of forensic genetics in the United States. One argument for permissive rules about the collection of CODIS genotypes is that the CODIS markers are thought to contain information relevant to identification only (such as a human fingerprint would), with little information about ancestry or traits. However, in the past 20 years, a quickly growing field has identified hundreds of thousands of genotype-trait associations. Here we conduct a survey of the landscape of such associations surrounding the CODIS loci as compared with non-CODIS STRs. We find that the regions around the CODIS markers are enriched for both known pathogenic variants (>90th percentile) and for SNPs identified as trait-associated in genome-wide association studies (GWAS) (≥95th percentile in 10kb and 100kb flanking regions), compared with other random sets of autosomal tetranucleotide-repeat STRs. Although it is not obvious how much phenotypic information CODIS would need to convey to strain the "DNA fingerprint" analogy, the CODIS markers, considered as a set, are in regions unusually dense with variants with known phenotypic associations.

Genome-wide contribution of common short-tandem repeats to Parkinson's disease genetic risk

Parkinson's disease is a complex neurodegenerative disorder with a strong genetic component, for which most known disease-associated variants are single nucleotide polymorphisms (SNPs) and small insertions and deletions (indels). DNA repetitive elements account for >50% of the human genome; however, little is known of their contribution to Parkinson's disease aetiology. While select short tandem repeats (STRs) within candidate genes have been studied in Parkinson's disease, their genome-wide contribution remains unknown. Here we present the first genome-wide association study of STRs in Parkinson's disease. Through a meta-analysis of 16 imputed genome-wide association study cohorts from the International Parkinson's Disease Genomic Consortium (IPDGC), totalling 39 087 individuals (16 642 cases and 22 445 controls of European ancestry), we identified 34 genome-wide significant STR loci (P < 5.34 × 10-6), with the strongest signal located in KANSL1 [chr17:44 205 351:[T]11, P = 3 × 10-39, odds ratio = 1.31 (95% confidence interval = 1.26-1.36)]. Conditional-joint analyses suggested that four significant STRs mapping nearby NDUFAF2, TRIML2, MIRNA-129-1 and NCOR1 were independent from known risk SNPs. Including STRs in heritability estimates increased the variance explained by SNPs alone. Gene expression analysis of STRs (eSTRs) in RNA sequencing data from 13 brain regions identified significant associations of STRs influencing the expression of multiple genes, including known Parkinson's disease genes. Further functional annotation of candidate STRs revealed that significant eSTRs within NUDFAF2 and ZSWIM7 overlap with regulatory features and are associated with change in the expression levels of nearby genes. Here, we show that STRs at known and novel candidate loci contribute to Parkinson's disease risk and have functional effects in disease-relevant tissues and pathways, supporting previously reported disease-associated genes and giving further evidence for their functional prioritization. These data represent a valuable resource for researchers currently dissecting Parkinson's disease risk loci.

Revisiting tandem repeats in psychiatric disorders from perspectives of genetics, physiology, and brain evolution

Genome-wide association studies (GWASs) have revealed substantial genetic components comprised of single nucleotide polymorphisms (SNPs) in the heritable risk of psychiatric disorders. However, genetic risk factors not covered by GWAS also play pivotal roles in these illnesses. Tandem repeats, which are likely functional but frequently overlooked by GWAS, may account for an important proportion in the "missing heritability" of psychiatric disorders. Despite difficulties in characterizing and quantifying tandem repeats in the genome, studies have been carried out in an attempt to describe impact of tandem repeats on gene regulation and human phenotypes. In this review, we have introduced recent research progress regarding the genomic distribution and regulatory mechanisms of tandem repeats. We have also summarized the current knowledge of the genetic architecture and biological underpinnings of psychiatric disorders brought by studies of tandem repeats. These findings suggest that tandem repeats, in candidate psychiatric risk genes or in different levels of linkage disequilibrium (LD) with psychiatric GWAS SNPs and haplotypes, may modulate biological phenotypes related to psychiatric disorders (e.g., cognitive function and brain physiology) through regulating alternative splicing, promoter activity, enhancer activity and so on. In addition, many tandem repeats undergo tight natural selection in the human lineage, and likely exert crucial roles in human brain evolution. Taken together, the putative roles of tandem repeats in the pathogenesis of psychiatric disorders is strongly implicated, and using examples from previous literatures, we wish to call for further attention to tandem repeats in the post-GWAS era of psychiatric disorders.

Further insight into the global variability of the OCA2-HERC2 locus for human pigmentation from multiallelic markers

The OCA2-HERC2 locus is responsible for the greatest proportion of eye color variation in humans. Numerous studies extensively described both functional SNPs and associated patterns of variation over this region. The goal of our study is to examine how these haplotype structures and allelic associations vary when highly variable markers such as microsatellites are used. Eleven microsatellites spanning 357 Kb of OCA2-HERC2 genes are analyzed in 3029 individuals from worldwide populations. We found that several markers display large differences in allele frequency (10% to 35% difference) among Europeans, East Asians and Africans. In Europe, the alleles showing increased frequency can also discriminate individuals with (IrisPlex) predicted blue and brown eyes. Distinct haplotypes are identified around the variants C and T of the functional SNP rs12913832 (associated to blue eyes), with linkage disequilibrium r² values significant up to 237 Kb. The haplotype carrying the allele rs12913832 C has high frequency (76%) in blue eye predicted individuals (30% in brown eye predicted individuals), while the haplotype associated to the allele rs12913832 T is restricted to brown eye predicted individuals. Finally, homozygosity values reach levels of 91% near rs12913832. Odds ratios show values of 4.2, 7.4 and 10.4 for four markers around rs12913832 and 7.1 for their core haplotype. Hence, this study provides an example on the informativeness of multiallelic markers that, despite their current limited potential contribution to forensic eye color prediction, supports the use of microsatellites for identifying causing variants showing similar genetic features and history.

Linkage and linkage disequilibrium among the markers in the forensic MPS panels

For the past two to three decades, forensic DNA evidence has been analyzed with a limited number of short tandem repeats (STRs), and these STRs are usually assumed to be independent for statistical calculations. With the development and implementation of the MPS technologies, more autosomal markers, both single nucleotide polymorphisms (SNPs) and STRs, can be analyzed. A number of these markers are physically very close to each other, and it may not be appropriate to assume all these markers are genetically unlinked or in linkage equilibrium. In this study, publicly accessible genomic data from five representative populations were used to evaluate the genetic linkage and linkage disequilibrium (LD) between autosomal markers represented in six major commercial panels (in total, 362 markers). Among the 3041 syntenic marker pairs, 1524 pairs had sex-average genetic distances <50 cM, and thus, these marker pairs can be considered as genetically linked. Among the 143 marker pairs with physical distances <1 Mb, 19 LD haplotype blocks (comprising 39 SNPs in total) were detected for at least one of the tested populations. Statistical methods for interpreting linked markers and/or markers in LD were suggested for various case scenarios.

Translational genomics and beyond in bipolar disorder

Genome-wide association studies (GWAS) have revealed multiple genomic loci conferring risk of bipolar disorder (BD), providing hints for its underlying pathobiology. However, there are still remaining questions to answer. For example, discordance exists between BD heritability estimated with earlier epidemiological evidence and that calculated based on common GWAS variations. Where is the "missing heritability"? How can we explain the biology of the disease based on genetic findings? In this review, we summarize the accomplishments and limitations of current BD GWAS, and discuss potential reasons for the "missing heritability." In addition, progresses of research for the biological mechanisms underlying BD genetic risk using brain tissues, reprogrammed cells, and model animals are reviewed. While our knowledge of BD genetic basis is significantly promoted by these efforts, the complexities of gene regulation in the genome, the spatial-temporal heterogeneity during brain development, and the limitations of different experimental models should always be considered. Notably, several genes have been widely studied given their relatively well-characterized involvement in BD (e.g., CACAN1C and ANK3), and findings of these genes are summarized to both outline possible biological mechanisms of BD and describe examples of translating GWAS discoveries into the pathophysiology.

Mathematical Properties of Linkage Disequilibrium Statistics Defined by Normalization of the Coefficient D = pAB - pApB

BACKGROUND

Many statistics for measuring linkage disequilibrium (LD) take the form of a normalization of the LD coefficient D. Different normalizations produce statistics with different ranges, interpretations, and arguments favoring their use.

METHODS

Here, to compare the mathematical properties of these normalizations, we consider 5 of these normalized statistics, describing their upper bounds, the mean values of their maxima over the set of possible allele frequency pairs, and the size of the allele frequency regions accessible given specified values of the statistics.

RESULTS

We produce detailed characterizations of these properties for the statistics d and ρ, analogous to computations previously performed for r2. We examine the relationships among the statistics, uncovering conditions under which some of them have close connections.

CONCLUSION

The results contribute insight into LD measurement, particularly the understanding of differences in the features of different LD measures when computed on the same data.

Modeling copy number variation in the genomic prediction of maize hybrids

Our study indicates that copy variants may play an essential role in the phenotypic variation of complex traits in maize hybrids. Moreover, predicting hybrid phenotypes by combining additive-dominance effects with copy variants has the potential to be a viable predictive model. Non-additive effects resulting from the actions of multiple loci may influence trait variation in single-cross hybrids. In addition, complementation of allelic variation could be a valuable contributor to hybrid genetic variation, especially when crossing inbred lines with higher contents of copy gains. With this in mind, we aimed (1) to study the association between copy number variation (CNV) and hybrid phenotype, and (2) to compare the predictive ability (PA) of additive and additive-dominance genomic best linear unbiased prediction model when combined with the effects of CNV in two datasets of maize hybrids (USP and HELIX). In the USP dataset, we observed a significant negative phenotypic correlation of low magnitude between copy number loss and plant height, revealing a tendency that more copy losses lead to lower plants. In the same set, when CNV was combined with the additive plus dominance effects, the PA significantly increased only for plant height under low nitrogen. In this case, CNV effects explicitly capture relatedness between individuals and add extra information to the model. In the HELIX dataset, we observed a pronounced difference in PA between additive (0.50) and additive-dominance (0.71) models for predicting grain yield, suggesting a significant contribution of dominance. We conclude that copy variants may play an essential role in the phenotypic variation of complex traits in maize hybrids, although the inclusion of CNVs into datasets does not return significant gains concerning PA.

A reference haplotype panel for genome-wide imputation of short tandem repeats

Short tandem repeats (STRs) are involved in dozens of Mendelian disorders and have been implicated in complex traits. However, genotyping arrays used in genome-wide association studies focus on single nucleotide polymorphisms (SNPs) and do not readily allow identification of STR associations. We leverage next-generation sequencing (NGS) from 479 families to create a SNP + STR reference haplotype panel. Our panel enables imputing STR genotypes into SNP array data when NGS is not available for directly genotyping STRs. Imputed genotypes achieve mean concordance of 97% with observed genotypes in an external dataset compared to 71% expected under a naive model. Performance varies widely across STRs, with near perfect concordance at bi-allelic STRs vs. 70% at highly polymorphic repeats. Imputation increases power over individual SNPs to detect STR associations with gene expression. Imputing STRs into existing SNP datasets will enable the first large-scale STR association studies across a range of complex traits.

Massive variation of short tandem repeats with functional consequences across strains of Arabidopsis thaliana

Short tandem repeat (STR) mutations may comprise more than half of the mutations in eukaryotic coding DNA, yet STR variation is rarely examined as a contributor to complex traits. We assessed this contribution across a collection of 96 strains of Arabidopsis thaliana, genotyping 2046 STR loci each, using highly parallel STR sequencing with molecular inversion probes. We found that 95% of examined STRs are polymorphic, with a median of six alleles per STR across these strains. STR expansions (large copy number increases) are found in most strains, several of which have evident functional effects. These include three of six intronic STR expansions we found to be associated with intron retention. Coding STRs were depleted of variation relative to noncoding STRs, and we detected a total of 56 coding STRs (11%) showing low variation consistent with the action of purifying selection. In contrast, some STRs show hypervariable patterns consistent with diversifying selection. Finally, we detected 133 novel STR-phenotype associations under stringent criteria, most of which could not be detected with SNPs alone, and validated some with follow-up experiments. Our results support the conclusion that STRs constitute a large, unascertained reservoir of functionally relevant genomic variation.

Linkage disequilibrium matches forensic genetic records to disjoint genomic marker sets

Combining genotypes across datasets is central in facilitating advances in genetics. Data aggregation efforts often face the challenge of record matching-the identification of dataset entries that represent the same individual. We show that records can be matched across genotype datasets that have no shared markers based on linkage disequilibrium between loci appearing in different datasets. Using two datasets for the same 872 people-one with 642,563 genome-wide SNPs and the other with 13 short tandem repeats (STRs) used in forensic applications-we find that 90-98% of forensic STR records can be connected to corresponding SNP records and vice versa. Accuracy increases to 99-100% when ∼30 STRs are used. Our method expands the potential of data aggregation, but it also suggests privacy risks intrinsic in maintenance of databases containing even small numbers of markers-including databases of forensic significance.

Linking short tandem repeat polymorphisms with cytosine modifications in human lymphoblastoid cell lines

Inter-individual variation in cytosine modifications has been linked to complex traits in humans. Cytosine modification variation is partially controlled by single nucleotide polymorphisms (SNPs), known as modified cytosine quantitative trait loci (mQTL). However, little is known about the role of short tandem repeat polymorphisms (STRPs), a class of structural genetic variants, in regulating cytosine modifications. Utilizing the published data on the International HapMap Project lymphoblastoid cell lines (LCLs), we assessed the relationships between 721 STRPs and the modification levels of 283,540 autosomal CpG sites. Our findings suggest that, in contrast to the predominant cis-acting mode for SNP-based mQTL, STRPs are associated with cytosine modification levels in both cis-acting (local) and trans-acting (distant) modes. In local scans within the ±1 Mb windows of target CpGs, 21, 9, and 21 cis-acting STRP-based mQTL were detected in CEU (Caucasian residents from Utah, USA), YRI (Yoruba people from Ibadan, Nigeria), and the combined samples, respectively. In contrast, 139,420, 76,817, and 121,866 trans-acting STRP-based mQTL were identified in CEU, YRI, and the combined samples, respectively. A substantial proportion of CpG sites detected with local STRP-based mQTL were not associated with SNP-based mQTL, suggesting that STRPs represent an independent class of mQTL. Functionally, genetic variants neighboring CpG-associated STRPs are enriched with genome-wide association study (GWAS) loci for a variety of complex traits and diseases, including cancers, based on the National Human Genome Research Institute (NHGRI) GWAS Catalog. Therefore, elucidating these STRP-based mQTL in addition to SNP-based mQTL can provide novel insights into the genetic architectures of complex traits.

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.

The landscape of human STR variation

Short tandem repeats are among the most polymorphic loci in the human genome. These loci play a role in the etiology of a range of genetic diseases and have been frequently utilized in forensics, population genetics, and genetic genealogy. Despite this plethora of applications, little is known about the variation of most STRs in the human population. Here, we report the largest-scale analysis of human STR variation to date. We collected information for nearly 700,000 STR loci across more than 1000 individuals in Phase 1 of the 1000 Genomes Project. Extensive quality controls show that reliable allelic spectra can be obtained for close to 90% of the STR loci in the genome. We utilize this call set to analyze determinants of STR variation, assess the human reference genome's representation of STR alleles, find STR loci with common loss-of-function alleles, and obtain initial estimates of the linkage disequilibrium between STRs and common SNPs. Overall, these analyses further elucidate the scale of genetic variation beyond classical point mutations.

Genome-wide association mapping of blood cell traits in mice

Genetic variations in blood cell parameters can impact clinical traits. We report here the mapping of blood cell traits in a panel of 100 inbred strains of mice of the Hybrid Mouse Diversity Panel (HMDP) using genome-wide association (GWA). We replicated a locus previously identified in using linkage analysis in several genetic crosses for mean corpuscular volume (MCV) and a number of other red blood cell traits on distal chromosome 7. Our peak for SNP association to MCV occurred in a linkage disequilibrium (LD) block spanning from 109.38 to 111.75 Mb that includes Hbb-b1, the likely causal gene. Altogether, we identified five loci controlling red blood cell traits (on chromosomes 1, 7, 11, 12, and 16), and four of these correspond to loci for red blood cell traits reported in a recent human GWA study. For white blood cells, including granulocytes, monocytes, and lymphocytes, a total of six significant loci were identified on chromosomes 1, 6, 8, 11, 12, and 15. An average of ten candidate genes were found at each locus and those were prioritized by examining functional variants in the HMDP such as missense and expression variants. These results provide intermediate phenotypes and candidate loci for genetic studies of atherosclerosis and cancer as well as inflammatory and immune disorders in mice.

A potential novel spontaneous preterm birth gene, AR, identified by linkage and association analysis of X chromosomal markers

Preterm birth is the major cause of neonatal mortality and morbidity. In many cases, it has severe life-long consequences for the health and neurological development of the newborn child. More than 50% of all preterm births are spontaneous, and currently there is no effective prevention. Several studies suggest that genetic factors play a role in spontaneous preterm birth (SPTB). However, its genetic background is insufficiently characterized. The aim of the present study was to perform a linkage analysis of X chromosomal markers in SPTB in large northern Finnish families with recurrent SPTBs. We found a significant linkage signal (HLOD  = 3.72) on chromosome locus Xq13.1 when the studied phenotype was being born preterm. There were no significant linkage signals when the studied phenotype was giving preterm deliveries. Two functional candidate genes, those encoding the androgen receptor (AR) and the interleukin-2 receptor gamma subunit (IL2RG), located near this locus were analyzed as candidates for SPTB in subsequent case-control association analyses. Nine single-nucleotide polymorphisms (SNPs) within these genes and an AR exon-1 CAG repeat, which was previously demonstrated to be functionally significant, were analyzed in mothers with preterm delivery (n = 272) and their offspring (n = 269), and in mothers with exclusively term deliveries (n = 201) and their offspring (n = 199), all originating from northern Finland. A replication study population consisting of individuals born preterm (n = 111) and term (n = 197) from southern Finland was also analyzed. Long AR CAG repeats (≥26) were overrepresented and short repeats (≤19) underrepresented in individuals born preterm compared to those born at term. Thus, our linkage and association results emphasize the role of the fetal genome in genetic predisposition to SPTB and implicate AR as a potential novel fetal susceptibility gene for SPTB.

Immunogenetics of hematopoietic stem cell transplantation: the contribution of microsatellite polymorphism studies

Polymorphisms of short tandem repeats of <10 nucleotides, or microsatellites (Msat), are largely used for post-transplant chimerism analyses in clinical hematopoietic stem cell transplantation (HSCT). Compared to single nucleotide polymorphisms (SNP), they have the advantage of a higher degree of allelic polymorphism and thus a potentially larger degree of informativity. Msat markers contribute to approximately 3% of the human genome and have been highly informative in disease association studies, population genetics, forensic medicine and organ and HSC transplantation. They allowed to expand our knowledge of the haplotypic structure of the HLA complex, including the noncoding sequences in the MHC, and to reach a better characterization of immunological phenotypes. Among the different immunogenetic studies in HSCT patients reviewed here, four Msat loci linked to cytokine genes have been analysed by a number of laboratories as potential candidates markers for HSCT outcome: IFNG, TNFd, IL-10(-1064) and IL-1RN. The low patient numbers and high diversity of clinical parameters account for some heterogeneity of the results. Among the trends starting to emerge from these studies, specific TNFd Msat alleles seem to be associated with acute graft-versus-host disease and mortality. Patient/donor Msat incompatibilities have also been used as surrogate markers to map biologically relevant polymorphisms, with a main focus on MHC-resident genetic variation. High throughput SNP typing and next-generation sequencing technologies will allow acquisition of large-scale genomic data and should allow refined analyses of clinically relevant genotypes in the transplantation settting, although the heterogeneity of the study cohorts will remain an issue. The analysis of Msat polymorphisms may still have a place in functional studies on the impact of Msat diversity in the control of immune response gene expression.

Including copy number variation in association studies to predict genotypic values

The objective of this study was to investigate, both empirically and deterministically, the ability to explain genetic variation resulting from a copy number polymorphism (CNP) by including the CNP, either by its genotype or by a continuous derivation thereof, alone or together with a nearby single nucleotide polymorphism (SNP) in the model. This continuous measure of a CNP genotype could be a raw hybridization measurement, or a predicted CNP genotype. Results from simulations showed that the linkage disequilibrium (LD) between an SNP and CNP was lower than LD between two SNPs, due to the higher mutation rate at the CNP loci. The model R(2) values from analysing the simulated data were very similar to the R(2) values predicted with the deterministic formulae. Under the assumption that x copies at a CNP locus lead to the effect of x times the effect of 1 copy, including a continuous measure of a CNP locus in the model together with the genotype of a nearby SNP increased power to explain variation at the CNP locus, even when the continuous measure explained only 15% of the variation at the CNP locus.

Interlocus non-random association of multiallelic polymorphisms spanning the coagulation factor VIII gene on human chromosome distalmost Xq28

The most common severe hereditary bleeding disorder phenotype in humans, the coagulation factor VIII (F8) deficiency haemophilia A (HEMA), maps on Xq28 band, a region that comprises 11.7% of genes and 14.2% of phenotypes on X chromosome. Information about the distribution and extent of gametic disequilibrium (GD) covering the F8 gene is scarce, despite its relevance for linkage and association studies. The aim of this study was to determine the patterns, by frequency and strength, of non-random multiallelic interallelic associations between two-locus combinations of seven microsatellite loci (REN90833, F8Int25.2, F8Int22, F8Int13.2, HEMA154311.3, TMLHEInt5 and HEMA154507.3, in that physical order) spanning 0.813 Mb on distalmost Xq28. We measured sign-based interallelic D' coefficients in 106 men and in 100 women drawn from a single unrelated Brazilian population. Significance and patterns of GD using haploid and phased diploid sample probabilities were close to conformity. Only 9.18% of the variance of D' could be accounted for by changes in length, indicating that GD is not a monotonically decreasing function of length. We defined two regions of overlapping long-range GD extending 698 735 base pairs (bp) (REN90833/TMLHEInt5 block) and 689 900 bp (F8Int13.2/HEMA154507.3 block) The extent of GD overlap is 575 637 bp (F8Int13.2/TMLHEInt5 interstice). Extended haplotype homozygosity analysis centred at the F8 intronic loci revealed that the most frequent core haplotypes decay the least in the flanking GD. The F8 intronic loci attend distinct non-random association forces; F8Int13.2 serves at maintenance of the long-range overlapping pattern of GD, whereas F8Int25.2 and F8Int22 serve at lessening it in force or effect.

DOGSET: pre-designed primer sets for fine-scale mapping and DNA sequence interrogation in the dog

DOGSET is an online resource that provides access to primer sequences that have been computationally mined from the reference genome using heuristic algorithms. The electronic repository includes PCR primers corresponding to 32,135 markers for genetic mapping and 334,657 sequence-tagged gene elements for targeted re-sequencing and mutation discovery. A customized report that tailors primer design to wet bench protocols can be exported for a region of interest by specifying genome coordinates in a graphical user interface.

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.

New genetic evidence for involvement of the dopamine system in migraine with aura

In order to systematically test the hypothesis that genetic variation in the dopamine system contributes to the susceptibility to migraine with aura (MA), we performed a comprehensive genetic association study of altogether ten genes from the dopaminergic system in a large German migraine with aura case-control sample. Based on the genotyping results of 53 variants across the ten genes in 270 MA cases and 272 controls, three genes-DBH, DRD2 and SLC6A3-were chosen to proceed to additional genotyping of 380 MA cases and 378 controls. Four of the 26 genotyped polymorphisms in these three genes displayed nominally significant allelic P-values in the sample of 650 MA patients and 650 controls. Three of these SNPs [rs2097629 in DBH (uncorrected allelic P value = 0.0012, OR = 0.77), rs7131056 in DRD2 (uncorrected allelic P value = 0.0018, OR = 1.28) and rs40184 in SLC6A3 (uncorrected allelic P value = 0.0082, OR = 0.81)] remained significant after gene-wide correction for multiple testing by permutation analysis. Further consideration of imputed genotype data from 2,937 British control individuals did not affirm the association with DRD2, but supported the associations with DBH and SLC6A3. Our data provide new evidence for an involvement of components of the dopaminergic system-in particular the dopamine-beta hydroxylase and dopamine transporter genes-to the pathogenesis of migraine with aura.

Variants of the elongator protein 3 (ELP3) gene are associated with motor neuron degeneration

Amyotrophic lateral sclerosis (ALS) is a spontaneous, relentlessly progressive motor neuron disease, usually resulting in death from respiratory failure within 3 years. Variation in the genes SOD1 and TARDBP accounts for a small percentage of cases, and other genes have shown association in both candidate gene and genome-wide studies, but the genetic causes remain largely unknown. We have performed two independent parallel studies, both implicating the RNA polymerase II component, ELP3, in axonal biology and neuronal degeneration. In the first, an association study of 1884 microsatellite markers, allelic variants of ELP3 were associated with ALS in three human populations comprising 1483 people (P = 1.96 × 10⁻⁹). In the second, an independent mutagenesis screen in Drosophila for genes important in neuronal communication and survival identified two different loss of function mutations, both in ELP3 (R475K and R456K). Furthermore, knock down of ELP3 protein levels using antisense morpholinos in zebrafish embryos resulted in dose-dependent motor axonal abnormalities [Pearson correlation: −0.49, P = 1.83 × 10⁻¹² (start codon morpholino) and −0.46, P = 4.05 × 10⁻⁹ (splice-site morpholino), and in humans, risk-associated ELP3 genotypes correlated with reduced brain ELP3 expression (P = 0.01). These findings add to the growing body of evidence implicating the RNA processing pathway in neurodegeneration and suggest a critical role for ELP3 in neuron biology and of ELP3 variants in ALS.

Cell-based Models for Discovery of Pharmacogenomic Markers of Anticancer Agent Toxicity

The field of pharmacogenomics is challenging because of the multigenic nature of drug response and toxicity. The candidate gene approach has been traditionally utilized to determine the contribution of genetic variation to a particular phenotype; however, the sequencing of the human genome and the genetic resource provided by the International HapMap Project has allowed researchers to perform genome-wide studies without a priori knowledge. Recent work has demonstrated the usefulness of cell-based models for pharmacogenomic discovery using the HapMap samples, which are a panel of well-genotyped, human lymphoblastoid cell lines (LCLs) derived from 90 Utah residents with ancestry from northern and western Europe (CEU), 90 Yoruba in Ibadan, Nigeria (YRI), 45 Japanese in Tokyo, Japan (JPT) and 45 Han Chinese in Beijing, China (CHB). Using these cell-based models, investigators are able to study not only individual variation in drug response, but also population differences in drug response. Finally, besides single nucleotide polymorphisms (SNPs) and gene expression, these cell-based models can also be used to investigate other genetic (e.g. copy number variants, CNVs), epigenetic or environmental factors responsible for drug response.