Hepatic gene expression variations in response to high-fat diet-induced impaired glucose tolerance using RNAseq analysis in collaborative cross mouse population.

Hepatic gene expression is known to differ between healthy and type 2 diabetes conditions. Identifying these variations will provide better knowledge to the development of gene-targeted therapies. The aim of this study is to assess diet-induced hepatic gene expression of susceptible versus resistant CC lines to T2D development. Next-generation RNA-sequencing was performed for 84 livers of diabetic and non-diabetic mice of 41 different CC lines (both sexes) following 12 weeks on high-fat diet (42% fat). Data analysis revealed significant variations of hepatic gene expression in diabetic versus non-diabetic mice with significant sex effect, where 601 genes were differentially expressed (DE) in overall population (males and females), 718 genes in female mice, and 599 genes in male mice. Top prioritized DE candidate genes were Lepr, Ins2, Mb, Ckm, Mrap2, and Ckmt2 for the overall population; for females-only group were Hdc, Serpina12, Socs1, Socs2, and Mb, while for males-only group were Serpine1, Mb, Ren1, Slc4a1, and Atp2a1. Data analysis for sex differences revealed 193 DE genes in health (Top: Lepr, Cav1, Socs2, Abcg2, and Col5a3), and 389 genes DE between diabetic females versus males (Top: Lepr, Clps, Ins2, Cav1, and Mrap2). Furthermore, integrating gene expression results with previously published QTL, we identified significant variants mapped at chromosomes at positions 36-49 Mb, 62-71 Mb, and 79-99 Mb, on chromosomes 9, 11, and 12, respectively. Our findings emphasize the complexity of T2D development and that significantly controlled by host complex genetic factors. As well, we demonstrate the significant sex differences between males and females during health and increasing to extent levels during disease/diabetes. Altogether, opening the venue for further studies targets the discovery of effective sex-specific and personalized preventions and therapies.

Crop wild relative populations of Beta vulgaris allow direct mapping of agronomically important genes.

Rapid identification of agronomically important genes is of pivotal interest for crop breeding. One source of such genes are crop wild relative (CWR) populations. Here we used a CWR population of <200 wild beets (B. vulgaris ssp. maritima), sampled in their natural habitat, to identify the sugar beet (Beta vulgaris ssp. vulgaris) resistance gene Rz2 with a modified version of mapping-by-sequencing (MBS). For that, we generated a draft genome sequence of the wild beet. Our results show the importance of preserving CWR in situ and demonstrate the great potential of CWR for rapid discovery of causal genes relevant for crop improvement. The candidate gene for Rz2 was identified by MBS and subsequently corroborated via RNA interference (RNAi). Rz2 encodes a CC-NB-LRR protein. Access to the DNA sequence of Rz2 opens the path to improvement of resistance towards rhizomania not only by marker-assisted breeding but also by genome editing.

Partial short-read sequencing of a highly inbred Iberian pig and genomics inference thereof.

Despite dramatic reduction in sequencing costs with the advent of next generation sequencing technologies, obtaining a complete mammalian genome sequence at sufficient depth is still costly. An alternative is partial sequencing. Here, we have sequenced a reduced representation library of an Iberian sow from the Guadyerbas strain, a highly inbred strain that has been used in numerous QTL studies because of its extreme phenotypic characteristics. Using the Illumina Genome Analyzer II (San Diego, CA, USA), we resequenced ∼ 1% of the genome with average 4 × depth, identifying 68,778 polymorphisms. Of these, 55,457 were putative fixed differences with respect to the assembly, based on the genome of a Duroc pig, and 13,321 were heterozygous positions within Guadyerbas. Despite being highly inbred, the estimate of heterozygosity within Guadyerbas was ∼ 0.78 kb(-1) in autosomes, after correcting for low depth. Nucleotide variability was consistently higher at the telomeric regions than on the rest of the chromosome, likely a result of increased recombination rates. Further, variability was 50% lower in the X-chromosome than in autosomes, which may be explained by a recent bottleneck or by selection. We divided the whole genome in 500 kb windows and we analyzed overrepresented gene ontology terms in regions of low and high variability. Multi organism process, pigmentation and cell killing were overrepresented in high variability regions and metabolic process ontology, within low variability regions. Further, a genome wide Hudson-Kreitman-Aguadé test was carried out per window; overall, variability was in agreement with neutral expectations.

