Here we report a high-quality draft genome sequence of the domestic dog (Canis familiaris), together with a dense map of single nucleotide polymorphisms (SNPs) across breeds. The dog is of particular interest because it provides important evolutionary information and because existing breeds show great phenotypic diversity for morphological, physiological and behavioural traits. We use sequence comparison with the primate and rodent lineages to shed light on the structure and evolution of genomes and genes. Notably, the majority of the most highly conserved non-coding sequences in mammalian genomes are clustered near a small subset of genes with important roles in development. Analysis of SNPs reveals long-range haplotypes across the entire dog genome, and defines the nature of genetic diversity within and across breeds. The current SNP map now makes it possible for genome-wide association studies to identify genes responsible for diseases and traits, with important consequences for human and companion animal health.

BACKGROUND The narrow host range of Mycobacterium leprae and the fact that it is refractory to growth in culture has limited research on and the biologic understanding of leprosy. Host genetic factors are thought to influence susceptibility to infection as well as disease progression. METHODS We performed a two-stage genomewide association study by genotyping 706 patients and 1225 controls using the Human610-Quad BeadChip (Illumina). We then tested three independent replication sets for an association between the presence of leprosy and 93 single-nucleotide polymorphisms (SNPs) that were most strongly associated with the disease in the genomewide association study. Together, these replication sets comprised 3254 patients and 5955 controls. We also carried out tests of heterogeneity of the associations (or lack thereof) between these 93 SNPs and disease, stratified according to clinical subtype (multibacillary vs. paucibacillary). RESULTS We observed a significant association (P<1.00×10 −10) between SNPs in the genes CCDC122, C13orf31, NOD2, TNFSF15, HLA-DR, and RIPK2 and a trend toward an association (P = 5.10×10 −5) with a SNP in LRRK2. The associations between the SNPs in C13orf31, LRRK2, NOD2, and RIPK2 and multibacillary leprosy were stronger than the associations between these SNPs and paucibacillary leprosy. CONCLUSIONS Variants of genes in the NOD2-mediated signaling pathway (which regulates the innate immune response) are associated with susceptibility to infection with M. leprae.

The hypoxic environment imposes severe selective pressure on species living at high altitude. To understand the genetic bases of adaptation to high altitude in dogs, we performed whole-genome sequencing of 60 dogs including five breeds living at continuous altitudes along the Tibetan Plateau from 800 to 5100 m as well as one European breed. More than 1503 sequencing coverage for each breed provides us with a comprehensive assessment of the genetic polymorphisms of the dogs, including Tibetan Mastiffs. Comparison of the breeds from different altitudes reveals strong signals of population differentiation at the locus of hypoxia-related genes including endothelial Per-Arnt-Sim (PAS) domain protein 1 (EPAS1) and beta hemoglobin cluster. Notably, four novel nonsynonymous mutations specific to high-altitude dogs are identified at EPAS1, one of which occurred at a quite conserved site in the PAS domain. The association testing between EPAS1 genotypes and blood-related phenotypes on additional high-altitude dogs reveals that the homozygous mutation is associated with decreased blood flow resistance, which may help to improve hemorheologic fitness. Interestingly, EPAS1 was also identified as a selective target in Tibetan highlanders, though no amino acid changes were found. Thus, our results not only indicate parallel evolution of humans and dogs in adaptation to high-altitude hypoxia, but also provide a new opportunity to study the role of EPAS1 in the adaptive processes.

Four pairs of XX and XY gonads from Nile tilapia were sequenced at four developmental stages, 5, 30, 90, and 180 days after hatching (dah) using Illumina HiseqTM technology. This produced 28 Gb sequences, which were mapped to 21,334 genes. Of these, 259 genes were found to be specifically expressed in XY gonads, and 69 were found to be specific to XX gonads. Totally, 187 XX- and 1,358 XY-enhanced genes were identified, and 2,978 genes were found to be co-expressed in XX and XY gonads. Almost all steroidogenic enzymes, including cyp19a1a, were up-regulated in XX gonads at 5 dah; but in XY gonads these enzymes, including cyp11b2, were significantly up-regulated at 90 dah, indicating that, at a time critical to sex determination, the XX fish produced estrogen and the XY fish did not produce androgens. The most pronounced expression of steroidogenic enzyme genes was observed at 30 and 90 dah for XX and XY gonads, corresponding to the initiation of germ cell meiosis in the female and male gonads, respectively. Both estrogen and androgen receptors were found to be expressed in XX gonads, but only estrogen receptors were expressed in XY gonads at 5 dah. This could explain why exogenous steroid treatment induced XX and XY sex reversal. The XX-enhanced expression of cyp19a1a and cyp19a1b at all stages suggests an important role for estrogen in female sex determination and maintenance of phenotypic sex. This work is the largest collection of gonadal transcriptome data in tilapia and lays the foundation for future studies into the molecular mechanisms of sex determination and maintenance of phenotypic sex in non-model teleosts.

The objective of this study was to evaluate changes in ruminal microorganisms and fermentation parameters due to dietary supplementation of soybean and linseed oil alone or in combination. Four dietary treatments were tested in a Latin square designed experiment using four primiparous rumen-cannulated dairy cows. Treatments were control (C, 60 : 40 forage to concentrate) or C with 4% soybean oil (S), 4% linseed oil (L) or 2% soybean oil plus 2% linseed oil (SL) in a 4 × 4 Latin square with four periods of 21 days. Forage and concentrate mixtures were fed at 0800 and 2000 h daily. Ruminal fluid was collected every 2 h over a 12-h period on day 19 of each experimental period and pH was measured immediately. Samples were prepared for analyses of concentrations of volatile fatty acids (VFA) by GLC and ammonia. Counts of total and individual bacterial groups (cellulolytic, proteolytic, amylolytic bacteria and total viable bacteria) were performed using the roll-tube technique, and protozoa counts were measured via microscopy in ruminal fluid collected at 0, 4 and 8 h after the morning feeding. Content of ruminal digesta was obtained via the rumen cannula before the morning feeding and used immediately for DNA extraction and quantity of specific bacterial species was obtained using real- time PCR. Ruminal pH did not differ but total VFA (110v. 105 mmol/l) were lower (P< 0.05) with oil supplementation compared with C. Concentration of ruminal NH3-N (4.4v. 5.6 mmol/l) was greater (P< 0.05) due to oil compared with C. Compared with C, oil supplementation resulted in lower (P< 0.05) cellulolytic bacteria (3.25 × 108v. 4.66 × 108colony-forming units (CFU)/ml) and protozoa (9.04 × 104v. 12.92 × 104cell/ml) colony counts. Proteolytic bacteria (7.01 × 108v. 6.08 × 108CFU/ml) counts, however, were greater in response to oil compared with C (P< 0.05). Among oil treatments, the amount ofButyrivibrio fibrisolvens,Fibrobacter succinogenesandRuminococcus flavefaciensin ruminal fluid was substantially lower (P< 0.05) when L was included. Compared to C, the amount ofRuminococcus albusdecreased by an average of 40% regardless of oil level or type. Overall, the results indicate that some ruminal microorganisms, except proteolytic bacteria, are highly susceptible to dietary unsaturated fatty acids supplementation, particularly when linolenic acid rich oils were fed. Dietary oil effects on ruminal fermentation parameters seemed associated with the profile of ruminal microorganisms.