Case Report: Whole-exome analysis of a child with polycystic kidney disease and ventriculomegaly.

Autosomal recessive polycystic kidney disease (ARPKD) is an inherited ciliopathy leading to progressive kidney and liver disease. Biallelic mutations in the PKHD1 gene underlie this condition. We describe a child with bilaterally enlarged cystic kidneys, portal hypertension, and cerebral ventriculomegaly. Molecular genetic investigations using whole-exome sequencing and confirmation using Sanger sequencing revealed a homozygous pathogenic mutation in PKHD1 underlying the clinical phenotype of ARPKD. Whole-exome data analysis was used to search for additional rare variants in additional ciliopathy genes that may have contributed to the unusual brain phenotype. Aside from a rare hypomorphic allele in MKS1, no other pathogenic variants were detected. We conclude that the homozygous pathogenic mutation in PKHD1 underlies the ciliopathy phenotype in this patient.

Bases adjacent to mononucleotide repeats show an increased single nucleotide polymorphism frequency in the human genome.

Mononucleotide repeats (MNRs) are abundant in eukaryotic genomes and exhibit a high degree of length variability due to insertion and deletion events. However, the relationship between these repeats and mutation rates in surrounding sequences has not been systematically investigated. We have analyzed the frequency of single nucleotide polymorphisms (SNPs) at positions close to and within MNRs in the human genome. Overall, we find a 2- to 4-fold increase in the SNP frequency at positions immediately adjacent to the boundaries of MNRs, relative to that at more distant bases. This relationship exhibits a strong asymmetry between 3' and 5' ends of repeat tracts and is dependent upon the repeat motif, length and orientation of surrounding repeats. Our analysis suggests that the incorporation or exclusion of bases adjacent to the boundary of the repeat through substitutions, in which these nucleotides mutate towards or away from the base present within the repeat, respectively, may be another mechanism by which MNRs expand and contract in the human genome.

Identification of the pathogenic pathways in osteoarthritic hip cartilage: commonality and discord between hip and knee OA.

OBJECTIVE: To define for the first time the transcriptomes of normal and end-stage osteoarthritis (OA) hip cartilage. MATERIALS AND METHODS: RNA was isolated from cartilage within 2h of joint replacement surgery. Gene expression was analyzed using Agilent GeneSpring GX 11 following hybridization to Illumina Human HT-12 V3 microarrays. Real-time reverse-transcription polymerase chain reaction (RT-PCR) was used to validate the expression of six genes identified by microarray as differentially expressed. Gene Set Enrichment Analysis (GSEA) and Ingenuity Pathway Analysis (IPA) were used to investigate enriched functions or canonical pathways amongst differentially expressed genes respectively. RESULTS: In total we identified 998 differentially expressed genes (fold change ≥ ±1.5, P-value ≤ 0.01) between neck of femur fracture (NOF) (n = 10) and OA hip (n = 9) patient cartilage. These differentially expressed genes were enriched within 71 canonical pathways. A comparison between a comparable knee dataset(20) only identified 229 genes similarly differentially expressed although remarkably 34 canonical pathways overlapped between experiments. CONCLUSIONS: This study is the first to report a comprehensive gene expression analysis of human hip OA cartilage compared to control (NOF) cartilage at the whole-genome level. Our differential gene expression dataset shows excellent correlation with similar defined studies using comparable tissue but reveals discord between hip and knee OA at the individual gene status but with commonality with regards the molecular pathways involved.

A de novo paradigm for male infertility.

De novo mutations are known to play a prominent role in sporadic disorders with reduced fitness. We hypothesize that de novo mutations play an important role in severe male infertility and explain a portion of the genetic causes of this understudied disorder. To test this hypothesis, we utilize trio-based exome sequencing in a cohort of 185 infertile males and their unaffected parents. Following a systematic analysis, 29 of 145 rare (MAF < 0.1%) protein-altering de novo mutations are classified as possibly causative of the male infertility phenotype. We observed a significant enrichment of loss-of-function de novo mutations in loss-of-function-intolerant genes (p-value = 1.00 × 10-5) in infertile men compared to controls. Additionally, we detected a significant increase in predicted pathogenic de novo missense mutations affecting missense-intolerant genes (p-value = 5.01 × 10-4) in contrast to predicted benign de novo mutations. One gene we identify, RBM5, is an essential regulator of male germ cell pre-mRNA splicing and has been previously implicated in male infertility in mice. In a follow-up study, 6 rare pathogenic missense mutations affecting this gene are observed in a cohort of 2,506 infertile patients, whilst we find no such mutations in a cohort of 5,784 fertile men (p-value = 0.03). Our results provide evidence for the role of de novo mutations in severe male infertility and point to new candidate genes affecting fertility.

