Search

Nonallelic homologous recombination (NAHR), occurring between low-copy repeats (LCRs) >10 kb in size and sharing >97% DNA sequence identity, is responsible for the majority of recurrent genomic rearrangements in the human genome. Recent studies have shown that transposable elements (TEs) can also mediate recurrent deletions and translocations, indicating the features of substrates that mediate NAHR may be significantly less stringent than previously believed. Using >4 kb length and >95% sequence identity criteria, we analyzed of the genome-wide distribution of long interspersed element (LINE) retrotransposon and their potential to mediate NAHR. We identified 17 005 directly oriented LINE pairs located

Summary In areas without newborn screening for severe combined immunodeficiency (SCID), disease-defining infections may lead to diagnosis, and in some cases, may not be identified prior to the first year of life. We describe a female infant who presented with disseminated vaccine-acquired varicella (VZV) and vaccine-acquired rubella infections at 13 months of age. Immunological evaluations demonstrated neutropenia, isolated CD4 lymphocytopenia, the presence of CD8+T cells, poor lymphocyte proliferation, hypergammaglobulinaemia and poor specific antibody production to VZV infection and routine immunizations. A combination of whole exome sequencing and custom-designed chromosomal microarray with exon coverage of primary immunodeficiency genes detected compound heterozygous mutations (one single nucleotide variant and one intragenic copy number variant involving one exon) within the IL7R gene. Mosaicism for wild-type allele (20–30%) was detected in pretransplant blood and buccal DNA and maternal engraftment (5–10%) demonstrated in pretransplant blood DNA. This may be responsible for the patient's unusual immunological phenotype compared to classical interleukin (IL)-7Rα deficiency. Disseminated VZV was controlled with anti-viral and immune-based therapy, and umbilical cord blood stem cell transplantation was successful. Retrospectively performed T cell receptor excision circle (TREC) analyses completed on neonatal Guthrie cards identified absent TREC. This case emphasizes the danger of live viral vaccination in severe combined immunodeficiency (SCID) patients and the importance of newborn screening to identify patients prior to high-risk exposures. It also illustrates the value of aggressive pathogen identification and treatment, the influence newborn screening can have on morbidity and mortality and the significant impact of newer genomic diagnostic tools in identifying the underlying genetic aetiology for SCID patients.

Summary: The Dirichlet-multinomial (DMN) distribution is a fundamental model for multicategory count data with overdispersion. This distribution has many uses in bioinformatics including applications to metagenomics data, transctriptomics and alternative splicing. The DMN distribution reduces to the multinomial distribution when the overdispersion parameter ψ is 0. Unfortunately, numerical computation of the DMN log-likelihood function by conventional methods results in instability in the neighborhood of . An alternative formulation circumvents this instability, but it leads to long runtimes that make it impractical for large count data common in bioinformatics. We have developed a new method for computation of the DMN log-likelihood to solve the instability problem without incurring long runtimes. The new approach is composed of a novel formula and an algorithm to extend its applicability. Our numerical experiments show that this new method both improves the accuracy of log-likelihood evaluation and the runtime by several orders of magnitude, especially in high-count data situations that are common in deep sequencing data. Using real metagenomic data, our method achieves manyfold runtime improvement. Our method increases the feasibility of using the DMN distribution to model many high-throughput problems in bioinformatics. We have included in our work an R package giving access to this method and a vingette applying this approach to metagenomic data. Availability and implementation: An implementation of the algorithm together with a vignette describing its use is available in Supplementary Data. Contact: pengyu.bio@gmail.com or cashaw@bcm.edu Supplementary information: Supplementary Data are available at Bioinformatics online.

