Rapid screening of 4000 individuals for germ-line variations in the BRAF gene

Background: The BRAF gene is frequently somatically altered in malignant melanoma. A majority of variations are at the valine 600 residue leading to a V600E substitution that constitutively activates the kinase. We screened 4000 patient and control DNAs for germ-line variations at the valine 600 residue. Methods: We developed a novel assay by adapting single-base variation assays and software for MALDI-TOF (matrix-assisted laser desorption/ionization time-of-flight) mass spectrometry to screen for all 5 reported variants at codon 600 of the BRAF gene. We screened a case-control collection comprising samples from 1082 melanoma patients and 154 of their unaffected relatives from 1278 families and from 2744 individuals from 659 unselected twin families with no history of melanoma. A panel of 66 melanoma cell lines was used for variation-positive controls. Results: All melanoma cell lines that we had found previously to carry a codon 600 variation were verified in this study. Three of the 4 possible variants (V600E n = 47, V600K n = 2, V600R n = 1) were detected, but no case of V600D was available. No germ-line variants were found in the samples from the 3980 melanoma patients or from the control individuals. Conclusions: This new assay is a high-throughput, automated alternative to standard sequencing and can be used as a rapid initial screen for somatic variants associated with melanoma. Germ-line variants at valine 600 are unlikely to exist and do not contribute to the reported role of the BRAF gene in melanoma predisposition. (c) 2006 American Association for Clinical Chemistry.

Staphylococcus aureus genotyping using novel real-time PCR formats

One approach to microbial genotyping is to make use of sets of single-nucleotide polymorphisms (SNPs) in combination with binary markers. Here we report the modification and automation of a SNP-plus-binary-marker-based approach to the genotyping of Staphylococcus aureus and its application to 391 S. aureus isolates from southeast Queensland, Australia. The SNPs used were arcC210, tpi243, arcC162, gmk318, pta294, tpi36, tpi241, and pta383. These provide a Simpson's index of diversity (D) of 0.95 with respect to the S. aureus multilocus sequence typing database and define 61 genotypes and the major clonal complexes. The binary markers used were pvl, cna, sdrE, pT181, and pUB110. Two novel real-time PCR formats for interrogating these markers were compared. One of these makes use of light upon extension (LUX) primers and biplexed reactions, while the other is a streamlined modification of kinetic PCR using SYBR green. The latter format proved to be more robust. In addition, automated methods for DNA template preparation, reaction setup, and data analysis were developed. A single SNP-based method for ST-93 (Queensland clone) identification was also devised. The genotyping revealed the numerical importance of the South West Pacific and Queensland community-acquired methicillin-resistant S. aureus (MRSA) clones and the clonal complex 239 Aus-1/Aus-2 hospital-associated MRSA. There was a strong association between the community-acquired clones and pvl.

Non-major-histocompatibility-complex genetics of ankylosing spondylitis

There is strong evidence from twin and family studies indicating that a substantial proportion of the heritability of susceptibility to ankylosing spondylitis (AS) and its clinical manifestations is encoded by non-major-histocompatibility-complex genes. Efforts to identify these genes have included genomewide linkage studies and candidate gene association studies. One region, the interleukin (IL)-I gene complex on chromosome 2, has been repeatedly associated with AS in both Caucasians and Asians. It is likely that more than one gene in this complex is involved in AS, with the strongest evidence to date implicating IL-IA. Identifying the genes underlying other linkage regions has been difficult due to the lack of obvious candidates and the low power of most studies to date to identify genes of the small to moderate magnitude that are likely to be involved. The field is moving towards genomewide association analysis, involving much larger datasets of unrelated cases and controls. Early successes using this approach in other diseases indicates that it is likely to identify genes in common diseases like AS, but there remains the risk that the common-variant, common-disease hypothesis will not hold true in AS. Nonetheless, it is appropriate for the field to be cautiously optimistic that the next few years will bring great advances in our understanding of the genetics of this condition.

Identification of the semaphorin receptor PLXNA2 as a candidate for susceptibility to schizophrenia

