Precision medicine in the age of big data: The present and future role of large-scale unbiased sequencing in drug discovery and development.

High throughput molecular and functional profiling of patients is a key driver of precision medicine. DNA and RNA characterization has been enabled at unprecedented cost and scale through rapid, disruptive progress in sequencing technology, but challenges persist in data management and interpretation. We analyze the state-of-the-art of large-scale unbiased sequencing in drug discovery and development, including technology, application, ethical, regulatory, policy and commercial considerations, and discuss issues of LUS implementation in clinical and regulatory practice.

Genome-wide scan of obesity in the Old Order Amish.

To identify the genetic determinants of typical obesity, we performed a genome-wide scan of obesity-related traits using data from the Amish. Multipoint linkage analysis was performed using a variance components procedure on body mass index (BMI), waist circumference, percentage of body fat, and serum leptin concentrations. All 672 individuals were genotyped for 357 markers in 22 autosomes. We observed modest evidence for linkage, with the maximum log odds (lod) scores for linkage for these traits occurring on chromosomes 3p (percentage of body fat: lod = 1.61, near the peroxisome proliferator-activated receptor-alpha gene), 14q (waist: lod = 1.80), and 16p (leptin: lod = 1.72; BMI: lod = 1.68). We also tested for linkage to BMI-adjusted leptin concentrations and observed suggestive evidence for linkage on chromosome 10p (lod = 2.73), approximately 10-20 cM telomeric from obesity loci previously reported in French and German Caucasians. Two additional linkage signals for this trait were observed on chromosomes 7q (lod = 1.77, approximately 20 cM from the leptin gene) and 14q (lod = 2.47). Follow-up studies may be warranted to pursue some of these linkage signals, especially those detected near known obesity candidate genes, and those in regions coinciding with linkage signals reported previously.

A set of canine interrepeat sequence PCR markers for high-throughput genotyping.

One hundred and sixteen interspersed repetitive DNA sequence (IRS)-PCR markers have been developed and characterized from Canis familiaris for high-throughput filter-based genotyping. We present a detailed analysis of markers produced by amplification using primers directed to the conserved regions of the C. familiaris short interspersed nuclear element (Can-SINE). The majority of IRS-PCR markers developed were moderately to highly polymorphic with mean heterozygosity (HET) and polymorphism information content (PIC) values of approximately 0.6. The HET value for 22.3% of the markers exceeded 0.7. We also demonstrate that sequence variation of Can-SINEs between breeds is significant and also represents a rich source of polymorphisms. Mapping of 73 of the markers to the existing integrated linkage-radiation hybrid map enriches the map as well as establishes the utility of the markers. The significance and utility of this new class of IRS-PCR Can-SINE-based markers for high-throughput genotyping is discussed. This method can also be extended to other species that are currently map-poor but have a sufficiently high density of SINEs to allow IRS-PCR.

Diabetes in the Old Order Amish: characterization and heritability analysis of the Amish Family Diabetes Study.

OBJECTIVE: The Old Order Amish (OOA) are a genetically well-defined closed Caucasian founder population. The Amish Family Diabetes Study was initiated to identify susceptibility genes for type 2 diabetes. This article describes the genetic epidemiology of type 2 diabetes and related traits in this unique population. RESEARCH DESIGN AND METHODS: The study cohort comprised Amish probands with diabetes who were diagnosed between 35 and 65 years of age and their extended adult family members. We recruited 953 adults who represented 45 multigenerational families. Phenotypic characterization included anthropometry, blood pressure, diabetes status, lipid profile, and leptin levels. RESULTS: The mean age of study participants was 46 years, and the mean BMI was 26.9 kg/m2. Subjects with type 2 diabetes were older, more obese, and had higher insulin levels. The prevalence of diabetes in the OOA was approximately half that of the Caucasian individuals who participated in the Third National Health and Nutrition Examination Survey (95% CI 0.23-0.84). The prevalence of diabetes in the siblings of the diabetic probands was 26.5% compared with a prevalence of 7.0% in spouses (lambdaS = 3.28, 95% CI 1.58-6.80). The heritability of diabetes-related quantitative traits was substantial (13-70% for obesity-related traits, 10-42% for glucose levels, and 11-24% for insulin levels during the oral glucose tolerance test; P = 0.01 to <0.0001). CONCLUSIONS: Type 2 diabetes in the Amish has similar phenotypic features to that of the overall Caucasian population, although the prevalence in the Amish community is lower than that of the Caucasian population. There is significant familial clustering of type 2 diabetes and related traits. This unique family collection will be an excellent resource for investigating the genetic underpinnings of type 2 diabetes.

