Daniel C. Ciobanu

Strain Differences in Stress Responsivity Are Associated with Divergent Amygdala Gene Expression and Glutamate-Mediated Neuronal Excitability

Stress is a major risk factor for numerous neuropsychiatric diseases. However, susceptibility to stress and the qualitative nature of stress effects on behavior differ markedly among individuals. This is partly because of the moderating influence of genetic factors. Inbred mouse strains provide a relatively stable and restricted range of genetic and environmental variability that is valuable for disentangling gene–stress interactions. Here, we screened a panel of inbred strains for anxiety- and depression-related phenotypes at baseline (trait) and after exposure to repeated restraint. Two strains, DBA/2J and C57BL/6J, differed in trait and restraint-induced anxiety-related behavior (dark/light exploration, elevated plus maze). Gene expression analysis of amygdala, medial prefrontal cortex, and hippocampus revealed divergent expression in DBA/2J and C57BL/6J both at baseline and after repeated restraint. Restraint produced strain-dependent expression alterations in various genes including glutamate receptors (e.g., Grin1, Grik1). To elucidate neuronal correlates of these strain differences, we performed ex vivo analysis of glutamate excitatory neurotransmission in amygdala principal neurons. Repeated restraint augmented amygdala excitatory postsynaptic signaling and altered metaplasticity (temporal summation of NMDA receptor currents) in DBA/2J but not C57BL/6J. Furthermore, we found that the C57BL/6J-like changes in anxiety-related behavior after restraint were absent in null mutants lacking the modulatory NMDA receptor subunit Grin2a, but not the AMPA receptor subunit Gria1. Grin2a null mutants exhibited significant (∼30%) loss of dendritic spines on amygdala principal neurons under nonrestraint conditions. Collectively, our data support a model in which genetic variation in glutamatergic neuroplasticity in corticolimbic circuitry underlies phenotypic variation in responsivity to stress.

Murine Gut Microbiota Is Defined by Host Genetics and Modulates Variation of Metabolic Traits

The gastrointestinal tract harbors a complex and diverse microbiota that has an important role in host metabolism. Microbial diversity is influenced by a combination of environmental and host genetic factors and is associated with several polygenic diseases. In this study we combined next-generation sequencing, genetic mapping, and a set of physiological traits of the BXD mouse population to explore genetic factors that explain differences in gut microbiota and its impact on metabolic traits. Molecular profiling of the gut microbiota revealed important quantitative differences in microbial composition among BXD strains. These differences in gut microbial composition are influenced by host-genetics, which is complex and involves many loci. Linkage analysis defined Quantitative Trait Loci (QTLs) restricted to a particular taxon, branch or that influenced the variation of taxa across phyla. Gene expression within the gastrointestinal tract and sequence analysis of the parental genomes in the QTL regions uncovered candidate genes with potential to alter gut immunological profiles and impact the balance between gut microbial communities. A QTL region on Chr 4 that overlaps several interferon genes modulates the population of Bacteroides, and potentially Bacteroidetes and Firmicutes–the predominant BXD gut phyla. Irak4, a signaling molecule in the Toll-like receptor pathways is a candidate for the QTL on Chr15 that modulates Rikenellaceae, whereas Tgfb3, a cytokine modulating the barrier function of the intestine and tolerance to commensal bacteria, overlaps a QTL on Chr 12 that influence Prevotellaceae. Relationships between gut microflora, morphological and metabolic traits were uncovered, some potentially a result of common genetic sources of variation.

Dissection of a QTL Hotspot on Mouse Distal Chromosome 1 that Modulates Neurobehavioral Phenotypes and Gene Expression

A remarkably diverse set of traits maps to a region on mouse distal chromosome 1 (Chr 1) that corresponds to human Chr 1q21–q23. This region is highly enriched in quantitative trait loci (QTLs) that control neural and behavioral phenotypes, including motor behavior, escape latency, emotionality, seizure susceptibility (Szs1), and responses to ethanol, caffeine, pentobarbital, and haloperidol. This region also controls the expression of a remarkably large number of genes, including genes that are associated with some of the classical traits that map to distal Chr 1 (e.g., seizure susceptibility). Here, we ask whether this QTL-rich region on Chr 1 (Qrr1) consists of a single master locus or a mixture of linked, but functionally unrelated, QTLs. To answer this question and to evaluate candidate genes, we generated and analyzed several gene expression, haplotype, and sequence datasets. We exploited six complementary mouse crosses, and combed through 18 expression datasets to determine class membership of genes modulated by Qrr1. Qrr1 can be broadly divided into a proximal part (Qrr1p) and a distal part (Qrr1d), each associated with the expression of distinct subsets of genes. Qrr1d controls RNA metabolism and protein synthesis, including the expression of ∼20 aminoacyl-tRNA synthetases. Qrr1d contains a tRNA cluster, and this is a functionally pertinent candidate for the tRNA synthetases. Rgs7 and Fmn2 are other strong candidates in Qrr1d. FMN2 protein has pronounced expression in neurons, including in the dendrites, and deletion of Fmn2 had a strong effect on the expression of few genes modulated by Qrr1d. Our analysis revealed a highly complex gene expression regulatory interval in Qrr1, composed of multiple loci modulating the expression of functionally cognate sets of genes.