The discovery of genetic factors that contribute to schizophrenia susceptibility is a key challenge in understanding the etiology of this disease. Here, we report the identification of a novel schizophrenia candidate gene on chromosome 1q32, plexin A2 (PLXNA2), in a genome-wide association study using 320 patients with schizophrenia of European descent and 325 matched controls. Over 25 000 single-nucleotide polymorphisms (SNPs) located within approximately 14 000 genes were tested. Out of 62 markers found to be associated with disease status, the most consistent finding was observed for a candidate locus on chromosome 1q32. The marker SNP rs752016 showed suggestive association with schizophrenia (odds ratio (OR) = 1.49, P = 0.006). This result was confirmed in an independent case control sample of European Americans (combined OR = 1.38, P = 0.035) and similar genetic effects were observed in smaller subsets of Latin Americans (OR = 1.26) and Asian Americans (OR = 1.37). Supporting evidence was also obtained from two family-based collections, one of which reached statistical significance (OR = 2.2, P = 0.02). High-density SNP mapping showed that the region of association spans approximately 60 kb of the PLXNA2 gene. Eight out of 14 SNPs genotyped showed statistically significant differences between cases and controls. These results are in accordance with previous genetic findings that identified chromosome 1q32 as a candidate region for schizophrenia. PLXNA2 is a member of the transmembrane semaphorin receptor family that is involved in axonal guidance during development and may modulate neuronal plasticity and regeneration. The PLXNA2 ligand semaphorin 3A has been shown to be upregulated in the cerebellum of individuals with schizophrenia. These observations, together with the genetic results, make PLXNA2 a likely candidate for the 1q32 schizophrenia susceptibility locus.

Further evidence of linkage at 7p22 with familial hyperaldosteronism type II

Primary aldosteronism (PAL) is caused by the autonomous over-production of aldosterone. Once thought rare, it is now reported to be responsible for 5–10% of hypertension. Familial hyperaldosteronism type II (FH-II), unlike familial hyperaldosteronism type I, is not glucocorticoid-remediable and not associated with the hybrid CYP11B1/CYP11B2 gene mutation. At least five times more common than FH-I, FH-II is clinically, biochemically and morphologically indistinguishable from apparently sporadic PAL, suggesting that its incidence maybe even higher. Studies performed in collaboration with C Stratakis (NIH, Bethesda) on our largest Australian FH-II family (eight affected members) demonstrated linkage at chromosome 7p22. Similar linkage at this region was also found in a South American FH-II family (DNA provided by MI New, Presbyterian Hospital, New York). Mutations in the exons and intron/exon boundaries of the PRKARIB gene (which resides at 7p22 and is closely related to PRKARIA gene mutated in Carney complex) have been excluded in our largest Australian FH-II family. Using more finely spaced markers, we have confirmed linkage at 7p22 in these 2 families, and identified a second Australian family with evidence of linkage at this locus. The combined multipoint LOD score for these 3 families is 4.87 (θ=0) with markers D7S462 and D7S2424, which exceeds the critical threshold for genome-wide significance suggested by Lander and Kruglyak (1995), providing strong support for this locus harbouring mutations responsible for FH-II. A newly identified recombination event in our largest Australian family has narrowed the region of linkage by 1.8 Mb, permitting exclusion of approximately half the genes residing in the original reported 5Mb linked locus. In addition, we have strongly excluded linkage to these key markers in two Australian families (maximum multipoint LOD scores −3.51 and −2.77), supporting the notion that FH-II may be genetically heterogeneous. In order to identify candidate genes at 7p22, more closely spaced markers will be used to refine the locus, as well as single nucleotide polymorphism analysis.

Increased prevalence of sex chromosome aneuploidies in specific language impairment and dyslexia

Aim: Sex chromosome aneuploidies increase the risk of spoken or written language disorders but individuals with specific language impairment (SLI) or dyslexia do not routinely undergo cytogenetic analysis. We assess the frequency of sex chromosome aneuploidies in individuals with language impairment or dyslexia. Method: Genome-wide single nucleotide polymorphism genotyping was performed in three sample sets: a clinical cohort of individuals with speech and language deficits (87 probands: 61 males, 26 females; age range 4 to 23 years), a replication cohort of individuals with SLI, from both clinical and epidemiological samples (209 probands: 139 males, 70 females; age range 4 to 17 years), and a set of individuals with dyslexia (314 probands: 224 males, 90 females; age range 7 to 18 years). Results: In the clinical language-impaired cohort, three abnormal karyotypic results were identified in probands (proband yield 3.4%). In the SLI replication cohort, six abnormalities were identified providing a consistent proband yield (2.9%). In the sample of individuals with dyslexia, two sex chromosome aneuploidies were found giving a lower proband yield of 0.6%. In total, two XYY, four XXY (Klinefelter syndrome), three XXX, one XO (Turner syndrome), and one unresolved karyotype were identified. Interpretation: The frequency of sex chromosome aneuploidies within each of the three cohorts was increased over the expected population frequency (approximately 0.25%) suggesting that genetic testing may prove worthwhile for individuals with language and literacy problems and normal non-verbal IQ. Early detection of these aneuploidies can provide information and direct the appropriate management for individuals. © 2013 The Authors. Developmental Medicine & Child Neurology published by John Wiley & Sons Ltd on behalf of Mac Keith Press.