Genomewide search for type 2 diabetes susceptibility genes in four American populations.

Type 2 diabetes is a serious, genetically influenced disease for which no fully effective treatments are available. Identification of biochemical or regulatory pathways involved in the disease syndrome could lead to innovative therapeutic interventions. One way to identify such pathways is the genetic analysis of families with multiple affected members where disease predisposing genes are likely to be segregating. We undertook a genomewide screen (389-395 microsatellite markers) in samples of 835 white, 591 Mexican American, 229 black, and 128 Japanese American individuals collected as part of the American Diabetes Association's GENNID study. Multipoint nonparametric linkage analyses were performed with diabetes, and diabetes or impaired glucose homeostasis (IH). Linkage to diabetes or IH was detected near markers D5S1404 (map position 77 cM, LOD = 2.80), D12S853 (map position 82 cM, LOD = 2.81) and GATA172D05 (X-chromosome map position 130 cM, LOD = 2.99) in whites, near marker D3S2432 (map position 51 cM, LOD = 3.91) in Mexican Americans, and near marker D10S1412 (map position 14 cM, LOD = 2.39) in African Americans mainly collected in phase 1 of the study. Further analyses showed evidence for interactions between the chromosome 5 locus and region on chromosome 12 containing the MODY 3 gene (map position 132 cM) and between the X-chromosome locus and region near D12S853 (map position 82 cM) in whites. Although these results were not replicated in samples collected in phase 2 of the GENNID study, the region on chromosome 12 was replicated in samples from whites described by Bektas et al. (1999).

A genomewide analysis provides evidence for novel linkages in inflammatory bowel disease in a large European cohort.

Inflammatory bowel disease (IBD) is characterized by a chronic relapsing intestinal inflammation, typically starting in early adulthood. IBD is subdivided into two subtypes, on the basis of clinical and histologic features: Crohn disease and ulcerative colitis (UC). Previous genomewide searches identified regions harboring susceptibility loci on chromosomes 1, 3, 4, 7, 12, and 16. To expand our understanding of the genetic risk profile, we performed a 9-cM genomewide search for susceptibility loci in 268 families containing 353 affected sibling pairs. Previous linkages on chromosomes 12 and 16 were replicated, and the chromosome 4 linkage was extended in this sample. New suggestive evidence for autosomal linkages was observed on chromosomes 1, 6, 10, and 22, with LOD scores of 2.08, 2.07, 2.30, and 1.52, respectively. A maximum LOD score of 1.76 was observed on the X chromosome, for UC, which is consistent with the clinical association of IBD with Ullrich-Turner syndrome. The linkage finding on chromosome 6p is of interest, given the possible contribution of human leukocyte antigen and tumor necrosis-factor genes in IBD. This genomewide linkage scan, done with a large family cohort, has confirmed three previous IBD linkages and has provided evidence for five additional regions that may harbor IBD predisposition genes.

Genome scan for human obesity and linkage to markers in 20q13.

Obesity is a highly prevalent, multigenic trait that predicts increased morbidity and mortality. Here we report results from a genome scan based on 354 markers in 513 members of 92 nuclear families ascertained through extreme obesity and normal body weight. The average marker interval was approximately 10 cM. We examined four correlated obesity phenotypes, including the body-mass index (BMI) (both as a quantitative trait and as a discrete trait with a threshold of BMI > or /=30 kg/m2) and percentage of fat (both as a quantitative trait and as a discrete trait with a threshold of 40%) as assessed by bioelectrical impedance. In the initial stage of the genome scan, four markers in 20q gave positive evidence for linkage, which was consistent across most obesity phenotypes and analytic methods. After saturating 20q with additional markers (25 markers total) in an augmented sample of 713 members from 124 families, we found linkage to several markers in a region, 20q13, previously implicated in both human and animal studies. Three markers (D20S107, D20S211, and D20S149) in 20q13 had empirical P values (based on Monte Carlo simulations, which controlled for multiple testing) < or /=. 01 for single-point analysis. In addition, the parametric, affecteds-only analysis for D20S476 yielded a LOD score of 3.06 (P=. 00009), and the affected-sib-pair test yielded a LOD score of 3.17 (P=.000067). Multipoint analyses further strengthened and localized these findings. This region includes several plausible candidate genes for obesity. Our results suggest that one or more genes affecting obesity are located in 20q13.

Genetic analysis of inflammatory bowel disease in a large European cohort supports linkage to chromosomes 12 and 16.

