A shared 336 kb haplotype associated with the belt pattern in three divergent cattle breeds

We recently mapped the belt mutation in Brown Swiss cattle to a 922 kb interval on BTA3. In this study, we analysed two additional cattle breeds with the belted phenotype: Galloway and Dutch Belted (Lakenvelder). By genotyping microsatellites in solid-coloured and belted Galloways, we confirmed that the belt mutation in Galloways is strongly associated with the same chromosomal locus as in Brown Swiss cattle. Subsequently, we analysed 36 SNPs in the belt interval in three breeds. We identified a single belt-associated haplotype for each of the analysed breeds. The three breed-specific belt haplotypes share alleles in four blocks. Three of these blocks comprise only one single or two consecutive markers, while the largest shared haplotype block encompasses nine consecutive SNPs in a 336 kb interval. The large shared haplotype across divergent breeds suggests a common mutation for the belt phenotype in all three breeds. We identified a potential candidate gene within this interval coding for the developmental transcription factor HES6. We re-sequenced the complete HES6 coding sequence in belted and solid-coloured cattle but did not find belt-associated polymorphisms. In conclusion, our data provide strong evidence in favour of a common founder for the belt phenotype in different cattle breeds and have resulted in an improved fine-mapping of the causative mutation.

Globally dispersed Y chromosomal haplotypes in wild and domestic sheep

To date, investigations of genetic diversity and the origins of domestication in sheep have utilised autosomal microsatellites and variation in the mitochondrial genome. We present the first analysis of both domestic and wild sheep using genetic markers residing on the ovine Y chromosome. Analysis of a single nucleotide polymorphism (oY1) in the SRY promoter region revealed that allele A-oY1 was present in all wild bighorn sheep (Ovis canadensis), two subspecies of thinhorn sheep (Ovis dalli), European Mouflon (Ovis musimon) and the Barbary (Ammontragis lervia). A-oY1 also had the highest frequency (71.4%) within 458 domestic sheep drawn from 65 breeds sampled from Africa, Asia, Australia, the Caribbean, Europe, the Middle East and Central Asia. Sequence analysis of a second locus, microsatellite SRYM18, revealed a compound repeat array displaying fixed differences, which identified bighorn and thinhorn sheep as distinct from the European Mouflon and domestic animals. Combined genotypic data identified 11 male-specific haplotypes that represented at least two separate lineages. Investigation of the geographical distribution of each haplotype revealed that one (H6) was both very common and widespread in the global sample of domestic breeds. The remaining haplotypes each displayed more restricted and informative distributions. For example, H5 was likely founded following the domestication of European breeds and was used to trace the recent transportation of animals to both the Caribbean and Australia. A high rate of Y chromosomal dispersal appears to have taken place during the development of domestic sheep as only 12.9% of the total observed variation was partitioned between major geographical regions.

Genome-wide search for markers associated with osteochondrosis in Hanoverian warmblood horses

A genome-wide scan was performed to detect quantitative trait loci (QTLs) for osteochondrosis (OC) and osteochondrosis dissecans (OCD) in horses. The marker set comprised 260 microsatellites. We collected data from 211 Hanoverian warmblood horses consisting of 14 paternal half-sib families. Traits used were OC (fetlock and/or hock joints affected), OCD (fetlock and/or hock joints affected), fetlock OC, fetlock OCD, hock OC, and hock OCD. The first genome scan included 172 microsatellite markers. In a second step 88 additional markers were chosen to refine putative QTLs found in the first scan. Genome-wide significant QTLs were located on equine chromosomes 2, 4, 5, and 16. QTLs for fetlock OC and hock OC partly overlapped on the same chromosomes, indicating that these traits may be genetically related. QTLs reached the chromosome-wide significance level on eight different equine chromosomes: 2, 3, 4, 5, 15, 16, 19, and 21. This whole-genome scan was a first step toward the identification of candidate genome regions harboring genes responsible for equine OC. Further investigations are necessary to refine the map positions of the QTLs already identified for OC.

Linkage mapping of ovine microphthalmia to chromosome 23, the sheep orthologue of human chromosome 18