Unique gene structure and paralogy define the 7D-cadherin family.

Cadherins are Ca2+-dependent transmembrane glycoproteins crucial for cell-cell adhesion in vertebrates and invertebrates. Classification of this superfamily due to their phylogenetic relationship is currently restricted to three major subfamilies: classical, desmosomal and protocadherins. Here we report evidence for a common phylogenetic origin of the kidney-specific Ksp- (Cdh16) and the intestine-specific LI-cadherin (Cdh17). Both genes consist of 18 exons and the positions of their exon-intron boundaries as well as their intron phases are perfectly conserved. We found an extensive paralogy of more than 40 megabases in mammals as well as teleost fish species encompassing the Ksp- and LI-cadherin genes. A comparable paralogy was not detected for other cadherin gene loci. These findings suggest that the Ksp- and LI-cadherin genes originated by chromosomal duplication early during vertebrate evolution and support our assumption that both proteins are paralogues within a separate cadherin family that we have termed 7D-cadherins.

Phylogenetic origin of LI-cadherin revealed by protein and gene structure analysis.

The intestine specific LI-cadherin differs in its overall structure from classical and desmosomal cadherins by the presence of seven instead of five cadherin repeats and a short cytoplasmic domain. Despite the low sequence similarity, a comparative protein structure analysis revealed that LI-cadherin may have originated from a five-repeat predecessor cadherin by a duplication of the first two aminoterminal repeats. To test this hypothesis, we cloned the murine LI-cadherin gene and compared its structure to that of other cadherins. The intron-exon organization, including the intron positions and phases, is perfectly conserved between repeats 3-7 of LI-cadherin and 1-5 of classical cadherins. Moreover, the genomic structure of the repeats 1-2 and 3-4 is identical for LI-cadherin and highly similar to that of the repeats 1-2 of classical cadherins. These findings strengthen our assumption that LI-cadherin originated from an ancestral cadherin with five domains by a partial gene duplication event.

Divergent genetic and epigenetic post-zygotic isolation mechanisms in Mus and Peromyscus.

Interspecific hybridization in the rodent genera Peromyscus and Mus results in abnormal placentation. In the Peromyscus interspecies hybrids, abnormal allelic interaction between an X-linked locus and the imprinted paternally expressed Peg3 locus was shown to cause the placental defects. In addition, loss-of-imprinting (LOI) of Peg3 was positively correlated with increased placental size. As in extreme cases this placental dysplasia constitutes a post-zygotic barrier against interspecies hybridization, this finding was the first direct proof that imprinted genes may be important in speciation and thus in evolution. In the Mus interspecies hybrids, a strong role of an X-linked locus in placental dysplasia has also been detected. However, here we show by backcross and allele specific expression analyses that neither LOI of Peg3 nor abnormal interactions between Peg3 and an X-linked locus are involved in generating placental dysplasia in Mus hybrids, although the placental phenotypes observed in the two genera seem to be identical. In contrast to this, another dysgenesis effect common to Peromyscus and Mus hybrids, altered foetal growth, is caused at least in part by the same X-chromosomal regions in both genera. These findings first underline the strong involvement of the X-chromosome in the genetics of speciation. Secondly, they indicate that disruption of epigenetic states, such as LOI, at specific loci may be involved in hybrid dysgenesis effects in one group, but not in another. Thus, we conclude that even in closely related groups divergent molecular mechanisms may be involved in the production of phenotypically similar post-zygotic barriers against hybridization.