BACKGROUND & AIMS: Inflammatory bowel disease (IBD) is a complex disorder of unknown etiology. Epidemiological investigations suggest a genetic basis for IBD. Recent genetic studies have identified several IBD linkages. The significance of these linkages will be determined by studies in large patient collections. The aim of this study was to replicate IBD linkages on chromosomes 12 and 16 in a large European cohort. METHODS: Three hundred fifty-nine affected sibling pairs from 274 kindreds were genotyped using microsatellite markers spanning chromosomes 12 and 16. Affection status of the sibling pairs was defined as Crohn's disease (CD) or ulcerative colitis (UC). RESULTS: Nonparametric statistical analyses showed linkage for both chromosomes. Two-point results for chromosome 12 peaked at D12S303 (logarithm of odds [LOD], 2.15; P = 0.003) for CD and at D12S75 (LOD, 0.92; P = 0.03) for UC. Multipoint analyses produced a peak LOD of 1.8 for CD. Chromosome 16 showed linkage for CD at marker D16S415 (LOD, 1.52; P = 0.007). Multipoint support peaked above markers D16S409 and D16S411 (LOD, 1.7). CONCLUSIONS: These data are consistent with linkage of IBD to chromosomes 12 and 16. The replication of genetic risk loci in a large independent family collection indicates important and common susceptibility genes in these regions and will facilitate identification of genes involved in IBD.

An autosomal genomic scan for loci linked to prediabetic phenotypes in Pima Indians.

Type 2 diabetes mellitus is a common chronic disease that is thought to have a substantial genetic basis. Identification of the genes responsible has been hampered by the complex nature of the syndrome. Abnormalities in insulin secretion and insulin action predict the development of type 2 diabetes and are, themselves, highly heritable traits. Since fewer genes may contribute to these precursors of type 2 diabetes than to the overall syndrome, such genes may be easier to identify. We, therefore, undertook an autosomal genomic scan to identify loci linked to prediabetic traits in Pima Indians, a population with a high prevalence of type 2 diabetes. 363 nondiabetic Pima Indians were genotyped at 516 polymorphic microsatellite markers on all 22 autosomes. Linkage analyses were performed using three methods (single-marker, nonparametric multipoint [MAPMAKER/SIBS], and variance components multipoint). These analyses provided evidence for linkage at several chromosomal regions, including 3q21-24 linked to fasting plasma insulin concentration and in vivo insulin action, 4p15-q12 linked to fasting plasma insulin concentration, 9q21 linked to 2-h insulin concentration during oral glucose tolerance testing, and 22q12-13 linked to fasting plasma glucose concentration. These results suggest loci that may harbor genes contributing to type 2 diabetes in Pima Indians. None of the linkages exceeded a LOD score of 3.6 (a 5% probability of occurring in a genome-wide scan). These findings must, therefore, be considered tentative until extended in this population or replicated in others.

Autosomal genomic scan for loci linked to obesity and energy metabolism in Pima Indians.

An autosomal genomic scan to search for linkage to obesity and energy metabolism was completed in Pima Indians, a population prone to obesity. Obesity was assessed by percent body fat (by hydrodensitometry) and fat distribution (the ratio of waist circumference to thigh circumference). Energy metabolism was measured in a respiratory chamber as 24-h metabolic rate, sleeping metabolic rate, and 24-h respiratory quotient (24RQ), an indicator of the ratio of carbohydrate oxidation to fat oxidation. Five hundred sixteen microsatellite markers with a median spacing of 6.4 cM were analyzed, in 362 siblings who had measurements of body composition and in 220 siblings who had measurements of energy metabolism. These comprised 451 sib pairs in 127 nuclear families, for linkage analysis to obesity, and 236 sib pairs in 82 nuclear families, for linkage analysis to energy metabolism. Pointwise and multipoint methods for regression of sib-pair differences in identity by descent, as well as a sibling-based variance-components method, were used to detect linkage. LOD scores >=2 were found at 11q21-q22, for percent body fat (LOD=2.1; P=.001), at 11q23-q24, for 24-h energy expenditure (LOD=2.0; P=.001), and at 1p31-p21 (LOD=2.0) and 20q11.2 (LOD=3.0; P=.0001), for 24RQ, by pointwise and multipoint analyses. With the variance-components method, the highest LOD score (LOD=2.3 P=.0006) was found at 18q21, for percent body fat, and at 1p31-p21 (LOD=2.8; P=.0003), for 24RQ. Possible candidate genes include LEPR (leptin receptor), at 1p31, and ASIP (agouti-signaling protein), at 20q11.2.