PURPOSE: To characterize the phenotype and map the locus responsible for autosomal recessive inherited ovine microphthalmia (OMO) in sheep. METHODS: Microphthalmia-affected lambs and their available relatives were collected in a field, and experimental matings were performed to obtain affected and normal lambs for detailed necropsy and histologic examinations. The matings resulted in 18 sheep families with 48 cases of microphthalmia. A comparative candidate gene approach was used to map the disease locus within the sheep genome. Initially, 27 loci responsible for the microphthalmia-anophthalmia phenotypes in humans or mice were selected to test for comparative linkage. Fifty flanking markers that were predicted from comparative genomic analysis to be closely linked to these genes were tested for linkage to the disease locus. After observation of statistical evidence for linkage, a confirmatory fine mapping strategy was applied by further genotyping of 43 microsatellites. RESULTS: The clinical and pathologic examinations showed slightly variable expressivity of isolated bilateral microphthalmia. The anterior eye chamber was small or absent, and a white mass admixed with cystic spaces extended from the papilla to the anterior eye chamber, while no recognizable vitreous body or lens was found within the affected eyes. Significant linkage to a single candidate region was identified at sheep chromosome 23. Fine mapping and haplotype analysis assigned the candidate region to a critical interval of 12.4 cM. This ovine chromosome segment encompasses an ancestral chromosomal breakpoint corresponding to two orthologue segments of human chromosomes 18, short and long arms. For the examined animals, we excluded the complete coding region and adjacent intronic regions of ovine TGIF1 to harbor disease-causing mutations. CONCLUSIONS: This is the first genetic localization for hereditary ovine isolated microphthalmia. It seems unlikely that a mutation in the TGIF1 gene is responsible for this disorder. The studied sheep represent a valuable large animal model for similar human ocular phenotypes.

Genomics of liver fibrosis and cirrhosis

Hepatic fibrosis is the response to chronic injury from viral, toxic, metabolic, cholestatic, or autoimmune liver injury. However, only a minority of affected individuals develop advanced fibrosis or cirrhosis, suggesting that modifiers determine the individual risk. Aside from well-established progression factors including gender, duration of exposure to the disease, and ethnicity, unknown host genetic factors are likely to influence disease progression and prognosis. Potential genetic susceptibility loci are single nucleotide polymorphisms in fibrosis-associated genes, point mutations, microsatellites, and haplotype blocks composed of genes pivotal for fibrosis development. However, the study of complex polygenetic diseases poses numerous pitfalls in contrast to the elucidation of monogenetic (i.e., Mendelian) diseases. Many publications on the role of certain genetic variants in modulating the progression of hepatic fibrosis have been limited by inadequate study design and low statistical power. At present, powerful research strategies are being developed that allow the application of knowledge derived from the successful sequencing of the human genome that will help to translate our newly acquired insight into practical improvements for research activities and medical practice.

Speciation with gene-flow in eastern North American red oaks (Quercus section Lobatae)

As foundational species, oaks (Quercus : Fagaceae) support the activities of both humans and wildlife. However, many oaks in North America are declining, a crisis exacerbated by the previous disappearance of other hard mast-producing trees. In addition, the economic demands placed on this drought-tolerant group may intensify if climate change extirpates other, relatively mesophytic species. Genetic tools can help address these management challenges. To this end, we developed a suite of 27 microsatellite markers, of which 22 are derived from expressed sequence tags (ESTs). Many of these markers bear significant homology to known genes and may be able to directly assay functional genetic variation. Markers obtained from enriched microsatellite libraries, on the other hand, are typically located in heterochromatic regions and should reflect demographic processes. Considered jointly, genic and genomic microsatellites can elucidate patterns of gene-flow and natural selection, which are fundamental to both an organism's evolutionary ecology and conservation biology. To this end, we employed the developed markers in an FST-based genome scan to detect the signature of divergent selection among the red oaks (Quercus section Lobatae). Three candidate genes with putative roles in stress responses demonstrated patterns of diversity consistent with adaptation to heterogeneous selective pressures. These genes may be important in both local genetic adaptation within species and divergence among them. Next, we used an isolation-with-migration model to quantify levels of gene-flow among four red oaks species during speciation. Both speciation in allopatry and speciation with gene-flow were found to be major drivers of red oak biodiversity. Loci playing a key role in speciation are also likely to be ecologically important within species

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.

