Genome-wide algorithm for detecting CNV associations with diseases

SNP genotyping arrays have been developed to characterize single-nucleotide polymorphisms (SNPs) and DNA copy number variations (CNVs). The quality of the inferences about copy number can be affected by many factors including batch effects, DNA sample preparation, signal processing, and analytical approach. Nonparametric and model-based statistical algorithms have been developed to detect CNVs from SNP genotyping data. However, these algorithms lack specificity to detect small CNVs due to the high false positive rate when calling CNVs based on the intensity values. Association tests based on detected CNVs therefore lack power even if the CNVs affecting disease risk are common. In this research, by combining an existing Hidden Markov Model (HMM) and the logistic regression model, a new genome-wide logistic regression algorithm was developed to detect CNV associations with diseases. We showed that the new algorithm is more sensitive and can be more powerful in detecting CNV associations with diseases than an existing popular algorithm, especially when the CNV association signal is weak and a limited number of SNPs are located in the CNV.^

Quantifying single gene copy number by measuring fluorescent probe lengths on combed genomic DNA

An approach was developed for the quantification of subtle gains and losses of genomic DNA. The approach relies on a process called molecular combing. Molecular combing consists of the extension and alignment of purified molecules of genomic DNA on a glass coverslip. It has the advantage that a large number of genomes can be combed per coverslip, which allows for a statistically adequate number of measurements to be made on the combed DNA. Consequently, a high-resolution approach to mapping and quantifying genomic alterations is possible. The approach consists of applying fluorescence hybridization to the combed DNA by using probes to identify the amplified region. Measurements then are made on the linear hybridization signals to ascertain the region's exact size. The reliability of the approach first was tested for low copy number amplifications by determining the copy number of chromosome 21 in a normal and trisomy 21 cell line. It then was tested for high copy number amplifications by quantifying the copy number of an oncogene amplified in the tumor cell line GTL-16. These results demonstrate that a wide range of amplifications can be accurately and reliably quantified. The sensitivity and resolution of the approach likewise was assessed by determining the copy number of a single allele (160 kb) alteration.

rrndb: the Ribosomal RNA Operon Copy Number Database

The Ribosomal RNA Operon Copy Number Database (rrndb) is an Internet-accessible database containing annotated information on rRNA operon copy number among prokaryotes. Gene redundancy is uncommon in prokaryotic genomes, yet the rRNA genes can vary from one to as many as 15 copies. Despite the widespread use of 16S rRNA gene sequences for identification of prokaryotes, information on the number and sequence of individual rRNA genes in a genome is not readily accessible. In an attempt to understand the evolutionary implications of rRNA operon redundancy, we have created a phylogenetically arranged report on rRNA gene copy number for a diverse collection of prokaryotic microorganisms. Each entry (organism) in the rrndb contains detailed information linked directly to external websites including the Ribosomal Database Project, GenBank, PubMed and several culture collections. Data contained in the rrndb will be valuable to researchers investigating microbial ecology and evolution using 16S rRNA gene sequences. The rrndb web site is directly accessible on the WWW at http://rrndb.cme.msu.edu.

High copy number I-Ab transgenes induce production of IgE through an interluekin 4-dependent mechanism.

To better understand the role of class II major histocompatibility complex molecules in both normal and autoimmune responses, we have produced a series of I-Ab transgenic mice. One of these transgenic constructs, designated NOD.PD, has the sequence of the NOD beta chain (Abeta(g7)) except at positions 56 and 57, where Pro-Asp replaces His-Ser. Several NOD.PD transgenic lines have been produced. One line of these mice carried a very high number of copies (>50) of the NOD.PD transgene. As has been described in other mice carrying high copy numbers of I-Ab transgenes, B-cell development was abnormal. The steady state numbers of mature B cells (IgM+/IgD(hi)) in the periphery were greatly reduced in transgenic mice compared to nontransgenic littermates. Surprisingly, rather than being accompanied by a generalized hypogammaglobulinemia, this B-cell deficiency was accompanied by elevated concentrations of IgG1 and IgE in the serum. Conversely, the levels of IgG2a were reduced in transgenic mice compared to nontransgenic littermates. Because this isotype pattern was characteristic of interleukin (IL)-4-induced class-switching, we then investigated the role of IL-4 in causing the observed phenotype. We crossed the high copy number transgenic mice with an IL-4-deficient strain of mice. As expected, the elevated levels of IgE in high copy number transgenic mice were eliminated when the IL-4 gene was inactivated. However, the reduction in the number of B cells was not ameliorated. These data indicate that the primary defect caused by the transgene was to reduce the number of B cells in these mice. This reduction was accompanied by a secondary increase in IL-4 production, which drove the remaining B cells toward the production of IgGl and IgE.