Germ-line mutations of TP53 in Li-Fraumeni families: an extended study of 39 families.

We have previously reported on the analysis of TP53 coding mutations in 12 classic Li-Fraumeni syndrome (LFS) families plus 9 families that were Li-Fraumeni-like (LFL) families (J. M. Birch et al., Cancer Res., 54: 1298-1304, 1994). Mutations were found in 6 of 12 LFS families and in 1 of 9 LFL families. We have now extended these studies to include an additional nine LFS and nine LFL families, and TP53 mutations have been detected in eight of nine LFS families and in three of nine LFL families. Six of the new mutations described here are the same as those previously identified in other Li-Fraumeni families and are missense mutations at codons 245, 248, and 273 (in two families); a nonsense mutation at codon 209; and a mutation at the splice donor site in exon 4. The other five mutations are novel germ-line mutations and include missense mutations at codons 136 and 344, a 2-bp deletion within codon 191, a splice acceptor mutation in intron 3, and a 167-bp deletion of part of exon 1 and intron 1. In addition, we have detected a codon 175 mutation in a family previously reported as TP53 negative. To summarize all of the data from the families we have studied in this and our previous report (J. M. Birch et al., Cancer Res., 54: 1298-1304, 1994), mutations have been detected in 15 of 21 LFS families (71%) and in 4 of 18 LFL families (22%). These figures are somewhat higher than those previously reported by us and others for the frequency of TP53 mutations in LFS and LFL families. This could reflect our analysis of all 11 exons of TP53, including noncoding regions, as well as the use of direct sequencing rather than other less-sensitive mutation detection methods.

The role of the RAS pathway in iAMP21-ALL.

Intrachromosomal amplification of chromosome 21 (iAMP21) identifies a high-risk subtype of acute lymphoblastic leukaemia (ALL), requiring intensive treatment to reduce their relapse risk. Improved understanding of the genomic landscape of iAMP21-ALL will ascertain whether these patients may benefit from targeted therapy. We performed whole-exome sequencing of eight iAMP21-ALL samples. The mutation rate was dramatically disparate between cases (average 24.9, range 5-51) and a large number of novel variants were identified, including frequent mutation of the RAS/MEK/ERK pathway. Targeted sequencing of a larger cohort revealed that 60% (25/42) of diagnostic iAMP21-ALL samples harboured 42 distinct RAS pathway mutations. High sequencing coverage demonstrated heterogeneity in the form of multiple RAS pathway mutations within the same sample and diverse variant allele frequencies (VAFs) (2-52%), similar to other subtypes of ALL. Constitutive RAS pathway activation was observed in iAMP21 samples that harboured mutations in the predominant clone (⩾35% VAF). Viable iAMP21 cells from primary xenografts showed reduced viability in response to the MEK1/2 inhibitor, selumetinib, in vitro. As clonal (⩾35% VAF) mutations were detected in 26% (11/42) of iAMP21-ALL, this evidence of response to RAS pathway inhibitors may offer the possibility to introduce targeted therapy to improve therapeutic efficacy in these high-risk patients.

Rare variant (Glu to Lys) in the leucine zipper region of CHOP (GADD153).

CHOP (GADD153) has been shown to be a dominant negative inhibitor of specific transcription factors. Direct sequencing of the gene, amplified from the DNA of a Li-Fraumeni family index case (liposarcoma, breast cancer) revealed a constitutional variant within the coding region. This alteration, though not responsible for the Li-Fraumeni phenotype, resulted in a glutamic acid to lysine switch within the leucine zipper domain, at a residue conserved between CHOP and its potential target molecules and between the human and hamster sequences. The variant created a Taq I restriction fragment length polymorphism (RFLP) facilitating screening. Analysis of 159 breast tumour DNA samples detected two encoding variant alleles (tumour and constitutional DNA).