Genetic diversity of Swiss sheep breeds in the focus of conservation research

There is constant pressure to improve evaluation of animal genetic resources in order to prevent their erosion. Maintaining the integrity of livestock species as well as their genetic diversity is of paramount interest for long-term agricultural policies. One major use of DNA techniques in conservation is to reveal genetic diversity within and between populations. Forty-one microsatellites were analysed to assess genetic diversity in nine Swiss sheep breeds and to measure the loss of the overall diversity when one breed would become extinct. The expected heterozygosities varied from 0.65 to 0.74 and 10.8% of the total genetic diversity can be explained by the variation among breeds. Based on the proportion of shared alleles, each of the nine breeds were clearly defined in their own cluster in the neighbour-joining tree describing the relationships among the breeds. Bayesian clustering methods assign individuals to groups based on their genetic similarity and infer the number of populations. In STRUCTURE, this approach pooled the Valais Blacknose and the Valais Red. With BAPS method the two Valais sheep breeds could be separated. Caballero & Toro approach (2002) was used to calculate the loss or gain of genetic diversity when each of the breeds would be removed from the set. The changes in diversity based on between-breed variation ranged from -12.2% (Valais Blacknose) to 0% (Swiss Black Brown Mountain and Mirror Sheep); based on within-breed diversity the removal of a breed could also produce an increase in diversity (-0.6% to + 0.6%). Allelic richness ranged from 4.9 (Valais Red) to 6.7 (Brown Headed Meat sheep and Red Engadine Sheep). Breed conservation decisions cannot be limited to genetic diversity alone. In Switzerland, conservation goals are embedded in the desire to carry the cultural legacy over to future generations.

Arachnomelia in Brown Swiss cattle maps to chromosome 5

Arachnomelia in Brown Swiss cattle is a monogenic autosomal recessive inherited congenital disorder of the skeletal system giving affected calves a spidery look (OMIA ID 000059). Over a period of 20 years 15 cases were sampled in the Swiss and Italian Brown cattle population. Pedigree data revealed that all affected individuals trace back to a single acknowledged carrier founder sire. A genome scan using 240 microsatellites spanning the 29 bovine autosomes showed homozygosity at three adjacent microsatellite markers on bovine Chr 5 in all cases. Linkage analysis confirmed the localization of the arachnomelia mutation in the region of the marker ETH10. Fine-mapping and haplotype analysis using a total of 34 markers in this region refined the critical region of the arachnomelia locus to a 7.19-Mb interval on bovine Chr 5. The disease-associated IBD haplotype was shared by 36 proven carrier animals and allows marker-assisted selection. As the corresponding human and mouse chromosome segments do not contain any clear functional candidate genes for this disorder, the mutation causing arachnomelia in the Brown Swiss cattle might help to identify an unknown gene in bone development.

Major histocompatibility complex and other allergy-related candidate genes associated with insect bite hypersensitivity in Icelandic horses

Insect bite hypersensitivity (IBH) is an allergic dermatitis of horses caused by bites of insects. IBH is a multifactorial disease with contribution of genetic and environmental factors. Candidate gene association analysis of IBH was performed in a group of 89 Icelandic horses all born in Iceland and imported to Europe. Horses were classified in IBH-affected and non-affected based on clinical signs and history of recurrent dermatitis, and on the results of an in vitro sulfidoleukotriene (sLT)-release assay with Culicoides nubeculosus and Simulium vittatum extract. Different genetic markers were tested for association with IBH by the Fisher's exact test. The effect of the major histocompatibility complex (MHC) gene region was studied by genotyping five microsatellites spanning the MHC region (COR112, COR113, COR114, UM011 and UMN-JH34-2), and exon 2 polymorphisms of the class II Eqca-DRA gene. Associations with Eqca-DRA and COR113 were identified (p < 0.05). In addition, a panel of 20 single nucleotide polymorphisms (SNPs) in 17 candidate allergy-related genes was tested. During the initial screen, no marker from the panel was significantly (p < 0.05) associated with IBH. Five SNPs associated with IBH at p < 0.10 were therefore used for analysis of combined genotypes. Out of them, SNPs located in the genes coding for the CD14 receptor (CD14), interleukin 23 receptor (IL23R), thymic stromal lymphopoietin (TSLP) and transforming growth factor beta 3 (TGFB3) molecules were associated with IBH as parts of complex genotypes. These results are supported by similar associations and by expression data from different horse populations and from human studies.

