Analysis of the pig genome for the identification of genomic regions affecting production traits

The aim of this work was to identify markers associated with production traits in the pig genome using different approaches. We focused the attention on Italian Large White pig breed using Genome Wide Association Studies (GWAS) and applying a selective genotyping approach to increase the power of the analyses. Furthermore, we searched the pig genome using Next Generation Sequencing (NSG) Ion Torrent Technology to combine selective genotyping approach and deep sequencing for SNP discovery. Other two studies were carried on with a different approach. Allele frequency changes for SNPs affecting candidate genes and at Genome Wide level were analysed to identify selection signatures driven by selection program during the last 20 years. This approach confirmed that a great number of markers may affect production traits and that they are captured by the classical selection programs. GWAS revealed 123 significant or suggestively significant SNP associated with Back Fat Thickenss and 229 associated with Average Daily Gain. 16 Copy Number Variant Regions resulted more frequent in lean or fat pigs and showed that different copies of those region could have a limited impact on fat. These often appear to be involved in food intake and behavior, beside affecting genes involved in metabolic pathways and their expression. By combining NGS sequencing with selective genotyping approach, new variants where discovered and at least 54 are worth to be analysed in association studies. The study of groups of pigs undergone to stringent selection showed that allele frequency of some loci can drastically change if they are close to traits that are interesting for selection schemes. These approaches could be, in future, integrated in genomic selection plans.

Bipolar disorder susceptibility region on chromosome 3q29 not confirmed in a case-control association study

We identified a bipolar disorder (BPD) susceptibility region on chromosome 3q29 in a genome-wide linkage scan (Bailer et al. 2002 (Biol Psychiatry 52: 40), NPL-score 4.09) and follow-up linkage analysis (Schosser et al. 2004 (J Psychiatr Res 38(3): 357), NPL-scores >3 with five markers). These findings were supported by further fine-mapping of this region (Schosser et al. 2007 (Eur Neuropsychopharmacol 17(6-7): 501)), finding NPL-scores >3.9 with SNPs (single nucleotide polymorphisms) spanning a region of 3.46 Mbp in BPD families. Since genetic association studies are more powerful than linkage studies for detecting susceptibility genes of small effect size, we aimed to replicate these findings in an independent case-control sample collected in London (UK) and Vienna (Austria).

An unusual splice defect in the mitofusin 2 gene (MFN2) is associated with degenerative axonopathy in Tyrolean Grey cattle

Tyrolean Grey cattle represent a local breed with a population size of approximately 5000 registered cows. In 2003, a previously unknown neurological disorder was recognized in Tyrolean Grey cattle. The clinical signs of the disorder are similar to those of bovine progressive degenerative myeloencephalopathy (weaver syndrome) in Brown Swiss cattle but occur much earlier in life. The neuropathological investigation of an affected calf showed axonal degeneration in the central nervous system (CNS) and femoral nerve. The pedigrees of the affected calves suggested a monogenic autosomal recessive inheritance. We localized the responsible mutation to a 1.9 Mb interval on chromosome 16 by genome-wide association and haplotype mapping. The MFN2 gene located in this interval encodes mitofusin 2, a mitochondrial membrane protein. A heritable human axonal neuropathy, Charcot-Marie-Tooth disease-2A2 (CMT2A2), is caused by MFN2 mutations. Therefore, we considered MFN2 a positional and functional candidate gene and performed mutation analysis in affected and control Tyrolean Grey cattle. We did not find any non-synonymous variants. However, we identified a perfectly associated silent SNP in the coding region of exon 20 of the MFN2 gene. This SNP is located within a putative exonic splice enhancer (ESE) and the variant allele leads to partial retention of the entire intron 19 and a premature stop codon in the aberrant MFN2 transcript. Thus we have identified a highly unusual splicing defect, where an exonic single base exchange leads to the retention of the preceding intron. This splicing defect represents a potential explanation for the observed degenerative axonopathy. Marker assisted selection can now be used to eliminate degenerative axonopathy from Tyrolean Grey cattle.

Microphthalmia in Texel sheep is associated with a missense mutation in the paired-like homeodomain 3 (PITX3) gene