An immediate early human gene encodes an Id-like helix-loop-helix protein and is regulated by protein kinase C activation in diverse cell types.

Transcription factors characterized by the presence of a helix-loop-helix (HLH) domain play a central role in the regulation of cell growth/differentiation and tumorigenesis. We report here the cDNA sequence of a human early-response gene, designated HLH 1R21, encoding a 15-kDa HLH protein that lacks a basic, DNA-binding domain and which by a number of criteria appears to be the human homologue of mouse HLH 462. Like its murine counterpart, HLH 1R21 protein functions as an Id (inhibitor of DNA binding) transcription factor by inhibiting the binding of E2A-containing protein complexes to muscle creatine kinase E-box enhancer oligonucleotide in vitro. However HLH 1R21 does not inhibit the binding of HLH Max protein to a Max-binding oligonucleotide in vitro, indicating that it has limited promiscuity in its ability to antagonize the function of other HLH transcription factors. In addition, HLH 1R21 mRNA transcripts are regulated by phorbol ester treatment of a diverse range of human cell lines and, when overexpressed in mouse NIH3T3 cells, HLH 1R21 induces a morphologically transformed phenotype.

A germline 2.35 kb deletion of p53 genomic DNA creating a specific loss of the oligomerization domain inherited in a Li-Fraumeni syndrome family.

The primary genetic cancer predisposing event in many Li-Fraumeni syndrome families is a germline mutation in the p53 gene. We describe an extended Li-Fraumeni family with a germline mutation in the p53 gene involving a deletion of exon 10. The mutation is a 2.35 kilobase intragenic deletion encompassing exon 10, which results in the specific loss of the entire p53 oligomerization domain. This mutation segregates with the cancer phenotype. A lymphoblastoid cell line developed from a mutation carrier shows accumulation of mutant p53 protein by immunoblotting. However, tumor tissues from two affected carriers are negative by immunohistochemical staining. A major structural alteration specifically involving the oligomerization domain of a germline p53 gene has not been previously described and occurs in a region rarely mutated in sporadic tumors. The oligomerization domain is dispensable for many wild-type p53 functions, including transactivation, sequence-specific DNA binding, and suppression of oncogenic transformation. However, the domain appears to be required for transcriptional repression, and DNA strand reassociation. The identification of this mutation in an LFS family may yield insights into the importance of the oligomerization domain for suppressor function of the p53 tumor suppressor gene.

Frequency of rare recessive mutations in unexplained late onset cerebellar ataxia.

Sporadic late onset cerebellar ataxia is a well-described clinical presentation with a broad differential diagnosis that adult neurologists should be familiar with. However, despite extensive clinical investigations, an acquired cause is identified in only a minority of cases. Thereafter, an underlying genetic basis is often considered, even in those without a family history. Here we apply whole exome sequencing to a cohort of 12 patients with late onset cerebellar ataxia. We show that 33% of 'idiopathic' cases harbor compound heterozygous mutations in known ataxia genes, including genes not included on multi-gene panels, or primarily associated with an ataxic presentation.

A phylogenetic approach to assessing the significance of missense mutations in disease genes.

The identification of deleterious mutations within candidate genes is a crucial step in the elucidation of the genetic bases of human disease. However, the significance of any base or amino acid change within a gene is unknown until detailed structural and functional analysis has been carried out. A potentially rapid way of identifying functionally important sites within a gene is to identify evolutionarily conserved regions. Mutations affecting such sites are assumed to be deleterious for the carrier. In this communication we generalize this approach and present a formal framework to assess whether a specific mutation is deleterious given sequence data from a set of homologues. We propose a score that takes into account the nature of the mutation, the conservation of the affected residue among the different species, and their phylogenetic relationships. Its performance is examined using published TP53 mutations and frequent polymorphic variants.

