Host gene expression signatures discriminate between ferrets infected with genetically similar H1N1 strains.

Different respiratory viruses induce virus-specific gene expression in the host. Recent evidence, including those presented here, suggests that genetically related isolates of influenza virus induce strain-specific host gene regulation in several animal models. Here, we identified systemic strain-specific gene expression signatures in ferrets infected with pandemic influenza A/California/07/2009, A/Mexico/4482/2009 or seasonal influenza A/Brisbane/59/2007. Using uncorrelated shrunken centroid classification, we were able to accurately identify the infecting influenza strain with a combined gene expression profile of 10 selected genes, independent of the severity of disease. Another gene signature, consisting of 7 genes, could classify samples based on lung pathology. Furthermore, we identified a gene expression profile consisting of 31 probes that could classify samples based on both strain and severity of disease. Thus, we show that expression-based analysis of non-infected tissue enables distinction between genetically related influenza viruses as well as lung pathology. These results open for development of alternative tools for influenza diagnostics.

De-novo transcriptome sequencing of a normalized cDNA pool from influenza infected ferrets.

The ferret is commonly used as a model for studies of infectious diseases. The genomic sequence of this animal model is not yet characterized, and only a limited number of fully annotated cDNAs are currently available in GenBank. The majority of genes involved in innate or adaptive immune response are still lacking, restricting molecular genetic analysis of host response in the ferret model. To enable de novo identification of transcriptionally active ferret genes in response to infection, we performed de-novo transcriptome sequencing of animals infected with H1N1 A/California/07/2009. We also included splenocytes induced with bacterial lipopolysaccharide to allow for identification of transcripts specifically induced by gram-negative bacteria. We pooled and normalized the cDNA library in order to delimit the risk of sequencing only highly expressed genes. While normalization of the cDNA library removes the possibility of assessing expression changes between individual animals, it has been shown to increase identification of low abundant transcripts. In this study, we identified more than 19,000 partial ferret transcripts, including more than 1000 gene orthologs known to be involved in the innate and the adaptive immune response.

Course of seasonal influenza A/Brisbane/59/07 H1N1 infection in the ferret.

Every year, influenza viruses infect approximately 5-20% of the population in the United States leading to over 200,000 hospitalizations and 36,000 deaths from flu-related complications. In this study, we characterized the immune and pathological progression of a seasonal strain of H1N1 influenza virus, A/Brisbane/59/2007 in a ferret model. The immune response of the animals showed a dose-dependent increase with increased virus challenge, as indicated by the presence of virus specific IgG, IgM, and neutralizing antibodies. Animals infected with higher doses of virus also experienced increasing severity of clinical symptoms and fever at 2 days post-infection (DPI). Interestingly, weight loss was more pronounced in animals infected with lower doses of virus compared to those infected with a higher dose; these results were consistent with viral titers of swabs collected from the nares, but not the throat. Analyzed specimens included nasal and throat swabs from 1, 3, 5, and 7 DPI as well as tissue samples from caudal lung and nasal turbinates. Viral titers of the swab samples in all groups were higher on 1 and 3 DPI and returned to baseline levels by 7 DPI. Analysis of nasal turbinates indicated presence of virus at 3 DPI in all infected groups, whereas virus was only detected in the lungs of animals in the two highest dose groups. Histological analysis of the lungs showed a range of pathology, such as chronic inflammation and bronchial epithelial hypertrophy. The results provided here offer important endpoints for preclinical testing of the efficacy of new antiviral compounds and experimental vaccines.

Transcriptome sequencing and development of an expression microarray platform for the domestic ferret.