Microphthalmia in sheep is an autosomal recessive inherited congenital anomaly found within the Texel breed. It is characterized by extremely small or absent eyes and affected lambs are absolutely blind. For the first time, we use a genome-wide ovine SNP array for positional cloning of a Mendelian trait in sheep. Genotyping 23 cases and 23 controls using Illumina's OvineSNP50 BeadChip allowed us to localize the causative mutation for microphthalmia to a 2.4 Mb interval on sheep chromosome 22 by association and homozygosity mapping. The PITX3 gene is located within this interval and encodes a homeodomain-containing transcription factor involved in vertebrate lens formation. An abnormal development of the lens vesicle was shown to be the primary event in ovine microphthalmia. Therefore, we considered PITX3 a positional and functional candidate gene. An ovine BAC clone was sequenced, and after full-length cDNA cloning the PITX3 gene was annotated. Here we show that the ovine microphthalmia phenotype is perfectly associated with a missense mutation (c.338G>C, p.R113P) in the evolutionary conserved homeodomain of PITX3. Selection against this candidate causative mutation can now be used to eliminate microphthalmia from Texel sheep in production systems. Furthermore, the identification of a naturally occurring PITX3 mutation offers the opportunity to use the Texel as a genetically characterized large animal model for human microphthalmia.

DEVELOPMENT AND APPLICATIONS OF GENOMIC RESOURCES FOR TWO NORTH AMERICAN HARDWOOD TAXA

