Discrete domains of gene expression in germinal layers distinguish the development of gyrencephaly.

Gyrencephalic species develop folds in the cerebral cortex in a stereotypic manner, but the genetic mechanisms underlying this patterning process are unknown. We present a large-scale transcriptomic analysis of individual germinal layers in the developing cortex of the gyrencephalic ferret, comparing between regions prospective of fold and fissure. We find unique transcriptional signatures in each germinal compartment, where thousands of genes are differentially expressed between regions, including ~80% of genes mutated in human cortical malformations. These regional differences emerge from the existence of discrete domains of gene expression, which occur at multiple locations across the developing cortex of ferret and human, but not the lissencephalic mouse. Complex expression patterns emerge late during development and map the eventual location of folds or fissures. Protomaps of gene expression within germinal layers may contribute to define cortical folds or functional areas, but our findings demonstrate that they distinguish the development of gyrencephalic cortices.

The murine model for Hantaan virus-induced lethal disease shows two distinct paths in viral evolutionary trajectory with and without ribavirin treatment.

In vitro, ribavirin acts as a lethal mutagen in Hantaan virus (HTNV)-infected Vero E6 cells, resulting in an increased mutation load and viral population extinction. In this study, we asked whether ribavirin treatment in the lethal, suckling mouse model of HTNV infection would act similarly. The HTNV genomic RNA (vRNA) copy number and infectious virus were measured in lungs of untreated and ribavirin-treated mice. In untreated, HTNV-infected mice, the vRNA copy number increased for 10 days postinfection (dpi) and thereafter remained constant through 26 dpi. Surprisingly, in ribavirin-treated, HTNV-infected mice, vRNA levels were similar to those in untreated mice between 10 and 26 dpi. Infectious virus levels, however, were different: in ribavirin-treated mice, the amount of infectious HTNV was significantly decreased relative to that in untreated mice, suggesting that ribavirin reduced the specific infectivity of the virus (amount of infectious virus produced per vRNA copy). Mutational analysis revealed a ribavirin-associated elevation in mutation frequency in HTNV vRNA similar to that previously reported in vitro. Codon-based analyses of rates of nonsynonymous (dN) and synonymous (dS) substitutions in the S segment revealed a positive selection for codons within the HTNV N protein gene in the ribavirin-treated vRNA population. In contrast, the vRNA population in untreated, HTNV-infected mice showed a lower level of diversity, reflecting purifying selection for the wild-type genome. In summary, these experiments show two different evolutionary paths that Hantavirus may take during infection in a lethal murine model of disease, as well as the importance of the in vivo host environment in the evolution of the virus, which was not apparent in our prior in vitro model system.

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.

Identification of novel spinal cholinergic genetic subtypes disclose Chodl and Pitx2 as markers for fast motor neurons and partition cells.

Spinal cholinergic neurons are critical for motor function in both the autonomic and somatic nervous systems and are affected in spinal cord injury and in diseases such as amyotrophic lateral sclerosis (ALS) and spinal muscular atrophy. Using two screening approaches and in situ hybridization, we identified 159 genes expressed in typical cholinergic patterns in the spinal cord. These include two general cholinergic neuron markers, one gene exclusively expressed in motor neurons, and nine genes expressed in unknown subtypes of somatic motor neurons. Further, we present evidence that chondrolectin (Chodl) is expressed by fast motor neurons and that estrogen-related receptor beta (ERRbeta) is a candidate marker for slow motor neurons. In addition, we suggest paired-like homeodomain transcription factor 2 (Pitx2) as a marker for cholinergic partition cells.

Recurrent genomic alterations in benign and malignant pheochromocytomas and paragangliomas revealed by whole-genome array comparative genomic hybridization analysis.

Pheochromocytomas and abdominal paragangliomas are adrenal and extra-adrenal catecholamine-producing tumours. They arise due to heritable cancer syndromes, or more frequently occur sporadically due to an unknown genetic cause. The majority of cases are benign, but malignant tumours are observed. Previous comparative genomic hybridization (CGH) and loss of heterozygosity studies have shown frequent deletions of chromosome arms 1p, 3q and 22q in pheochromocytomas. We applied high-resolution whole-genome array CGH on 53 benign and malignant pheochromocytomas and paragangliomas to narrow down candidate regions as well as to identify chromosomal alterations more specific to malignant tumours. Minimal overlapping regions (MORs) were identified on 16 chromosomes, with the most frequent MORs of deletion (> or = 32%) occurring on chromosome arms 1p, 3q, 11p/q, 17p and 22q, while the chromosome arms 1q, 7p, 12q and 19p harboured the most common MORs of gain (> or = 14%). The most frequent MORs (61-75%) in the pheochromocytomas were identified at 1p, and the four regions of common losses encompassed 1p36, 1p32-31, 1p22-21 and 1p13. Tumours that did not show 1p loss generally demonstrated aberrations on chromosome 11. Gain of chromosomal material was significantly more frequent among the malignant cases. Moreover, gain at 19q, trisomy 12 and loss at 11q were positively associated with malignant pheochromocytomas, while 1q gain was commonly observed in the malignant paragangliomas. Our study revealed novel and narrow recurrent chromosomal regions of loss and gain at several autosomes, a prerequisite for identifying candidate tumour suppressor genes and oncogenes involved in the development of adrenal and extra-adrenal catecholamine-producing tumours.

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.

Frequent genetic differences between matched primary and metastatic breast cancer provide an approach to identification of biomarkers for disease progression.

Breast cancer is a major cause of morbidity and mortality in women and its metastatic spread is the principal reason behind the fatal outcome. Metastasis-related research of breast cancer is however underdeveloped when compared with the abundant literature on primary tumors. We applied an unexplored approach comparing at high resolution the genomic profiles of primary tumors and synchronous axillary lymph node metastases from 13 patients with breast cancer. Overall, primary tumors displayed 20% higher number of aberrations than metastases. In all but two patients, we detected in total 157 statistically significant differences between primary lesions and matched metastases. We further observed differences that can be linked to metastatic disease and there was also an overlapping pattern of changes between different patients. Many of the differences described here have been previously linked to poor patient survival, suggesting that this is a viable approach toward finding biomarkers for disease progression and definition of new targets useful for development of anticancer drugs. Frequent genetic differences between primary tumors and metastases in breast cancer also question, at least to some extent, the role of primary tumors as a surrogate subject of study for the systemic disease.

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.