ABERRATIONS OF A PUTATIVE TUMOR SUPPRESSOR GENE SEL1L IN PANCREATIC DUCTAL ADENOCARCINOMA

Introduction: Pancreatic cancer is the fourth leading cause of cancer-related death among males and females in the United States. Sel-1-like (SEL1L) is a putative tumor suppressor gene that is downregulated in a significant proportion of human pancreatic ductal adenocarcinoma (PDAC). It was hypothesized that SEL1L expression could be down-modulated by somatic mutation, loss of heterozygosity (LOH), CpG island hypermethylation and/or aberrantly expressed microRNAs (miRNAs). Material and methods: In 42 PDAC tumors, the SEL1L coding region was amplified using reverse transcription polymerase chain reaction (RT-PCR), and analyzed by agarose gel electrophoresis and sequenced to search for mutations. Using fluorescent fragment analysis, two intragenic microsatellites in the SEL1L gene region were examined to detect LOH in a total of 73 pairs of PDAC tumors and normal-appearing adjacent tissues. Bisulfite DNA sequencing was performed to determine the methylation status of the SEL1L promoter in 41 PDAC tumors and 6 PDAC cell lines. Using real-time quantitative PCR, the expression levels of SEL1L mRNA and 7 aberrantly upregulated miRNAs that potentially target SEL1L were assessed in 42 PDAC tumor and normal pairs. Statistical methods were applied to evaluate the correlation between SEL1L mRNA and the miRNAs. Further the interaction was determined by functional analysis using a molecular biological approach. Results: No mutations were detected in the SEL1L coding region. More than 50% of the samples displayed abnormally alternate or aberrant spliced transcripts of SEL1L. About 14.5% of the tumors displayed LOH at the CAR/CAL microsatellite locus and 10.7% at the RepIN20 microsatellite locus. However, the presence of LOH did not show significant association with SEL1L downregulation. No methylation was observed in the SEL1L promoter. Statistical analysis showed that SEL1L mRNA expression levels significantly and inversely correlated with the expression of hsa-mir-143, hsa-mir-155, and hsa-mir-223. Functional analysis indicated that hsa-mir-155 acted as a suppressor of SEL1L in PL18 and MDAPanc3 PDAC cell lines. Discussion: Evidence from these studies suggested that SEL1L was possibly downregulated by aberrantly upregulated miRNAs in PDAC. Future studies should be directed towards developing a better understanding of the mechanisms for generation of aberrant SEL1L transcripts, and further analysis of miRNAs that may downregulate SEL1L.

Genetics of recurrent airway obstruction (RAO)

Recurrent airway obstruction (RAO) is a multifactorial and polygenic disease. Affected horses are typically 7 years of age or older and show exercise intolerance, increased breathing effort, coughing, airway neutrophilia, mucus accumulation and hyperreactivity as well as cholinergic bronchospasm. The environmental factors responsible are predominantly allergens and irritants in haydust, but the immunological mechanisms underlying RAO are still unclear. Several studies have demonstrated a familiar predisposition for RAO and it is now proven that the disease has a genetic basis. In offspring, the risk of developing RAO is 3-fold increased when one parent is affected and increases to almost 5-fold when both parents have RAO. Segregation analysis in two high-prevalence families demonstrated a high heritability and a complex inheritance with several major genes. A whole genomescan showed chromosome-wide significant linkage of seven chromosomal regions with RAO. Of the microsatellites, which were located near atopy candidate genes, those in a region of chromosome 13 harboring the IL4R gene were strongly associated with the RAO phenotype in the offspring of one RAO-affected stallion. Furthermore, IgE-levels are influenced by hereditary factors in the horse, and we have evidence that RAO-affected offspring of the same stallion have increased levels of specific IgE against moldspore allergens. The identification of genetic markers and ultimately of the responsible genes will not only allow for an improved prophylaxis, i.e. early identification of susceptible individuals and avoidance of high-risk matings, but also improve our ability to find new therapeutic targets and to optimize existing treatments.