Advanced integrated mouse YAC map including BAC framework.

Functional characterization of the mouse genome requires the availability of a comprehensive physical map to obtain molecular access to chromosomal regions of interest. Positional cloning remains a crucial way of linking phenotype with particular genes. A key step and frequent stumbling block in positional cloning is making a contig of a genetically defined candidate region. The most efficient first step is isolating YAC (Yeast Artificial Chromosome) clones. A robust, detailed YAC contig map is thus an important tool. Employing Interspersed Repetitive Sequence (IRS)-PCR genomics, we have generated an advanced second-generation YAC contig map of the mouse genome that doubles both the depth of clones and the density of markers available. In addition to the primarily YAC-based map, we located 1942 BAC (Bacterial Artificial Chromosome) clones. This allows us to present for the first time a dense framework of BACs spanning the genome of the mouse, which, for instance, can serve as a nucleus for genomic sequencing. Four large-insert mouse YAC libraries from three different strains are included in our data, and our analysis incorporates the data of Hunter et al. and Nusbaum et al. There is a total of 20,205 markers on the final map, 12,033 from our own data, and a total of 56,093 YACs, of which 44,401 are positive for more than one marker.

UniGene cDNA array-based monitoring of transcriptome changes during mouse placental development.

The placenta is a highly specialized organ essential for embryonic growth and development. Here, we have applied cDNA subtraction between extraembryonic tissues of early- (day 7.5 of gestation) and late-stage embryos (day 17.5) to generate stage-specific cDNA pools that were used for screening of high-density mouse UniGene cDNA arrays containing 25,000 clones. A total of 638 clones were identified, 488 with the e7.5-specific probe and 150 with the e17.5-specific probe. Importantly, 363/638 (56.9%) of the hybridizing clones were not known to be expressed during placental development before. Differential regulation was confirmed by Northern blot and in situ hybridization for a total of 44/44 of positive clones. Thus, this combination of cDNA subtraction and array hybridization was highly successful for identification of genes expressed and regulated during placental development. These included growth factors and receptors, components of the transcriptional and translational machinery, cell cycle regulators, molecular chaperones, and cytoskeletal elements. The extensive in situ hybridization analysis revealed extraembryonic structures with a high density of differentially expressed genes, most strikingly the ectoplacental cone and the spongiotrophoblast. This large-scale identification of genes regulated during placentogenesis is extremely useful to further elucidate the molecular basis of extraembryonic development.

High-throughput scanning of the rat genome using interspersed repetitive sequence-PCR markers.

We report the establishment of a hybridization-based marker system for the rat genome based on the PCR amplification of interspersed repetitive sequences (IRS). Overall, 351 IRS markers were mapped within the rat genome. The IRS marker panel consists of 210 nonpolymorphic and 141 polymorphic markers that were screened for presence/absence polymorphism patterns in 38 different rat strains and substrains that are commonly used in biomedical research. The IRS marker panel was demonstrated to be useful for rapid genome screening in experimental rat crosses and high-throughput characterization of large-insert genomic library clones. Information on corresponding YAC clones is made available for this IRS marker set distributed over the whole rat genome. The two existing rat radiation hybrid maps were integrated by placing the IRS markers in both maps. The genetic and physical mapping data presented provide substantial information for ongoing positional cloning projects in the rat.

cDNA subtraction cloning reveals novel genes whose temporal and spatial expression indicates association with trophoblast invasion.