Hardwoods comprise about half of the biomass of forestlands in North America and present many uses including economic, ecological and aesthetic functions. Forest trees rely on the genetic variation within tree populations to overcome the many biotic, abiotic, anthropogenic factors which are further worsened by climate change, that threaten their continued survival and functionality. To harness these inherent genetic variations of tree populations, informed knowledge of the genomic resources and techniques, which are currently lacking or very limited, are imperative for forest managers. The current study therefore aimed to develop genomic microsatellite markers for the leguminous tree species, honey locust, Gleditsia triacanthos L. and test their applicability in assessing genetic variation, estimation of gene flow patterns and identification of a full-sib mapping population. We also aimed to test the usefulness of already developed nuclear and gene-based microsatellite markers in delineation of species and taxonomic relationships between four of the taxonomically difficult Section Lobatae species (Quercus coccinea, Q. ellipsoidalis, Q. rubra and Q. velutina. We recorded 100% amplification of G. triacanthos genomic microsatellites developed using Illumina sequencing techniques in a panel of seven unrelated individuals with 14 of these showing high polymorphism and reproducibility. When characterized in 36 natural population samples, we recorded 20 alleles per locus with no indication for null alleles at 13 of the 14 microsatellites. This is the first report of genomic microsatellites for this species. Honey locust trees occur in fragmented populations of abandoned farmlands and pastures and is described as essentially dioecious. Pollen dispersal if the main source of gene flow within and between populations with the ability to offset the effects of random genetic drift. Factors known to influence gene include fragmentation and degree of isolation, which make the patterns gene flow in fragmented populations of honey locust a necessity for their sustainable management. In this follow-up study, we used a subset of nine of the 14 developed gSSRs to estimate gene flow and identify a full-sib mapping population in two isolated fragments of honey locust. Our analyses indicated that the majority of the seedlings (65-100% - at both strict and relaxed assignment thresholds) were sired by pollen from outside the two fragment populations. Only one selfing event was recorded confirming the functional dioeciousness of honey locust and that the seed parents are almost completely outcrossed. From the Butternut Valley, TN population, pollen donor genotypes were reconstructed and used in paternity assignment analyses to identify a relatively large full-sib family comprised of 149 individuals, proving the usefulness of isolated forest fragments in identification of full-sib families. In the Ames Plantation stand, contemporary pollen dispersal followed a fat-tailed exponential-power distribution, an indication of effective gene flow. Our estimate of δ was 4,282.28 m, suggesting that insect pollinators of honey locust disperse pollen over very long distances. The high proportion of pollen influx into our sampled population implies that our fragment population forms part of a large effectively reproducing population. The high tendency of oak species to hybridize while still maintaining their species identity make it difficult to resolve their taxonomic relationships. Oaks of the section Lobatae are famous in this regard and remain unresolved at both morphological and genetic markers. We applied 28 microsatellite markers including outlier loci with potential roles in reproductive isolation and adaptive divergence between species to natural populations of four known interfertile red oaks, Q. coccinea, Q. ellpsoidalis, Q. rubra and Q. velutina. To better resolve the taxonomic relationships in this difficult clade, we assigned individual samples to species, identified hybrids and introgressive forms and reconstructed phylogenetic relationships among the four species after exclusion of genetically intermediate individuals. Genetic assignment analyses identified four distinct species clusters, with Q. rubra most differentiated from the three other species, but also with a comparatively large number of misclassified individuals (7.14%), hybrids (7.14%) and introgressive forms (18.83%) between Q. ellipsoidalis and Q. velutina. After the exclusion of genetically intermediate individuals, Q. ellipsoidalis grouped as sister species to the largely parapatric Q. coccinea with high bootstrap support (91 %). Genetically intermediate forms in a mixed species stand were located proximate to both potential parental species, which supports recent hybridization of Q. velutina with both Q. ellipsoidalis and Q. rubra. Analyses of genome-wide patterns of interspecific differentiation can provide a better understanding of speciation processes and taxonomic relationships in this taxonomically difficult group of red oak species.

Accumulating mutations in series of haplotypes at the KIT and MITF loci are major determinants of white markings in Franches-Montagnes horses.

Coat color and pattern variations in domestic animals are frequently inherited as simple monogenic traits, but a number are known to have a complex genetic basis. While the analysis of complex trait data remains a challenge in all species, we can use the reduced haplotypic diversity in domestic animal populations to gain insight into the genomic interactions underlying complex phenotypes. White face and leg markings are examples of complex traits in horses where little is known of the underlying genetics. In this study, Franches-Montagnes (FM) horses were scored for the occurrence of white facial and leg markings using a standardized scoring system. A genome-wide association study (GWAS) was performed for several white patterning traits in 1,077 FM horses. Seven quantitative trait loci (QTL) affecting the white marking score with p-values p≤10(-4) were identified. Three loci, MC1R and the known white spotting genes, KIT and MITF, were identified as the major loci underlying the extent of white patterning in this breed. Together, the seven loci explain 54% of the genetic variance in total white marking score, while MITF and KIT alone account for 26%. Although MITF and KIT are the major loci controlling white patterning, their influence varies according to the basic coat color of the horse and the specific body location of the white patterning. Fine mapping across the MITF and KIT loci was used to characterize haplotypes present. Phylogenetic relationships among haplotypes were calculated to assess their selective and evolutionary influences on the extent of white patterning. This novel approach shows that KIT and MITF act in an additive manner and that accumulating mutations at these loci progressively increase the extent of white markings.

A frameshift mutation in the cubilin gene (CUBN) in Border Collies with Imerslund-Gräsbeck syndrome (selective cobalamin malabsorption).

Imerslund-Gräsbeck syndrome (IGS) or selective cobalamin malabsorption has been described in humans and dogs. IGS occurs in Border Collies and is inherited as a monogenic autosomal recessive trait in this breed. Using 7 IGS cases and 7 non-affected controls we mapped the causative mutation by genome-wide association and homozygosity mapping to a 3.53 Mb interval on chromosome 2. We re-sequenced the genome of one affected dog at ∼10× coverage and detected 17 non-synonymous variants in the critical interval. Two of these non-synonymous variants were in the cubilin gene (CUBN), which is known to play an essential role in cobalamin uptake from the ileum. We tested these two CUBN variants for association with IGS in larger cohorts of dogs and found that only one of them was perfectly associated with the phenotype. This variant, a single base pair deletion (c.8392delC), is predicted to cause a frameshift and premature stop codon in the CUBN gene. The resulting mutant open reading frame is 821 codons shorter than the wildtype open reading frame (p.Q2798Rfs*3). Interestingly, we observed an additional nonsense mutation in the MRC1 gene encoding the mannose receptor, C type 1, which was in perfect linkage disequilibrium with the CUBN frameshift mutation. Based on our genetic data and the known role of CUBN for cobalamin uptake we conclude that the identified CUBN frameshift mutation is most likely causative for IGS in Border Collies.

A COL11A2 mutation in Labrador retrievers with mild disproportionate dwarfism.

We describe a mild form of disproportionate dwarfism in Labrador Retrievers, which is not associated with any obvious health problems such as secondary arthrosis. We designate this phenotype as skeletal dysplasia 2 (SD2). It is inherited as a monogenic autosomal recessive trait with incomplete penetrance primarily in working lines of the Labrador Retriever breed. Using 23 cases and 37 controls we mapped the causative mutation by genome-wide association and homozygosity mapping to a 4.44 Mb interval on chromosome 12. We re-sequenced the genome of one affected dog at 30x coverage and detected 92 non-synonymous variants in the critical interval. Only two of these variants, located in the lymphotoxin A (LTA) and collagen alpha-2(XI) chain gene (COL11A2), respectively, were perfectly associated with the trait. Previously described COL11A2 variants in humans or mice lead to skeletal dysplasias and/or deafness. The dog variant associated with disproportionate dwarfism, COL11A2:c.143G>C or p.R48P, probably has only a minor effect on collagen XI function, which might explain the comparatively mild phenotype seen in our study. The identification of this candidate causative mutation thus widens the known phenotypic spectrum of COL11A2 mutations. We speculate that non-pathogenic COL11A2 variants might even contribute to the heritable variation in height.

A genome-to-genome analysis of associations between human genetic variation, HIV-1 sequence diversity, and viral control

HIV-1 sequence diversity is affected by selection pressures arising from host genomic factors. Using paired human and viral data from 1071 individuals, we ran >3000 genome-wide scans, testing for associations between host DNA polymorphisms, HIV-1 sequence variation and plasma viral load (VL), while considering human and viral population structure. We observed significant human SNP associations to a total of 48 HIV-1 amino acid variants (p<2.4 × 10−12). All associated SNPs mapped to the HLA class I region. Clinical relevance of host and pathogen variation was assessed using VL results. We identified two critical advantages to the use of viral variation for identifying host factors: (1) association signals are much stronger for HIV-1 sequence variants than VL, reflecting the ‘intermediate phenotype’ nature of viral variation; (2) association testing can be run without any clinical data. The proposed genome-to-genome approach highlights sites of genomic conflict and is a strategy generally applicable to studies of host–pathogen interaction.

Mouse models and quantitative trait loci of susceptibility to bleomycin- and radiation-induced pulmonary fibrosis

Radiotherapy involving the thoracic cavity and chemotherapy with the drug bleomycin are both dose limited by the development of pulmonary fibrosis. From evidence that there is variation in the population in susceptibility to pulmonary fibrosis, and animal data, it was hypothesized that individual variation in susceptibility to bleomycin-induced, or radiation-induced, pulmonary fibrosis is, in part, genetically controlled. In this thesis a three generation mouse genetic model of C57BL/6J (fibrosis prone) and C3Hf/Kam (fibrosis resistant) mouse strains and F1 and F2 (F1 intercross) progeny derived from the parental strains was developed to investigate the genetic basis of susceptibility to fibrosis. In the bleomycin studies the mice received 100 mg/kg (125 for females) of bleomycin, via mini osmotic pump. The animals were sacrificed at eight weeks following treatment or when their breathing rate indicated respiratory distress. In the radiation studies the mice were given a single dose of 14 or 16 Gy (Co$\sp{60})$ to the whole thorax and were sacrificed when moribund. The phenotype was defined as the percent of fibrosis area in the left lung as quantified with image analysis of histological sections. Quantitative trait loci (QTL) mapping was used to identify the chromosomal location of genes which contribute to susceptibility to bleomycin-induced pulmonary fibrosis in C57BL/6J mice compared to C3Hf/Kam mice and to determine if the QTL's which influence susceptibility to bleomycin-induced lung fibrosis in these progenitor strains could be implicated in susceptibility to radiation-induced lung fibrosis. For bleomycin, a genome wide scan revealed QTL's on chromosome 17, at the MHC, (LOD = 11.7 for males and 7.2 for females) accounting for approximately 21% of the phenotypic variance, and on chromosome 11 (LOD = 4.9), in male mice only, adding 8% of phenotypic variance. The bleomycin QTL on chromosome 17 was also implicated for susceptibility to radiation-induced fibrosis (LOD = 5.0) and contributes 7% of the phenotypic variance in the radiation study. In conclusion, susceptibility to both bleomycin-induced and radiation-induced pulmonary fibrosis are heritable traits, and are influenced by a genetic factor which maps to a genomic region containing the MHC. ^

DQB1 locus alone explains most of the risk and protection in narcolepsy with cataplexy in Europe

STUDY OBJECTIVE Prior research has identified five common genetic variants associated with narcolepsy with cataplexy in Caucasian patients. To replicate and/or extend these findings, we have tested HLA-DQB1, the previously identified 5 variants, and 10 other potential variants in a large European sample of narcolepsy with cataplexy subjects. DESIGN Retrospective case-control study. SETTING A recent study showed that over 76% of significant genome-wide association variants lie within DNase I hypersensitive sites (DHSs). From our previous GWAS, we identified 30 single nucleotide polymorphisms (SNPs) with P < 10(-4) mapping to DHSs. Ten SNPs tagging these sites, HLADQB1, and all previously reported SNPs significantly associated with narcolepsy were tested for replication. PATIENTS AND PARTICIPANTS For GWAS, 1,261 narcolepsy patients and 1,422 HLA-DQB1*06:02-matched controls were included. For HLA study, 1,218 patients and 3,541 controls were included. MEASUREMENTS AND RESULTS None of the top variants within DHSs were replicated. Out of the five previously reported SNPs, only rs2858884 within the HLA region (P < 2x10(-9)) and rs1154155 within the TRA locus (P < 2x10(-8)) replicated. DQB1 typing confirmed that DQB1*06:02 confers an extraordinary risk (odds ratio 251). Four protective alleles (DQB1*06:03, odds ratio 0.17, DQB1*05:01, odds ratio 0.56, DQB1*06:09 odds ratio 0.21, DQB1*02 odds ratio 0.76) were also identified. CONCLUSION An overwhelming portion of genetic risk for narcolepsy with cataplexy is found at DQB1 locus. Since DQB1*06:02 positive subjects are at 251-fold increase in risk for narcolepsy, and all recent cases of narcolepsy after H1N1 vaccination are positive for this allele, DQB1 genotyping may be relevant to public health policy.

Comparative genetic analyses point to HCP5 as susceptibility locus for HCV-associated hepatocellular carcinoma

BACKGROUND & AIMS: Recently, genetic variations in MICA (lead single nucleotide polymorphism [SNP] rs2596542) were identified by a genome-wide association study (GWAS) to be associated with hepatitis C virus (HCV)-related hepatocellular carcinoma (HCC) in Japanese patients. In the present study, we sought to determine whether this SNP is predictive of HCC development in the Caucasian population as well. METHODS: An extended region around rs2596542 was genotyped in 1924 HCV-infected patients from the Swiss Hepatitis C Cohort Study (SCCS). Pair-wise correlation between key SNPs was calculated both in the Japanese and European populations (HapMap3: CEU and JPT). RESULTS: To our surprise, the minor allele A of rs2596542 in proximity of MICA appeared to have a protective impact on HCC development in Caucasians, which represents an inverse association as compared to the one observed in the Japanese population. Detailed fine-mapping analyses revealed a new SNP in HCP5 (rs2244546) upstream of MICA as strong predictor of HCV-related HCC in the SCCS (univariable p=0.027; multivariable p=0.0002, odds ratio=3.96, 95% confidence interval=1.90-8.27). This newly identified SNP had a similarly directed effect on HCC in both Caucasian and Japanese populations, suggesting that rs2244546 may better tag a putative true variant than the originally identified SNPs. CONCLUSIONS: Our data confirms the MICA/HCP5 region as susceptibility locus for HCV-related HCC and identifies rs2244546 in HCP5 as a novel tagging SNP. In addition, our data exemplify the need for conducting meta-analyses of cohorts of different ethnicities in order to fine map GWAS signals.

Congenital Hepatic Fibrosis in the Franches-Montagnes Horse Is Associated with the Polycystic Kidney and Hepatic Disease 1 (PKHD1) Gene

Congenital hepatic fibrosis has been described as a lethal disease with monogenic autosomal recessive inheritance in the Swiss Franches-Montagnes horse breed. We performed a genome-wide association study with 5 cases and 12 controls and detected an association on chromosome 20. Subsequent homozygosity mapping defined a critical interval of 952 kb harboring 10 annotated genes and loci including the polycystic kidney and hepatic disease 1 (autosomal recessive) gene (PKHD1). PKHD1 represents an excellent functional candidate as variants in this gene were identified in human patients with autosomal recessive polycystic kidney and hepatic disease (ARPKD) as well as several mouse and rat mutants. Whereas most pathogenic PKHD1 variants lead to polycystic defects in kidney and liver, a small subset of the human ARPKD patients have only liver symptoms, similar to our horses with congenital hepatic fibrosis. The PKHD1 gene is one of the largest genes in the genome with multiple alternative transcripts that have not yet been fully characterized. We sequenced the genomes of an affected foal and 46 control horses to establish a comprehensive list of variants in the critical interval. We identified two missense variants in the PKHD1 gene which were strongly, but not perfectly associated with congenital hepatic fibrosis. We speculate that reduced penetrance and/or potential epistatic interactions with hypothetical modifier genes may explain the imperfect association of the detected PKHD1 variants. Our data thus indicate that horses with congenital hepatic fibrosis represent an interesting large animal model for the liver-restricted subtype of human ARPKD.

Deletion in the EVC2 gene causes chondrodysplastic dwarfism in Tyrolean Grey cattle.

During the summer of 2013 seven Italian Tyrolean Grey calves were born with abnormally short limbs. Detailed clinical and pathological examination revealed similarities to chondrodysplastic dwarfism. Pedigree analysis showed a common founder, assuming autosomal monogenic recessive transmission of the defective allele. A positional cloning approach combining genome wide association and homozygosity mapping identified a single 1.6 Mb genomic region on BTA 6 that was associated with the disease. Whole genome re-sequencing of an affected calf revealed a single candidate causal mutation in the Ellis van Creveld syndrome 2 (EVC2) gene. This gene is known to be associated with chondrodysplastic dwarfism in Japanese Brown cattle, and dwarfism, abnormal nails and teeth, and dysostosis in humans with Ellis-van Creveld syndrome. Sanger sequencing confirmed the presence of a 2 bp deletion in exon 19 (c.2993_2994ACdel) that led to a premature stop codon in the coding sequence of bovine EVC2, and was concordant with the recessive pattern of inheritance in affected and carrier animals. This loss of function mutation confirms the important role of EVC2 in bone development. Genetic testing can now be used to eliminate this form of chondrodysplastic dwarfism from Tyrolean Grey cattle.

Looking the cow in the eye: deletion in the NID1 gene is associated with recessive inherited cataract in Romagnola cattle

Cataract is a known condition leading to opacification of the eye lens causing partial or total blindness. Mutations are known to cause autosomal dominant or recessive inherited forms of cataracts in humans, mice, rats, guinea pigs and dogs. The use of large-sized animal models instead of those using mice for the study of this condition has been discussed due to the small size of rodent lenses. Four juvenile-onset cases of bilateral incomplete immature nuclear cataract were recently observed in Romagnola cattle. Pedigree analysis suggested a monogenic autosomal recessive inheritance. In addition to the cataract, one of the cases displayed abnormal head movements. Genome-wide association and homozygosity mapping and subsequent whole genome sequencing of a single case identified two perfectly associated sequence variants in a critical interval of 7.2 Mb on cattle chromosome 28: a missense point mutation located in an uncharacterized locus and an 855 bp deletion across the exon 19/intron 19 border of the bovine nidogen 1 (NID1) gene (c.3579_3604+829del). RT-PCR showed that NID1 is expressed in bovine lenses while the transcript of the second locus was absent. The NID1 deletion leads to the skipping of exon 19 during transcription and is therefore predicted to cause a frameshift and premature stop codon (p.1164fs27X). The truncated protein lacks a C-terminal domain essential for binding with matrix assembly complexes. Nidogen 1 deficient mice show neurological abnormalities and highly irregular crystal lens alterations. This study adds NID1 to the list of candidate genes for inherited cataract in humans and is the first report of a naturally occurring mutation leading to non-syndromic catarct in cattle provides a potential large animal model for human cataract.