Rapid and accurate G M1-gangliosidosis diagnosis using a parentage testing microsatellite

The G M1-gangliosidosis is an autosomal recessive lysosomal storage disease caused by structural defects of the beta-galactosidase gene (GLB1) which lead to a severe phenotypical impairment in homozygous individuals, whereas heterozygous carriers remain clinically normal. Currently employed DNA parentage tests include the analysis of microsatellites, which also have a diagnostic predictive value. The aim of this study was to provide a reliable tool for genotyping the canine GLB1 which can be effectively integrated in parentage testing investigations. For this purpose the association between the GLB1 gene and the AHT K253 microsatellite was analyzed in 30 Alaskan huskies (11 GLB1+/+, 17 GLB1+/- and 2 GLB1-/- dogs). The 143 bp AHT K253 microsatellite allele was identified only in GLB1+/- and GLB1-/- animals and was in strong linkage disequilibrium with the causative mutation for G M1-gangliosidosis, a 19 bp duplication within exon 15 of the GLB1 gene. The results of the present study revealed a 100% concordance between the previous established genotypes and those obtained after the analysis of the AHT K253 microsatellite. Thus, the genotype of the AHT K253 microsatellite, which is routinely determined during dog parentage testing, has a high predictive value for the G M1-gangliosidosis carrier status.

A comparative radiation hybrid map of sheep chromosome 10

Comparative radiation hybrid (RH) maps of individual ovine chromosomes are essential to identify genes governing traits of economic importance in sheep, a livestock species for which whole genome sequence data are not yet available. The USUoRH5000 radiation hybrid panel was used to generate a RH map of sheep chromosome 10 (OAR10) with 59 markers that span 1,422 cR over an estimated 92 Mb of the chromosome, thus providing markers every 2 Mb (equivalent to every 24 cR). The markers were derived from 46 BAC end sequences (BESs), a single EST, and 12 microsatellites. Comparative analysis showed that OAR10 shares remarkable conservation of gene order along the entire length of cattle chromosome 12 and that OAR10 contains four major homologous synteny blocks, each related to segments of the homologous human chromosome 13. Extending the comparison to the horse, dog, mouse, and chicken genome showed that these blocks share conserved synteny across species.

Hybrid ‘superswarm’ leads to rapid divergence and establishment of populations during a biological invasion

Understanding the genetic background of invading species can be crucial information clarifying why they become invasive. Intraspecific genetic admixture among lineages separated in the native ranges may promote the rate and extent of an invasion by substantially increasing standing genetic variation. Here we examine the genetic relationships among threespine stickleback that recently colonized Switzerland. This invasion results from several distinct genetic lineages that colonized multiple locations and have since undergone range expansions, where they coexist and admix in parts of their range. Using 17 microsatellites genotyped for 634 individuals collected from 17 Swiss and two non-Swiss European sites, we reconstruct the invasion of stickleback and investigate the potential and extent of admixture and hybridization among the colonizing lineages from a population genetic perspective. Specifically we test for an increase in standing genetic variation in populations where multiple lineages coexist. We find strong evidence of massive hybridization early on, followed by what appears to be recent increased genetic isolation and the formation of several new genetically distinguishable populations, consistent with a hybrid ‘superswarm’. This massive hybridization and population formation event(s) occurred over approximately 140 years and likely fuelled the successful invasion of a diverse range of habitats. The implications are that multiple colonizations coupled with hybridization can lead to the formation of new stable genetic populations potentially kick-starting speciation and adaptive radiation over a very short time.