Organisation of the equine immunoglobulin constant heavy chain genes. II. Equine cgamma genes.

The number of immunoglobulin G constant heavy chain genes (cgamma genes) varies broadly among mammalian species, reflecting structural and functional differences between expressed immunoglobulin G (IgG) isotypes and allotypes. Up to now equine IgG isotypes have been defined only at the biochemical and serological level. It is still not clear how many IgG isotypes exist in horses and whether there are any allotypes. Here, we describe the isolation and characterisation of equine cgamma genes. An equine genomic lambda phage library was screened with a human cgamma4 probe. Cross-hybridising equine cgamma sequences were cloned twice and characterised by restriction mapping with the human cgamma4 and a murine sgamma1 probe. Genomic equine DNA probes for both, cgamma genes and corresponding switch regions (sgamma), were isolated and used for a more detailed BamHI restriction analysis, comparing genomic DNA of various horses. This analysis reveals the existence of at least five, or probably six cgamma genes in the equine haploid genome. Beside the porcine system, this is the highest number of cgamma genes described for any mammalian species. Moreover, for two of these cgamma genes, BamHI restriction fragment length polymorphism became evident.

Proteins encoded in genomic regions associated with immune-mediated disease physically interact and suggest underlying biology

Genome-wide association studies (GWAS) have defined over 150 genomic regions unequivocally containing variation predisposing to immune-mediated disease. Inferring disease biology from these observations, however, hinges on our ability to discover the molecular processes being perturbed by these risk variants. It has previously been observed that different genes harboring causal mutations for the same Mendelian disease often physically interact. We sought to evaluate the degree to which this is true of genes within strongly associated loci in complex disease. Using sets of loci defined in rheumatoid arthritis (RA) and Crohn's disease (CD) GWAS, we build protein-protein interaction (PPI) networks for genes within associated loci and find abundant physical interactions between protein products of associated genes. We apply multiple permutation approaches to show that these networks are more densely connected than chance expectation. To confirm biological relevance, we show that the components of the networks tend to be expressed in similar tissues relevant to the phenotypes in question, suggesting the network indicates common underlying processes perturbed by risk loci. Furthermore, we show that the RA and CD networks have predictive power by demonstrating that proteins in these networks, not encoded in the confirmed list of disease associated loci, are significantly enriched for association to the phenotypes in question in extended GWAS analysis. Finally, we test our method in 3 non-immune traits to assess its applicability to complex traits in general. We find that genes in loci associated to height and lipid levels assemble into significantly connected networks but did not detect excess connectivity among Type 2 Diabetes (T2D) loci beyond chance. Taken together, our results constitute evidence that, for many of the complex diseases studied here, common genetic associations implicate regions encoding proteins that physically interact in a preferential manner, in line with observations in Mendelian disease.

The genomic Echinococcus microsatellite EmsB sequences: from a molecular marker to the epidemiological tool

In the field of molecular and epidemiological parasitology, characterization of fast evolving genetic markers appears as an important challenge to consider the diversity and genetic structure of parasites. The study of respective populations can help us to understand their adaptive strategies to survive and perpetuate the species within different host populations, all trying to resist infection. In the past, the relative monomorphic features of Echinococcus multilocularis, the causative agent of alveolar echinococcosis and a severe human parasitic disease, did not stimulate studies dealing with the genetic variability of Echinococcus species or respective populations. A recently developed, characterized and validated original multilocus microsatellite, named EmsB, tandemly repeated in the genome, offered an additional opportunity for this line of investigation. We have compiled in this review new insights brought by this molecular tracker on the transmission activity of Echinococcus among different hosts and at different geographical scales.

Genomic Data Reveal a Complex Making of Humans

In the last few years, two paradigms underlying human evolution have crumbled. Modern humans have not totally replaced previous hominins without any admixture, and the expected signatures of adaptations to new environments are surprisingly lacking at the genomic level. Here we review current evidence about archaic admixture and lack of strong selective sweeps in humans. We underline the need to properly model differential admixture in various populations to correctly reconstruct past demography. We also stress the importance of taking into account the spatial dimension of human evolution, which proceeded by a series of range expansions that could have promoted both the introgression of archaic genes and background selection.

Evolutionary forces shaping genomic islands of population differentiation in humans

