Germline mutation rates in young adults predict longevity and reproductive lifespan.

Ageing may be due to mutation accumulation across the lifespan, leading to tissue dysfunction, disease, and death. We tested whether germline autosomal mutation rates in young adults predict their remaining survival, and, for women, their reproductive lifespans. Age-adjusted mutation rates (AAMRs) in 61 women and 61 men from the Utah CEPH (Centre d'Etude du Polymorphisme Humain) families were determined. Age at death, cause of death, all-site cancer incidence, and reproductive histories were provided by the Utah Population Database, Utah Cancer Registry, and Utah Genetic Reference Project. Higher AAMRs were significantly associated with higher all-cause mortality in both sexes combined. Subjects in the top quartile of AAMRs experienced more than twice the mortality of bottom quartile subjects (hazard ratio [HR], 2.07; 95% confidence interval [CI], 1.21-3.56; p = 0.008; median survival difference = 4.7 years). Fertility analyses were restricted to women whose age at last birth (ALB) was ≥ 30 years, the age when fertility begins to decline. Women with higher AAMRs had significantly fewer live births and a younger ALB. Adult germline mutation accumulation rates are established in adolescence, and later menarche in women is associated with delayed mutation accumulation. We conclude that germline mutation rates in healthy young adults may provide a measure of both reproductive and systemic ageing. Puberty may induce the establishment of adult mutation accumulation rates, just when DNA repair systems begin their lifelong decline.

A hyperactivating proinflammatory RIPK2 allele associated with early-onset osteoarthritis.

Osteoarthritis (OA) is a common debilitating disease characterized by abnormal remodeling of the cartilage and bone of the articular joint. Ameliorating therapeutics are lacking due to limited understanding of the molecular pathways affecting disease initiation and progression. Notably, although a link between inflammation and overt OA is well established, the role of inflammation as a driver of disease occurrence is highly disputed. We analyzed a family with dominant inheritance of early-onset OA and found that affected individuals harbored a rare variant allele encoding a significant amino acid change (p.Asn104Asp) in the kinase domain of receptor interacting protein kinase 2 (RIPK2), which transduces signals from activated bacterial peptidoglycan sensors through the NF-κB pathway to generate a proinflammatory immune response. Functional analyses of RIPK2 activity in zebrafish embryos indicated that the variant RIPK2104Asp protein is hyperactive in its signaling capacity, with augmented ability to activate the innate immune response and the NF-κB pathway and to promote upregulation of OA-associated genes. Further we show a second allele of RIPK2 linked to an inflammatory disease associated with arthritis also has enhanced activity stimulating the NF-κB pathway. Our studies reveal for the first time the inflammatory response can function as a gatekeeper risk factor for OA.

Germline but not somatic de novo mutations are common in human congenital diaphragmatic hernia.

OBJECTIVES: Congenital diaphragmatic hernia (CDH) is a developmental defect of the diaphragm that causes high newborn morbidity and mortality. CDH is considered to be a multifactorial disease, with strong evidence implicating genetic factors. Although recent studies suggest the biological role of deleterious germline de novo variants, the effect of gene variants specific to the diaphragm remains unclear, and few single genes have been definitively implicated in human disease. METHODS: We performed genome sequencing on 16 individuals with CDH and their unaffected parents, including 10 diaphragmatic samples. RESULTS: We did not detect damaging somatic mutations in diaphragms, but identified germline heterozygous de novo functional mutations of 14 genes in nine patients. Although the majority of these genes are not known to be associated with CDH, one patient with CDH and cardiac anomalies harbored a frameshift mutation in NR2F2 (aka COUP-TFII), generating a premature truncation of the protein. This patient also carried a missense variant predicted to be damaging in XIRP2 (aka Myomaxin), a transcriptional target of MEF2A. Both NR2F2 and MEF2A map to chromosome 15q26, where recurring de novo deletions and unbalanced translocations have been observed in CDH. CONCLUSIONS: Somatic variants are not common in CDH. To our knowledge, this is the second case of a germline de novo frameshift mutation in NR2F2 in CDH. Since NR2F2 null mice exhibit a diaphragmatic defect, and XIRP2 is implicated in cardiac development, our data suggest the role of these two variants in the etiology of CDH, and possibly cardiac anomalies.

A defect in myoblast fusion underlies Carey-Fineman-Ziter syndrome.

Multinucleate cellular syncytial formation is a hallmark of skeletal muscle differentiation. Myomaker, encoded by Mymk (Tmem8c), is a well-conserved plasma membrane protein required for myoblast fusion to form multinucleated myotubes in mouse, chick, and zebrafish. Here, we report that autosomal recessive mutations in MYMK (OMIM 615345) cause Carey-Fineman-Ziter syndrome in humans (CFZS; OMIM 254940) by reducing but not eliminating MYMK function. We characterize MYMK-CFZS as a congenital myopathy with marked facial weakness and additional clinical and pathologic features that distinguish it from other congenital neuromuscular syndromes. We show that a heterologous cell fusion assay in vitro and allelic complementation experiments in mymk knockdown and mymkinsT/insT zebrafish in vivo can differentiate between MYMK wild type, hypomorphic and null alleles. Collectively, these data establish that MYMK activity is necessary for normal muscle development and maintenance in humans, and expand the spectrum of congenital myopathies to include cell-cell fusion deficits.

Trans-ethnic fine-mapping of genetic loci for body mass index in the diverse ancestral populations of the Population Architecture using Genomics and Epidemiology (PAGE) Study reveals evidence for multiple signals at established loci.

Most body mass index (BMI) genetic loci have been identified in studies of primarily European ancestries. The effect of these loci in other racial/ethnic groups is less clear. Thus, we aimed to characterize the generalizability of 170 established BMI variants, or their proxies, to diverse US populations and trans-ethnically fine-map 36 BMI loci using a sample of >102,000 adults of African, Hispanic/Latino, Asian, European and American Indian/Alaskan Native descent from the Population Architecture using Genomics and Epidemiology Study. We performed linear regression of the natural log of BMI (18.5-70 kg/m2) on the additive single nucleotide polymorphisms (SNPs) at BMI loci on the MetaboChip (Illumina, Inc.), adjusting for age, sex, population stratification, study site, or relatedness. We then performed fixed-effect meta-analyses and a Bayesian trans-ethnic meta-analysis to empirically cluster by allele frequency differences. Finally, we approximated conditional and joint associations to test for the presence of secondary signals. We noted directional consistency with the previously reported risk alleles beyond what would have been expected by chance (binomial p < 0.05). Nearly, a quarter of the previously described BMI index SNPs and 29 of 36 densely-genotyped BMI loci on the MetaboChip replicated/generalized in trans-ethnic analyses. We observed multiple signals at nine loci, including the description of seven loci with novel multiple signals. This study supports the generalization of most common genetic loci to diverse ancestral populations and emphasizes the importance of dense multiethnic genomic data in refining the functional variation at genetic loci of interest and describing several loci with multiple underlying genetic variants.

Genome-wide analyses identify common variants associated with macular telangiectasia type 2.

Idiopathic juxtafoveal retinal telangiectasis type 2 (macular telangiectasia type 2; MacTel) is a rare neurovascular degenerative retinal disease. To identify genetic susceptibility loci for MacTel, we performed a genome-wide association study (GWAS) with 476 cases and 1,733 controls of European ancestry. Genome-wide significant associations (P < 5 × 10-8) were identified at three independent loci (rs73171800 at 5q14.3, P = 7.74 × 10-17; rs715 at 2q34, P = 9.97 × 10-14; rs477992 at 1p12, P = 2.60 × 10-12) and then replicated (P < 0.01) in an independent cohort of 172 cases and 1,134 controls. The 5q14.3 locus is known to associate with variation in retinal vascular diameter, and the 2q34 and 1p12 loci have been implicated in the glycine/serine metabolic pathway. We subsequently found significant differences in blood serum levels of glycine (P = 4.04 × 10-6) and serine (P = 2.48 × 10-4) between MacTel cases and controls.

A Rare Variant in CACNA1D Segregates with 7 Bipolar I Disorder Cases in a Large Pedigree.

Whole-genome sequencing was performed on 3 bipolar I disorder (BPI) cases from a multiplex pedigree of European ancestry with 7 BPI cases. Within CACNA1D, a gene implicated by genome-wide association studies, a G to C nucleotide transversion at 53,835,340 base pairs (bps) was found predicting the substitution of proline for alanine at amino acid position 1751 (A1751P). Using Sanger sequencing, the DNA variant was shown to co-segregate with the remaining 4 BPI cases within the pedigree. A high-resolution DNA denaturing curve method was then used to screen for the presence of the A1751P change in 4,150 BPI cases from the NIMH Genetics Initiative. The A1751P variant was found in 4 BPI cases. A second variant within exon 43, a C to T nucleotide transition, was found in 1 case at 53,835,355 bps, predicting the substitution of tryptophan for arginine at amino acid position 1771 (R1771W). In the NHLBI Exome Sequencing Project database, the heterozygous A1751P variant was present in 3 of 4,300 subjects of European ancestry, and the R1771W change was not present in any subject. Given the rarity of these variants, large-scale case/control rare variant sequencing studies will be required for definitive conclusions.

Demonstration of Protein-Based Human Identification Using the Hair Shaft Proteome.

Human identification from biological material is largely dependent on the ability to characterize genetic polymorphisms in DNA. Unfortunately, DNA can degrade in the environment, sometimes below the level at which it can be amplified by PCR. Protein however is chemically more robust than DNA and can persist for longer periods. Protein also contains genetic variation in the form of single amino acid polymorphisms. These can be used to infer the status of non-synonymous single nucleotide polymorphism alleles. To demonstrate this, we used mass spectrometry-based shotgun proteomics to characterize hair shaft proteins in 66 European-American subjects. A total of 596 single nucleotide polymorphism alleles were correctly imputed in 32 loci from 22 genes of subjects' DNA and directly validated using Sanger sequencing. Estimates of the probability of resulting individual non-synonymous single nucleotide polymorphism allelic profiles in the European population, using the product rule, resulted in a maximum power of discrimination of 1 in 12,500. Imputed non-synonymous single nucleotide polymorphism profiles from European-American subjects were considerably less frequent in the African population (maximum likelihood ratio = 11,000). The converse was true for hair shafts collected from an additional 10 subjects with African ancestry, where some profiles were more frequent in the African population. Genetically variant peptides were also identified in hair shaft datasets from six archaeological skeletal remains (up to 260 years old). This study demonstrates that quantifiable measures of identity discrimination and biogeographic background can be obtained from detecting genetically variant peptides in hair shaft protein, including hair from bioarchaeological contexts.

Exome analysis of a family with Wolff-Parkinson-White syndrome identifies a novel disease locus.

Wolff-Parkinson-White (WPW) syndrome is a common cause of supraventricular tachycardia that carries a risk of sudden cardiac death. To date, mutations in only one gene, PRKAG2, which encodes the 5'-AMP-activated protein kinase subunit γ-2, have been identified as causative for WPW. DNA samples from five members of a family with WPW were analyzed by exome sequencing. We applied recently designed prioritization strategies (VAAST/pedigree VAAST) coupled with an ontology-based algorithm (Phevor) that reduced the number of potentially damaging variants to 10: a variant in KCNE2 previously associated with Long QT syndrome was also identified. Of these 11 variants, only MYH6 p.E1885K segregated with the WPW phenotype in all affected individuals and was absent in 10 unaffected family members. This variant was predicted to be damaging by in silico methods and is not present in the 1,000 genome and NHLBI exome sequencing project databases. Screening of a replication cohort of 47 unrelated WPW patients did not identify other likely causative variants in PRKAG2 or MYH6. MYH6 variants have been identified in patients with atrial septal defects, cardiomyopathies, and sick sinus syndrome. Our data highlight the pleiotropic nature of phenotypes associated with defects in this gene.

Identification of rare DNA sequence variants in high-risk autism families and their prevalence in a large case/control population.

BACKGROUND: Genetics clearly plays a major role in the etiology of autism spectrum disorders (ASDs), but studies to date are only beginning to characterize the causal genetic variants responsible. Until recently, studies using multiple extended multi-generation families to identify ASD risk genes had not been undertaken. METHODS: We identified haplotypes shared among individuals with ASDs in large multiplex families, followed by targeted DNA capture and sequencing to identify potential causal variants. We also assayed the prevalence of the identified variants in a large ASD case/control population. RESULTS: We identified 584 non-conservative missense, nonsense, frameshift and splice site variants that might predispose to autism in our high-risk families. Eleven of these variants were observed to have odds ratios greater than 1.5 in a set of 1,541 unrelated children with autism and 5,785 controls. Three variants, in the RAB11FIP5, ABP1, and JMJD7-PLA2G4B genes, each were observed in a single case and not in any controls. These variants also were not seen in public sequence databases, suggesting that they may be rare causal ASD variants. Twenty-eight additional rare variants were observed only in high-risk ASD families. Collectively, these 39 variants identify 36 genes as ASD risk genes. Segregation of sequence variants and of copy number variants previously detected in these families reveals a complex pattern, with only a RAB11FIP5 variant segregating to all affected individuals in one two-generation pedigree. Some affected individuals were found to have multiple potential risk alleles, including sequence variants and copy number variants (CNVs), suggesting that the high incidence of autism in these families could be best explained by variants at multiple loci. CONCLUSIONS: Our study is the first to use haplotype sharing to identify familial ASD risk loci. In total, we identified 39 variants in 36 genes that may confer a genetic risk of developing autism. The observation of 11 of these variants in unrelated ASD cases further supports their role as ASD risk variants.

Fine Mapping and Identification of BMI Loci in African Americans.

Genome-wide association studies (GWASs) primarily performed in European-ancestry (EA) populations have identified numerous loci associated with body mass index (BMI). However, it is still unclear whether these GWAS loci can be generalized to other ethnic groups, such as African Americans (AAs). Furthermore, the putative functional variant or variants in these loci mostly remain under investigation. The overall lower linkage disequilibrium in AA compared to EA populations provides the opportunity to narrow in or fine-map these BMI-related loci. Therefore, we used the Metabochip to densely genotype and evaluate 21 BMI GWAS loci identified in EA studies in 29,151 AAs from the Population Architecture using Genomics and Epidemiology (PAGE) study. Eight of the 21 loci (SEC16B, TMEM18, ETV5, GNPDA2, TFAP2B, BDNF, FTO, and MC4R) were found to be associated with BMI in AAs at 5.8 × 10(-5). Within seven out of these eight loci, we found that, on average, a substantially smaller number of variants was correlated (r(2) > 0.5) with the most significant SNP in AA than in EA populations (16 versus 55). Conditional analyses revealed GNPDA2 harboring a potential additional independent signal. Moreover, Metabochip-wide discovery analyses revealed two BMI-related loci, BRE (rs116612809, p = 3.6 × 10(-8)) and DHX34 (rs4802349, p = 1.2 × 10(-7)), which were significant when adjustment was made for the total number of SNPs tested across the chip. These results demonstrate that fine mapping in AAs is a powerful approach for both narrowing in on the underlying causal variants in known loci and discovering BMI-related loci.

A systematic mapping approach of 16q12.2/FTO and BMI in more than 20,000 African Americans narrows in on the underlying functional variation: results from the Population Architecture using Genomics and Epidemiology (PAGE) study.

Genetic variants in intron 1 of the fat mass- and obesity-associated (FTO) gene have been consistently associated with body mass index (BMI) in Europeans. However, follow-up studies in African Americans (AA) have shown no support for some of the most consistently BMI-associated FTO index single nucleotide polymorphisms (SNPs). This is most likely explained by different race-specific linkage disequilibrium (LD) patterns and lower correlation overall in AA, which provides the opportunity to fine-map this region and narrow in on the functional variant. To comprehensively explore the 16q12.2/FTO locus and to search for second independent signals in the broader region, we fine-mapped a 646-kb region, encompassing the large FTO gene and the flanking gene RPGRIP1L by investigating a total of 3,756 variants (1,529 genotyped and 2,227 imputed variants) in 20,488 AAs across five studies. We observed associations between BMI and variants in the known FTO intron 1 locus: the SNP with the most significant p-value, rs56137030 (8.3 × 10(-6)) had not been highlighted in previous studies. While rs56137030was correlated at r(2)>0.5 with 103 SNPs in Europeans (including the GWAS index SNPs), this number was reduced to 28 SNPs in AA. Among rs56137030 and the 28 correlated SNPs, six were located within candidate intronic regulatory elements, including rs1421085, for which we predicted allele-specific binding affinity for the transcription factor CUX1, which has recently been implicated in the regulation of FTO. We did not find strong evidence for a second independent signal in the broader region. In summary, this large fine-mapping study in AA has substantially reduced the number of common alleles that are likely to be functional candidates of the known FTO locus. Importantly our study demonstrated that comprehensive fine-mapping in AA provides a powerful approach to narrow in on the functional candidate(s) underlying the initial GWAS findings in European populations.

Identification of rare recurrent copy number variants in high-risk autism families and their prevalence in a large ASD population.

Structural variation is thought to play a major etiological role in the development of autism spectrum disorders (ASDs), and numerous studies documenting the relevance of copy number variants (CNVs) in ASD have been published since 2006. To determine if large ASD families harbor high-impact CNVs that may have broader impact in the general ASD population, we used the Affymetrix genome-wide human SNP array 6.0 to identify 153 putative autism-specific CNVs present in 55 individuals with ASD from 9 multiplex ASD pedigrees. To evaluate the actual prevalence of these CNVs as well as 185 CNVs reportedly associated with ASD from published studies many of which are insufficiently powered, we designed a custom Illumina array and used it to interrogate these CNVs in 3,000 ASD cases and 6,000 controls. Additional single nucleotide variants (SNVs) on the array identified 25 CNVs that we did not detect in our family studies at the standard SNP array resolution. After molecular validation, our results demonstrated that 15 CNVs identified in high-risk ASD families also were found in two or more ASD cases with odds ratios greater than 2.0, strengthening their support as ASD risk variants. In addition, of the 25 CNVs identified using SNV probes on our custom array, 9 also had odds ratios greater than 2.0, suggesting that these CNVs also are ASD risk variants. Eighteen of the validated CNVs have not been reported previously in individuals with ASD and three have only been observed once. Finally, we confirmed the association of 31 of 185 published ASD-associated CNVs in our dataset with odds ratios greater than 2.0, suggesting they may be of clinical relevance in the evaluation of children with ASDs. Taken together, these data provide strong support for the existence and application of high-impact CNVs in the clinical genetic evaluation of children with ASD.

A family-based paradigm to identify candidate chromosomal regions for isolated congenital diaphragmatic hernia.

Congenital diaphragmatic hernia (CDH) is a developmental defect of the diaphragm that causes high newborn mortality. Isolated or non-syndromic CDH is considered a multifactorial disease, with strong evidence implicating genetic factors. As low heritability has been reported in isolated CDH, family-based genetic methods have yet to identify the genetic factors associated with the defect. Using the Utah Population Database, we identified distantly related patients from several extended families with a high incidence of isolated CDH. Using high-density genotyping, seven patients were analyzed by homozygosity exclusion rare allele mapping (HERAM) and phased haplotype sharing (HapShare), two methods we developed to map shared chromosome regions. Our patient cohort shared three regions not previously associated with CDH, that is, 2q11.2-q12.1, 4p13 and 7q11.2, and two regions previously involved in CDH, that is, 8p23.1 and 15q26.2. The latter regions contain GATA4 and NR2F2, two genes implicated in diaphragm formation in mice. Interestingly, three patients shared the 8p23.1 locus and one of them also harbored the 15q26.2 segment. No coding variants were identified in GATA4 or NR2F2, but a rare shared variant was found in intron 1 of GATA4. This work shows the role of heritability in isolated CDH. Our family-based strategy uncovers new chromosomal regions possibly associated with disease, and suggests that non-coding variants of GATA4 and NR2F2 may contribute to the development of isolated CDH. This approach could speed up the discovery of the genes and regulatory elements causing multifactorial diseases, such as isolated CDH.

Exome analysis of a family with pleiotropic congenital heart disease.

BACKGROUND: A number of single gene defects have been identified in patients with isolated or nonsyndromic congenital heart defects (CHDs). However, due to significant genetic heterogeneity, candidate gene approaches have had limited success in finding high-risk alleles in most cases. The purpose of this study was to use exome sequencing to identify high-risk gene variants in a family with highly penetrant pleiotropic CHD. METHODS AND RESULTS: DNA samples from 2 members of a family with diverse CHD were analyzed by exome sequencing. Variants were filtered to eliminate common variants and sequencing artifacts and then prioritized based on the predicted effect of the variant and on gene function. The remainder of the family was screened using polymerase chain reaction, high-resolution melting analysis, and DNA sequencing to evaluate variant segregation. After filtering, >2000 rare variants (including single nucleotide substitutions and indels) were shared by the 2 individuals. Of these, 46 were nonsynonymous, 3 were predicted to alter splicing, and 6 resulted in a frameshift. Prioritization reduced the number of variants potentially involved in CHD to 18. None of the variants completely segregated with CHD in the kindred. However, 1 variant, Myh6 Ala290Pro, was identified in all but 1 affected individual. This variant was previously identified in a patient with tricuspid atresia and large secundum atrial septal defect. CONCLUSIONS: It is likely that next-generation sequencing will become the method of choice for unraveling the complex genetics of CHD, but information gained by analysis of transmission through families will be crucial.

C21orf91 genotypes correlate with herpes simplex labialis (cold sore) frequency: description of a cold sore susceptibility gene.

BACKGROUND: Herpes simplex virus type 1 (HSV-1) infects >70% of the United States population. We identified a 3-megabase region on human chromosome 21 containing 6 candidate genes associated with herpes simplex labialis (HSL, "cold sores"). METHODS: We conducted single nucleotide polymorphism (SNP) scans of the chromosome 21 region to define which of 6 possible candidate genes were associated with cold sore frequency. We obtained the annual HSL frequency for 355 HSV-1 seropositive individuals and determined the individual genotypes by SNPlex for linkage analysis and parental transmission disequilibrium testing (ParenTDT). RESULTS: Two-point linkage analysis showed positive linkage between cold sore frequency and 2 SNPs within the C21orf91 region, 1 of which is nonsynonymous. ParenTDT analysis revealed a strong association between another C21orf91 SNP, predicted to lie in the 3' untranslated region, and frequent HSL (P = .0047). C21orf 91 is a predicted open reading frame of unknown function that encodes a cytosolic protein. CONCLUSIONS: We evaluated candidate genes in the cold sore susceptibility region using fine mapping with 45 SNP markers. 2 complementary techniques identified C21orf91 as a gene of interest for susceptibility to HSL. We propose that C21orf91 be designated the Cold Sore Susceptibility Gene-1 (CSSG1).

Carriers of rare missense variants in IFIH1 are protected from psoriasis.

Testing of ∼25,000 putative functional single-nucleotide polymorphisms (SNPs) across the human genome in a genetic association study has identified three psoriasis genes, IL12B, IL23R, and IL13. We now report evidence for the association of psoriasis risk with missense SNPs in the interferon induced with helicase C domain 1 gene (IFIH1). The rare alleles of two independent SNPs were associated with decreased risk of psoriasis--rs35667974 (Ile923Val): odds ratio (OR) for minor allele carriers is 0.43, P=2.36 × 10(-5) (2,098 cases vs. 1,748 controls); and rs10930046 (His460Arg): OR for minor allele carriers is 0.51, P=6.47 × 10(-4) (2,098 cases vs. 1,744 controls). Compared to noncarriers, carriers of the 923Val and/or 460Arg variants were protected from psoriasis (OR=0.46, P=5.56 × 10(-8)). To our knowledge, these results suggest that IFIH1 is a previously unreported psoriasis gene.

Colorectal adenomas and cancer link to chromosome 13q22.1-13q31.3 in a large family with excess colorectal cancer.

BACKGROUND: Colorectal cancer is the fourth most common type of cancer and the second most common cause of cancer death. Fewer than 5% of colon cancers arise in the presence of a clear hereditary cancer condition; however, current estimates suggest that an additional 15-25% of colorectal cancers arise on the basis of unknown inherited factors. AIM: To identify additional genetic factors responsible for colon cancer. METHODS: A large kindred with excess colorectal cancer was identified through the Utah Population Database and evaluated clinically and genetically for inherited susceptibility. RESULTS: A major genetic locus segregating with colonic polyps and cancer in this kindred was identified on chromosome 13q with a non-parametric linkage score of 24 (LOD score of 2.99 and p=0.001). The genetic region spans 21 Mbp and contains 27 RefSeq genes. Sequencing of all candidate genes in this region failed to identify a clearly deleterious mutation; however, polymorphisms segregating with the phenotype were identified. Chromosome 13q is commonly gained and overexpressed in colon cancers and correlates with metastasis, suggesting the presence of an important cancer progression gene. Evaluation of tumours from the kindred revealed a gain of 13q as well. CONCLUSIONS: This identified region may contain a novel gene responsible for colon cancer progression in a significant proportion of sporadic cancers. Identification of the precise gene and causative genetic change in the kindred will be an important next step to understanding cancer progression and metastasis.

(1)H nuclear magnetic resonance metabolomics analysis identifies novel urinary biomarkers for lung function.

Chronic obstructive pulmonary disease (COPD), characterized by chronic airflow limitation, is a serious public health concern. In this study, we used proton nuclear magnetic resonance ((1)H NMR) spectroscopy to identify and quantify metabolites associated with lung function in COPD. Plasma and urine were collected from 197 adults with COPD and from 195 without COPD. Samples were assayed using a 600 MHz NMR spectrometer, and the resulting spectra were analyzed against quantitative spirometric measures of lung function. After correcting for false discoveries and adjusting for covariates (sex, age, smoking) several spectral regions in urine were found to be significantly associated with baseline lung function. These regions correspond to the metabolites trigonelline, hippurate and formate. Concentrations of each metabolite, standardized to urinary creatinine, were associated with baseline lung function (minimum p-value = 0.0002 for trigonelline). No significant associations were found with plasma metabolites. Urinary hippurate and formate are often related to gut microflora. This could suggest that the microbiome varies between individuals with different lung function. Alternatively, the associated metabolites may reflect lifestyle differences affecting overall health. Our results will require replication and validation, but demonstrate the utility of NMR metabolomics as a screening tool for identifying novel biomarkers of pulmonary outcomes.