Structural organisation and chromosomal mapping of the human Id-3 gene.

The helix-loop-helix (HLH) family of transcription factors plays a central role in the regulation of cell growth, differentiation and tumourigenesis. Members of the Id (inhibitor of DNA binding) class of these nuclear proteins are able to heterodimerise with and thereby antagonise the functions of other transcription factors of this family. We report here on the genomic organisation of the human Id3 (HLH 1R21/heir1) gene. Comparison with the two other mammalian Id genes, Id1 and Id2, reveals a highly conserved protein coding gene organisation consistent with evolution from a common, ancestral Id-like gene. In addition, by using a yeast artificial chromosome (YAC) clone of Id3, we have fine-scale mapped the gene to chromosome band 1p36.1 by fluorescence in situ hybridisation (FISH) and, using the same FISH technique, we have detected heterogeneity in tumour-associated 1p36 chromosome translocations.

The expression of aphidicolin-induced fragile sites in familial breast cancer patients.

The expression frequency of aphidicolin-induced fragile sites was examined in familial breast cancer patients to determine whether this parameter could be used as a marker of genetic susceptibility in at-risk individuals. No difference was found in expression frequency between the breast cancer patients and a group of normal individuals (p = 0.61). This indicates that the expression frequency of aphidicolin-induced fragile sites is not a suitable marker for assessing genetic susceptibility in familial breast cancer.

Increased expression of fatty acid binding protein 4 and leptin in resident macrophages characterises atherosclerotic plaque rupture.

OBJECTIVE: Resident macrophages play an important role in atheromatous plaque rupture. The macrophage gene expression signature associated with plaque rupture is incompletely defined due to the complex cellular heterogeneity in the plaque. We aimed to characterise differential gene expression in resident plaque macrophages from ruptured and stable human atheromatous lesions. METHODS AND RESULTS: We performed genome-wide expression analyses of isolated macrophage-rich regions of stable and ruptured human atherosclerotic plaques. Plaques present in carotid endarterectomy specimens were designated as stable or ruptured using clinical, radiological and histopathological criteria. Macrophage-rich regions were excised from 5 ruptured and 6 stable plaques by laser micro-dissection. Transcriptional profiling was performed using Affymetrix microarrays. The profiles were characteristic of activated macrophages. At a false discovery rate of 10%, 914 genes were differentially expressed between stable and ruptured plaques. The findings were confirmed in fourteen further stable and ruptured samples for a subset of eleven genes with the highest expression differences (p < 0.05). Pathway analysis revealed that components of the PPAR/Adipocytokine signaling pathway were the most significantly upregulated in ruptured compared to stable plaques (p = 5.4 × 10(-7)). Two key components of the pathway, fatty-acid binding-protein 4 (FABP4) and leptin, showed nine-fold (p = 0.0086) and five-fold (p = 0.0012) greater expression respectively in macrophages from ruptured plaques. CONCLUSIONS: We found differences in gene expression signatures between macrophages isolated from stable and ruptured human atheromatous plaques. Our findings indicate the involvement of FABP4 and leptin in the progression of atherosclerosis and plaque rupture, and suggest that down-regulation of PPAR/adipocytokine signaling within plaques may have therapeutic potential.

Identification of candidate genes involved in neuroblastoma progression by combining genomic and expression microarrays with survival data.

Identifying genes, whose expression is consistently altered by chromosomal gains or losses, is an important step in defining genes of biological relevance in a wide variety of tumour types. However, additional criteria are needed to discriminate further among the large number of candidate genes identified. This is particularly true for neuroblastoma, where multiple genomic copy number changes of proven prognostic value exist. We have used Affymetrix microarrays and a combination of fluorescent in situ hybridization and single nucleotide polymorphism (SNP) microarrays to establish expression profiles and delineate copy number alterations in 30 primary neuroblastomas. Correlation of microarray data with patient survival and analysis of expression within rodent neuroblastoma cell lines were then used to define further genes likely to be involved in the disease process. Using this approach, we identify >1000 genes within eight recurrent genomic alterations (loss of 1p, 3p, 4p, 10q and 11q, 2p gain, 17q gain, and the MYCN amplicon) whose expression is consistently altered by copy number change. Of these, 84 correlate with patient survival, with the minimal regions of 17q gain and 4p loss being enriched significantly for such genes. These include genes involved in RNA and DNA metabolism, and apoptosis. Orthologues of all but one of these genes on 17q are overexpressed in rodent neuroblastoma cell lines. A significant excess of SNPs whose copy number correlates with survival is also observed on proximal 4p in stage 4 tumours, and we find that deletion of 4p is associated with improved outcome in an extended cohort of tumours. These results define the major impact of genomic copy number alterations upon transcription within neuroblastoma, and highlight genes on distal 17q and proximal 4p for downstream analyses. They also suggest that integration of discriminators, such as survival and comparative gene expression, with microarray data may be useful in the identification of critical genes within regions of loss or gain in many human cancers.