Prenatal environmental exposure and neurodevelopmentally important gene expression in malformed brain tissue from pediatric intractable epilepsy patients

The primary objective of this proposal was to determine whether mitochondrial oxidative stress and variation in a particular mtDNA lineage contribute to the risk of developing cortical dysplasia and are potential contributing factors in epileptogenesis in children. The occurrence of epilepsy in children is highly associated with malformations of cortical development (MCD). It appears that MCD might arise from developmental errors due to environmental exposures in combination with inherited variation in response to environmental exposures and mitochondrial function. Therefore, it is postulated that variation in a particular mtDNA lineage of children contributes to the effects of mitochondrial DNA damage on MCD phenotype. Quantitative PCR and dot blot were used to examine mitochondrial oxidative damage and single nucleotide polymorphism (SNP) in the mitochondrial genome in brain tissue from 48 pediatric intractable epilepsy patients from Miami Children’s Hospital and 11 control samples from NICHD Brain and Tissue Bank for Developmental Disorders. Epilepsy patients showed higher mtDNA copy number compared to normal health subjects (controls). Oxidative mtDNA damage was lower in non-neoplastic but higher in neoplastic epilepsy patients compared to controls. There was a trend of lower mtDNA oxidative damage in the non-neoplastic (MCD) patients compared to controls, yet, the reverse was observed in neoplastic (MCD and Non-MCD) epilepsy patients. The presence of mtDNA SNP and haplogroups did not show any statistically significant relationships with epilepsy phenotypes. However, SNPs G9804A and G9952A were found in higher frequencies in epilepsy samples. Logistic regression analysis showed no relationship between mtDNA oxidative stress, mtDNA copy number, mitochondrial haplogroups and SNP variations with epilepsy in pediatric patients. The levels of mtDNA copy number and oxidative mtDNA damage and the SNPs G9952A and T10010C predicted neoplastic epilepsy, however, this was not significant due to a small sample size of pediatric subjects. Findings of this study indicate that an increase in mtDNA content may be compensatory mechanisms for defective mitochondria in intractable epilepsy and brain tumor. Further validation of these findings related to mitochondrial genotypes and mitochondrial dysfunction in pediatric epilepsy and MCD may lay the ground for the development of new therapies and prevention strategies during embryogenesis.

A Case-Only Genome-wide Association Study of Gender- and Age-specific Risk Markers for Childhood Leukemia

Males and age group 1 to 5 years show a much higher risk for childhood acute lymphoblastic leukemia (ALL). We performed a case-only genome-wide association study (GWAS), using the Illumina Infinium HumanCoreExome Chip, to unmask gender- and age-specific risk variants in 240 non-Hispanic white children with ALL recruited at Texas Children’s Cancer Center, Houston, Texas. Besides statistically most significant results, we also considered results that yielded the highest effect sizes. Existing experimental data and bioinformatic predictions were used to complement results, and to examine the biological significance of statistical results. Our study identified novel risk variants for childhood ALL. The SNP, rs4813720 (RASSF2), showed the statistically most significant gender-specific associations (P < 2 x 10-6). Likewise, rs10505918 (SOX5) yielded the lowest P value (P < 1 x 10-5) for age-specific associations, and also showed the statistically most significant association with age-at-onset (P < 1 x 10-4). Two SNPs, rs12722042 and 12722039, from the HLA-DQA1 region yielded the highest effect sizes (odds ratio (OR) = 15.7; P = 0.002) for gender-specific results, and the SNP, rs17109582 (OR = 12.5; P = 0.006), showed the highest effect size for age-specific results. Sex chromosome variants did not appear to be involved in gender-specific associations. The HLA-DQA1 SNPs belong to DQA1*01:07and confirmed previously reported male-specific association with DQA1*01:07. Twenty one of the SNPs identified as risk markers for gender- or age-specific associations were located in the transcription factor binding sites and 56 SNPs were non-synonymous variants, likely to alter protein function. Although bioinformatic analysis did not implicate a particular mechanism for gender- and age-specific associations, RASSF2 has an estrogen receptor-alpha binding site in its promoter. The unknown mechanisms may be due to lack of interest in gender- and age-specificity in associations. These results provide a foundation for further studies to examine the gender- and age-differential in childhood ALL risk. Following replication and mechanistic studies, risk factors for one gender or age group may have a potential to be used as biomarkers for targeted intervention for prevention and maybe also for treatment.