Calculus of variations and supplementary notes and exercises, 1945-1946.

Reproduced from type-written copy.

Interpreting Human Genetic Variation through Genomics and In Vivo Models

Improvements in genomic technology, both in the increased speed and reduced cost of sequencing, have expanded the appreciation of the abundance of human genetic variation. However the sheer amount of variation, as well as the varying type and genomic content of variation, poses a challenge in understanding the clinical consequence of a single mutation. This work uses several methodologies to interpret the observed variation in the human genome, and presents novel strategies for the prediction of allele pathogenicity.

Using the zebrafish model system as an in vivo assay of allele function, we identified a novel driver of Bardet-Biedl Syndrome (BBS) in CEP76. A combination of targeted sequencing of 785 cilia-associated genes in a cohort of BBS patients and subsequent in vivo functional assays recapitulating the human phenotype gave strong evidence for the role of CEP76 mutations in the pathology of an affected family. This portion of the work demonstrated the necessity of functional testing in validating disease-associated mutations, and added to the catalogue of known BBS disease genes.

Further study into the role of copy-number variations (CNVs) in a cohort of BBS patients showed the significant contribution of CNVs to disease pathology. Using high-density array comparative genomic hybridization (aCGH) we were able to identify pathogenic CNVs as small as several hundred bp. Dissection of constituent gene and in vivo experiments investigating epistatic interactions between affected genes allowed for an appreciation of several paradigms by which CNVs can contribute to disease. This study revealed that the contribution of CNVs to disease in BBS patients is much higher than previously expected, and demonstrated the necessity of consideration of CNV contribution in future (and retrospective) investigations of human genetic disease.

Finally, we used a combination of comparative genomics and in vivo complementation assays to identify second-site compensatory modification of pathogenic alleles. These pathogenic alleles, which are found compensated in other species (termed compensated pathogenic deviations [CPDs]), represent a significant fraction (from 3 – 10%) of human disease-associated alleles. In silico pathogenicity prediction algorithms, a valuable method of allele prioritization, often misrepresent these alleles as benign, leading to omission of possibly informative variants in studies of human genetic disease. We created a mathematical model that was able to predict CPDs and putative compensatory sites, and functionally showed in vivo that second-site mutation can mitigate the pathogenicity of disease alleles. Additionally, we made publically available an in silico module for the prediction of CPDs and modifier sites.

These studies have advanced the ability to interpret the pathogenicity of multiple types of human variation, as well as made available tools for others to do so as well.

Behavior Genetics and Post Genomics

The science of genetics is undergoing a paradigm shift. Recent discoveries, including the activity of retrotransposons, the extent of copy number variations, somatic and chromosomal mosaicism, and the nature of the epigenome as a regulator of DNA expressivity, are challenging a series of dogmas concerning the nature of the genome and the relationship between genotype and phenotype. DNA, once held to be the unchanging template of heredity, now appears subject to a good deal of environmental change; considered to be identical in all cells and tissues of the body, there is growing evidence that somatic mosaicism is the normal human condition; and treated as the sole biological agent of heritability, we now know that the epigenome, which regulates gene expressivity, can be inherited via the germline. These developments are particularly significant for behavior genetics for at least three reasons: First, these phenomena appear to be particularly prevalent in the human brain, and likely are involved in much of human behavior; second, they have important implications for the validity of heritability and gene association studies, the methodologies that largely define the discipline of behavior genetics; and third, they appear to play a critical role in development during the perinatal period, and in enabling phenotypic plasticity in offspring in particular. I examine one of the central claims to emerge from the use of heritability studies in the behavioral sciences, the principle of “minimal shared maternal effects,” in light of the growing awareness that the maternal perinatal environment is a critical venue for the exercise of adaptive phenotypic plasticity. This consideration has important implications for both developmental and evolutionary biology

