Improved Characterization of Complex ß-Globin Gene Cluster Structural Variants Using Long-Read Sequencing.

Complex insertion-deletion (indel) events in the globin genes manifest in widely variable clinical phenotypes. Many are incompletely characterized because of a historic lack of efficient methods. A more complete assessment enables improved prediction of clinical impact, which guides emerging therapeutic choices. Current methods have limited capacity for breakpoint assignment and accurate assessment of mutation extent, especially in cases containing duplications or multiple deletions and insertions. Technology, such as long-read sequencing, holds promise for significant impact in the characterization of indel events because of read lengths that span large regions, resulting in improved resolution. Four known complex ß-globin gene cluster indel types were assessed using single-molecule, real-time sequencing technology and showed high correlation with previous reports, including the Caribbean locus control deletion (g.5,305,478_5,310,336del), a large ß-gene duplication containing the Hb S mutation (g.4,640,335_5,290,171dup with g.5,248,232T>A, c.20A>T; variant allele fraction, 64%), and two nested variants (double deletions with intervening inversion): the Indian Gγ(Aγδß)0-thalassemia (g.5,246,804-5,254,275del, g.5,254,276_5,269,600inv, and g.5,269,601_5,270,442del) and the Turkish/Macedonian (δß)0 thalassemia (g.5,235,064_5,236,652del, g.5,236,653_5,244,280inv, and g.5,244,281_5,255,766del). Our data confirm long-read sequencing as an efficient and accurate method to identify these clinically significant complex events. Limitations include high-complexity sample preparation requirements, which hinder routine use in clinical laboratories. Continued improvements in sample and data workflow processes are needed to accommodate volumes in a tertiary clinical laboratory.

DNA methylation profiling: comparison of genome-wide sequencing methods and the Infinium Human Methylation 450 Bead Chip.

AIMS: To compare the performance of four sequence-based and one microarray methods for DNA methylation profiling. METHODS: DNA from two cell lines were profiled by reduced representation bisulfite sequencing, methyl capture sequencing (SS-Meth Seq), NimbleGen SeqCapEpi CpGiant(Nimblegen MethSeq), methylated DNA immunoprecipitation (MeDIP) and the Human Methylation 450 Bead Chip (Meth450K). RESULTS & CONCLUSION: Despite differences in genome-wide coverage, high correlation and concordance were observed between different methods. Significant overlap of differentially methylated regions was identified between sequenced-based platforms. MeDIP provided the best coverage for the whole genome and gene body regions, while RRBS and Nimblegen MethSeq were superior for CpGs in CpG islands and promoters. Methylation analyses can be achieved by any of the five methods but understanding their differences may better address the research question being posed.

Linkage analysis of obesity phenotypes in pre- and post-menopausal women from a United States mid-western population.

BACKGROUND: Obesity has a strong genetic influence, with some variants showing stronger associations among women than men. Women are also more likely to distribute weight in the abdomen following menopause. We investigated whether genetic loci link with obesity-related phenotypes differently by menopausal status. METHODS: We performed univariate and bivariate linkage analysis for the phenotypes of body mass index (BMI), waist (W) and hip (H) circumferences (WC, HC), and WH ratio (WHR) separately among 172 pre-menopausal and 405 post-menopausal women from 90 multigenerational families using a genome scan with 403 microsatellite markers. Bivariate analysis used pair-wise combinations of obesity phenotypes to detect linkage at loci with pleiotropic effects for genetically correlated traits. BMI was adjusted in models of WC, HC and WHR. RESULTS: Pre-menopausal women, compared to post-menopausal women, had higher heritability for BMI (h2 = 94% versus h2 = 39%, respectively) and for HC (h2 = 99% versus h2 = 43%, respectively), and lower heritability for WC (h2 = 29% versus h2 = 61%, respectively) and for WHR (h2 = 39% versus h2 = 57%, respectively). Among pre-menopausal women, the strongest evidence for linkage was for the combination of BMI and HC traits at 3p26 (bivariate LOD = 3.65) and at 13q13-q14 (bivariate LOD = 3.59). Among post-menopausal women, the highest level of evidence for genetic linkage was for HC at 4p15.3 (univariate LOD = 2.70) and 14q13 (univariate LOD = 2.51). WC was not clearly linked to any locus. CONCLUSIONS: These results support a genetic basis for fat deposition that differs by menopausal status, and suggest that the same loci encode genes that influence general obesity (BMI) and HC, specifically, among pre-menopausal women. However, lower heritability among pre-menopausal women for WC and WHR suggests that pre-menopausal waist girth may be influenced to a greater extent by controllable environmental factors than post-menopausal waist girth. Possibly, targeted interventions for weight control among pre-menopausal women may prevent or attenuate post-menopausal abdominal weight deposition.