Trophoblast invasion is a critical process in development of most mammals that shares similarities with the invasive behavior of tumor cells. In the present investigation, a cDNA subtraction library was constructed between invasive trophoblast at day 8 of murine development and mature noninvasive placenta at day 18 of gestation. One of the differentially expressed clones, Epcs26, was mapped to the X chromosome and revealed no homology to any known gene. It was predominantly expressed in parietal endoderm, undifferentiated cells of the ectoplacental cone, and a few trophoblast giant cells. Another gene, designated Epcs50, was mapped to chromosome 19. It exhibited homologies to the mouse Mps1 gene and, like Mps1, may have a distant relationship to the lytic protein perforin. High expression was detected in parietal endoderm cells and in a subset of secondary trophoblast giant cells. Two sequences, Epcs24 and Epcs68, exhibited an extensive open reading frame that shared the common features of the cysteine proteinase cathepsin L. Expression was confined to an undefined subpopulation of trophoblast giant cells. Both genes were mapped to chromosome 13 in close proximity to cathepsins L and J. The known functions of MPS1 and cathepsin L proteins indicate that the related proteins EPCS50, EPCS24, and EPCS68 participate in conferring invasive properties to the mouse trophoblast.

The mouse Psma1 gene coding for the alpha-type C2 proteasome subunit: structural and functional analysis, mapping, and colocalization with Pde3b on mouse chromosome 7.

We have isolated and functionally characterized the mouse gene for the C2 subunit of the 20S proteasome. The gene contains 10 exons distributed over a region of 12 kb on the distal end of mouse chromosome 7. Its exon-intron structure differs from those of the other few known proteasome genes. Transfection assays revealed that 1.5 kb of 5' flanking sequence is active as promoter in cultured myoblasts. Deletion reporter constructs narrowed this presumptive promoter region to within 450 bp upstream of the translation initiation site. Several consensus motifs for transcription factor binding sites were identified in this upstream region of the gene. Psma1 was mapped to mouse chromosome 7 using the interspecific backcross DNA panels from The Jackson Laboratory. Additional mapping studies showed that the mouse genes Psma1 and Pde3b are closely linked, residing between cM 53 and 53.3 in a region syntenic to human chromosome 11p15. Our results extend the structural and functional analysis of genes encoding the 20S proteasome subunits and provide the basis for the study of their regulation.

Molecular cloning and characterization of murine ICOS and identification of B7h as ICOS ligand.

Human ICOS (huICOS) is a T cell-specific molecule structurally related to CD28 and CTLA-4 with potent co-stimulatory activities on T cell proliferation, cytokine induction and T cell help for B cells. We have now cloned and characterized murine ICOS (muICOS). muICOS mRNA of 1.5 kb and 3.3 kb encodes a protein with a deduced molecular mass of 20.3 kDa, which is 71.7 % identical to huICOS. On the cell surface, muICOS is expressed as a disulfide-linked, glycosylated homodimer of 47-57 kDa, with subunits of approximately 26 kDa. With a panel of monoclonal antibodies we have determined the expression of muICOS in vitro and in vivo. Following activation of splenic T cells via CD3, muICOS became detectable at 12 h and reached a maximum of expression at around 48 h, thus exhibiting expression kinetics similar to huICOS. In vivo, muICOS was found to be substantially expressed in the thymic medulla and in the germinal centers and T cell zones of lymph nodes and Peyer's patches. Non-lymphoid tissue was ICOS negative. The muICOS gene was mapped to a region of chromosome 1 also harboring the CD28 and CTLA-4 genes. Using recombinant chimeric muICOS-Ig we determined that B7h, a recently cloned B7-like molecule, is a ligand for muICOS.

Interspersed repetitive sequence (IRS)-PCR for typing of whole genome radiation hybrid panels.

The typing of a radiation hybrid (RH) panel is generally achieved using a unique primer pair for each marker. We here describe a complementing approach utilizing IRS-PCR. Advantages of this technology include the use of a single universal primer to specify any locus, the rapid typing of RH lines by hybridization, and the conservative use of hybrid DNA. The technology allows the mapping of a clone without the requirement for STS generation. To test the technique, we have mapped 48 BAC clones derived from mouse chromosome 12 which we mostly identified using complex probes. As mammalian genomes are repeat-rich, the technology can easily be adapted to species other than mouse.

