Mutation scanning by meltMADGE: validations using BRCA1 and LDLR, and demonstration of the potential to identify severe, moderate, silent, rare, and paucimorphic mutations in the general population

Family-specific aggregation of lipid GWAS variants confers the susceptibility to familial hypercholesterolemia in a large Austrian family

BACKGROUND AND AIMS

Hypercholesterolemia confers susceptibility to cardiovascular disease (CVD). Both serum total cholesterol (TC) and LDL-cholesterol (LDL-C) exhibit a strong genetic component (heritability estimates 0.41-0.50). However, a large part of this heritability cannot be explained by the variants identified in recent extensive genome-wide association studies (GWAS) on lipids. Our aim was to find genetic causes leading to high LDL-C levels and ultimately CVD in a large Austrian family presenting with what appears to be autosomal dominant inheritance for familial hypercholesterolemia (FH).

METHODS

We utilized linkage analysis followed by whole-exome sequencing and genetic risk score analysis using an Austrian multi-generational family with various dyslipidemias, including elevated TC and LDL-C, and one family branch with elevated lipoprotein (a) (Lp(a)).

RESULTS

We did not find evidence for genome-wide significant linkage for LDL-C or apparent causative variants in the known FH genes rather, we discovered a particular family-specific combination of nine GWAS LDL-C SNPs (p = 0.02 by permutation), and putative less severe familial hypercholesterolemia mutations in the LDLR and APOB genes in a subset of the affected family members. Separately, high Lp(a) levels observed in one branch of the family were explained primarily by the LPA locus, including short (<23) Kringle IV repeats and rs3798220.

CONCLUSIONS

Taken together, some forms of FH may be explained by family-specific combinations of LDL-C GWAS SNPs.

Systematic analysis of variants related to familial hypercholesterolemia in families with premature myocardial infarction

Familial hypercholesterolemia (FH) is an oligogenic disorder characterized by markedly elevated low-density lipoprotein cholesterol (LDLC) levels. Variants in four genes have been reported to cause the classical autosomal-dominant form of the disease. FH is largely under-diagnosed in European countries. As FH increases the risk for coronary artery disease (CAD) and myocardial infarction (MI), it might be specifically overlooked in the large number of such patients. Here, we systematically examined the frequency of potential FH-causing variants by exome sequencing in 255 German patients with premature MI and a positive family history for CAD. We further performed co-segregation analyses in an average of 5.5 family members per MI patient. In total, we identified 11 potential disease-causing variants that co-segregate within the families, that is, 5% of patients with premature MI and positive CAD family history had FH. Eight variants were previously reported as disease-causing and three are novel (LDLR.c.811G>A p.(V271I)), PCSK9.c.610G>A (p.(D204N)) and STAP1.c.139A>G (p.(T47A))). Co-segregation analyses identified multiple additional family members carrying one of these FH variants and the clinical phenotype of either FH (n=2) or FH and premature CAD (n=15). However, exome sequencing also revealed that some variants in FH genes, which have been reported to cause FH, do not co-segregate with FH. The data reveal that a large proportion of FH patients escape the diagnosis, even when they have premature MI. Hence, systematic molecular-genetic screening for FH in such patients may reveal a substantial number of cases and thereby allow a timely LDLC-lowering in both FH/MI patients as well as their variant-carrying family members.

Two-dye based arrayed primer extension for simultaneous multigene detection in lipid metabolism

BACKGROUND

Cardiovascular disease (CVD) is one of the major causes of death worldwide. Numerous genetic risk factors in lipid metabolism, including mutations of LDLR, APOB, and PCSK9, as well as polymorphisms of CETP and APOE, have been found to associate with CVD.

METHODS

In this study, a two-dye based arrayed primer extension (APEX) microarray assay for simultaneous multigene (LDLR, APOB, PCSK9, CETP, and APOE) detection was developed. The DNA templates, originating from 1 DNA sample of known genotype and 7 blind DNA samples, were amplified by uniplex PCR.

RESULTS

Optimized conditions for the APEX reaction were determined to include a hybridization temperature of 55°C and a DNA template size of 50-150bp. The total assay including PCR, purification, fragmentation, APEX reaction, and image analysis could be performed in 6h. In total, 48 genotypes were identified among 8 individual DNA samples by APEX analysis.

CONCLUSIONS

The data suggest that this APEX microarray offers a robust, fast, and versatile option for screening these genotypes in hypercholesterolemia patients.

Population mutation scanning of human GHR by meltMADGE and identification of a paucimorphic variant