Resolving patterns of population genetic and phylogeographic structure to inform control and eradication initiatives for brown rats Rattus norvegicus on South Georgia

Funded by UK Government's Overseas Territories Environmental Programme (OTEP)

Cell Fate Specification and the Regulation of RNA-dependent DNA Methylation in the Arabidopsis Root Meristem

The Arabidopsis root apical meristem (RAM) is a complex tissue capable of generating all the cell types that ultimately make up the root. The work presented in this thesis takes advantage of the versatility of high-throughput sequencing to address two independent questions about the root meristem. Although a lot of information is known regarding the cell fate decisions that occur at the RAM, cortex specification and differentiation remain poorly understood. In the first part of this thesis, I used an ethylmethanesulfonate (EMS) mutagenized marker line to perform a forward genetics screen. The goal of this screen was to identify novel genes involved in the specification and differentiation of the cortex tissue. Mapping analysis from the results obtained in this screen revealed a new allele of BRASSINOSTEROID4 with abnormal marker expression in the cortex tissue. Although this allele proved to be non-cortex specific, this project highlights new technology that allows mapping of EMS-generated mutations without the need to map-cross or back-cross. In the second part of this thesis, using fluorescence activated cell sorting (FACS) coupled with high throughput sequencing, my collaborators and I generated single-base resolution whole genome DNA methylomes, mRNA transcriptomes, and smallRNA transcriptomes for six different populations of cell types in the Arabidopsis root meristem. We were able to discover that the columella is hypermethylated in the CHH context within transposable elements. This hypermethylation is accompanied by upregulation of the RNA-dependent DNA methylation pathway (RdDM), including higher levels of 24-nt silencing RNAs (siRNAs). In summary, our studies demonstrate the versatility of high-throughput sequencing as a method for identifying single mutations or to perform complex comparative genomic analyses.

Identification of independent association signals and putative functional variants for breast cancer risk through fine-scale mapping of the 12p11 locus