Autism is a neurodevelopmental disorder with increasing evidence of heterogeneous genetic etiology including de novo and inherited copy number variants (CNVs). We performed array comparative genomic hybridization using a custom Agilent 1 M oligonucleotide array intended to cover 197 332 unique exons in RefSeq genes; 98% were covered by at least one probe and 95% were covered by three or more probes with the focus on detecting relatively small CNVs that would implicate a single protein-coding gene. The study group included 99 trios from the Simons Simplex Collection. The analysis identified and validated 55 potentially pathogenic CNVs, categorized as de novo autosomal heterozygous, inherited homozygous autosomal, complex autosomal and hemizygous deletions on the X chromosome of probands. Twenty percent (11 of 55) of these CNV calls were rare when compared with the Database of Genomic Variants. Thirty-six percent (20 of 55) of the CNVs were also detected in the same samples in an independent analysis using the 1 M Illumina single-nucleotide polymorphism array. Findings of note included a common and sometimes homozygous 61 bp exonic deletion in SLC38A10, three CNVs found in lymphoblast-derived DNA but not present in whole-blood derived DNA and, most importantly, in a male proband, an exonic deletion of the TMLHE (trimethyllysine hydroxylase epsilon) that encodes the first enzyme in the biosynthesis of carnitine. Data for CNVs present in lymphoblasts but absent in fresh blood DNA suggest that these represent clonal outgrowth of individual B cells with pre-existing somatic mutations rather than artifacts arising in cell culture. GEO accession number GSE23765 (http://www.ncbi.nlm.nih.gov/geo/, date last accessed on 30 August 2011). Genboree accession: http://genboree.org/java-bin/gbrowser.jsp?refSeqId=1868&entryPointId=chr17&from=53496072&to=53694382&isPublic=yes, date last accessed on 30 August 2011.

Background. Annual vaccination is the primary means for preventing influenza. However, great interindividual variability exists in vaccine responses, the cellular events that take place in vivo after vaccination are poorly understood, and appropriate biomarkers for vaccine responsiveness have not been developed. Methods. We immunized a cohort of healthy male adults with a licensed trivalent influenza vaccine and performed a timed assessment of global gene expression before and after vaccination. We analyzed the relationship between gene expression patterns and the humoral immune response to vaccination. Results. Marked up regulation of expression of genes involved in interferon signaling, positive IL-6 regulation, and antigen processing and presentation, were detected within 24 hours of immunization. The late vaccine response showed a transcriptional pattern suggestive of increased protein biosynthesis and cellular proliferation. Integrative analyses revealed a 494-gene expression signature—including STAT1, CD74, and E2F2—which strongly correlates with the magnitude of the antibody response. High vaccine responder status correlates with increased early expression of interferon signaling and antigen processing and presentation genes. Conclusions. The results highlight the role of a systems biology approach in understanding the molecular events that take place in vivo after influenza vaccination and in the development of better predictors of vaccine responsiveness.

Genomic architecture, higher order structural features of the human genome, can provide molecular substrates for recurrent sub-microscopic chromosomal rearrangements, or may result in genomic instability by forming structures susceptible to DNA double-strand breaks. Pelizaeus-Merzbacher disease (PMD) is a genomic disorder most commonly arising from genomic duplications of the dosage-sensitive proteolipid protein gene (PLP1). Unlike many other genomic disorders that result from non-allelic homologous recombination utilizing flanking low-copy repeats (LCRs) as substrates, generating a common and recurrent rearrangement, the breakpoints of PLP1 duplications have been reported not to cluster, yielding duplicated genomic segments of varying lengths. This suggests a distinct molecular mechanism underlying PLP1 duplication events. To determine whether structural features of the genome also facilitate PLP1 duplication events, we analyzed extensively the genomic architecture of the PLP1 region and defined several novel LCRs (LCR-PMDs). Array comparative genomic hybridization showed that PLP1 duplication sizes differed, but revealed a subgroup of patients with apparently similar PLP1 duplication breakpoints. Pulsed-field gel electrophoresis analysis using probes adjacent to the LCR-PMDs detected unique recombination-specific junction fragments in 12 patients, enabled us to associate the LCR-PMDs with breakpoint regions, and revealed rearrangements inconsistent with simple tandem duplications in four patients. Two-color fluorescence in situ hybridization was consistent with directly oriented duplications. Our study provides evidence that PLP1 duplication events may be stimulated by LCRs, possibly non-homologous pairs at both the proximal and distal breakpoints in some cases, and further supports an alternative role of genomic architecture in rearrangements responsible for genomic disorders.

Cartilage-hair hypoplasia (CHH), also known as metaphyseal chondrodysplasia McKusick type (OMIM no. 250250), is an autosomal recessive, multi-systemic disease characterized by disproportionate short stature, fine and sparse hair, deficient cellular immunity and a predisposition to malignancy. It is caused by mutations in RMRP, the RNA component of the ribonucleoprotein complex RNase MRP, and, thus, CHH represents one of few Mendelian disorders caused by mutations in a nuclear encoded, non-coding RNA. While studies in yeast indicate that RMRP contributes to diverse cellular functions, the pathogenesis of the human condition is unknown. Studies of our CHH patient cohort revealed mutations in both the promoter and the transcribed region of RMRP. While mutations in the promoter abolished transcription in vitro, RMRP RNA levels in patients with transcribed mutations were also decreased suggesting an unstable RNA. RMRP mutations introduced into the yeast ortholog, NME1, exhibited normal mitochondrial function, chromosomal segregation and cell cycle progression, while a CHH fibroblast cell line exhibited normal mitochondrial content. However, the most commonly found mutation in CHH patients, 70A>G, caused an alteration in ribosomal processing by altering the ratio of the short versus the long form of the 5.8S rRNA in yeast. Transcriptional profiling of CHH patient RNAs showed upregulation of several cytokines and cell cycle regulatory genes, one of which has been implicated in chondrocyte hypertrophy. These data suggest that alteration of ribosomal processing in CHH is associated with altered cytokine signalling and cell cycle progression in terminally differentiating cells in the lymphocytic and chondrocytic cell lineages.

Terminal deletions of 1p36 result in a mental retardation syndrome that is presumably caused byhaploinsufficiency of a number of genes. Although monosomy 1p36 is the most commonly observedterminal deletion syndrome in humans, the molecular mechanism(s) that generates and stabilizes terminaldeletions of 1p36 is not completely understood. Our previous molecular analysis of a large cohort ofmonosomy 1p36 subjects demonstrated that deletion sizes vary widely from 1Mbto>10.5Mb in the mostdistal portion of 1p36 with no single common breakpoint. In this report, we have identified the precisebreakpoint junctions in three subjects with apparently pure terminal deletions of 1p36 ranging from 2.5 to4.25Mb.These junctionsrevealedone deletion tobestabilizedbytelomericrepeat sequences andtwoto haveterminal deletions associated with cryptic interrupted inverted duplications at the ends of the chromosomes.These interrupted inverted duplication/deletion breakpoints are reminiscent of those seen in tumor cell linesthat have undergone breakage–fusion–bridge (BFB) cycles leading to gene amplification. We propose a pre-meiotic model for the formation of these deletions in which a terminally deleted chromosome is generated inthe germ line and passes through at least one BFB cycle to produce gametes with terminal deletionsassociated with interrupted inverted duplications. These data suggest that, on a molecular level, seeminglypure terminal deletions visualized cytogenetically may be more complex, and BFB cycles may play animportant role in generating terminal deletions associated with genetic disease in humans.INTRODUCTIONTerminal deletions are one of the most commonly observedstructural chromosome abnormalities detected by routinechromosome analysis. Telomeric deficiencies result in severalwell-known mental retardation and multiple congenital anomalysyndromes including monosomy 1p36 (1p-), Wolf–Hirschhorn(4p-), Cri-du-chat (5p-), Miller–Dieker (17p-), 18q- and 22q-.In addition, cryptic telomeric abnormalities undetectable at thelevel of the light microscope have recently emerged as a major,yet previously under-recognized, cause of mental retardation(reviewedin1).Although terminaldeletionshavebeendescribedfor every human chromosome, the molecular mechanism(s) thatgenerates and stabilizes terminal deletions is not understood.Telomeres are specialized protein–DNA complexes that playcritical roles in the replication of linear chromosome ends andthe prevention of chromosome fusions (reviewed in 2,3).Chromosomes that lose a telomere due to improper telomeremaintenance or by a double-strand break near the telomere mustacquire a new telomeric cap to be structurally stable. Over 50years ago, McClintock first demonstrated in maize that chromo-somes that have lost their telomeres form end-to-end fusions thataresubsequentlybrokenbypassageofthenewlyformeddicentricchromosome through mitosis (4,5). This chromosome instabilityis perpetuated as a result of repeated breakage–fusion–bridge(BFB)cyclesduringeachcelldivisionuntilastabletelomericcapis obtained. Recently, telomere loss followed by BFB cycles hasbeen shown to play an important role in chromosome instabilityand gene amplification during cancer progression (6–10).However, to our knowledge, BFB cycles have not yet beendemonstrated to play a role in generating constitutionalchromosome rearrangements resulting in human genetic disease.Terminal deletions of 1p36 occur in 1 in 5000 live births,making it the most frequently observed terminal deletion andone of the most commonly occurring mental retardationsyndromes in humans (11). Monosomy 1p36 results in a

We developed an algorithm, HMZDelFinder, that uses whole exome sequencing (WES) data to identify rare and intragenic homozygous and hemizygous (HMZ) deletions that may represent complete loss-of-function of the indicated gene. HMZDelFinder was applied to 4866 samples in the Baylor–Hopkins Center for Mendelian Genomics (BHCMG) cohort and detected 773 HMZ deletion calls (567 homozygous or 206 hemizygous) with an estimated sensitivity of 86.5% (82% for single-exonic and 88% for multi-exonic calls) and precision of 78% (53% single-exonic and 96% for multi-exonic calls). Out of 773 HMZDelFinder-detected deletion calls, 82 were subjected to array comparative genomic hybridization (aCGH) and/or breakpoint PCR and 64 were confirmed. These include 18 single-exon deletions out of which 8 were exclusively detected by HMZDelFinder and not by any of seven other CNV detection tools examined. Further investigation of the 64 validated deletion calls revealed at least 15 pathogenic HMZ deletions. Of those, 7 accounted for 17–50% of pathogenic CNVs in different disease cohorts where 7.1–11% of the molecular diagnosis solved rate was attributed to CNVs. In summary, we present an algorithm to detect rare, intragenic, single-exon deletion CNVs using WES data; this tool can be useful for disease gene discovery efforts and clinical WES analyses.

Summary: Gene Ontology (GO) annotations have become a major tool for analysis of genome-scale experiments. We have created OntologyTraverser—an R package for GO analysis of gene lists. Our system is a major advance over previous work because (1) the system can be installed as an R package, (2) the system uses Java to instantiate the GO structure and the SJava system to integrate R and Java and (3) the system is also deployed as a publicly available web tool.Availability: Our software is academically available through http://franklin.imgen.bcm.tmc.edu/OntologyTraverser/. Both the R package and the web tool are accessible.Contact: cashaw@bcm.tmc.edu