Current studies of human genetic diversity are focused in two areas: first, detection of rare mutations in highly selected clinical cases; and second, in common single-nucleotide polymorphism (SNP) and haplotype effects in the general population. Less frequent SNPs and "paucimorphisms" remain underexplored, although lower frequency coding SNPs are more likely to have functional impact. We have developed a cost-efficient mutation scanning technology, meltMADGE, for population mutation scanning. Previous research in GHR has explored its role in extreme (-3 SD) growth retardation and, subsequently, "moderate" (-2 SD) growth retardation cases. Here, we describe meltMADGE assays for the entire coding region of GHR. As a first step we have established long polymerase chain reaction subbanks for GHR from 2423 unselected subjects and have applied meltMADGE scanning assays of exons 4 and 5 to these subbanks. A novel paucimorphism present at 439+30A>C (allele frequency: 0.0021) in intron 5 (location chr5:42,695,221 in GRCh37/hg19) was identified in 10 individuals, confirmed by sequencing and analysis made for major phenotypic effects. This approach is relevant to the deep sampling of populations for less frequent sequence diversity, some of which is expected to exert significant phenotypic effects.

Genetic screening for homozygous and heterozygous familial hypercholesterolemia

Familial hypercholesterolemia (FH) is a common inherited disorder that results in premature atherosclerosis. Diagnosis of FH is suspected on the basis of clinical criteria, but confirmation requires genetic testing. In the era of statins, early diagnosis and initiation of treatment can modify disease progression and outcomes. Therefore, cascade screening with a combination of lipid concentration measurements and DNA testing should be used to identify relatives of index cases with a clinical diagnosis of FH. Autosomal dominant FH is related to mutations in the low-density lipoprotein receptor (LDLR), apolipoprotein B-100 (APOB), or proprotein convertase subtilisin/kexin type 9 (PCSK9) genes. Genetic screening of the LDLR gene is challenging to achieve at a feasible cost, especially in people who do not have a founder effect. Nucleotide sequencing of all exons and flanking splicing regions in combination with multiplex ligation probe amplification to detect large insertions or deletions is considered the gold-standard approach to screen for LDLR mutations. Alternatively, the cDNA can be sequenced; however, this procedure is not suitable for use in large populations, because of the need of RNA extraction. Multiplex analysis can be appropriate for population with founder effects or a low number of different mutations. Finally, there are many techniques for a mutation scanning approach, which have some benefits over sequencing, and also with the potential for detecting known and novel mutations. Familial defective Apo B is amenable to genetic diagnosis by screening for a few mutations. Recently, gain-of-function mutations in PCSK9 gene have been demonstrated to cause FH phenotype. Strategies for population screening, cost-effectiveness of genetic screening, ethical aspects, and insurance policies are discussed and need implementation worldwide.

MeltMADGE for mutation scanning of specific genes in population studies

MeltMADGE reconfigures the mutation scanning process of denaturing gradient gel electrophoresis so that the independent variable is time rather than space and the dependent (denaturing) variable is temperature rather than concentration of chemical denaturant. Use of a thermal ramp enables the use of a homogeneous gel and therefore of high-density arrays of wells such as those of microplate array diagonal gel electrophoresis (MADGE). In this configuration, electrophoresis of products on 10-12 96-well meltMADGE gels can be conducted in a 1- to 2-liter tank in a 1- to 2-h run, enabling the scanning of a target amplicon in over 1,000 subjects simultaneously. Gels are read by imaging the fluorescence of UV-excited ethidium bromide, giving a simple, economical system for identifying rarer sequence variants in target genes; it is suitable for large-scale case-control or population studies and other comparable applications. Different amplicons with similar melting characteristics can also be combined in the same run.

Haplotype analyses, mechanism and evolution of common double mutants in the human LDL receptor gene

Familial hypercholesterolemia (FH), an autosomal dominant inherited disorder resulting in increased levels of circulating plasma low-density lipoprotein (LDL), tendon xanthomas and premature coronary artery disease (CAD), is caused by defects in the LDL receptor gene (LDLR). Three widespread LDLR alterations not causing FH (c.1061-8T>C, c.2177C>T and c.829G>A) and one mutation (c.12G>A) with narrow geographical distribution and thought to cause disease were investigated. In an attempt to improve knowledge on their origin, spread and possible selective effects, estimations of the ages of these variants (t generations) and haplotype analysis were performed by genotyping 86 healthy individuals and 98 FH patients in Spain for five LDLR SNPs: c.81T>C, c.1413G>A, c.1725C>T, c.1959T>C, and c.2232G>A; most patients carried two of these LDLR variants simultaneously. It was found that both the c.1061-8T>C (t = 54) and c.2177C>T alterations (t = 62) arose at about the same time (54 and 62 generations ago, respectively) in the CGCTG haplotype, while the c.12G>A mutation (t = 70) appeared in a CGCCG haplotype carrying an earlier c.829G>A alteration (t = 83). The estimated ages of selectively neutral alterations could explain their distribution by migrations. The origin of the c.12G>A mutation could be in the Iberian Peninsula; despite its estimated age, a low selective pressure could explain its conservation in Spain from where it could have spread to China and Mexico, since the sixteenth century through the Spanish/Portuguese colonial expeditions.