A gene-based SNP linkage map for Pacific white shrimp, Litopenaeus vannamei.

Pacific white shrimp (Litopenaeus vannamei) are of particular economic importance to the global shrimp aquaculture industry. However, limited genomics information is available for the penaeid species. We utilized the limited public information available, mainly single nucleotide polymorphisms (SNPs) and expressed sequence tags, to discover markers for the construction of the first SNP genetic map for Pacific white shrimp. In total, 1344 putative SNPs were discovered, and out of 825 SNPs genotyped, 418 SNP markers from 347 contigs were mapped onto 45 sex-averaged linkage groups, with approximate coverage of 2071 and 2130 cm for the female and male maps, respectively. The average-squared correlation coefficient (r(2)), a measure of linkage disequilibrium, for markers located more than 50 cm apart on the same linkage group, was 0.15. Levels of r(2) increased with decreasing inter-marker distance from approximately 80 cm, and increased more rapidly from approximately 30 cm. A QTL for shrimp gender was mapped on linkage group 13. Comparative mapping to model organisms, Daphnia pulex and Drosophila melanogaster, revealed extensive rearrangement of genome architecture for L. vannamei, and that L. vannamei was more related to Daphnia pulex. This SNP genetic map lays the foundation for future shrimp genomics studies, especially the identification of genetic markers or regions for economically important traits.

Detection, Validation, and Downstream Analysis of Allelic Variation in Gene Expression

Common sequence variants within a gene often generate important differences in expression of corresponding mRNAs. This high level of local (allelic) control—or cis modulation—rivals that produced by gene targeting, but expression is titrated finely over a range of levels. We are interested in exploiting this allelic variation to study gene function and downstream consequences of differences in expression dosage. We have used several bioinformatics and molecular approaches to estimate error rates in the discovery of cis modulation and to analyze some of the biological and technical confounds that contribute to the variation in gene expression profiling. Our analysis of SNPs and alternative transcripts, combined with eQTL maps and selective gene resequencing, revealed that between 17 and 25% of apparent cis modulation is caused by SNPs that overlap probes rather than by genuine quantitative differences in mRNA levels. This estimate climbs to 40–50% when qualitative differences between isoform variants are included. We have developed an analytical approach to filter differences in expression and improve the yield of genuine cis-modulated transcripts to ∼80%. This improvement is important because the resulting variation can be successfully used to study downstream consequences of altered expression on higher-order phenotypes. Using a systems genetics approach we show that two validated cis-modulated genes, Stk25 and Rasd2, are likely to control expression of downstream targets and affect disease susceptibility.

Joint mouse-human phenome-wide association to test gene function and disease risk

Phenome-wide association is a novel reverse genetic strategy to analyze genome-to-phenome relations in human clinical cohorts. Here we test this approach using a large murine population segregating for ∼5 million sequence variants, and we compare our results to those extracted from a matched analysis of gene variants in a large human cohort. For the mouse cohort, we amassed a deep and broad open-access phenome consisting of ∼4,500 metabolic, physiological, pharmacological and behavioural traits, and more than 90 independent expression quantitative trait locus (QTL), transcriptome, proteome, metagenome and metabolome data sets—by far the largest coherent phenome for any experimental cohort (www.genenetwork.org). We tested downstream effects of subsets of variants and discovered several novel associations, including a missense mutation in fumarate hydratase that controls variation in the mitochondrial unfolded protein response in both mouse and Caenorhabditis elegans, and missense mutations in Col6a5 that underlies variation in bone mineral density in both mouse and human.

A major SNP resource for dissection of phenotypic and genetic variation in Pacific white shrimp (Litopenaeus vannamei)

Bioinformatics and re-sequencing approaches were used for the discovery of sequence polymorphisms in Litopenaeus vannamei. A total of 1221 putative single nucleotide polymorphisms (SNPs) were identified in a pool of individuals from various commercial populations. A set of 211 SNPs were selected for further molecular validation and 88% showed variation in 637 samples representing three commercial breeding lines. An association analysis was performed between these markers and several traits of economic importance for shrimp producers including resistance to three major viral diseases. A small number of SNPs showed associations with test weekly gain, grow-out survival and resistance to Taura Syndrome Virus. Very low levels of linkage disequilibrium were revealed between most SNP pairs, with only 11% of SNPs showing an r(2)-value above 0.10 with at least one other SNP. Comparison of allele frequencies showed small changes over three generations of the breeding programme in one of the commercial breeding populations. This unique SNP resource has the potential to catalyse future studies of genetic dissection of complex traits, tracing relationships in breeding programmes, and monitoring genetic diversity in commercial and wild populations of L. vannamei.