Familiality of physical and metabolic characteristics that predict the development of non-insulin-dependent diabetes mellitus in Pima Indians.

Susceptibility to non-insulin-dependent diabetes mellitus (NIDDM) is largely genetically determined. In Pima Indians, obesity, insulin resistance, and a low acute insulin response (AIR) to an intravenous glucose infusion are each predictors of the disease. To ascertain whether these phenotypes are genetically determined, we estimated their familiality in nondiabetic Pima Indians with a maximum-likelihood method. Percentage body fat (PFAT) was highly familial (h2 =.76), whereas waist/ thigh circumference ratio (W/T ratio) was not significantly familial after controlling for PFAT (h2 = .16). AIR was also highly familial (h2 = .80 at 10 min), even after controlling for PFAT and insulin action (h2 = .70). Insulin action at physiologic plasma insulin concentrations was familial (h2 = .61) but less so after controlling for PFAT and W/T ratio (h2 = .38). At maximally stimulating insulin concentrations, insulin action was familial (h2 = .45) and was less influenced by controlling for PFAT and W/T ratio (h2 = .49). We conclude that in Pima Indians (1) PFAT and AIR are highly familial traits, (2) central distribution of fat is not a familial trait when controlled for PFAT, (3) 38%-49% of the variance in insulin action, independent of the effect of obesity, is familial, and (4) PFAT, AIR, and insulin action are useful traits to study genetic susceptibility to NIDDM. Because genetic parameter estimates are applicable only to the populations from which they were estimated, it is important to determine whether these estimates of familialities in Pima Indians can be confirmed in other populations before the utility of these traits in searching for NIDDM susceptibility genes in those populations can be fully advocated.

Cryopreservation of embryos as a means of germ plasm conservation in sheep.

Embryos were collected from 4 lines of Targhee sheep between 1986 and 1990. The lines were selected for preweaning growth rate (Lines DH and HW) or for multiple births (Line HT); Line C served as an unselected control group. Estrus was synchronized using fluorogestone acetate-impregnated vaginal pessaries, and ewes were superovulated with FSH. Embryos at the morula or blastocyst stage were surgically recovered from mature ewes at Days 5 to 6 and were frozen following morphological evaluation. The overall average number of freezable embryos per collection was 2.9, and did not differ significantly among years or among lines. Of the embryos collected between 1986 and 1988, 92 were transferred to 53 recipients in 1989, producing 53 lambs. Survival rates were 60.9 and 47.8%, respectively, for embryos evaluated as good and fair after thawing. Good-quality blastocysts yielded the highest survival rate (64.4%). Analyses indicated no significant effects of line, developmental stage or embryo evaluation on the incidence of lambing. It was concluded that embryos of morula or blastocyst stage can be successfully frozen for extended periods. The data on embryo yield and survival following cryopreservation were used to calculate numbers of donors needed to preserve and reconstitute a population of specified size.

Evaluation of Australian merino and U.S. sheep breeds for growth and carcass traits.

One hundred twenty 4-mo-old wether lambs born to Targhee ewes and sired by six rams each from Merino (Finewool, FM and Strongwool, SM), Rambouillet (Dubois, DR and Texas, TR), and Targhee (T) breeds were randomly assigned to predetermined slaughter weight groups of 43, 48, 52, and 57 kg and evaluated for growth and carcass traits. Overall mean ADG and feed conversion rate (FC, kilograms of feed/kilogram of gain) were .28 and 6.4 kg, respectively; T grew the fastest (.31 kg) and FM grew the slowest (.23 kg). Targhee and SM had the best FC (6.2), whereas FM (6.8) had the poorest FC (P < .05). Overall mean backfat thickness (BT) and carcass fat (CF) were 4.8 mm and 25.6%, respectively. Targhee had the lowest (24%) and FM the highest (27.8%) CF percentage (P < .05). Differences (P < .05) were observed for BT and CF among slaughter weight groups; overall means for both traits gradually increased from the 43- to the 57-kg groups. Separate analysis of the 43- and the 48-kg groups indicated nonsignificant breed differences for feed traits, whereas significant differences still existed for CF. It was concluded that Merino strains grew more slowly, were less efficient in postweaning growth, and had higher carcass fat content than U.S. breeds at a constant slaughter weight; SM were more comparable to U.S. breeds than were FM. Merino-cross lambs should be slaughtered at lighter BW to avoid excessive carcass fat.