Evaluation of a new high-dimensional miRNA profiling platform.

BACKGROUND: MicroRNAs (miRNAs) are a class of approximately 22 nucleotide long, widely expressed RNA molecules that play important regulatory roles in eukaryotes. To investigate miRNA function, it is essential that methods to quantify their expression levels be available. METHODS: We evaluated a new miRNA profiling platform that utilizes Illumina's existing robust DASL chemistry as the basis for the assay. Using total RNA from five colon cancer patients and four cell lines, we evaluated the reproducibility of miRNA expression levels across replicates and with varying amounts of input RNA. The beta test version was comprised of 735 miRNA targets of Illumina's miRNA profiling application. RESULTS: Reproducibility between sample replicates within a plate was good (Spearman's correlation 0.91 to 0.98) as was the plate-to-plate reproducibility replicates run on different days (Spearman's correlation 0.84 to 0.98). To determine whether quality data could be obtained from a broad range of input RNA, data obtained from amounts ranging from 25 ng to 800 ng were compared to those obtained at 200 ng. No effect across the range of RNA input was observed. CONCLUSION: These results indicate that very small amounts of starting material are sufficient to allow sensitive miRNA profiling using the Illumina miRNA high-dimensional platform. Nonlinear biases were observed between replicates, indicating the need for abundance-dependent normalization. Overall, the performance characteristics of the Illumina miRNA profiling system were excellent.

Strong evidence of a genetic determinant for mammographic density, a major risk factor for breast cancer.

Increased mammographic density (MD), the proportion of dense tissue visible on a mammogram, is a strong risk factor for breast cancer, common in the population and clusters in families. We conducted the first genome-wide linkage scan to identify genes influencing MD. DNA was obtained from 889 relatives (756 women, 133 men) from 89 families. Percent MD was estimated on 618 (82%) female family members using a validated computer-assisted thresholding method. The genome-wide scan included 403 microsatellite DNA markers with an average spacing of 9 cM. Fine mapping of a region of chromosome 5p (5p13.1-5p15.1) was done using 21 additional closely spaced DNA markers. Linkage analyses were conducted to quantify the evidence for a gene responsible for MD across the genome. The maximum log odds for linkage (LOD) score from the genome-wide scan was on chromosome 5p (LOD = 2.9, supporting linkage by a factor of 10(2.9) or 794 to 1) with a 1-LOD interval spanning 28.6 cM. Two suggestive regions for linkage were also identified on chromosome 12 (LOD = 2.6, 1-LOD interval of 14.8 cM; and LOD = 2.5, 1-LOD interval of 17.2 cM). Finer mapping of the region surrounding the maximum LOD on chromosome 5p resulted in stronger and statistically significant evidence for linkage (LOD = 4.2) and a narrowed 1-LOD interval (13.4 cM). The putative locus on chromosome 5p is likely to account for up to 22% of variation in MD. Hence, 1 or more of the 45 candidate genes in this region could explain a large proportion of MD and, potentially, breast cancer.

Comparison of microsatellites versus single-nucleotide polymorphisms in a genome linkage screen for prostate cancer-susceptibility Loci.

Prostate cancer is one of the most common cancers among men and has long been recognized to occur in familial clusters. Brothers and sons of affected men have a 2-3-fold increased risk of developing prostate cancer. However, identification of genetic susceptibility loci for prostate cancer has been extremely difficult. Although the suggestion of linkage has been reported for many chromosomes, the most promising regions have been difficult to replicate. In this study, we compare genome linkage scans using microsatellites with those using single-nucleotide polymorphisms (SNPs), performed in 467 men with prostate cancer from 167 families. For the microsatellites, the ABI Prism Linkage Mapping Set version 2, with 402 microsatellite markers, was used, and, for the SNPs, the Early Access Affymetrix Mapping 10K array was used. Our results show that the presence of linkage disequilibrium (LD) among SNPs can lead to inflated LOD scores, and this seems to be an artifact due to the assumption of linkage equilibrium that is required by the current genetic-linkage software. After excluding SNPs with high LD, we found a number of new LOD-score peaks with values of at least 2.0 that were not found by the microsatellite markers: chromosome 8, with a maximum model-free LOD score of 2.2; chromosome 2, with a LOD score of 2.1; chromosome 6, with a LOD score of 4.2; and chromosome 12, with a LOD score of 3.9. The LOD scores for chromosomes 6 and 12 are difficult to interpret, because they occurred only at the extreme ends of the chromosomes. The greatest gain provided by the SNP markers was a large increase in the linkage information content, with an average information content of 61% for the SNPs, versus an average of 41% for the microsatellite markers. The strengths and weaknesses of microsatellite versus SNP markers are illustrated by the results of our genome linkage scans.