BACKGROUND: The ferret (Mustela putorius furo) represents an attractive animal model for the study of respiratory diseases, including influenza. Despite its importance for biomedical research, the number of reagents for molecular and immunological analysis is restricted. We present here a parallel sequencing effort to produce an extensive EST (expressed sequence tags) dataset derived from a normalized ferret cDNA library made from mRNA from ferret blood, liver, lung, spleen and brain. RESULTS: We produced more than 500000 sequence reads that were assembled into 16000 partial ferret genes. These genes were combined with the available ferret sequences in the GenBank to develop a ferret specific microarray platform. Using this array, we detected tissue specific expression patterns which were confirmed by quantitative real time PCR assays. We also present a set of 41 ferret genes with even transcription profiles across the tested tissues, indicating their usefulness as housekeeping genes. CONCLUSION: The tools developed in this study allow for functional genomic analysis and make further development of reagents for the ferret model possible.

Characterization of novel and complex genomic aberrations in glioblastoma using a 32K BAC array.

Glioblastomas (GBs) are malignant CNS tumors often associated with devastating symptoms. Patients with GB have a very poor prognosis, and despite treatment, most of them die within 12 months from diagnosis. Several pathways, such as the RAS, tumor protein 53 (TP53), and phosphoinositide kinase 3 (PIK3) pathways, as well as the cell cycle control pathway, have been identified to be disrupted in this tumor. However, emerging data suggest that these aberrations represent only a fraction of the genetic changes involved in gliomagenesis. In this study, we have applied a 32K clone-based genomic array, covering 99% of the current assembly of the human genome, to the detailed genetic profiling of a set of 78 GBs. Complex patterns of aberrations, including high and narrow copy number amplicons, as well as a number of homozygously deleted loci, were identified. Amplicons that varied both in number (three on average) and in size (1.4 Mb on average) were frequently detected (81% of the samples). The loci encompassed not only previously reported oncogenes (EGFR, PDGFRA, MDM2, and CDK4) but also numerous novel oncogenes as GRB10, MKLN1, PPARGC1A, HGF, NAV3, CNTN1, SYT1, and ADAMTSL3. BNC2, PTPLAD2, and PTPRE, on the other hand, represent novel candidate tumor suppressor genes encompassed within homozygously deleted loci. Many of these genes are already linked to several forms of cancer; others represent new candidate genes that may serve as prognostic markers or even as therapeutic targets in the future. The large individual variation observed between the samples demonstrates the underlying complexity of the disease and strengthens the demand for an individualized therapy based on the genetic profile of the patient.

Distal 22q11.2 microduplication encompassing the BCR gene.

Chromosome 22 band q11.2 has been recognized to be highly susceptible to subtle microdeletions and microduplications, which have been attributed to the presence of several large segmental duplications; also known as low copy repeats (LCRs). These LCRs function as mediators of non-allelic homologous recombination (NAHR), which results in these chromosomal rearrangements as a result of unequal crossover. The four centromeric LCRs at proximal 22q11.2 have been previously implicated in recurrent chromosomal rearrangements including the DiGeorge/Velocardiofacial syndrome (DG/VCFs) microdeletion and its reciprocal microduplication. Recently, we and others have demonstrated that the four telomeric LCRs at distal 22q11.2 are causally implicated in a newly recognized recurrent distal 22q11.2 microdeletion syndrome in the region immediately telomeric to the DG/VCFs typically deleted region. Here we report on the clinical, cytogenetic, and array CGH studies of a 4.5-year-old girl with history of failure to thrive, developmental delay (DD), and relative macrocephaly. She carries a paternally inherited approximately 2.1 Mb microduplication at distal 22q11.2, which spans approximately 34 annotated genes, and is flanked by two of the four telomeric 22q11.2 LCRs. We conclude that the four telomeric LCRs at distal 22q11.2 can mediate both deletions and duplications in this genomic region. Both deletions and duplication of this region present with subtle clinical features including mild to moderate mental retardation, DD, and mild dysmorphic features.

Somatic mosaicism for copy number variation in differentiated human tissues.