An expectation-maximization algorithm for the analysis of allelic expression imbalance.

A significant proportion of the variation between individuals in gene expression levels is genetic, and it is likely that these differences correlate with phenotypic differences or with risk of disease. Cis-acting polymorphisms are important in determining interindividual differences in gene expression that lead to allelic expression imbalance, which is the unequal expression of homologous alleles in individuals heterozygous for such a polymorphism. This expression imbalance can be detected using a transcribed polymorphism, and, once it is established, the next step is to identify the polymorphisms that are responsible for or predictive of allelic expression levels. We present an expectation-maximization algorithm for such analyses, providing a formal statistical framework to test whether a candidate polymorphism is associated with allelic expression differences.

Overexpression of aurora B kinase (AURKB) in primary non-small cell lung carcinoma is frequent, generally driven from one allele, and correlates with the level of genetic instability.

Aurora kinases are key regulators of chromosome segregation during mitosis. We have previously shown by microarray analysis of primary lung carcinomas and matched normal tissue that AURKB (22 out of 37) and AURKA (15 out of 37) transcripts are frequently over-represented in these tumours. We now confirm these observations in a second series of 44 carcinomas and also show that aurora B kinase protein levels are raised in the tumours compared to normal tissue. Elevated levels of expression in tumours are not a consequence of high-level amplification of the AURKB gene. Using a coding sequence polymorphism we show that in most cases (seven out of nine) tumour expression is predominantly driven from one AURKB allele. Given the function of aurora B kinase, we examined whether there was an association between expression levels and genetic instability. We defined two groups of high and low AURKB expression. Using a panel of 10 microsatellite markers, we found that the group showing the higher level of expression had a higher frequency of allelic imbalance (P=0.0012). Analysis of a number of other genes that are strongly and specifically expressed in tumour over normal lung, including SERPINB5, TERT and PRAME, showed marked allelic expression imbalances in the tumour tissue in the context of balanced or only marginally imbalanced relative allelic copy numbers. Our data support a model of early carcinogenesis wherein defects in the process of inactivation of lung stem-cell associated genes during differentiation, contributes to the development of carcinogenesis.

Dinucleotide frequencies in different reading frame positions of coding mammalian DNA sequences.

A statistical model for the assessment of suppressions or preferences of 16 dinucleotides in DNA sequences was developed. It is based on the description by a hypergeometric distribution of the doublet frequencies in randomly "scrambled" DNA sequences. The statistical test is sequential and extracts one after another dinucleotides that differ significantly from their expected values. It is shown that in mammalian DNA only TA and CG are consistently depressed in all three reading frame positions. The deviations of other dinucleotides are either restricted to one frame position or not significant. The possibility that the coding commitments of the DNA sequences may be the causes of the non-random distribution was studied. Only in position 1/2 of the reading frame is the frequency behavior of TA adequately explained by the amino acid sequence coded for. It is concluded that TA and CG are avoided wherever possible for reasons that do not reside in the coding function of mammalian DNA sequences.

Dinucleotide frequencies in different reading frame positions of coding bacterial DNA sequences.

Results of the analysis of the dinucleotide frequencies in different frame positions of coding bacterial sequences are presented. They are compared with those obtained in mammalian sequences. It is concluded that the dinucleotide frequencies in both types of sequences are caused by different influence factors.