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.

Phylogeny of horse chromosome 5q in the genus Equus and centromere repositioning

Horses, asses and zebras belong to the genus Equus and are the only extant species of the family Equidae in the order Perissodactyla. In a previous work we demonstrated that a key factor in the rapid karyotypic evolution of this genus was evolutionary centromere repositioning, that is, the shift of the centromeric function to a new position without alteration of the order of markers along the chromosome. In search of previously undiscovered evolutionarily new centromeres, we traced the phylogeny of horse chromosome 5, analyzing the order of BAC markers, derived from a horse genomic library, in 7 Equus species (E. caballus, E. hemionus onager, E. kiang, E. asinus, E. grevyi, E. burchelli and E. zebra hartmannae). This analysis showed that repositioned centromeres are present in E. asinus (domestic donkey, EAS) chromosome 16 and in E. burchelli (Burchell's zebra, EBU) chromosome 17, confirming that centromere repositioning is a strikingly frequent phenomenon in this genus. The observation that the neocentromeres in EAS16 and EBU17 are in the same chromosomal position suggests that they may derive from the same event and therefore, E. asinus and E. burchelli may be more closely related than previously proposed; alternatively, 2 centromere repositioning events, involving the same chromosomal region, may have occurred independently in different lineages, pointing to the possible existence of hot spots for neocentromere formation. Our comparative analysis also showed that, while E. caballus chromosome 5 seems to represent the ancestral configuration, centric fission followed by independent fusion events gave rise to 3 different submetacentric chromosomes in other Equus lineages.

Cloning of a full-length cDNA encoding bovine interleukin 4 by the polymerase chain reaction

A human interleukin 4 (hIL-4)-encoding cDNA (hIL4) probe was used to screen a bovine genomic library, and three clones containing sequences with homology to the human and mouse IL4 cDNAs were isolated. Sequence information obtained from one of these genomic clones was used to design an oligodeoxyribonucleotide primer corresponding to the transcription start point region for use in the polymerase chain reaction (PCR). The PCR-RACE protocol, designed for the rapid amplification of cDNA ends, was successfully used to generate a full-length bovine IL4 (bIL4) cDNA clone from polyadenylated RNA isolated from concanavalin A-stimulated bovine lymph node cells. The bIL4 cDNA is 570 bp in length and contains an open reading frame of 405 nucleotides (nt), coding for a 15.1-kDa precursor of 135 amino acids (aa), which should be reduced to 12.6 kDa for unglycosylated bIL4 after cleavage of a putative hydrophobic leader sequence of 24 aa. The aa sequence contains one possible Asn-linked glycosylation site. Bovine IL4 is shorter than mouse (mIL4) and hIL4, because of a 51-nt deletion in the coding region. Comparison of the overall nt and deduced aa sequences shows a greater homology of bIL4 with hIL4 than with mIL4. This homology is not evenly distributed, however, with the nt sequences 5' and 3' of the coding region showing a much greater homology between all three species than the coding sequence.

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.

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.

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.

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.

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.