Defining the regions of homology between human chromosome 16p12--13 and mouse chromosomes

In this study, the evolutionary relationship between human chromosome 16p12-p13 and mouse chromosomes was investigated by determining the order of marker loci in the region and then identifying the chromosomal locations of the homologous loci in mice. Eighteen genes from human 16 were mapped to fifteen subchromosomal regions by a variety of mapping approaches.^ Thirteen of the genes were mapped in the mouse. Linkage analysis with backcross mice and segregation analysis in a mouse - Chinese Hamster Ovary (CHO) somatic cell hybrid panel informative for different regions of mouse genome were used. The results assigned the thirteen genes to three different mouse chromosomes.^ A group of six genes on mouse 16 was found to be closely linked to Scid. The order of Myh11 and Mrp remains ambiguous since no recombination was detected in backcross analysis. Their relative position in human is also uncertain since they were shown to be very close to each other. For the other mouse loci, an unambiguous gene order could be determined and was found to be identical to that in human. Therefore, they comprise a new conserved linkage group between the two species. The orientation of the group was inverted relative to the centromeres, i.e. the proximal loci in one species become distal in another. The size of the group was estimated to be from 4.4 to 8 Mb and 10 to 32 cM in human. In mouse, it was about 21 cM in the backcross analysis. The two boundaries of the conserved linkage were defined within a 1 Mb range. It is now possible to predict the locations of mouse homologs for some human disease genes based on their locations on human 16p.^ The six human 16p genes that map to MMU7 showed a different gene order in mouse than in human. No recombination was found between Crym and Umod while Crym was distal to D16S79A and proximal to D16S92. The location of Stp and Cdr2 with respect to the above four loci was not determined since they were not mapped in the same set of backcross mice. These genes greatly expanded an existing conserved synteny group between the human 16p12-p13 region and the MMU7. It now consists of eleven loci that span a region of probably more than 10 Mb in human. The gene order derived from this study provided further evidence for chromosomal rearrangements within the conserved synteny. (Abstract shortened by UMI.) ^

Human TNF receptor p75/80 gene structure and promoter characterization