Maternal Copy Number Variants (CNV) transmission to their Autism Spectrum Disorder (ASD) sons correlates with phenotypic traits

Background: Common neurodevelopmental disorder, global prevalence ~1 %; Persistent deficits in social communication and social interaction; restricted and repetitive behavior, interests, or activities; Highly heterogeneous clinical presentation; Male to female ratio ~4:1.

From genes to radioresistance in head and neck squamous cell carcinoma

Head and Neck Cancers (HNC) are a group of tumours located in the upper aero-digestive tract. Head and Neck Squamous Cell Carcinoma (HNSCC) represent about 90% of all HNC cases. It has been considered the sixth most malignant tumour worldwide and, despite clinical and technological advances, the five-year survival rate has not improved much in the last years. Nowadays, HNSCC is well established as a heterogeneous disease and that its development is due to accumulation of genetic events. Apart from the majority of the patients being diagnosed in an advanced stage, HNSCC is also a disease with poor therapeutic outcome. One of the therapeutic approaches is radiotherapy. However, this approach has different drawbacks like the radioresistance acquired by some tumour cells, leading to a worse prognosis. A major knowledge in radiation biology is imperative to improve this type of treatment and avoid late toxicities, maintaining patient quality of life in the subsequent years after treatment. Then, identification of genetic markers associated to radiotherapy response in patients and possible alterations in cells after radiotherapy are essential steps towards an improved diagnosis, higher survival rate and a better life quality. Not much is known about the radiation effects on cells, so, the principal aim of this study was to contribute to a more extensive knowledge about radiation treatment in HNSCC. For this, two commercial cell lines, HSC-3 and BICR-10, were used and characterized resorting to karyotyping, aCGH and MS-MLPA. These cell lines were submitted to different doses of irradiation and the resulting genetic and methylation alterations were evaluated. Our results showed a great difference in radiation response between the two cell lines, allowing the conclusion that HSC-3 was much more radiosensitive than BICR-10. Bearing this in mind, analysis of cell death, cell cycle and DNA damages was performed to try to elucidate the motifs behind this difference. The characterization of both cell lines allowed the confirmation that HSC-3 was derived from a metastatic tumour and the hypothesis that BICR-10 was derived from a dysplasia. Furthermore, this pilot study enabled the suggestion of some genetic and epigenetic alterations that cells suffer after radiation treatment. Additionally, it also allowed the association of some genetic characteristics that could be related to the differences in radiation response observable in this two cell lines. Taken together all of our results contribute to a better understanding of radiation effects on HNSCC allowing one further step towards the prediction of patients’ outcome, better choice of treatment approaches and ultimately a better quality of life.

Identification of outcome-oriented cut-offs for copy number alterations in multiple myeloma: predictive biomarkers with improved prognostic accuracy