Two major types of genetic variation are known: single nucleotide polymorphisms (SNPs), and a more recently discovered structural variation, involving changes in copy number (CNVs) of kilobase- to megabase-sized chromosomal segments. It is unknown whether CNVs arise in somatic cells, but it is, however, generally assumed that normal cells are genetically identical. We tested 34 tissue samples from three subjects and, having analyzed for each tissue < or =10(-6) of all cells expected in an adult human, we observed at least six CNVs, affecting a single organ or one or more tissues of the same subject. The CNVs ranged from 82 to 176 kb, often encompassing known genes, potentially affecting gene function. Our results indicate that humans are commonly affected by somatic mosaicism for stochastic CNVs, which occur in a substantial fraction of cells. The majority of described CNVs were previously shown to be polymorphic between unrelated subjects, suggesting that some CNVs previously reported as germline might represent somatic events, since in most studies of this kind, only one tissue is typically examined and analysis of parents for the studied subjects is not routinely performed. A considerable number of human phenotypes are a consequence of a somatic process. Thus, our conclusions will be important for the delineation of genetic factors behind these phenotypes. Consequently, biobanks should consider sampling multiple tissues to better address mosaicism in the studies of somatic disorders.

Autoantibodies to glutathione S-transferase theta 1 in patients with primary sclerosing cholangitis and other autoimmune diseases.

Primary sclerosing cholangitis (PSC) is an enigmatic disorder with a suggested autoimmune basis. A variety of autoantigens have been suggested but no specific or highly directed epitope has been identified. To address this issue, we constructed a cDNA library from normal human choledochus and screened expressing clones with serum from a patient with PSC and inflammatory bowel disease (IBD). Based on this screening, glutathione S-transferase theta 1 (GSTT1) was identified as a potential autoantigenic target. To study the specificity of GSTT1, we determined immunoreactivity using a panel of 58 patients with PSC, with and without IBD, 57 patients with IBD, 31 patients with Hashimoto's thyroiditis, 30 patients with primary biliary cirrhosis (PBC), 20 patients with insulin dependent diabetes mellitus, 22 patients with autoimmune polyendocrine syndrome type I, 10 patients with systemic lupus erythematosus (SLE), 20 patients with Sjögren's syndrome, 12 patients with autoimmune pancreatitis, 28 patients with Addison's disease, 27 patients with Grave's disease, 17 with myasthenia gravis, and 118 healthy controls. Reactivity against GSTT1 was found with PSC and IBD as well as some patients with other autoimmune pathology, indicating that this population of antibodies is neither specific nor a sensitive serologic marker for PSC, but the frequency was clearly higher in autoimmune patients than controls. GSTT1-antibodies have been described in persons with GSTT1-null genotype and are suggested to develop as an alloimmune response to blood transfusions from GSTT1-positive donors or pregnancies with GSTT1-positive children. Therefore, two IBD patients with and 15 PSC patients without GSTT1-antibodies were genotyped for GSTT1 to investigate if the presence of GSTT1-antibodies was associated with the GSTT1-null genotype and possibly caused by an alloimmune response. Both IBD patients and three of the PSC patients were of the GSTT1-null genotype. We note that the frequency of GSTT1-antibodies in this study is more than 100-fold higher than the frequency described earlier in patients with autoimmune diseases. We also observe an increased frequency of GSTT1-antibodies in patients with autoimmune diseases compared to healthy controls. This increased frequency can be explained by an autoimmune phenotype which increases susceptibility to such autoantibodies, or by a high frequency of the GSTT1-null genotype in autoimmune disease.

Phenotypically concordant and discordant monozygotic twins display different DNA copy-number-variation profiles.

The exploration of copy-number variation (CNV), notably of somatic cells, is an understudied aspect of genome biology. Any differences in the genetic makeup between twins derived from the same zygote represent an irrefutable example of somatic mosaicism. We studied 19 pairs of monozygotic twins with either concordant or discordant phenotype by using two platforms for genome-wide CNV analyses and showed that CNVs exist within pairs in both groups. These findings have an impact on our views of genotypic and phenotypic diversity in monozygotic twins and suggest that CNV analysis in phenotypically discordant monozygotic twins may provide a powerful tool for identifying disease-predisposition loci. Our results also imply that caution should be exercised when interpreting disease causality of de novo CNVs found in patients based on analysis of a single tissue in routine disease-related DNA diagnostics.