Tumor necrosis factor receptor p75/80 ((TNF-R p75/80) is a 75 kDa type 1 transmembrane protein expressed predominately on cells of hematopoietic lineage. TNF-R p75/80 belongs to the TNF receptor superfamily characterized by cysteine-rich extracellular regions composed of three to six disulfide-linked domains. In the present report, we have characterized, for the first time, the complete gene structure for human TNF-R p75/80 which spans approximately 43 kbp. The gene consists of 10 exons (ranging from 34 bp to 2.5 kbp) and 9 introns (343 bp to 19 kbp). Consensus elements for transcription factors involved in T cell development and activation were noted in the 5$\sp\prime$ flanking region including TCF-1, Ikaros, AP-1, CK-2, IL-6RE, ISRE, GAS, NF-$\kappa$B and SP1, as well as an unusually high GC content and CpG frequency that appears characteristic of some TNF-R family members. The unusual (GATA)$\sb{\rm n}$ and (GAA)(GGA) repeats found within intron 1 may prove useful for further genome analysis within the 1p36 chromosomal locus. The human TNF-R p75/80 gene structure will permit further assessment of its involvement in normal hematopoietic cell development and function, autoimmune disease, and non-random translocations in hematopoietic malignancies. The region 1.8 kb 5$\sp\prime$ of the ATG was able to drive luciferase expression when transfected into cell lines expressing TNF-R p75/80. Further characterization of the 5$\sp\prime$-regulatory region will aid in determining factors and signal transduction pathways involved in regulating TNF-R p75/80 expression. ^

Frederiksenia canicola gen. nov., sp. nov. isolated from dogs and human dog-bite wounds.

Polyphasic analysis was done on 24 strains of Bisgaard taxon 16 from five European countries and mainly isolated from dogs and human dog-bite wounds. The isolates represented a phenotypically and genetically homogenous group within the family Pasteurellaceae. Their phenotypic profile was similar to members of the genus Pasteurella. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry clearly identified taxon 16 and separated it from all other genera of Pasteurellaceae showing a characteristic peak combination. Taxon 16 can be further separated and identified by a RecN protein signature sequence detectable by a specific PCR. In all phylogenetic analyses based on 16S rRNA, rpoB, infB and recN genes, taxon 16 formed a monophyletic branch with intraspecies sequence similarity of at least 99.1, 90.8, 96.8 and 97.2 %, respectively. Taxon 16 showed closest genetic relationship with Bibersteinia trehalosi as to the 16S rRNA gene (95.9 %), the rpoB (89.8 %) and the recN (74.4 %), and with Actinobacillus lignieresii for infB (84.9 %). Predicted genome similarity values based on the recN gene sequences between taxon 16 isolates and the type strains of known genera of Pasteurellaceae were below the genus level. Major whole cell fatty acids for the strain HPA 21(T) are C14:0, C16:0, C18:0 and C16:1 ω7c/C15:0 iso 2OH. Major respiratory quinones are menaquinone-8, ubiquinone-8 and demethylmenaquinone-8. We propose to classify these organisms as a novel genus and species within the family of Pasteurellaceae named Frederiksenia canicola gen. nov., sp. nov. The type strain is HPA 21(T) (= CCUG 62410(T) = DSM 25797(T)).

Comparative methodology to investigate the presence of Escherichia coli K-12 strains in environmental and human stool samples

This study was undertaken to evaluate the specificity and efficiency of different methods to detect Escherichia coli K-12 strains. Another aim was to determine the frequency of E. coli K-12 strains among wild-type E. coli isolates from different sources. The detection of K-12 strains was performed both genotypically by K-12 specific polymerase chain reaction (PCR) and on the basis of phenotypical tests. In addition, the genome structures of E. coli strains were characterized by pulsed-field gel electrophoresis (PFGE). The most specific results could be obtained by the genotypical tests PCR and PFGE as well as by the K-12 specific phage assay. In total, 131 stool and 95 water isolates as well as 14 K-12 derivatives were examined by the different methods. No E. coli K-12 strains were detected among the wild-type isolates.

Heterologous expression from the human D-Loop in organello

We report the expression of a linear reporter construct in isolated human mitochondria. The reporter construct contained the entire human D-Loop with adjacent tRNA (MTT) genes (mt.15956-647), the human ND1 gene with an in frame GFP gene and adjacent endogenous MTT genes and heterologous rat MTT genes. Natural competence of isolated human mitochondria of HepG2 cells was used to import reporter constructs. The import efficiency of various fluorescently labelled PCR-generated import substrates in the range of 250bp up to 3.5kb was assessed by quantitative PCR and evaluated by confocal microscopy. Heterologous expression of the imported construct was confirmed at RNA level by a circular RNA (cRNA)-RT-PCR assay for the expression of tRNAs and by in organello [α-(32)P]-UTP labelling and subsequent hybridisation to reporter-specific sequences for monitoring mRNA expression. Heterologous expression of rat mitochondrial tRNA(Leu(UUR)) (rMT-TL1) was confirmed by co-/post-transcriptional trinucleotide (CCA) addition. Interestingly, the rat-specific MT-TL1 was correctly processed in isolated human mitochondria at the 3' end, but showed an aberrant 5' end processing. Correct 3' end processing of the heterologous expressed mitochondrial rat tRNA(Ser2) (MT-TS2) was detected. These findings demonstrate the feasibility of genetic manipulation of human mitochondria, providing a tool for characterisation of cis-acting elements of the human mitochondrial genome and for the study of human mitochondrial tRNA processing in organello.

Evaluation of genome-wide loci of iron metabolism in hereditary hemochromatosis identifies PCSK7 as a host risk factor of liver cirrhosis

Genome-wide association studies (GWAS) have revealed genetic determinants of iron metabolism, but correlation of these with clinical phenotypes is pending. Homozygosity for HFE C282Y is the predominant genetic risk factor for hereditary hemochromatosis (HH) and may cause liver cirrhosis. However, this genotype has a low penetrance. Thus, detection of yet unknown genetic markers that identify patients at risk of developing severe liver disease is necessary for better prevention. Genetic loci associated with iron metabolism (TF, TMPRSS6, PCSK7, TFR2 and Chr2p14) in recent GWAS and liver fibrosis (PNPLA3) in recent meta-analysis were analyzed for association with either liver cirrhosis or advanced fibrosis in 148 German HFE C282Y homozygotes. Replication of associations was sought in additional 499 Austrian/Swiss and 112 HFE C282Y homozygotes from Sweden. Only variant rs236918 in the PCSK7 gene (proprotein convertase subtilisin/kexin type 7) was associated with cirrhosis or advanced fibrosis (P = 1.02 × 10(-5)) in the German cohort with genotypic odds ratios of 3.56 (95% CI 1.29-9.77) for CG heterozygotes and 5.38 (95% CI 2.39-12.10) for C allele carriers. Association between rs236918 and cirrhosis was confirmed in Austrian/Swiss HFE C282Y homozygotes (P = 0.014; ORallelic = 1.82 (95% CI 1.12-2.95) but not in Swedish patients. Post hoc combined analyses of German/Swiss/Austrian patients with available liver histology (N = 244, P = 0.00014, ORallelic = 2.84) and of males only (N = 431, P = 2.17 × 10(-5), ORallelic = 2.54) were consistent with the premier finding. Association between rs236918 and cirrhosis was not confirmed in alcoholic cirrhotics, suggesting specificity of this genetic risk factor for HH. PCSK7 variant rs236918 is a risk factor for cirrhosis in HH patients homozygous for the HFE C282Y mutation.

A Genome-Wide Association Study Reveals Loci Influencing Height and Other Conformation Traits in Horses

The molecular analysis of genes influencing human height has been notoriously difficult. Genome-wide association studies (GWAS) for height in humans based on tens of thousands to hundreds of thousands of samples so far revealed ∼200 loci for human height explaining only 20% of the heritability. In domestic animals isolated populations with a greatly reduced genetic heterogeneity facilitate a more efficient analysis of complex traits. We performed a genome-wide association study on 1,077 Franches-Montagnes (FM) horses using ∼40,000 SNPs. Our study revealed two QTL for height at withers on chromosomes 3 and 9. The association signal on chromosome 3 is close to the LCORL/NCAPG genes. The association signal on chromosome 9 is close to the ZFAT gene. Both loci have already been shown to influence height in humans. Interestingly, there are very large intergenic regions at the association signals. The two detected QTL together explain ∼18.2% of the heritable variation of height in horses. However, another large fraction of the variance for height in horses results from ECA 1 (11.0%), although the association analysis did not reveal significantly associated SNPs on this chromosome. The QTL region on ECA 3 associated with height at withers was also significantly associated with wither height, conformation of legs, ventral border of mandible, correctness of gaits, and expression of the head. The region on ECA 9 associated with height at withers was also associated with wither height, length of croup and length of back. In addition to these two QTL regions on ECA 3 and ECA 9 we detected another QTL on ECA 6 for correctness of gaits. Our study highlights the value of domestic animal populations for the genetic analysis of complex traits.

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.

Genome profiling of sterol synthesis shows convergent evolution in parasites and guides chemotherapeutic attack

Sterols are an essential class of lipids in eukaryotes, where they serve as structural components of membranes and play important roles as signaling molecules. Sterols are also of high pharmacological significance: cholesterol-lowering drugs are blockbusters in human health, and inhibitors of ergosterol biosynthesis are widely used as antifungals. Inhibitors of ergosterol synthesis are also being developed for Chagas's disease, caused by Trypanosoma cruzi. Here we develop an in silico pipeline to globally evaluate sterol metabolism and perform comparative genomics. We generate a library of hidden Markov model-based profiles for 42 sterol biosynthetic enzymes, which allows expressing the genomic makeup of a given species as a numerical vector. Hierarchical clustering of these vectors functionally groups eukaryote proteomes and reveals convergent evolution, in particular metabolic reduction in obligate endoparasites. We experimentally explore sterol metabolism by testing a set of sterol biosynthesis inhibitors against trypanosomatids, Plasmodium falciparum, Giardia, and mammalian cells, and by quantifying the expression levels of sterol biosynthetic genes during the different life stages of T. cruzi and Trypanosoma brucei. The phenotypic data correlate with genomic makeup for simvastatin, which showed activity against trypanosomatids. Other findings, such as the activity of terbinafine against Giardia, are not in agreement with the genotypic profile.

A Multi-Breed Genome-Wide Association Analysis for Canine Hypothyroidism Identifies a Shared Major Risk Locus on CFA12.

Hypothyroidism is a complex clinical condition found in both humans and dogs, thought to be caused by a combination of genetic and environmental factors. In this study we present a multi-breed analysis of predisposing genetic risk factors for hypothyroidism in dogs using three high-risk breeds-the Gordon Setter, Hovawart and the Rhodesian Ridgeback. Using a genome-wide association approach and meta-analysis, we identified a major hypothyroidism risk locus shared by these breeds on chromosome 12 (p = 2.1x10-11). Further characterisation of the candidate region revealed a shared ~167 kb risk haplotype (4,915,018-5,081,823 bp), tagged by two SNPs in almost complete linkage disequilibrium. This breed-shared risk haplotype includes three genes (LHFPL5, SRPK1 and SLC26A8) and does not extend to the dog leukocyte antigen (DLA) class II gene cluster located in the vicinity. These three genes have not been identified as candidate genes for hypothyroid disease previously, but have functions that could potentially contribute to the development of the disease. Our results implicate the potential involvement of novel genes and pathways for the development of canine hypothyroidism, raising new possibilities for screening, breeding programmes and treatments in dogs. This study may also contribute to our understanding of the genetic etiology of human hypothyroid disease, which is one of the most common endocrine disorders in humans.

Models of hybridization during range expansions and their application to recent human evolution

Several lines of genetic, archeological and paleontological evidence suggest that anatomically modern humans (Homo sapiens) colonized the world in the last 60,000 years by a series of migrations originating from Africa (e.g. Liu et al., 2006; Handley et al., 2007; Prugnolle, Manica, and Balloux, 2005; Ramachandran et al. 2005; Li et al. 2008; Deshpande et al. 2009; Mellars, 2006a, b; Lahr and Foley, 1998; Gravel et al., 2011; Rasmussen et al., 2011). With the progress of ancient DNA analysis, it has been shown that archaic humans hybridized with modern humans outside Africa. Recent direct analyses of fossil nuclear DNA have revealed that 1–4 percent of the genome of Eurasian has been likely introgressed by Neanderthal genes (Green et al., 2010; Reich et al., 2010; Vernot and Akey, 2014; Sankararaman et al., 2014; Prufer et al., 2014; Wall et al., 2013), with Papua New Guineans and Australians showing even larger levels of admixture with Denisovans (Reich et al., 2010; Skoglund and Jakobsson, 2011; Reich et al., 2011; Rasmussen et al., 2011). It thus appears that the past history of our species has been more complex than previously anticipated (Alves et al., 2012), and that modern humans hybridized several times with local hominins during their expansion out of Africa, but the exact mode, time and location of these hybridizations remain to be clarifi ed (Ibid.; Wall et al., 2013). In this context, we review here a general model of admixture during range expansion, which lead to some predictions about expected patterns of introgression that are relevant to modern human evolution.

A RAB3GAP1 SINE Insertion in Alaskan Huskies with Polyneuropathy, Ocular Abnormalities and Neuronal Vacuolation (POANV) Resembling Human Warburg Micro Syndrome 1 (WARBM1)

We observed a hereditary phenotype in Alaskan Huskies, which was characterized by polyneuropathy with ocular abnormalities and neuronal vacuolation (POANV). The affected dogs developed a progressive severe ataxia, which led to euthanasia between 8 and 16 months of age. The pedigrees were consistent with a monogenic autosomal recessive inheritance. We localized the causative genetic defect to a 4 Mb interval on chromosome 19 by a combined linkage and homozygosity mapping approach. Whole genome sequencing of one affected dog, an obligate carrier and an unrelated control revealed a 218 bp SINE insertion into exon 7 of the RAB3GAP1 gene. The SINE insertion was perfectly associated with the disease phenotype in a cohort of 43 Alaskan Huskies and it was absent from 541 control dogs of diverse other breeds. The SINE insertion induced aberrant splicing and led to a transcript with a greatly altered exon 7. RAB3GAP1 loss-of-function variants in humans cause Warburg Micro Syndrome 1 (WARBM1), which is characterized by additional developmental defects compared to canine POANV, whereas Rab3gap1 deficient mice have a much milder phenotype than either humans or dogs. Thus the RAB3GAP1 mutant Alaskan Huskies provide an interesting intermediate phenotype that may help to better understand the function of RAB3GAP1 in development. Furthermore, the identification of the presumed causative genetic variant will enable genetic testing to avoid the non-intentional breeding of affected dogs.

Detection of convergent genome-wide signals of adaptation to tropical forests in humans

Tropical forests are believed to be very harsh environments for human life. It is unclear whether human beings would have ever subsisted in those environments without external resources. It is therefore possible that humans have developed recent biological adaptations in response to specific selective pressures to cope with this challenge. To understand such biological adaptations we analyzed genome-wide SNP data under a Bayesian statistics framework, looking for outlier markers with an overly large extent of differentiation between populations living in a tropical forest, as compared to genetically related populations living outside the forest in Africa and the Americas. The most significant positive selection signals were found in genes related to lipid metabolism, the immune system, body development, and RNA Polymerase III transcription initiation. The results are discussed in the light of putative tropical forest selective pressures, namely food scarcity, high prevalence of pathogens, difficulty to move, and inefficient thermoregulation. Agreement between our results and previous studies on the pygmy phenotype, a putative prototype of forest adaptation, were found, suggesting that a few genetic regions previously described as associated with short stature may be evolving under similar positive selection in Africa and the Americas. In general, convergent evolution was less pervasive than local adaptation in one single continent, suggesting that Africans and Amerindians may have followed different routes to adapt to similar environmental selective pressures.

Inference of evolutionary forces acting on human biological pathways

Because natural selection is likely to act on multiple genes underlying a given phenotypic trait, we study here the potential effect of ongoing and past selection on the genetic diversity of human biological pathways. We first show that genes included in gene sets are generally under stronger selective constraints than other genes and that their evolutionary response is correlated. We then introduce a new procedure to detect selection at the pathway level based on a decomposition of the classical McDonald–Kreitman test extended to multiple genes. This new test, called 2DNS, detects outlier gene sets and takes into account past demographic effects and evolutionary constraints specific to gene sets. Selective forces acting on gene sets can be easily identified by a mere visual inspection of the position of the gene sets relative to their two-dimensional null distribution. We thus find several outlier gene sets that show signals of positive, balancing, or purifying selection but also others showing an ancient relaxation of selective constraints. The principle of the 2DNS test can also be applied to other genomic contrasts. For instance, the comparison of patterns of polymorphisms private to African and non-African populations reveals that most pathways show a higher proportion of nonsynonymous mutations in non-Africans than in Africans, potentially due to different demographic histories and selective pressures.

New insights into the performance of human whole-exome capture platforms.

Whole exome sequencing (WES) is increasingly used in research and diagnostics. WES users expect coverage of the entire coding region of known genes as well as sufficient read depth for the covered regions. It is, however, unknown which recent WES platform is most suitable to meet these expectations. We present insights into the performance of the most recent standard exome enrichment platforms from Agilent, NimbleGen and Illumina applied to six different DNA samples by two sequencing vendors per platform. Our results suggest that both Agilent and NimbleGen overall perform better than Illumina and that the high enrichment performance of Agilent is stable among samples and between vendors, whereas NimbleGen is only able to achieve vendor- and sample-specific best exome coverage. Moreover, the recent Agilent platform overall captures more coding exons with sufficient read depth than NimbleGen and Illumina. Due to considerable gaps in effective exome coverage, however, the three platforms cannot capture all known coding exons alone or in combination, requiring improvement. Our data emphasize the importance of evaluation of updated platform versions and suggest that enrichment-free whole genome sequencing can overcome the limitations of WES in sufficiently covering coding exons, especially GC-rich regions, and in characterizing structural variants.