Background Levels of differentiation among populations depend both on demographic and selective factors: genetic drift and local adaptation increase population differentiation, which is eroded by gene flow and balancing selection. We describe here the genomic distribution and the properties of genomic regions with unusually high and low levels of population differentiation in humans to assess the influence of selective and neutral processes on human genetic structure. Methods Individual SNPs of the Human Genome Diversity Panel (HGDP) showing significantly high or low levels of population differentiation were detected under a hierarchical-island model (HIM). A Hidden Markov Model allowed us to detect genomic regions or islands of high or low population differentiation. Results Under the HIM, only 1.5% of all SNPs are significant at the 1% level, but their genomic spatial distribution is significantly non-random. We find evidence that local adaptation shaped high-differentiation islands, as they are enriched for non-synonymous SNPs and overlap with previously identified candidate regions for positive selection. Moreover there is a negative relationship between the size of islands and recombination rate, which is stronger for islands overlapping with genes. Gene ontology analysis supports the role of diet as a major selective pressure in those highly differentiated islands. Low-differentiation islands are also enriched for non-synonymous SNPs, and contain an overly high proportion of genes belonging to the 'Oncogenesis' biological process. Conclusions Even though selection seems to be acting in shaping islands of high population differentiation, neutral demographic processes might have promoted the appearance of some genomic islands since i) as much as 20% of islands are in non-genic regions ii) these non-genic islands are on average two times shorter than genic islands, suggesting a more rapid erosion by recombination, and iii) most loci are strongly differentiated between Africans and non-Africans, a result consistent with known human demographic history.

Multi-locus microsatellite analysis supports the hypothesis of an autochthonous focus of Echinococcus multilocularis in northern Italy

Echinococcus multilocularis is characterised by a wide geographical distribution, encompassing three continents (North America, Asia and Europe) yet very low genetic variability is documented. Recently, this parasite has been detected in red foxes (Vulpes vulpes) circulating in an Alpine region of Italy, close to Austria. This finding raised the question as to whether an autochthonous cycle exists in Italy or whether the infected foxes originated from the neighbouring regions of Austria. Studies have shown that multi-locus microsatellite analysis can identify genomic regions carrying mutations that result in a local adaptation. We used a tandem repeated multi-locus microsatellite (EmsB) to evaluate the genetic differences amongst adult worms of E. multilocularis collected in Italy, worms from neighbouring Austria and from other European and extra-European countries. Fluorescent PCR was performed on a panel of E. multilocularis samples to assess intra-specific polymorphism. The analysis revealed four closed genotypes for Italian samples of E. multilocularis which were unique compared with the other 25 genotypes from Europe and the five genotypes from Alaska. An analysis in the Alpine watershed, comparing Italian adult worms with those from neighbouring areas in Austria, showed a unique cluster for Italian samples. This result supports the hypothesis of the presence of an autochthonous cycle of E. multilocularis in Italy. EmsB can be useful for 'tracking' the source of infection of this zoonotic parasite and developing appropriate measures for preventing or reducing the risk of human alveolar echinococcosis.

A non-coding genomic duplication at the HMX1 locus is associated with crop ears in highland cattle.

Highland cattle with congenital crop ears have notches of variable size on the tips of both ears. In some cases, cartilage deformation can be seen and occasionally the external ears are shortened. We collected 40 cases and 80 controls across Switzerland. Pedigree data analysis confirmed a monogenic autosomal dominant mode of inheritance with variable expressivity. All affected animals could be traced back to a single common ancestor. A genome-wide association study was performed and the causative mutation was mapped to a 4 Mb interval on bovine chromosome 6. The H6 family homeobox 1 (HMX1) gene was selected as a positional and functional candidate gene. By whole genome re-sequencing of an affected Highland cattle, we detected 6 non-synonymous coding sequence variants and two variants in an ultra-conserved element at the HMX1 locus with respect to the reference genome. Of these 8 variants, only a non-coding 76 bp genomic duplication (g.106720058_106720133dup) located in the conserved region was perfectly associated with crop ears. The identified copy number variation probably results in HMX1 misregulation and possible gain-of-function. Our findings confirm the role of HMX1 during the development of the external ear. As it is sometimes difficult to phenotypically diagnose Highland cattle with slight ear notches, genetic testing can now be used to improve selection against this undesired trait.

Genomic cloning and expression of three murine UDP-galactose: beta-N-acetylglucosamine beta1,3-galactosyltransferase genes.