Towards identification of individual homologous chromosomes: comparative genomic hybridization and spectral karyotyping discriminate between paternal and maternal euchromatin in Mus musculus x M. spretus interspecific hybrids.

We have developed an in situ technique to label individual euchromatic chromosome arms in interspecific crosses between Mus musculus (MMU) and M. spretus (MSP). The MMU and MSP genomes diverged 2-3 million years ago and show an overall sequence divergence of approximately 1%. Comparative hybridization of MMU versus MSP DNA and subsequent spectral analysis of the euchromatic hybridization profiles discriminated between maternal (MMU) and paternal (MSP) chromosomes in F(1) hybrids. Dispersed repetitive DNA elements were the preferred hybridization target of MMU DNA on maternal chromosomes and of MSP DNA on paternal chromosomes. Differences in centromeric satellite DNAs were detected by conventional fluorescence in situ hybridization and served as internal controls. Our experiments suggest that it is possible, in principle, to discriminate between paternal and maternal chromosomes on the basis of sequence differences.

Panel of microsatellite markers for whole-genome scans and radiation hybrid mapping and a mouse family tree.

To facilitate whole-genome scan experiments, we selected a panel of 128 microsatellite markers on the basis of spacing and polymorphism in the strains DBA/2, BALB/c, AKR, C57BL/6, C57BL/10, A/J, C3H, 129/J, SJL/J, JF1, and PWB. Many of the primer pairs were redesigned for better performance. The last four strains were not characterized previously using these markers. JF1 and PWB are particularly interesting for intersubspecific crosses offering high polymorphism. We provide allele size data for the markers on these strains and add them to the emerging radiation hybrid framework map, which is not continuous except for chromosome 17 and 13. Information on the interrelationships of strains is useful both because of the importance of polymorphism in designing crosses and the background in assessing phenotypes. Microsatellites offer a widely dispersed, selectively neutral set of characters that lends itself conceptually to parsimony methods of analysis. The microsatellite allele size data were recoded as binary discrete characters in such a way that adjacent sizes differ by one step. Trees were generated using a Wagner parsimony method. As expected, the non-Mus domesticus strains, PWB (musculus) and JF1 (molossinus), are excluded from the domesticus strains. Among the domesticus strains, C57BL/6 and C57BL/10 (derived from the same founding pair) form a strongly supported group, as do C3H, A/J, and BALB/c (derived from the Bagg albino stock). No unique branching order for SJL/J, AKR, and DBA/2 is strongly supported, which may reflect a complicated history. Strain 129/J is clearly placed as the most deeply diverged of the domesticus strains represented.

Identification and characterization of G90, a novel mouse RNA that lacks an extensive open reading frame.

We describe the cloning and characterization of the murine G90 gene, identified by subtractive hybridization based on the differential presence of its transcript in large and small intestine. The full-length cDNA and genomic sequences were cloned and found to produce a 1.5kb transcript that is polyadenylated but has no open reading frame larger than 249bp. The G90 gene was mapped to the proximal region of mouse chromosome 6. Expression analysis by Northern blotting showed that G90 is transcribed at very high levels in the small intestine and at lower levels in large intestine, testis and kidney of the mouse. In situ hybridization analysis on sections of small and large intestine and testis showed that G90 transcripts are present only in post-mitotic cells.

Technology development at the interface of proteome research and genomics: mapping nonpolymorphic proteins on the physical map of mouse chromosomes.