A segmental maximum a posteriori approach to genome-wide copy number profiling.

MOTIVATION: Copy number profiling methods aim at assigning DNA copy numbers to chromosomal regions using measurements from microarray-based comparative genomic hybridizations. Among the proposed methods to this end, Hidden Markov Model (HMM)-based approaches seem promising since DNA copy number transitions are naturally captured in the model. Current discrete-index HMM-based approaches do not, however, take into account heterogeneous information regarding the genomic overlap between clones. Moreover, the majority of existing methods are restricted to chromosome-wise analysis. RESULTS: We introduce a novel Segmental Maximum A Posteriori approach, SMAP, for DNA copy number profiling. Our method is based on discrete-index Hidden Markov Modeling and incorporates genomic distance and overlap between clones. We exploit a priori information through user-controllable parameterization that enables the identification of copy number deviations of various lengths and amplitudes. The model parameters may be inferred at a genome-wide scale to avoid overfitting of model parameters often resulting from chromosome-wise model inference. We report superior performances of SMAP on synthetic data when compared with two recent methods. When applied on our new experimental data, SMAP readily recognizes already known genetic aberrations including both large-scale regions with aberrant DNA copy number and changes affecting only single features on the array. We highlight the differences between the prediction of SMAP and the compared methods and show that SMAP accurately determines copy number changes and benefits from overlap consideration.

Statistical issues in the analysis of DNA Copy Number Variations.

Approaches to assess copy number variation have advanced rapidly and are being incorporated into genetic studies. While the technology exists for CNV genotyping, a further understanding and discussion of how to use the CNV data for association analyses is warranted. We present the options available for processing and analysing CNV data. We break these steps down into choice of genotyping platform, normalisation of the array data, calling algorithm, and statistical analysis.

Profiling of copy number variations (CNVs) in healthy individuals from three ethnic groups using a human genome 32 K BAC-clone-based array.

To further explore the extent of structural large-scale variation in the human genome, we assessed copy number variations (CNVs) in a series of 71 healthy subjects from three ethnic groups. CNVs were analyzed using comparative genomic hybridization (CGH) to a BAC array covering the human genome, using DNA extracted from peripheral blood, thus avoiding any culture-induced rearrangements. By applying a newly developed computational algorithm based on Hidden Markov modeling, we identified 1,078 autosomal CNVs, including at least two neighboring/overlapping BACs, which represent 315 distinct regions. The average size of the sequence polymorphisms was approximately 350 kb and involved in total approximately 117 Mb or approximately 3.5% of the genome. Gains were about four times more common than deletions, and segmental duplications (SDs) were overrepresented, especially in larger deletion variants. This strengthens the notion that SDs often define hotspots of chromosomal rearrangements. Over 60% of the identified autosomal rearrangements match previously reported CNVs, recognized with various platforms. However, results from chromosome X do not agree well with the previously annotated CNVs. Furthermore, data from single BACs deviating in copy number suggest that our above estimate of total variation is conservative. This report contributes to the establishment of the common baseline for CNV, which is an important resource in human genetics.

A patient with 22q11.2 deletion and Opitz syndrome-like phenotype has the same deletion as velocardiofacial patients.

Five patients were previously described with the Opitz (GBBB) syndrome (OMIM 145410) phenotype and 22q11.2 deletion determined by FISH but the precise limits of their deletions have not been determined. Since one locus for Opitz syndrome maps to 22q11.2 and chromosomal arrangements are frequently complex and could inactivate such a locus, we performed high-resolution array-based comparative genomic hybridization (CGH) on a new Opitz syndrome-like phenotype patient with a 22q11.2 deletion. He shares the same deletion as patients with velocardiofacial and DiGeorge syndrome.

A previously unrecognized microdeletion syndrome on chromosome 22 band q11.2 encompassing the BCR gene.