dbSNP in the detail and copy number complexities

dbSNP is a general catalog of genetic polymorphism maintained by NCBI, mainly collating information for single nucleotide variations, many of which will be single nucleotide polymorphisms (SNPs), but also including small indels. It takes submissions from many sources, now also including large numbers of sequence variants identified by next-generation sequencing. A number of differently designed studies have attempted to estimate the error rates in data archived in dbSNP. Most recently, a study added to earlier studies identifying specific issues for duplicons and copy number variations (CNVs); earlier analyses have focused on stop codons, splice sites, and the general content of dbSNP. This article overviews dbSNP itself, these studies, and their implications.

Development of a universal chemiluminometric genotyping method for high-throughput detection of 7 LDLR gene mutations in Greek population

OBJECTIVES

Familial hypercholesterolemia (FH) is caused by mutations in the LDL receptor (LDLR) gene. We report the application of a universal method with high allele discrimination properties to the simultaneous genotyping of 7 LDLR mutations in Greeks, in dry-reagent format.

DESIGN AND METHODS

We genotyped mutations C858A, C939A, G1285A, T1352C, G1646A, G1775A, C/T81G. Unpurified amplicons from a multiplex PCR that produced fragments encompassing all 7 mutations were subjected to probe extension reactions in the presence of fluorescein-modified dCTP, and a microtiter well-based assay of extension products with a peroxidase-antifluorescein conjugate and a chemiluminogenic substrate. We used lyophilized dry reagents and assigned genotypes by the signal ratio of normal-to-mutant-specific probe.

RESULTS

We standardized the method and optimised all steps for specificity. The method was validated by genotyping blindly 119 (833 genotypings). Results were fully concordant with other methods used as standards.

CONCLUSIONS

This method is accurate, simple, rapid and robust. The microtiter well format allows genotyping of a large number of samples in parallel for several mutations.

Cubic exact solutions for the estimation of pairwise haplotype frequencies: implications for linkage disequilibrium analyses and a web tool 'CubeX'

BACKGROUND

The frequency of a haplotype comprising one allele at each of two loci can be expressed as a cubic equation (the 'Hill equation'), the solution of which gives that frequency. Most haplotype and linkage disequilibrium analysis programs use iteration-based algorithms which substitute an estimate of haplotype frequency into the equation, producing a new estimate which is repeatedly fed back into the equation until the values converge to a maximum likelihood estimate (expectation-maximisation).

RESULTS

We present a program, "CubeX", which calculates the biologically possible exact solution(s) and provides estimated haplotype frequencies, D', r2 and chi2 values for each. CubeX provides a "complete" analysis of haplotype frequencies and linkage disequilibrium for a pair of biallelic markers under situations where sampling variation and genotyping errors distort sample Hardy-Weinberg equilibrium, potentially causing more than one biologically possible solution. We also present an analysis of simulations and real data using the algebraically exact solution, which indicates that under perfect sample Hardy-Weinberg equilibrium there is only one biologically possible solution, but that under other conditions there may be more.

CONCLUSION

Our analyses demonstrate that lower allele frequencies, lower sample numbers, population stratification and a possible |D'| value of 1 are particularly susceptible to distortion of sample Hardy-Weinberg equilibrium, which has significant implications for calculation of linkage disequilibrium in small sample sizes (eg HapMap) and rarer alleles (eg paucimorphisms, q < 0.05) that may have particular disease relevance and require improved approaches for meaningful evaluation.

Molecular genetics of human growth hormone, insulin-like growth factors and their pathways in common disease