Genetic variation in two conserved local Romanian pig breeds using type 1 DNA markers

Analysis of the genetic variation of an endangered population is an important component for the success of conservation. Animals from two local Romanian pig breeds, the Mangalitsa and Bazna, were analysed for variation at a number of genetic loci using PCR-based DNA tests. Polymorphism was assessed at loci which 1) are known to cause phenotypic variation, 2) are potentially involved in trait differences or 3) are putative candidate genes. The traits considered are disease resistance, growth, coat colour, meat quality and prolificacy. Even though the populations are small and the markers are limited to specific genes, we found significant differences in five of the ten characterised loci. In some cases the observed allele frequencies were interesting in relation to gene function and the phenotype of the breed. These breeds are part of a conservation programme in Romania and marker information may be useful in preserving a representative gene pool in the populations. The use of polymorphisms in type 1 (gene) markers may be a useful complement to analysis based on anonymous markers.

Multiple Genes on Chromosome 7 Regulate Dopaminergic Amacrine Cell Number in the Mouse Retina

PURPOSE The size of neuronal populations is modulated by gene variants that influence cell production and survival, in turn influencing neuronal connectivity, function, and disease risk. The size of the dopaminergic amacrine (DA) cell population is a highly heritable trait exhibiting sixfold variation among inbred strains of mice and is used here to identify genes that modulate the number of DA cells. METHODS The entire population was counted in retinal wholemounts from 37 genetically defined lines of mice, including six standard inbred strains, 25 recombinant inbred strains (AXB/BXA), reciprocal F1 hybrids, a chromosome (Chr) 7 consomic line, and three additional genetically modified lines. RESULTS Much of this variation was mapped to a broad locus on Chr 7 (Dopaminergic amacrine cell number control, Chr 7 [Dacnc7]). The Dacnc7 locus is flanked by two candidate genes known to modulate the number of other types of retinal neuron-the proapoptotic gene, Bax, and tyrosinase. The Tyr mutation was shown to modulate DA cell number modestly, though in the direction opposite that predicted. In contrast, Bax deficiency increased the population fourfold. Bax expression was significantly greater in the A/J than in the C57BL/6J strain, an effect that may be attributed to an SNP in a p53 consensus binding site known to modulate transcription. Finally, we note a strong candidate situated at the peak of the Dacnc7 locus, Lrrk1, a Parkinsons disease gene exhibiting missense mutations segregating within the AXB/BXA cross. CONCLUSIONS Multiple polymorphic genes on Chr 7 modulate the size of the population of DA cells.

High-throughput sequencing of the DBA/2J mouse genome

BackgroundThe DBA/2J mouse is not only the oldest inbred strain,but also one of the most widely used strains. DBA/2Jexhibits many unique anatomical, physiological, andbehavior traits. In addition, DBA/2J is one parent of thelarge BXD family of recombinant inbred strains [1]. Thegenome of the other parent of this BXD family—C57BL/6J—has been sequenced and serves as themouse reference genome [2]. We sequenced the gen-ome of DBA/2J using SOLiD and Illumina highthroughput short read protocols to generate a compre-hensive set of ~5 million sequence variants segregatingin the BXD family that ultimately cause developmental,anatomical, functional and behavioral differences amongthese 80+ strains.ResultsWe generated approximately 13.2 and 38.9× whole-gen-ome short reads of DBA/2J females using Illumina GA2and ABI SOLiD massively parallel DNA sequencingplatforms. Comparing to the C57BL/6J reference gen-ome sequence, we identified over 4.5 million singlenucleotide polymorphisms (SNPs), including 84 non-sense and ~11,000 missense mutations, 78% of whichare novel. We also detected ~568,000 insertions anddeletions (indels) within single short reads and ~9,400between mate-paired reads. Approximately 300 inver-sions were detected by SOLiD mate-pair reads, 46 ofwhich span at least one exon. In addition, we identified~22,000 copy number variants (CNVs) in the range of1 Kb to 100 Kb (Figure 1).ConclusionOur study generates the first consensus sequence for theDBA/2J and creates a compendium of sequence andstructural variations that will be used by the communityof researchers who study complex traits in mouse mod-els.Thesequencedataprovideanovelresourcewithwhich to initiate reverse genetic analysis of complex