Susceptibility of the chromosome 22q11.2 region to rearrangements has been recognized on the basis of common clinical disorders such as the DiGeorge/velocardiofacial syndrome (DG/VCFs). Recent evidence has implicated low-copy repeats (LCRs); also known as segmental duplications; on 22q as mediators of nonallelic homologous recombination (NAHR) that result in rearrangements of 22q11.2. It has been shown that both deletion and duplication events can occur as a result of NAHR caused by unequal crossover of LCRs. Here we report on the clinical, cytogenetic and array CGH studies of a 15-year-old Hispanic boy with history of learning and behavior problems. We suggest that he represents a previously unrecognized microdeletion syndrome on chromosome 22 band q11.2 just telomeric to the DG/VCFs typically deleted region and encompassing the BCR gene. Using a 32K BAC array CGH chip we were able to refine and precisely narrow the breakpoints of this microdeletion, which was estimated to be 1.55-1.92 Mb in size and to span approximately 20 genes. This microdeletion region is flanked by LCR clusters containing several modules with a very high degree of sequence homology (>95%), and therefore could play a causal role in its origin.

Overlapping phenotype of Wolf-Hirschhorn and Beckwith-Wiedemann syndromes in a girl with der(4)t(4;11)(pter;pter).

We report on an 8-month-old girl with a novel unbalanced chromosomal rearrangement, consisting of a terminal deletion of 4p and a paternal duplication of terminal 11p. Each of these is associated with the well-known clinical phenotypes of Wolf-Hirschhorn syndrome (WHS) and Beckwith-Wiedemann syndrome (BWS), respectively. She presented for clinical evaluation of dysmorphic facial features, developmental delay, atrial septal defect (ASD), and left hydronephrosis. High-resolution cytogenetic analysis revealed a normal female karyotype, but subtelomeric fluorescence in situ hybridization (FISH) analysis revealed a der(4)t(4;11)(pter;pter). Both FISH and microarray CGH studies clearly demonstrated that the WHS critical regions 1 and 2 were deleted, and that the BWS imprinted domains (ID) 1 and 2 were duplicated on the der(4). Parental chromosome analysis revealed that the father carried a cryptic balanced t(4;11)(pter;pter). As expected, our patient manifests findings of both WHS (a growth retardation syndrome) and BWS (an overgrowth syndrome). We compare her unique phenotypic features with those that have been reported for both syndromes.

Growth hormone overexpression in the central nervous system results in hyperphagia-induced obesity associated with insulin resistance and dyslipidemia.

It is well known that peripherally administered growth hormone (GH) results in decreased body fat mass. However, GH-deficient patients increase their food intake when substituted with GH, suggesting that GH also has an appetite stimulating effect. Transgenic mice with an overexpression of bovine GH in the central nervous system (CNS) were created to investigate the role of GH in CNS. This study shows that overexpression of GH in the CNS differentiates the effect of GH on body fat mass from that on appetite. The transgenic mice were not GH-deficient but were obese and showed increased food intake as well as increased hypothalamic expression of agouti-related protein and neuropeptide Y. GH also had an acute effect on food intake following intracerebroventricular injection of C57BL/6 mice. The transgenic mice were severely hyperinsulinemic and showed a marked hyperplasia of the islets of Langerhans. In addition, the transgenic mice displayed alterations in serum lipid and lipoprotein levels and hepatic gene expression. In conclusion, GH overexpression in the CNS results in hyperphagia-induced obesity indicating a dual effect of GH with a central stimulation of appetite and a peripheral lipolytic effect.

Gene expression profiling shows that macrophages derived from mouse embryonic stem cells is an improved in vitro model for studies of vascular disease.