Data obtained from protein spots by peptide mass fingerprinting are used to identify the corresponding genes in sequence databases. The relevant cDNAs are obtained as clones from the Integrated Molecular Analysis of Genome Expression (I.M.A.G.E.) consortium. Mapping of I.M.A.G.E. clones is performed in two steps: first, cDNA clones are hybridized against a 10-hit genomic mouse bacterial artificial chromosome (BAC) library. Second, interspersed repetitive sequence polymerase chain reaction (IRS-PCR) using a single primer directed against the mouse B1 repeat element is performed on BACs. As each cDNA detects several BACs, and each individual BAC has a 50% chance to recover an IRS-PCR fragment, the majority of cDNAs produce at least a single IRS-PCR fragment. Individual IRS fragments are hybridized against high-density spotted filter grids containing the three-dimensional permutated pools of yeast artificial chromosome (YAC) library resources that are currently being used to construct a physical map of the mouse genome. IRS fragments that hybridize to YAC clones already placed into contigs immediately provide highly precise map positions. This technology therefore is able to draw links between proteins detected by 2-D gel electrophoresis and the corresponding gene loci in the mouse genome.

Comparative mapping of mouse and rat chromosomes by fluorescence in situ hybridization.

Comparative fluorescence in situ hybridization mapping using DNA libraries from flow-sorted mouse chromosomes and region-specific mouse BAC clones on rat chromosomes reveals chromosomal homologies between mouse (Mus musculus, MMU) and rat (Rattus norvegicus, RNO). Each of the MMU 2, 3, 4, 6, 7, 9, 12, 14, 15, 16, 18, 19, and X chromosomes paints only a single rat chromosome or chromosome segment and, thus, the chromosomes are largely conserved between the two species. In contrast, the painting probes for MMU chromosomes 1, 5, 8, 10, 11, 13, and 17 produce split hybridization signals in the rat, disclosing evolutionary chromosome rearrangements. Comparative mapping data delineate several large linkage groups on RNO 1, 2, 4, 7, and 14 that are conserved in human but diverged in the mouse. On the other hand, there are linkage groups in the mouse, i.e., on MMU 1, 8, 10, and 11, that are disrupted in both rat and human. In addition, we have hybridized probes for Nap2, p57, Igf2, H19, and Sh3d2c from MMU 7 to RNO 1q and found the orientation of the imprinting gene cluster and Sh3d2c to be the same in mouse and rat. Hybridization of rat genomic DNA shows blocks of (rat-specific) repetitive sequences in the pericentromeric region of RNO chromosomes 3-5, 7-13, and 20; on the short arms of RNO chromosomes 3, 12, and 13; and on the entire Y chromosome.

Expression of the von Hippel-Lindau-binding protein-1 (Vbp1) in fetal and adult mouse tissues.

The von Hippel-Lindau (VHL) tumour suppressorgene product is believed to be involved in the down-regulation of transcriptional elongation by preventing the association of elongin B and C with the catalytic subunit elongin A. Alterations in the human VHL gene lead to VHL disease which is associated with various rare neoplasias, including haemangioblastoma of the central nervous system, retinal angioma, clear cell renal carcinoma and pheochromocytoma. Recently, a protein (VBP1) was isolated that was found to bind to the VHL protein in vivo. We have used the murine Vbp1 homologous cDNA to investigate the expression of the Vbp1 mRNA in the mouse by in situ hybridization and northern blot analysis. In fetal stages between days 9 and 18 of gestation, Vbp1 was expressed mainly in the central nervous system, retina and liver. In addition, at day 12, high expression was observed in the labyrinthine region of the placenta. In later stage placentas, Vbp1 expression was, however, considerably reduced. Northern blot analysis of adult mouse tissues showed that Vbp1 was ubiquitously expressed. In situ analysis on several adult tissues showed that in most tissues, transcripts were evenly distributed. In brain, eye, kidney and intestine, however, Vbp1 was expressed in specific cell types. Moreover, expression of the human VBP1 gene was investigated in cerebellum and in various tumours of VHL patients encompassinghaemangioblastomas, renal cell carcinomas and pheochromocytomas. In all of these tissues, VBP1 was ubiquitously expressed at low levels. However, no consistent differences in VBP1 expression levels could be detected between tumours and normal tissue. Mapping of the murine Vbp1 gene revealed conserved chromosomal localization between mouse and human in a region homologous to human Xq28.