Based on the detection of expressed sequence tags that are similar to known galactosyltransferase sequences, we have isolated three novel UDP-galactose:beta-N-acetylglucosamine beta1, 3-galactosyltransferase (beta3GalT) genes from a mouse genomic library. The three genes, named beta3GalT-I, -II, and -III, encode type II transmembrane proteins of 326, 422, and 331 amino acids, respectively. The three proteins constitute a distinct subfamily as they do not share any sequence identity with other eucaryotic galactosyltransferases. Also, the entire protein-coding region of the three beta3GalT genes was contained in a single exon, which contrasts with the genomic organization of the beta1,4- and alpha1, 3-galactosyltransferase genes. The three beta3GalT genes were mainly expressed in brain tissue. The expression of the full-length murine genes as recombinant baculoviruses in insect cells revealed that the beta3GalT enzymes share the same acceptor specificity for beta-linked GlcNAc, although they differ in their Km for this acceptor and the donor UDP-Gal. The identification of beta3GalT genes emphasizes the structural diversity present in the galactosyltransferase gene family.

Extensive Copy-Number Variation of Young Genes across Stickleback Populations

Duplicate genes emerge as copy-number variations (CNVs) at the population level, and remain copy-number polymorphic until they are fixed or lost. The successful establishment of such structural polymorphisms in the genome plays an important role in evolution by promoting genetic diversity, complexity and innovation. To characterize the early evolutionary stages of duplicate genes and their potential adaptive benefits, we combine comparative genomics with population genomics analyses to evaluate the distribution and impact of CNVs across natural populations of an eco-genomic model, the three-spined stickleback. With whole genome sequences of 66 individuals from populations inhabiting three distinct habitats, we find that CNVs generally occur at low frequencies and are often only found in one of the 11 populations surveyed. A subset of CNVs, however, displays copy-number differentiation between populations, showing elevated within-population frequencies consistent with local adaptation. By comparing teleost genomes to identify lineage-specific genes and duplications in sticklebacks, we highlight rampant gene content differences among individuals in which over 30% of young duplicate genes are CNVs. These CNV genes are evolving rapidly at the molecular level and are enriched with functional categories associated with environmental interactions, depicting the dynamic early copy-number polymorphic stage of genes during population differentiation.

Genome-Wide Analyses Suggest Mechanisms Involving Early B-Cell Development in Canine IgA Deficiency.

Immunoglobulin A deficiency (IgAD) is the most common primary immune deficiency disorder in both humans and dogs, characterized by recurrent mucosal tract infections and a predisposition for allergic and other immune mediated diseases. In several dog breeds, low IgA levels have been observed at a high frequency and with a clinical resemblance to human IgAD. In this study, we used genome-wide association studies (GWAS) to identify genomic regions associated with low IgA levels in dogs as a comparative model for human IgAD. We used a novel percentile groups-approach to establish breed-specific cut-offs and to perform analyses in a close to continuous manner. GWAS performed in four breeds prone to low IgA levels (German shepherd, Golden retriever, Labrador retriever and Shar-Pei) identified 35 genomic loci suggestively associated (p <0.0005) to IgA levels. In German shepherd, three genomic regions (candidate genes include KIRREL3 and SERPINA9) were genome-wide significantly associated (p <0.0002) with IgA levels. A ~20kb long haplotype on CFA28, significantly associated (p = 0.0005) to IgA levels in Shar-Pei, was positioned within the first intron of the gene SLIT1. Both KIRREL3 and SLIT1 are highly expressed in the central nervous system and in bone marrow and are potentially important during B-cell development. SERPINA9 expression is restricted to B-cells and peaks at the time-point when B-cells proliferate into antibody-producing plasma cells. The suggestively associated regions were enriched for genes in Gene Ontology gene sets involving inflammation and early immune cell development.

A 4 Mb high resolution BAC contig on bovine chromosome 1q12 and comparative analysis with human chromosome 21q22

The bovine RPCI-42 BAC library was screened to construct a sequence-ready ~4 Mb single contig of 92 BAC clones on BTA 1q12. The contig covers the region between the genes KRTAP8P1 and CLIC6. This genomic segment in cattle is of special interest as it contains the dominant gene responsible for the hornless or polled phenotype in cattle. The construction of the BAC contig was initiated by screening the bovine BAC library with heterologous cDNA probes derived from 12 human genes of the syntenic region on HSA 21q22. Contig building was facilitated by BAC end sequencing and chromosome walking. During the construction of the contig, 165 BAC end sequences and 109 single-copy STS markers were generated. For comparative mapping of 25 HSA 21q22 genes, genomic PCR primers were designed from bovine EST sequences and the gene-associated STSs mapped on the contig. Furthermore, bovine BAC end sequence comparisons against the human genome sequence revealed significant matches to HSA 21q22 and allowed the in silico mapping of two new genes in cattle. In total, 31 orthologues of human genes located on HSA 21q22 were directly mapped within the bovine BAC contig, of which 16 genes have been cloned and mapped for the first time in cattle. In contrast to the existing comparative bovine-human RH maps of this region, these results provide a better alignment and reveal a completely conserved gene order in this 4 Mb segment between cattle, human and mouse. The mapping of known polled linked BTA 1q12 microsatellite markers allowed the integration of the physical contig map with existing linkage maps of this region and also determined the exact order of these markers for the first time. Our physical map and transcript map may be useful for positional cloning of the putative polled gene in cattle.