BACKGROUND: Multiple recent genome-wide association studies (GWAS) have identified a single nucleotide polymorphism (SNP), rs10771399, at 12p11 that is associated with breast cancer risk. METHOD: We performed a fine-scale mapping study of a 700 kb region including 441 genotyped and more than 1300 imputed genetic variants in 48,155 cases and 43,612 controls of European descent, 6269 cases and 6624 controls of East Asian descent and 1116 cases and 932 controls of African descent in the Breast Cancer Association Consortium (BCAC; http://bcac.ccge.medschl.cam.ac.uk/ ), and in 15,252 BRCA1 mutation carriers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Stepwise regression analyses were performed to identify independent association signals. Data from the Encyclopedia of DNA Elements project (ENCODE) and the Cancer Genome Atlas (TCGA) were used for functional annotation. RESULTS: Analysis of data from European descendants found evidence for four independent association signals at 12p11, represented by rs7297051 (odds ratio (OR) = 1.09, 95 % confidence interval (CI) = 1.06-1.12; P = 3 × 10(-9)), rs805510 (OR = 1.08, 95 % CI = 1.04-1.12, P = 2 × 10(-5)), and rs1871152 (OR = 1.04, 95 % CI = 1.02-1.06; P = 2 × 10(-4)) identified in the general populations, and rs113824616 (P = 7 × 10(-5)) identified in the meta-analysis of BCAC ER-negative cases and BRCA1 mutation carriers. SNPs rs7297051, rs805510 and rs113824616 were also associated with breast cancer risk at P < 0.05 in East Asians, but none of the associations were statistically significant in African descendants. Multiple candidate functional variants are located in putative enhancer sequences. Chromatin interaction data suggested that PTHLH was the likely target gene of these enhancers. Of the six variants with the strongest evidence of potential functionality, rs11049453 was statistically significantly associated with the expression of PTHLH and its nearby gene CCDC91 at P < 0.05. CONCLUSION: This study identified four independent association signals at 12p11 and revealed potentially functional variants, providing additional insights into the underlying biological mechanism(s) for the association observed between variants at 12p11 and breast cancer risk

Electrochemical RNA genosensors for toxic algal species: enhancing selectivity and sensitivity

Harmful algal blooms (HABs) are becoming more frequent as climate changes, with tropical species moving northward. Monitoring programs detecting the presence of toxic algae before they bloom are of paramount importance to protect aquatic ecosystems, aquaculture, human health and local economies. Rapid and reliable species identification methods using molecular barcodes coupled to biosensor detection tools have received increasing attention over the past decade as an alternative to the impractical standard microscopic counting-based techniques. This work reports on a PCR amplification-free electrochemical genosensor for the enhanced selective and sensitive detection of RNA from multiple Mediterranean toxic algal species. For a sandwich hybridization (SHA), we designed longer capture and signal probes for more specific target discrimination against a single base-pair mismatch from closely related species and for reproducible signals. We optimized experimental conditions, viz., minimal probe concentration in the SHA on a screen-printed gold electrode and selected the best electrochemical mediator. Probes from 13 Mediterranean dinoflagellate species were tested under optimized conditions and the format further tested for quantification of RNA from environmental samples. We not only enhanced the selectivity and sensitivity of the state-of-the-art toxic algal genosensors but also increased the repertoire of toxic algal biosensors in the Mediterranean, towards an integral and automatic monitoring system.

Electrochemical RNA genosensors for toxic algal species: enhancing selectivity and sensitivity

Harmful algal blooms (HABs) are becoming more frequent as climate changes, with tropical species moving northward. Monitoring programs detecting the presence of toxic algae before they bloom are of paramount importance to protect aquatic ecosystems, aquaculture, human health and local economies. Rapid and reliable species identification methods using molecular barcodes coupled to biosensor detection tools have received increasing attention over the past decade as an alternative to the impractical standard microscopic counting-based techniques. This work reports on a PCR amplification-free electrochemical genosensor for the enhanced selective and sensitive detection of RNA from multiple Mediterranean toxic algal species. For a sandwich hybridization (SHA), we designed longer capture and signal probes for more specific target discrimination against a single base-pair mismatch from closely related species and for reproducible signals. We optimized experimental conditions, viz., minimal probe concentration in the SHA on a screen-printed gold electrode and selected the best electrochemical mediator. Probes from 13 Mediterranean dinoflagellate species were tested under optimized conditions and the format further tested for quantification of RNA from environmental samples. We not only enhanced the selectivity and sensitivity of the state-of-the-art toxic algal genosensors but also increased the repertoire of toxic algal biosensors in the Mediterranean, towards an integral and automatic monitoring system.

A compendium of myeloma-associated chromosomal copy number abnormalities and their prognostic value.

To obtain a comprehensive genomic profile of presenting multiple myeloma cases we performed high-resolution single nucleotide polymorphism mapping array analysis in 114 samples alongside 258 samples analyzed by U133 Plus 2.0 expression array (Affymetrix). We examined DNA copy number alterations and loss of heterozygosity (LOH) to define the spectrum of minimally deleted regions in which relevant genes of interest can be found. The most frequent deletions are located at 1p (30%), 6q (33%), 8p (25%), 12p (15%), 13q (59%), 14q (39%), 16q (35%), 17p (7%), 20 (12%), and 22 (18%). In addition, copy number-neutral LOH, or uniparental disomy, was also prevalent on 1q (8%), 16q (9%), and X (20%), and was associated with regions of gain and loss. Based on fluorescence in situ hybridization and expression quartile analysis, genes of prognostic importance were found to be located at 1p (FAF1, CDKN2C), 1q (ANP32E), and 17p (TP53). In addition, we identified common homozygously deleted genes that have functions relevant to myeloma biology. Taken together, these analyses indicate that the crucial pathways in myeloma pathogenesis include the nuclear factor-κB pathway, apoptosis, cell-cycle regulation, Wnt signaling, and histone modifications. This study was registered at http://isrctn.org as ISRCTN68454111.

Homozygous deletion mapping in myeloma samples identifies genes and an expression signature relevant to pathogenesis and outcome.

PURPOSE: Myeloma is a clonal malignancy of plasma cells. Poor-prognosis risk is currently identified by clinical and cytogenetic features. However, these indicators do not capture all prognostic information. Gene expression analysis can be used to identify poor-prognosis patients and this can be improved by combination with information about DNA-level changes. EXPERIMENTAL DESIGN: Using single nucleotide polymorphism-based gene mapping in combination with global gene expression analysis, we have identified homozygous deletions in genes and networks that are relevant to myeloma pathogenesis and outcome. RESULTS: We identified 170 genes with homozygous deletions and corresponding loss of expression. Deletion within the "cell death" network was overrepresented and cases with these deletions had impaired overall survival. From further analysis of these events, we have generated an expression-based signature associated with shorter survival in 258 patients and confirmed this signature in data from two independent groups totaling 800 patients. We defined a gene expression signature of 97 cell death genes that reflects prognosis and confirmed this in two independent data sets. CONCLUSIONS: We developed a simple 6-gene expression signature from the 97-gene signature that can be used to identify poor-prognosis myeloma in the clinical environment. This signature could form the basis of future trials aimed at improving the outcome of poor-prognosis myeloma.