Aim of the present study was to develop a statistical approach to define the best cut-off Copy number alterations (CNAs) calling from genomic data provided by high throughput experiments, able to predict a specific clinical end-point (early relapse, 18 months) in the context of Multiple Myeloma (MM). 743 newly diagnosed MM patients with SNPs array-derived genomic and clinical data were included in the study. CNAs were called both by a conventional (classic, CL) and an outcome-oriented (OO) method, and Progression Free Survival (PFS) hazard ratios of CNAs called by the two approaches were compared. The OO approach successfully identified patients at higher risk of relapse and the univariate survival analysis showed stronger prognostic effects for OO-defined high-risk alterations, as compared to that defined by CL approach, statistically significant for 12 CNAs. Overall, 155/743 patients relapsed within 18 months from the therapy start. A small number of OO-defined CNAs were significantly recurrent in early-relapsed patients (ER-CNAs) - amp1q, amp2p, del2p, del12p, del17p, del19p -. Two groups of patients were identified either carrying or not ≥1 ER-CNAs (249 vs. 494, respectively), the first one with significantly shorter PFS and overall survivals (OS) (PFS HR 2.15, p<0001; OS HR 2.37, p<0.0001). The risk of relapse defined by the presence of ≥1 ER-CNAs was independent from those conferred both by R-IIS 3 (HR=1.51; p=0.01) and by low quality (< stable disease) clinical response (HR=2.59 p=0.004). Notably, the type of induction therapy was not descriptive, suggesting that ER is strongly related to patients’ baseline genomic architecture. In conclusion, the OO- approach employed allowed to define CNAs-specific dynamic clonality cut-offs, improving the CNAs calls’ accuracy to identify MM patients with the highest probability to ER. As being outcome-dependent, the OO-approach is dynamic and might be adjusted according to the selected outcome variable of interest.

Distal 22q11.2 Microduplication Combined With Typical 22q11.2 Proximal Deletion: A Case Report.

The 22q11 chromosomal region contains low copy repeats (LCRs) sequences that mediate non-allelic homologous recombination, which predisposes to copy number variations (CNVs) at this locus. Hemizygous deletions of the proximal 22q11.2 region result in the 22q11.2 deletion syndrome (22q11.2 DS). In addition, 22q11.2 duplications involving the distal LCR22s have been reported. This article describes a patient presenting a 2.5-Mb de novo deletion at proximal 22q11.21 region (between LCRs A-D), combined with a 1.3-Mb maternally inherited duplication at distal 22q11.23 region (between LCRs F-H). The presence of concomitant chromosomal imbalances found in this patient has not been reported previously. Clinical and molecular data were compared with literature, in order to contribute to genotype-phenotype correlation. These findings exemplify the complexity and genetic heterogeneity observed in 22q11.2 deletion syndrome and highlights the difficulty to make genetic counseling and predict phenotypic consequences in these situations.

Polar organic marker compounds in atmospheric aerosols during the LBA-SMOCC 2002 biomass burning experiment in Rondonia, Brazil: sources and source processes, time series, diel variations and size distributions

Measurements of polar organic marker compounds were performed on aerosols that were collected at a pasture site in the Amazon basin (Rondonia, Brazil) using a high-volume dichotomous sampler (HVDS) and a Micro-Orifice Uniform Deposit Impactor (MOUDI) within the framework of the 2002 LBA-SMOCC (Large-Scale Biosphere Atmosphere Experiment in Amazonia - Smoke Aerosols, Clouds, Rainfall, and Climate: Aerosols From Biomass Burning Perturb Global and Regional Climate) campaign. The campaign spanned the late dry season (biomass burning), a transition period, and the onset of the wet season (clean conditions). In the present study a more detailed discussion is presented compared to previous reports on the behavior of selected polar marker compounds, including levoglucosan, malic acid, isoprene secondary organic aerosol (SOA) tracers and tracers for fungal spores. The tracer data are discussed taking into account new insights that recently became available into their stability and/or aerosol formation processes. During all three periods, levoglucosan was the most dominant identified organic species in the PM(2.5) size fraction of the HVDS samples. In the dry period levoglucosan reached concentrations of up to 7.5 mu g m(-3) and exhibited diel variations with a nighttime prevalence. It was closely associated with the PM mass in the size-segregated samples and was mainly present in the fine mode, except during the wet period where it peaked in the coarse mode. Isoprene SOA tracers showed an average concentration of 250 ng m(-3) during the dry period versus 157 ng m(-3) during the transition period and 52 ng m(-3) during the wet period. Malic acid and the 2-methyltetrols exhibited a different size distribution pattern, which is consistent with different aerosol formation processes (i.e., gas-to-particle partitioning in the case of malic acid and heterogeneous formation from gas-phase precursors in the case of the 2-methyltetrols). The 2-methyltetrols were mainly associated with the fine mode during all periods, while malic acid was prevalent in the fine mode only during the dry and transition periods, and dominant in the coarse mode during the wet period. The sum of the fungal spore tracers arabitol, mannitol, and erythritol in the PM(2.5) fraction of the HVDS samples during the dry, transition, and wet periods was, on average, 54 ng m(-3), 34 ng m(-3), and 27 ng m(-3), respectively, and revealed minor day/night variation. The mass size distributions of arabitol and mannitol during all periods showed similar patterns and an association with the coarse mode, consistent with their primary origin. The results show that even under the heavy smoke conditions of the dry period a natural background with contributions from bioaerosols and isoprene SOA can be revealed. The enhancement in isoprene SOA in the dry season is mainly attributed to an increased acidity of the aerosols, increased NO(x) concentrations and a decreased wet deposition.