Macrophages (Mphi) play an important role in the initiation and progress of the atherogenic process. They contribute to the growth of atherosclerotic plaque by affecting lipoprotein metabolism, matrix homeostasis, lipoprotein modification and cholesterol accumulation. Access to in vitro Mphi models is therefore important for understanding the mechanisms involved in the transition of the relatively simple fatty streak to more complex type of lesions. The aim of the present work was to compare the expression profile of macrophages differentiated from the hematopoietic lineage to peritoneal mouse macrophages and two commonly used mouse macrophage cell lines (J774.A1 and RAW264.7). Our results showed that the embryonic stem cell-derived macrophages (ES Mphi) had a more similar expression profile to peritoneal macrophages than the two mouse macrophage cell lines. The ES Mphi had unchanged expression of the majority of cholesterol efflux mediators when compared to mouse peritoneal macrophages, whereas the cell lines showed altered expression of several of these genes. A key gene in this process is apolipoprotein E, which is expressed in ES Mphi but not in macrophage cell lines. In conclusion, ES Mphi are likely to provide a better in vitro model than mouse Mphi cell lines to study macrophage involvement in atherosclerosis.

Premature ageing in mice expressing defective mitochondrial DNA polymerase.

Point mutations and deletions of mitochondrial DNA (mtDNA) accumulate in a variety of tissues during ageing in humans, monkeys and rodents. These mutations are unevenly distributed and can accumulate clonally in certain cells, causing a mosaic pattern of respiratory chain deficiency in tissues such as heart, skeletal muscle and brain. In terms of the ageing process, their possible causative effects have been intensely debated because of their low abundance and purely correlative connection with ageing. We have now addressed this question experimentally by creating homozygous knock-in mice that express a proof-reading-deficient version of PolgA, the nucleus-encoded catalytic subunit of mtDNA polymerase. Here we show that the knock-in mice develop an mtDNA mutator phenotype with a threefold to fivefold increase in the levels of point mutations, as well as increased amounts of deleted mtDNA. This increase in somatic mtDNA mutations is associated with reduced lifespan and premature onset of ageing-related phenotypes such as weight loss, reduced subcutaneous fat, alopecia (hair loss), kyphosis (curvature of the spine), osteoporosis, anaemia, reduced fertility and heart enlargement. Our results thus provide a causative link between mtDNA mutations and ageing phenotypes in mammals.

Strong conservation of the human NF2 locus based on sequence comparison in five species.

We analyzed 137 kb covering human neurofibromatosis 2 ( NF2) tumor suppressor locus and orthologous loci from baboon, mouse, rat, and pufferfish Takifugu rubripes. A predominant feature of human-rodent conservation is a very similar distribution of conserved islands, regarding length, position, and degree of identity. By use of a threshold of 75% identity over > or =100 bp of gap-free alignment, comparisons of human-mouse sequences resulted in 3.58% for extra-exonic conservation, which can be compared to 4.5% of exonic sequence content within the human locus. We identified a duplication of neurofibromin 2 in pufferfish, which resulted in two putative proteins with 74% and 76% identity to the human protein. One distinct island (called inter 1), conserved between all analyzed species, was located between promoters of the NIPSNAP1 and NF2 genes. Inter 1 might represent a novel regulatory element, important for the function of this locus. The high level of intronic conservation in the NF2 locus suggests that a number of unknown regulatory elements might exist within this gene. These elements could be affected by disease-causing mutations in NF2 patients and NF2-associated tumors. Alternatively, this conservation might be explained by presence of not yet characterized transcriptional unit(s) within this locus.

Refined mapping of the Welander distal myopathy region on chromosome 2p13 positions the new candidate region telomeric of the DYSF locus.

Welander distal myopathy (WDM) is a late adult-onset autosomal dominant disorder, characterized by a slow progression and distal limb weakness of the extremity muscles. The WDM locus has been mapped to chromosome 2p13. Within this region a common shared haplotype co-segregates in all affected patients, indicating a founder effect. By undertaking an extended linkage analysis we have significantly reduced the WDM locus to a critical interval of approximately 1.2 Mb flanked by markers D2S358 and PAC3-H52. The dysferlin gene, a strong candidate gene responsible for two other distal myopathies in the same region, is located centromeric to PAC3-H52 and can thereby formally be excluded as cause for WDM.