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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

NF2 tumor suppressor gene: a comprehensive and efficient detection of somatic mutations by denaturing HPLC and microarray-CGH.

The NF2 tumor suppressor gene, located in chromosome 22q12, is involved in the development of multiple tumors of the nervous system, either associated with neurofibromatosis 2 or sporadic ones, mainly schwannomas and meningiomas. In order to evaluate the role of the NF2 gene in sporadic central nervous system (CNS) tumors, we analyzed NF2 mutations in 26 specimens: 14 meningiomas, 4 schwannomas, 4 metastases, and 4 other histopathological types of neoplasms. Denaturing high performance liquid chromatography (denaturing HPLC) and comparative genomic hybridization on a DNA microarray (microarray- CGH) were used as scanning methods for small mutations and gross rearrangements respectively. Small mutations were identified in six out of seventeen meningiomas and schwannomas, one mutation was novel. Large deletions were detected in six meningiomas. All mutations were predicted to result in truncated protein or in the absence of a large protein domain. No NF2 mutations were found in other histopathological types of CNS tumors. These results provide additional evidence that mutations in the NF2 gene play an important role in the development of sporadic meningiomas and schwannomas. Denaturing HPLC analysis of small mutations and microarray-CGH of large deletions are complementary, fast, and efficient methods for the detection of mutations in tumor tissues.

High-resolution profiling of an 11 Mb segment of human chromosome 22 in sporadic schwannoma using array-CGH.

Previous low-resolution schwannoma studies have reported diverse frequencies (30-80%) of 22q deletions, involving the neurofibromatosis-2 tumor suppressor (NF2) gene. We constructed an array spanning 11 million base pairs of 22q encompassing the NF2 gene, with 100% coverage and an average resolution of 58 kb. Moreover, the 220 kb genomic sequence encompassing the NF2 gene was covered by 13 cosmids to further enhance the resolution of analysis. The rationale of this array-CGH study was to map and size 22q deletions around the NF2 gene in sporadic schwannoma using a reliable method with maximal resolution. We studied tumor and constitutional DNA from 47 patients and detected heterozygous deletions in 21 (45%) tumors, which could be classified into three profiles. The predominant profile (12/21) was a continuous deletion of the 11 Mb segment, consistent with monosomy 22. The second profile, comprising five schwannomas, was also in agreement with a continuous 11 Mb heterozygous deletion. However, these displayed a distinctly different level of deletion when compared to the first profile, suggesting a considerable amount of normal tissue in the tumor samples. This is the first report demonstrating the sensitivity of array-CGH to discriminate such samples. The third profile was composed of four cases displaying interstitial deletions of various sizes. Two of these did not encompass the NF2 locus, which further emphasize the importance of other loci in schwannoma development. This is the first high-resolution study performed on a large series of tumors, using an array continuously covering 1/3 of a human chromosome. Our findings warrant further studies of an extended tumor series on a full 22q genomic array, to better define additional, putative 22q-located loci important for schwannoma development. Our array also provides a new diagnostic tool for analysis of NF2 gene deletions in patients affected with neurofibromatosis-2.

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.

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.

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.

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.

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.