Identification of a new genomic hot spot of evolutionary diversification of protein function.

Establishment of phylogenetic relationships remains a challenging task because it is based on computational analysis of genomic hot spots that display species-specific sequence variations. Here, we identify a species-specific thymine-to-guanine sequence variation in the Glrb gene which gives rise to species-specific splice donor sites in the Glrb genes of mouse and bushbaby. The resulting splice insert in the receptor for the inhibitory neurotransmitter glycine (GlyR) conveys synaptic receptor clustering and specific association with a particular synaptic plasticity-related splice variant of the postsynaptic scaffold protein gephyrin. This study identifies a new genomic hot spot which contributes to phylogenetic diversification of protein function and advances our understanding of phylogenetic relationships.

Genomic knock-out of rancRNAs in the archaeon H. volcanii

As translation is the final step in gene expression it is particularly important to understand the processes involved in translation regulation. It was shown in the last years that a class of RNA, the non-protein-coding RNAs (ncRNAs), is involved in regulation of gene expression via various mechanisms [1]. Herein included is the prominent example of gene silencing caused by micro RNAs (miRNAs) and small interfering RNAs (siRNAs). Almost all of these ncRNA discovered so far target the mRNA in order to modulate protein biosynthesis, this is rather unexpected considering the crucial role of the ribosome during gene expression. However, recent data from our laboratory showed that there is a new class of RNAs among the well-studied ncRNAs that target the ribosome itself [2,3]. These so called ribosome-associated ncRNAs (rancRNAs) have an impact on translation regulation, mainly by interfering / modulating the rate of protein biosynthesis. Recent studies show the presence of small regulatory RNAs (sRNAs) in archaea which are involved in many biological processes including stress response and metabolic regulation [4]. To date the biological function and the targets of these archaeal sRNAs are only described for a few examples. There are reports of sRNAs binding to the 5’ as well as to the 3’ of mRNAs [5,6]. In addition to these findings, a tRNA derived fragment (tRF) of Valine tRNA was found in a genomic screen of RNAs associated with the ribosome in H. volcanii in our laboratory [3]. This Valine tRF seems to be processed in a stress-dependent manner and showed in vitro binding to the ribosome and inhibited in vitro translation. These results showed that Valine tRF is capable to regulate translation in H. volcanii by targeting the ribosome. The main goal of this project is to identify and describe novel potential regulatory rancRNAs in H. volcanii with the focus on intergenic candidates. Northern blot analyses already revealed interactions with the ribosome and showed differential expression patterns in response to stress conditions. To investigate the biological relevance of some of the ribosome-associated ncRNA candidates, knock-out and phenotypic characterization studies are done. The genomic knock out of a hypothetical ORF (198nt), where one putative rancRNA candidate (46nt) named IG33 was detected in the library at the beginning of the ORF, showed interesting growth phenotype under specific stress conditions. Furthermore a strain with an introduced start to stop codon mutation in this hypothetical ORF still shows the same phenotype indicating that rather the missing protein than the missing sRNA causes this growth phenotype.

Profiling invasiveness in head and neck cancer: recent contributions of genomic and transcriptomic approaches.

High-throughput molecular profiling approaches have emerged as precious research tools in the field of head and neck translational oncology. Such approaches have identified and/or confirmed the role of several genes or pathways in the acquisition/maintenance of an invasive phenotype and the execution of cellular programs related to cell invasion. Recently published new-generation sequencing studies in head and neck squamous cell carcinoma (HNSCC) have unveiled prominent roles in carcinogenesis and cell invasion of mutations involving NOTCH1 and PI3K-patwhay components. Gene-expression profiling studies combined with systems biology approaches have allowed identifying and gaining further mechanistic understanding into pathways commonly enriched in invasive HNSCC. These pathways include antigen-presenting and leucocyte adhesion molecules, as well as genes involved in cell-extracellular matrix interactions. Here we review the major insights into invasiveness in head and neck cancer provided by high-throughput molecular profiling approaches.