The human growth hormone gene (GH1) and the insulin-like growth factor 1 and 2 genes (IGF1 and IGF2) encode the central elements of a key pathway influencing growth in humans. This "growth pathway" also includes transcription factors, agonists, antagonists, receptors, binding proteins, and endocrine factors that constitute an intrincate network of feedback loops. GH1 is evolutionarily coupled with other genes in linkage disequilibrium in 17q24.2, and the same applies to IGF2 in 11p15.5. In contrast, IGF1 in 12q22-24.1 is not in strong linkage disequilibrium with neighbouring genes. Knowledge of the functional architecture of these regions is important for the understanding of the combined evolution and function of GH1, IGF2 and IGF1 in relation to complex diseases. A number of mutations accounting for rare Mendelian disorders have been described in GH-IGF elements. The constellation of genes in this key pathway contains potential candidates in a number of complex diseases, including growth disorders, metabolic syndrome, diabetes (notably IGF2BP2) cardiovascular disease, and central nervous system diseases, and in longevity, aging and cancer. We review these genes and their associations with disease phenotypes, with special attention to metabolic risk traits.

Prevalence and functionality of paucimorphic and privateMC4Rmutations in a large, unselected European British population, scanned by meltMADGE

Identification of unknown mutations has remained laborious, expensive, and only viable for studies of selected cases. Population-based "reference ranges" of rarer sequence diversity are not available. However, the research and diagnostic interpretation of sequence variants depends on such information. Additionally, this is the only way to determine prevalence of severe, moderate, and silent mutations and is also relevant to the development of screening programs. We previously described a system, meltMADGE, suitable for mutation scanning at the population level. Here we describe its application to a population-based study of MC4R (melanocortin 4 receptor) mutations, which are associated with obesity. We developed nine assays representing MC4R and examined a population sample of 1,100 subjects. Two "paucimorphisms" were identified (c.307G>A/p.Val103Ile in 27 subjects and c.-178A>C in 22 subjects). Neither exhibited any anthropometric effects, whereas there would have been >90% power to detect a body mass index (BMI) effect of 0.5 kg/m(2) at P=0.01. Two "private" variants were also identified. c.335C>T/p.Thr112Met has been previously described and appears to be silent. A novel variant, c.260C>A/p.Ala87Asp, was observed in a subject with a BMI of 31.5 kg/m(2) (i.e., clinically obese) but not on direct assay of a further 3,525 subjects. This mutation was predicted to be deleterious and analysis using a cyclic AMP (cAMP) responsive luciferase reporter assay showed substantial loss of function of the mutant receptor. This population-based mutation scan of MC4R suggests that there is no severe MC4R mutation with high prevalence in the United Kingdom, but that obesity-causing MC4R mutation at 1 in 1,100 might represent one of the commonest autosomal dominant disorders in man.

Analysis of LDLR mutations in familial hypercholesterolemia patients in Greece by use of the NanoChip® Microelectronic Array Technology

BACKGROUND

Three mutations in the low density lipoprotein receptor (LDLR) gene account for 49% of familial hypercholesterolemia (FH) cases in Greece.

METHODS

We used the microelectronic array technology of the NanoChip Molecular Biology Workstation to develop a multiplex method to analyze these single-nucleotide polymorphisms (SNPs). Primer pairs amplified the region encompassing each SNP. The biotinylated PCR amplicon was electronically addressed to streptavidin-coated microarray sites. Allele-specific fluorescently labeled oligonucleotide reporters were designed and used for detection of wild-type and SNP sequences. Genotypes were compared to PCR-restriction fragment length polymorphism (PCR-RFLP).

RESULTS

We developed three monoplex assays (1 SNP/site) and an optimized multiplex assay (3SNPs/site). We performed 92 Greece II, 100 Genoa, and 98 Afrikaner-2 NanoChip monoplex assays (addressed to duplicate sites and analyzed separately). Of the 580 monoplex genotypings (290 samples), 579 agreed with RFLP. Duplicate sites of one sample were not in agreement with each other. Of the 580 multiplex genotypings, 576 agreed with the monoplex results. Duplicate sites of three samples were not in agreement with each other, indicating requirement for repetition upon which discrepancies were resolved.

CONCLUSIONS

The multiplex assay detects common LDLR mutations in Greek FH patients and can be extended to accommodate additional mutations.

TILLING: practical single-nucleotide mutation discovery

In the post-genomic sequencing era, an expanding portfolio of genomic technologies has been applied to the study of gene function. Reverse genetics approaches that provide targeted inactivation of genes identified by sequence analysis include TILLING (for Targeting Local Lesions IN Genomes). TILLING searches the genomes of mutagenized organisms for mutations in a chosen gene, typically single base-pair substitutions. This review covers practical aspects of the technology, ranging from building the mutagenized population to mutation discovery, and discusses possible improvements to current protocols and the impact of new genomic methods for mutation discovery in relation to the future of the TILLING approach.