European ancestry and polymorphisms in DNA repair genes modify the risk of melanoma: A case-control study in a high UV index region in Brazil

Background: UV radiation is the major environmental factor related to development of cutaneous melanoma. Besides sun exposure and the influence of latitude, some host characteristics such as skin phototype and hair and eye color are also risk factors for melanoma. Polymorphisms in DNA repair genes could be good candidates for susceptibility genes, mainly in geographical regions exposed to high solar radiation. Objective: Evaluate the role of host characteristic.; and DNA repair polymorphism in melanoma risk in Brazil. Methods: We carried out a hospital-based case-control study in Brazil to evaluate the contribution of host factors and polymorphisms in DNA repair to melanoma risk. A total of 412 patients (202 with melanoma and 210 controls) were analyzed regarding host characteristics for melanoma risk as well as for 11 polymorphisms in DNA repair genes. Results: We found an association of host characteristics with melanoma development, such as eye and hair color, fair skin, history of pigmented lesions removed, sunburns in childhood and adolescence, and also European ancestry. Regarding DNA repair gene polymorphisms, we found protection for the XPG 1104 His/His genotype (OR 0.32; 95% CI 0.13-0.75), and increased risk for three polymorphisms in the XPC gene (PAT+; IV-6A and 939Gln), which represent a haplotype for XPC. Melanoma risk was higher in individuals carrying the complete XPC haplotype than each individual polymorphism (OR 3.64; 95% CI 1.77-7.48). Conclusions: Our data indicate that the host factors European ancestry and XPC polymorphisms contributed to melanoma risk in a region exposed to high sun radiation. (C) 2011 Japanese Society for Investigative Dermatology. Published by Elsevier Ireland Ltd. All rights reserved.

Efficiency of the DNA repair and polymorphisms of the XRCC1, XRCC3 and XRCC4 DNA repair genes in systemic lupus erythematosus

Impaired DNA repair efficiency in systematic lupus erythematosus (SLE) patients has been reported ill some studies, mainly regarding the repair of oxidative damage, but little is known about repair kinetics towards primarily single-stranded DNA breaks. In the present study, we aimed to investigate: (a) the efficiency of SLE peripheral blood leucocytes in repairing DNA damage induced by ionizing radiation and (b) the association of DNA repair gene (XRCC1 Arg399Gln, XRCC3 Thr241Met and XRCC4 Ile401Thr) polymorphisms in SLE patients, considering the whole group, or stratified sub-groups according to clinical and laboratory features. A total of 163 SLE patients and 125 healthy control were studied. The kinetics of DNA strand break repair was evaluated by the comet assay, and genotyping for DNA repair genes was performed by PCR-RFLP. Compared with controls. SLE leucocytes exhibited decreased efficiency of DNA repair evaluated at 30 min following irradiation. A significant association with DNA repair gene polymorphisms was not observed for the whole group of SLE patients; however, the XRCC1Arg399Gln polymorphism was associated with the presence of anti-dsDNA antibody. The concomitance of two DNA repair polymorphic sites was associated with the presence of neuropsychiatric manifestations and antiphospholipid antibody syndrome. Taken together, these results indicated that SLE leucocytes repair less efficiently the radiation-induced DNA damage, and DNA repair polymorphic sites may predispose to the development of particular clinical and laboratory features. Lupus (2008) 17, 988-995.