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

Genome-wide association with footrot in hair and wool sheep

Ovine footrot is an infectious disease with important contributions from Dichelobacter nodosus and Fusobacterium necrophorum. Footrot is characterized by separation of the hoof from underlying tissue, and this causes severe lameness that negatively impacts animal wellbeing, growth, and profitability. Large economic losses result from lost production as well as treatment costs, and improved genetic tools to address footrot are a valuable long-term goal. Prior genetic studies had examined European wool sheep, but hair sheep breeds such as Katahdin and Blackbelly have been reported to have increased resistance to footrot, as well as to intestinal parasites. Thus, footrot condition scores were collected from 251 U.S. sheep including Katahdin, Blackbelly, and European-influenced crossbred sheep with direct and imputed genotypes at OvineHD array (>500,000 single nucleotide polymorphism) density. Genome-wide association was performed with a mixed model accounting for farm and principal components derived from animal genotypes, as well as a random term for the genomic relationship matrix. We identified three genome-wide significant associations, including SNPs in or near GBP6 and TCHH. We also identified 33 additional associated SNPs with genome-wide suggestive evidence, including a cluster of 6 SNPs in a peak near the genome-wide significance threshold located near the glutamine transporter gene SLC38A1. These findings suggest genetic susceptibility to footrot may be influenced by genes involved in divergent biological processes such as immune responses, nutrient availability, and hoof growth and integrity. This is the first genome-wide study to investigate susceptibility to footrot by including hair sheep and also the first study of any kind to identify multiple genome-wide significant associations with ovine footrot. These results provide a foundation for developing genetic tests for marker-assisted selection to improve resistance to ovine footrot once additional steps like fine mapping and validation are complete.

An inherited night blindness in Wiltshire sheep

Twelve cases of adult-onset blindness were identified in a flock of 130 polled Wiltshire sheep in New Zealand over a 3-year period. Affected sheep developed night blindness between 2 and 3 years of age, which progressed to complete blindness by 4 to 5 years of age. Fundic examination findings included progressive tapetal hyperreflectivity and attenuation of retinal blood vessels. Histologically, the retinas had a selective loss of rod photoreceptors with initial preservation of cone photoreceptors. Retinal degeneration was not accompanied by any other ocular or central nervous system abnormalities, and pedigree analysis suggested an inherited basis for the disease. Mating an affected Wiltshire ram to 2 affected Wiltshire ewes resulted in 6 progeny that all developed retinal degeneration by 2 years of age, while mating of the same affected ram to 6 unaffected ewes resulted in 8 unaffected progeny, consistent with autosomal recessive inheritance. Homozygosity mapping of 5 affected Wiltshire sheep and 1 unaffected Wiltshire sheep using an OvineSNP50 Genotyping BeadChip revealed an identical-by-descent region on chromosome 5, but none of the genes within this region were considered plausible candidate genes. Whole-genome sequencing of 2 affected sheep did not reveal any significant mutations in any of the genes associated with retinitis pigmentosa in humans or progressive retinal atrophy in dogs. Inherited progressive retinal degeneration affecting rod photoreceptors has not been previously reported in sheep, but this disease has several similarities to inherited retinal dystrophies in other species.

A Nonsense Variant in CCDC65 Gene Causes Respiratory Failure Associated with Increased Lamb Mortality in French Lacaune Dairy Sheep

We recently demonstrated that the Lacaune deficient homozygous haplotype 6 (LDHH6) potentially hosts a recessive perinatal lethal mutation in Lacaune dairy sheep mapped on OAR3. In the present study, we have analyzed the whole-genome sequences of two Lacaune ram heterozygous carriers of LDHH6. After variant calling and filtering against the variants of 86 non-carrier rams, we have identified a single nucleotide variant (SNV) in the two LDHH6 carriers whose variant allele induced a premature stop codon (p.Glu111*) in the Coiled-Coil Domain Containing 65 (CCDC65) gene. CCDC65 is involved in the assembly of the nexin-dynein regulatory complex for the formation of microtubules in ciliated cells. In order to identify the phenotype in homozygous sheep, we generated at-risk matings (n = 17) between rams and ewes heterozygous for the candidate variant in CCDC65. A total of 16 lambs were born alive with five genotyped as homozygous carriers. The homozygous lambs suffered from respiratory problems, and four of them died within the first month of life. At necropsy, we observed a broad hepatization of lung lobes possibly induced by infectious pneumonia. The management of this lethal recessive allele (frequency of 0.06) through reasoned mating in the Lacaune sheep selection schemes could reduce lamb mortality by 2%.

Direct and indirect contributions of molecular genetics to farm animal welfare: a review

Since domestication, farm animals have played a key role to increase the prosperity of humankind, while animal welfare (AW) is debated even today. This paper aims to comprehensively review the contributions of developing molecular genetics to farm animal welfare (FAW) and to raise awareness among both scientists and farmers about AW. Welfare is a complex trait affected by genetic structure and environmental factors. Therefore, the best welfare status can be achieved not only to enhance environmental factors such as management and feeding practices, but also the genetic structure of animals must be improved. In this regard, advances in molecular genetics have made great contributions to improve the genetic structure of farm animals, which has increased AW. Today, by sequencing and/or molecular markers, genetic diseases may be detected and eliminated in local herds. Additionally, genes related to diseases or adaptations are investigated by molecular techniques, and the frequencies of desired genotypes are increased in farm animals to keep welfare at an optimized level. Furthermore, stress on animals can be reduced with DNA extraction from stool and feather samples which reduces physical contact between animals and veterinarians. Together with molecular genetics, advances in genome editing tools and biotechnology are promising to improve FAW in the future.

Genes involved in immune, gene translation and chromatin organization pathways associated with Mycoplasma ovipneumoniae presence in nasal secretions of domestic sheep

Mycoplasma ovipneumoniae contributes to polymicrobial pneumonia in domestic sheep. Elucidation of host genetic influences of M. ovipneumoniae nasal detection has the potential to reduce the incidence of polymicrobial pneumonia in sheep through implementation of selective breeding strategies. Nasal mucosal secretions were collected from 647 sheep from a large US sheep flock. Ewes of three breeds (Polypay n = 222, Rambouillet n = 321, and Suffolk n = 104) ranging in age from one to seven years, were sampled at three different times in the production cycle (February, April, and September/October) over four years (2015 to 2018). The presence and DNA copy number of M. ovipneumoniae was determined using a newly developed species-specific qPCR. Breed (P<0.001), age (P<0.024), sampling time (P<0.001), and year (P<0.001) of collection affected log10 transformed M. ovipneumoniae DNA copy number, where Rambouillet had the lowest (P<0.0001) compared with both Polypay and Suffolk demonstrating a possible genetic component to detection. Samples from yearlings, April, and 2018 had the highest (P<0.046) detected DNA copy number mean. Sheep genomic DNA was genotyped with the Illumina OvineHD BeadChip. Principal component analysis identified most of the variation in the dataset was associated with breed. Therefore, genome wide association analysis was conducted with a mixed model (EMMAX), with principal components 1 to 6 as fixed and a kinship matrix as random effects. Genome-wide significant (P<9x10-8) SNPs were identified on chromosomes 6 and 7 in the all-breed analysis. Individual breed analysis had genome-wide significant (P<9x10-8) SNPs on chromosomes 3, 4, 7, 9, 10, 15, 17, and 22. Annotated genes near these SNPs are part of immune (ANAPC7, CUL5, TMEM229B, PTPN13), gene translation (PIWIL4), and chromatin organization (KDM2B) pathways. Immune genes are expected to have increased expression when leukocytes encounter M. ovipneumoniae which would lead to chromatin reorganization. Work is underway to narrow the range of these associated regions to identify the underlying causal mutations.

Deletion of Lrp4 increases the incidence of microphthalmia

Microphthalmia is a malformation that reduces the size of the ocular globe. The etiologies of this anomaly are various, but the genetic background appears to have a predominant influence on its development through mutations of genes controlling ocular developmental processes. LRP4 is a type I single transmembrane protein that is essential for the formation of neuromuscular junctions. We created and experimented on homozygous Lrp4-deficient mice and found the microphthalmia phenotype in their eyes. The loss of Lrp4 resulted in an elevated incidence of microphthalmia and affected the mRNA expression of the members of bone morphogenetic protein, fibroblast growth factor, Sonic hedgehog, and WNT signaling pathways and of several pathogenic genes for microphthalmia. Moreover, the loss of Lrp4 enhanced the incidence of aberrant retinal folds, which appeared pleated and corrugated in the eyeball.

Congenital Ocular Abnormalities in Free-Ranging White-Tailed Deer

Congenital ocular abnormalities in cervids have been previously reported as individual cases from various regions of the United States and include microphthalmia, anophthalmia, congenital cataracts, dermoids, and colobomata. A common underlying cause for these abnormalities, such as nutritional deficiencies, environmental toxin exposures, or genetic mutations, has not been established. This retrospective study summarized and compared cases of suspected congenital ocular abnormalities in free-ranging white-tailed deer ( Odocoileus virginianus) submitted to the Southeastern Cooperative Wildlife Disease Study (SCWDS) in Athens, Georgia, to the preexisting literature. Of 3645 accessions of white-tailed deer submitted to SCWDS, 15 qualifying case records were found. An additional 15 cases were reported previously in the literature. Conditions described in SCWDS cases included microphthalmia (8/15), congenital cataracts (3/15), anophthalmia (2/15), colobomata (1/15), anterior segment dysgenesis (1/15), ectopic lacrimal gland tissue (1/15), and congenital blindness with corneal opacity (1/15). Most (11/15; 73%) of the SCWDS cases were male fawns with an average age of 4 months at presentation, consistent with previously described cases. Most animals had bilateral abnormalities with few extraocular congenital abnormalities, also consistent with existing reports. Cases were variably tested for various infectious agents at the time of submission; 2 cases were seropositive for bluetongue virus. Spatiotemporal clustering of cases was not evident. This study provided a concise and systematic summary of known existing cases of congenital ocular defects in fawns but did not identify a cause.

Unveiling genomic regions that underlie differences between Afec-Assaf sheep and its parental Awassi breed

BACKGROUND

Sheep production in Israel has improved by crossing the fat-tailed local Awassi breed with the East Friesian and later, with the Booroola Merino breed, which led to the formation of the highly prolific Afec-Assaf strain. This strain differs from its parental Awassi breed in morphological traits such as tail and horn size, coat pigmentation and wool characteristics, as well as in production, reproductive and health traits. To identify major genes associated with the formation of the Afec-Assaf strain, we genotyped 41 Awassi and 141 Afec-Assaf sheep using the Illumina Ovine SNP50 BeadChip array, and analyzed the results with PLINK and EMMAX software. The detected variable genomic regions that differed between Awassi and Afec-Assaf sheep (variable genomic regions; VGR) were compared to selection signatures that were reported in 48 published genome-wide association studies in sheep. Because the Afec-Assaf strain, but not the Awassi breed, carries the Booroola mutation, association analysis of BMPR1B used as the test gene was performed to evaluate the ability of this study to identify a VGR that includes such a major gene.

RESULTS

Of the 20 detected VGR, 12 were novel to this study. A ~7-Mb VGR was identified on Ovies aries chromosome OAR6 where the Booroola mutation is located. Similar to other studies, the most significant VGR was detected on OAR10, in a region that contains candidate genes affecting horn type (RXFP2), climate adaptation (ALOX5AP), fiber diameter (KATNAl1), coat pigmentation (FRY) and genes associated with fat distribution. The VGR on OAR2 included BNC2, which is also involved in controlling coat pigmentation in sheep. Six other VGR contained genes that were shown to be involved in coat pigmentation by analyzing their mammalian orthologues. Genes associated with fat distribution in humans, including GRB14 and COBLL1, were located in additional VGR. Sequencing DNA from Awassi and Afec-Assaf individuals revealed non-synonymous mutations in some of these candidate genes.

CONCLUSIONS

Our results highlight VGR that differentiate the Awassi breed from the Afec-Assaf strain, some of which may include genes that confer an advantage to Afec-Assaf and Assaf over Awassi sheep with respect to intensive sheep production under Mediterranean conditions.

Genomic prediction and genome-wide association study for dagginess and host internal parasite resistance in New Zealand sheep

Background

Dagginess (faecal soiling of the perineum region) and host nematode parasite resistance are important animal welfare traits in New Zealand sheep. Genomic prediction (GP) estimates the genetic merit, as a molecular breeding value (mBV), for each trait based on many SNPs. The additional information the mBV provides (as determined by its accuracy) has led to its incorporation into breeding schemes. Some GP methods give SNP effects, which provide additional information to identify genome-wide associations (GWAS) for a trait of interest. Here we report results from a GP and GWAS study for dagginess and host nematode parasite resistance in a New Zealand sheep industry resource.

Results

Genomic prediction analysis was performed using 50K SNP chip data and parent average-removed, de-regressed BVs for five traits, from a resource of 8705 pedigree recorded animals. The five traits were dag score at three and eight months (DAG3, DAG8) and nematode faecal egg count in summer (FEC1), autumn (FEC2) and as an adult (AFEC). The resource consisted of Romney, Coopworth, Perendale, Texel and various breed crosses (designated: CompRCP, CompRCPT and CompCRP). The pure breeds, apart from Texel, plus CompRCP were used to develop the GP. The resulting SNP effects were used to identify genetic regions associated with dagginess and parasite resistance. Accuracies of the weighted correlation between mBV and true BV ranged between −0.07 (Texel) and 0.56 (Coopworth) for DAG3 and DAG8. For FEC1, FEC2 and AFEC accuracies ranged between −0.22 (CompRCPT) and 0.69 (Coopworth). The weighted average individual accuracy (calculated from theory) ranges were 0.13 (Texel) to 0.52 (Coopworth) and 0.11 (Texel) to 0.55 (Coopworth) respectively, for dagginess and parasite traits. There was one SNP for DAG8 that reached Bonferroni significance threshold (P < 1 × 10⁻⁶) on OAR15, the same two SNPs for each of the parasite traits (OAR26) and none for DAG3. A notable peak was also observed on OAR7 for all the parasite traits, however, it did not reach the Bonferroni significance threshold.

Conclusions

This study presents the first results of a GWAS on dagginess and faecal egg count traits in New Zealand sheep. The results suggest that there are quantitative trait loci on OAR 15 for dagginess and on OAR26 and seven for faecal egg count.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2148-2) contains supplementary material, which is available to authorized users.

Genome-Wide Association Identifies SLC2A9 and NLN Gene Regions as Associated with Entropion in Domestic Sheep

Entropion is an inward rolling of the eyelid allowing contact between the eyelashes and cornea that may lead to blindness if not corrected. Although many mammalian species, including humans and dogs, are afflicted by congenital entropion, no specific genes or gene regions related to development of entropion have been reported in any mammalian species to date. Entropion in domestic sheep is known to have a genetic component therefore, we used domestic sheep as a model system to identify genomic regions containing genes associated with entropion. A genome-wide association was conducted with congenital entropion in 998 Columbia, Polypay, and Rambouillet sheep genotyped with 50,000 SNP markers. Prevalence of entropion was 6.01%, with all breeds represented. Logistic regression was performed in PLINK with additive allelic, recessive, dominant, and genotypic inheritance models. Two genome-wide significant (empirical P<0.05) SNP were identified, specifically markers in SLC2A9 (empirical P = 0.007; genotypic model) and near NLN (empirical P = 0.026; dominance model). Six additional genome-wide suggestive SNP (nominal P<1x10^-5) were identified including markers in or near PIK3CB (P = 2.22x10^-6; additive model), KCNB1 (P = 2.93x10^-6; dominance model), ZC3H12C (P = 3.25x10^-6; genotypic model), JPH1 (P = 4.68x20^-6; genotypic model), and MYO3B (P = 5.74x10^-6; recessive model). This is the first report of specific gene regions associated with congenital entropion in any mammalian species, to our knowledge. Further, none of these genes have previously been associated with any eyelid traits. These results represent the first genome-wide analysis of gene regions associated with entropion and provide target regions for the development of sheep genetic markers for marker-assisted selection.

Conserved genetic pathways associated with microphthalmia, anophthalmia, and coloboma

The human eye is a complex organ whose development requires extraordinary coordination of developmental processes. The conservation of ocular developmental steps in vertebrates suggests possible common genetic mechanisms. Genetic diseases involving the eye represent a leading cause of blindness in children and adults. During the last decades, there has been an exponential increase in genetic studies of ocular disorders. In this review, we summarize current success in identification of genes responsible for microphthalmia, anophthalmia, and coloboma (MAC) phenotypes, which are associated with early defects in embryonic eye development. Studies in animal models for the orthologous genes identified overlapping phenotypes for most factors, confirming the conservation of their function in vertebrate development. These animal models allow for further investigation of the mechanisms of MAC, integration of various identified genes into common developmental pathways and finally, provide an avenue for the development and testing of therapeutic interventions.

Selection signatures in worldwide sheep populations

The diversity of populations in domestic species offers great opportunities to study genome response to selection. The recently published Sheep HapMap dataset is a great example of characterization of the world wide genetic diversity in sheep. In this study, we re-analyzed the Sheep HapMap dataset to identify selection signatures in worldwide sheep populations. Compared to previous analyses, we made use of statistical methods that (i) take account of the hierarchical structure of sheep populations, (ii) make use of linkage disequilibrium information and (iii) focus specifically on either recent or older selection signatures. We show that this allows pinpointing several new selection signatures in the sheep genome and distinguishing those related to modern breeding objectives and to earlier post-domestication constraints. The newly identified regions, together with the ones previously identified, reveal the extensive genome response to selection on morphology, color and adaptation to new environments.

Signatures of selection in sheep bred for resistance or susceptibility to gastrointestinal nematodes

Background

Gastrointestinal nematodes are one of the most serious causes of disease in domestic ruminants worldwide. There is considerable variation in resistance to gastrointestinal nematodes within and between sheep breeds, which appears to be due to underlying genetic diversity. Through selection of resistant animals, rapid genetic progress has been demonstrated in both research and commercial flocks. Recent advances in genome sequencing and genomic technologies provide new opportunities to understand the ovine host response to gastrointestinal nematodes at the molecular level, and to identify polymorphisms conferring nematode resistance.

Results

Divergent lines of Romney and Perendale sheep, selectively bred for high and low faecal nematode egg count, were genotyped using the Illumina® Ovine SNP50 BeadChip. The resulting genome-wide SNP data were analysed for selective sweeps on loci associated with resistance or susceptibility to gastrointestinal nematode infection. Population differentiation using F_ST and Peddrift revealed sixteen regions, which included candidate genes involved in chitinase activity and the cytokine response. Two of the sixteen regions identified were contained within previously identified QTLs associated with nematode resistance.

Conclusions

In this study we identified fourteen novel regions associated with resistance or susceptibility to gastrointestinal nematodes. Results from this study support the hypothesis that host resistance to internal nematode parasites is likely to be controlled by a number of loci of moderate to small effects.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-637) contains supplementary material, which is available to authorized users.

Genomic variation and population structure detected by single nucleotide polymorphism arrays in Corriedale, Merino and Creole sheep

The aim of this study was to investigate the genetic diversity within and among three breeds of sheep: Corriedale, Merino and Creole. Sheep from the three breeds (Merino n = 110, Corriedale n = 108 and Creole n = 10) were genotyped using the Illumina Ovine SNP50 beadchip^®. Genetic diversity was evaluated by comparing the minor allele frequency (MAF) among breeds. Population structure and genetic differentiation were assessed using STRUCTURE software, principal component analysis (PCA) and fixation index (F_ST). Fixed markers (MAF = 0) that were different among breeds were identified as specific breed markers. Using a subset of 18,181 single nucleotide polymorphisms (SNPs), PCA and STUCTURE analysis were able to explain population stratification within breeds. Merino and Corriedale divergent lines showed high levels of polymorphism (89.4% and 86% of polymorphic SNPs, respectively) and moderate genetic differentiation (F_ST = 0.08) between them. In contrast, Creole had only 69% polymorphic SNPs and showed greater genetic differentiation from the other two breeds (F_ST = 0.17 for both breeds). Hence, a subset of molecular markers present in the OvineSNP50 is informative enough for breed assignment and population structure analysis of commercial and Creole breeds.

Development of the catfish 250K SNP array for genome-wide association studies

BACKGROUND

Quantitative traits, such as disease resistance, are most often controlled by a set of genes involving a complex array of regulation. The dissection of genetic basis of quantitative traits requires large numbers of genetic markers with good genome coverage. The application of next-generation sequencing technologies has allowed discovery of over eight million SNPs in catfish, but the challenge remains as to how to efficiently and economically use such SNP resources for genetic analysis.

RESULTS

In this work, we developed a catfish 250K SNP array using Affymetrix Axiom genotyping technology. The SNPs were obtained from multiple sources including gene-associated SNPs, anonymous genomic SNPs, and inter-specific SNPs. A set of 640K high-quality SNPs obtained following specific requirements of array design were submitted. A panel of 250,113 SNPs was finalized for inclusion on the array. The performance evaluated by genotyping individuals from wild populations and backcross families suggested the good utility of the catfish 250K SNP array.

CONCLUSIONS

This is the first high-density SNP array for catfish. The array should be a valuable resource for genome-wide association studies (GWAS), fine QTL mapping, high-density linkage map construction, haplotype analysis, and whole genome-based selection.

Novel and recurrent PITX3 mutations in Belgian families with autosomal dominant congenital cataract and anterior segment dysgenesis have similar phenotypic and functional characteristics

Background

Congenital cataracts are clinically and genetically heterogeneous with more than 45 known loci and 38 identified genes. They can occur as isolated defects or in association with anterior segment developmental anomalies. One of the disease genes for congenital cataract with or without anterior segment dysgenesis (ASD) is PITX3, encoding a transcription factor with a crucial role in lens and anterior segment development. Only five unique PITX3 mutations have been described, of which the 17-bp duplication c.640_656dup, p.(Gly220Profs*95), is the most common one and the only one known to cause cataract with ASD. The aim of this study was to perform a genetic study of the PITX3 gene in five probands with autosomal dominant congenital cataract (ADCC) and ASD, to compare their clinical presentations to previously reported PITX3-associated phenotypes and to functionally evaluate the PITX3 mutations found.

Methods

Sanger sequencing of the coding region and targeted exons of PITX3 was performed in probands and family members respectively. Transactivation, DNA-binding and subcellular localization assays were performed for the PITX3 mutations found. Ophthalmological examinations included visual acuity measurement, slit-lamp biomicroscopy, tonometry and fundoscopy.

Results

In four Belgian families with ADCC and ASD the recurrent 17-bp duplication c.640_656dup, p.(Gly220Profs*95), was found in a heterozygous state. A novel PITX3 mutation c.573del, p.(Ser192Alafs*117), was identified in heterozygous state in a Belgo-Romanian family with a similar phenotype. Functional assays showed that this novel mutation retains its nuclear localization but results in decreased DNA-binding and transactivation activity, similar to the recurrent duplication.

Conclusions

Our study identified a second PITX3 mutation leading to congenital cataract with ASD. The similarity in phenotypic expression was substantiated by our in vitro functional studies which demonstrated comparable molecular consequences for the novel p.(Ser192Alafs*117) and the recurrent p.(Gly220Profs*95) mutations.

Identification of a 31-bp deletion in the RELN gene causing lissencephaly with cerebellar hypoplasia in sheep

Lissencephaly is an inherited developmental disorder in which neuronal migration is impaired. A type of lissencephaly associated with cerebellar hypoplasia (LCH) was diagnosed in a commercial flock of Spanish Churra sheep. The genotyping of 7 affected animals and 33 controls with the OvineSNP50 BeadChip enabled the localization of the causative mutation for ovine LCH to a 4.8-Mb interval on sheep chromosome 4 using genome-wide association and homozygosity mapping. The RELN gene, which is located within this interval, was considered a strong positional and functional candidate because it plays critical roles in neuronal migration and layer formation. By performing a sequencing analysis of this gene’s specific mRNA in a control lamb, we obtained the complete CDS of the ovine RELN gene. The cDNA sequence from an LCH-affected lamb revealed a deletion of 31 bp (c.5410_5440del) in predicted exon 36 of RELN, resulting in a premature termination codon. A functional analysis of this mutation revealed decreased levels of RELN mRNA and a lack of reelin protein in the brain cortex and blood of affected lambs. This mutation showed a complete concordance with the Mendelian recessive pattern of inheritance observed for the disease. The identification of the causal mutation of LCH in Churra sheep will facilitate the implementation of gene-assisted selection to detect heterozygous mutants, which will help breeders avoid at-risk matings in their flocks. Moreover, the identification of this naturally occurring RELN mutation provides an opportunity to use Churra sheep as a genetically characterized large animal model for the study of reelin functions in the developing and mature brain.

A genome-wide scan study identifies a single nucleotide substitution in ASIP associated with white versus non-white coat-colour variation in sheep (Ovis aries)

In sheep, coat colour (and pattern) is one of the important traits of great biological, economic and social importance. However, the genetics of sheep coat colour has not yet been fully clarified. We conducted a genome-wide association study of sheep coat colours by genotyping 47 303 single-nucleotide polymorphisms (SNPs) in the Finnsheep population in Finland. We identified 35 SNPs associated with all the coat colours studied, which cover genomic regions encompassing three known pigmentation genes (TYRP1, ASIP and MITF) in sheep. Eighteen of these associations were confirmed in further tests between white versus non-white individuals, but none of the 35 associations were significant in the analysis of only non-white colours. Across the tests, the s66432.1 in ASIP showed significant association (P=4.2 × 10(-11) for all the colours; P=2.3 × 10(-11) for white versus non-white colours) with the variation in coat colours and strong linkage disequilibrium with other significant variants surrounding the ASIP gene. The signals detected around the ASIP gene were explained by differences in white versus non-white alleles. Further, a genome scan for selection for white coat pigmentation identified a strong and striking selection signal spanning ASIP. Our study identified the main candidate gene for the coat colour variation between white and non-white as ASIP, an autosomal gene that has been directly implicated in the pathway regulating melanogenesis. Together with ASIP, the two other newly identified genes (TYRP1 and MITF) in the Finnsheep, bordering associated SNPs, represent a new resource for enriching sheep coat-colour genetics and breeding.

Genomewide association for a dominant pigmentation gene in sheep

Most published genomewide association studies (GWAS) in sheep have investigated recessively inherited monogenic traits. The objective here was to assess the feasibility of performing GWAS for a dominant trait for which the genetic basis was already known. A total of 42 Manchega and Rasa Aragonesa sheep that segregate solid black or white coat pigmentation were genotyped using the SNP50 BeadChip. Previous analysis in Manchegas demonstrated a complete association between the pigmentation trait and alleles of the MC1R gene, setting an a priori expectation for GWAS. Multiple methods were used to identify and quantify the strength of population substructure between black and white animals, before allelic association testing was performed for 49,034 SNPs. Following correction for substructure, GWAS identified the most strongly associated SNP (s26449) was also the closest to the MC1R gene. The finding was strongly supported by the permutation tree-based random forest (RF) analysis. Importantly, GWAS identified unlinked SNP with only slightly lower p-values than for s26449. Random forest analysis indicated these were false positives, suggesting interpretation based on both approaches was beneficial. The results indicate that a combined analytical approach can be successful in studies where a modest number of animals are available and substantial population stratification exists.

Retroviral infections in sheep and goats: small ruminant lentiviruses and host interaction

Small ruminant lentiviruses (SRLV) are members of the Retrovirus family comprising the closely related Visna/Maedi Virus (VMV) and the Caprine Arthritis-Encephalitis Virus (CAEV), which infect sheep and goats. Both infect cells of the monocyte/macrophage lineage and cause lifelong infections. Infection by VMV and CAEV can lead to Visna/Maedi (VM) and Caprine Arthritis-Encephalitis (CAE) respectively, slow progressive inflammatory diseases primarily affecting the lungs, nervous system, joints and mammary glands. VM and CAE are distributed worldwide and develop over a period of months or years, always leading to the death of the host, with the consequent economic and welfare implications. Currently, the control of VM and CAE relies on the control of transmission and culling of infected animals. However, there is evidence that host genetics play an important role in determining Susceptibility/Resistance to SRLV infection and disease progression, but little work has been performed in small ruminants. More research is necessary to understand the host-SRLV interaction.

Genome-wide genetic diversity and differentially selected regions among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep

Sheep are among the major economically important livestock species worldwide because the animals produce milk, wool, skin, and meat. In the present study, the Illumina OvineSNP50 BeadChip was used to investigate genetic diversity and genome selection among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds from the United States. After quality-control filtering of SNPs (single nucleotide polymorphisms), we used 48,026 SNPs, including 46,850 SNPs on autosomes that were in Hardy-Weinberg equilibrium and 1,176 SNPs on chromosome × for analysis. Phylogenetic analysis based on all 46,850 SNPs clearly separated Suffolk from Rambouillet, Columbia, Polypay, and Targhee, which was not surprising as Rambouillet contributed to the synthesis of the later three breeds. Based on pair-wise estimates of F_ST, significant genetic differentiation appeared between Suffolk and Rambouillet (F_ST = 0.1621), while Rambouillet and Targhee had the closest relationship (F_ST = 0.0681). A scan of the genome revealed 45 and 41 differentially selected regions (DSRs) between Suffolk and Rambouillet and among Rambouillet-related breed populations, respectively. Our data indicated that regions 13 and 24 between Suffolk and Rambouillet might be good candidates for evaluating breed differences. Furthermore, ovine genome v3.1 assembly was used as reference to link functionally known homologous genes to economically important traits covered by these differentially selected regions. In brief, our present study provides a comprehensive genome-wide view on within- and between-breed genetic differentiation, biodiversity, and evolution among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds. These results may provide new guidance for the synthesis of new breeds with different breeding objectives.

Occurrence of congenital disorders in Swiss sheep

BACKGROUND

The rates of congenital disorders in Swiss sheep were determined by a questionnaire which was sent to 3,183 members of the Swiss Sheep Breeders' Association.

FINDINGS

A total of 993 questionnaires were returned, giving a response rate of 31.2%. Of these, 862 questionnaires originated from farms keeping one of the predominant Swiss sheep breeds: Swiss White Alpine sheep, Brown-Headed Meat sheep, Swiss Black Brown Mountain sheep and Valais Blacknose sheep. During a 10-year-period, entropion was reported in 33.6% of the farms, brachygnathia inferior in 29.5%, abdominal/umbilical hernia in 15.9%, cryptorchidism in 10.5% and torticollis in 10.5%. The most significant difference between the four breeds (P<0.001) occurred for entropion in Swiss White Alpine sheep and Brown-Headed Meat sheep, brachygnathia inferior in Swiss Black Brown Mountain sheep, and scrotal/inguinal hernia in Valais Blacknose sheep. The Swiss White Alpine breed showed a significantly higher animal prevalence of entropion (6.2% in 2011 and 5.5% in 2012) than other breeds (P<0.001).

CONCLUSIONS

These findings indicate a breed-specific necessity for action, particularly regarding Swiss animal welfare legislation, especially entropion in Swiss White Alpine sheep is concerned. In general, careful selection of breeding stock is to be recommended.

Genome-wide association and regional heritability mapping to identify loci underlying variation in nematode resistance and body weight in Scottish Blackface lambs

The genetic architecture underlying nematode resistance and body weight in Blackface lambs was evaluated comparing genome-wide association (GWA) and regional heritability mapping (RHM) approaches. The traits analysed were faecal egg count (FEC) and immunoglobulin A activity against third-stage larvae from Teladorsagia circumcincta, as indicators of nematode resistance, and body weight in a population of 752 Scottish Blackface lambs, genotyped with the 50k single-nucleotide polymorphism (SNP) chip. FEC for both Nematodirus and Strongyles nematodes (excluding Nematodirus), as well as body weight were collected at approximately 16, 20 and 24 weeks of age. In addition, a weighted average animal effect was estimated for both FEC and body weight traits. After quality control, 44 388 SNPs were available for the GWA analysis and 42 841 for the RHM, which utilises only mapped SNPs. The same fixed effects were used in both analyses: sex, year, management group, litter size and age of dam, with day of birth as covariate. Some genomic regions of interest for both nematode resistance and body weight traits were identified, using both GWA and RHM approaches. For both methods, strong evidence for association was found on chromosome 14 for Nematodirus average animal effect, chromosome 6 for Strongyles FEC at 16 weeks and chromosome 6 for body weight at 16 weeks. Across the entire data set, RHM identified more regions reaching the suggestive level than GWA, suggesting that RHM is capable of capturing some of the variation not detected by GWA analyses.

Genetic diversity and investigation of polledness in divergent goat populations using 52 088 SNPs

The recent availability of a genome-wide SNP array for the goat genome dramatically increases the power to investigate aspects of genetic diversity and to conduct genome-wide association studies in this important domestic species. We collected and analysed genotypes from 52 088 SNPs in Boer, Cashmere and Rangeland goats that had both polled and horned individuals. Principal components analysis revealed a clear genetic division between animals for each population, and model-based clustering successfully detected evidence of admixture that matched aspects of their recorded history. For example, shared co-ancestry was detected, suggesting Boer goats have been introgressed into the Rangeland population. Further, allele frequency data successfully tracked the altered genetic profile that has taken place after 40 years of breeding Australian Cashmere goats using the Rangeland animals as the founding population. Genome-wide association mapping of the POLL locus revealed a strong signal on goat chromosome 1. The 769-kb critical interval contained the polled intersex syndrome locus, confirming the genetic basis in non-European animals is the same as identified previously in Saanen goats. Interestingly, analysis of the haplotypes carried by a small set of sex-reversed animals, known to be associated with polledness, revealed some animals carried the wild-type chromosome associated with the presence of horns. This suggests a more complex basis for the relationship between polledness and the intersex condition than initially thought while validating the application of the goat SNP50 BeadChip for fine-mapping traits in goat.

GWA analysis for milk production traits in dairy sheep and genetic support for a QTN influencing milk protein percentage in the LALBA gene

In this study, we used the Illumina OvineSNP50 BeadChip to conduct a genome-wide association (GWA) analysis for milk production traits in dairy sheep by analyzing a commercial population of Spanish Churra sheep. The studied population consisted of a total of 1,681 Churra ewes belonging to 16 half-sib families with available records for milk yield (MY), milk protein and fat yields (PY and FY) and milk protein and fat contents (PP and FP). The most significant association identified reached experiment-wise significance for PP and FP and was located on chromosome 3 (OAR3). These results confirm the population-level segregation of a previously reported QTL affecting PP and suggest that this QTL has a significant pleiotropic effect on FP. Further associations were detected at the chromosome-wise significance level on 14 other chromosomal regions. The marker on OAR3 showing the highest significant association was located at the third intron of the alpha-lactalbumin (LALBA) gene, which is a functional and positional candidate underlying this association. Sequencing this gene in the 16 Churra rams of the studied resource population identified additional polymorphisms. One out of the 31 polymorphisms identified was located within the coding gene sequence (LALBA_g.242T>C) and was predicted to cause an amino acid change in the protein (Val27Ala). Different approaches, including GWA analysis, a combined linkage and linkage disequilibrium study and a concordance test with the QTL segregating status of the sires, were utilized to assess the role of this mutation as a putative QTN for the genetic effects detected on OAR3. Our results strongly support the polymorphism LALBA_g.242T>C as the most likely causal mutation of the studied OAR3 QTL affecting PP and FP, although we cannot rule out the possibility that this SNP is in perfect linkage disequilibrium with the true causal polymorphism.

Transcriptome analysis and SSR/SNP markers information of the blunt snout bream (Megalobrama amblycephala)

Background

Blunt snout bream (Megalobrama amblycephala) is an herbivorous freshwater fish species native to China and has been recognized as a main aquaculture species in the Chinese freshwater polyculture system with high economic value. Right now, only limited EST resources were available for M. amblycephala. Recent advances in large-scale RNA sequencing provide a fast, cost-effective, and reliable approach to generate large expression datasets for functional genomic analysis, which is especially suitable for non-model species with un-sequenced genomes.

Methodology and Principal Findings

Using 454 pyrosequencing, a total of 1,409,706 high quality reads (total length 577 Mbp) were generated from the normalized cDNA of pooled M. amblycephala individuals. These sequences were assembled into 26,802 contigs and 73,675 singletons. After BLAST searches against the NCBI non-redundant (NR) and UniProt databases with an arbitrary expectation value of E⁻¹⁰, over 40,000 unigenes were functionally annotated and classified using the FunCat functional annotation scheme. A comparative genomics approach revealed a substantial proportion of genes expressed in M. amblycephala tanscriptome to be shared across the genomes of zebrafish, medaka, tetraodon, fugu, stickleback, human, mouse, and chicken, and identified a substantial number of potentially novel M. amblycephala genes. A total number of 4,952 SSRs were found and 116 polymorphic loci have been characterized. A significant number of SNPs (25,697) and indels (23,287) were identified based on specific filter criteria in the M. amblycephala.

Conclusions

This study is the first comprehensive transcriptome analysis for a fish species belonging to the genus Megalobrama. These large EST resources are expected to be valuable for the development of molecular markers, construction of gene-based linkage map, and large-scale expression analysis of M. amblycephala, as well as comparative genome analysis for the genus Megalobrama fish species. The identified SSR and SNP markers will greatly benefit its breeding program and whole genome association studies.

Brachygnathia, cardiomegaly and renal hypoplasia syndrome (BCRHS) in Merino sheep maps to a 1.1-megabase region on ovine chromosome OAR2

A genome scan was conducted to map the autosomal recessive lethal disorder brachygnathia, cardiomegaly and renal hypoplasia syndrome (BCRHS) in Poll Merino sheep. The scan involved 10 affected and 27 unaffected animals from a single Poll Merino/Merino sheep flock, which were genotyped with the Illumina Ovine SNP50 BeadChip. Association and homozygosity mapping analyses located the disorder in a region comprising 20 consecutive SNPs spanning 1.1 Mb towards the distal end of chromosome OAR2. All affected animals and none of the unaffected animals were homozygous for the associated haplotype in this region. These results provide the basis for identifying the causative mutation(s) and should enable the development of a DNA test to identify carriers in the Poll Merino sheep population. Understanding the molecular control of BCRHS may provide insight into the fundamental genetic control and regulation of the affected organ systems.

Genetic testing for phenotype-causing variants in sheep and goats

This review gives an overview on ovine and caprine defects/disorders, disease predispositions, production traits and coat colours for which causal gene variants are known. Most phenotypes are inherited autosomal-recessive or dominant and in the majority are caused by single nucleotide substitutions or deletions. Causative sequence variants mainly were identified by sequencing candidate genes in the past, and recently also by whole genome analysis using the ovine 50k SNP chip. While PCR-fragment length polymorphism analyses were developed for the majority of causative sequence variants, other low- to medium-throughput PCR-based methods as PCR-single strand conformation analysis and allele-specific PCR were also established frequently. For processing large sample numbers, high-throughput methods as MALDI-ToF MS or real-time PCR are available for some gene variants. Further progress in development of ovine and caprine genome sequences and SNP chips will be beneficial for the discovery of additional causative variants in these two species.

Genomic scan of selective sweeps in thin and fat tail sheep breeds for identifying of candidate regions associated with fat deposition

BACKGROUND

Identification of genomic regions that have been targets of selection for phenotypic traits is one of the most important and challenging areas of research in animal genetics. However, currently there are relatively few genomic regions identified that have been subject to positive selection. In this study, a genome-wide scan using ~50,000 Single Nucleotide Polymorphisms (SNPs) was performed in an attempt to identify genomic regions associated with fat deposition in fat-tail breeds. This trait and its modification are very important in those countries grazing these breeds.

RESULTS

Two independent experiments using either Iranian or Ovine HapMap genotyping data contrasted thin and fat tail breeds. Population differentiation using FST in Iranian thin and fat tail breeds revealed seven genomic regions. Almost all of these regions overlapped with QTLs that had previously been identified as affecting fat and carcass yield traits in beef and dairy cattle. Study of selection sweep signatures using FST in thin and fat tail breeds sampled from the Ovine HapMap project confirmed three of these regions located on Chromosomes 5, 7 and X. We found increased homozygosity in these regions in favour of fat tail breeds on chromosome 5 and X and in favour of thin tail breeds on chromosome 7.

CONCLUSIONS

In this study, we were able to identify three novel regions associated with fat deposition in thin and fat tail sheep breeds. Two of these were associated with an increase of homozygosity in the fat tail breeds which would be consistent with selection for mutations affecting fat tail size several thousand years after domestication.

Lethal genetic disorder in Poll Merino/Merino sheep in Australia

OBJECTIVES

Characterise a lethal genetic disorder in Poll Merino/Merino sheep

DESIGN

Pathological description of a new congenital multisystem disorder in a commercial sheep flock, and analysis of breeding data collected each lambing season between 2004 and mid-lambing season 2010.

PROCEDURE

Necropsies were conducted on six affected lambs and the mode of inheritance of the disorder was determined by pedigree and segregation analyses.

RESULTS

The affected lambs were dwarfs with multiple defects in several organs, including skeleton, heart, liver and kidneys. The disorder has been named brachygnathia, cardiomegaly and renal hypoplasia syndrome (BCRHS). Segregation analysis suggests the disorder is transmitted as an autosomal trait with a recessive mode of inheritance. An annual incidence of the disorder in the discovery flock of up to 2.5% was recorded.

CONCLUSIONS

As a lethal disorder, the occurrence of BCRHS raises potential ethical and economic concerns for Merino breeders. The development of a DNA test would be useful to investigate its distribution in the Australian wool-sheep population. As the disorder affects both the skeleton and several critical organs, including the heart, it may provide a potential animal model for investigating key developmental processes in humans and other animals.

Using regulatory and epistatic networks to extend the findings of a genome scan: identifying the gene drivers of pigmentation in merino sheep

Extending genome wide association analysis by the inclusion of gene expression data may assist in the dissection of complex traits. We examined piebald, a pigmentation phenotype in both human and Merino sheep, by analysing multiple data types using a systems approach. First, a case control analysis of 49,034 ovine SNP was performed which confirmed a multigenic basis for the condition. We combined these results with gene expression data from five tissue types analysed with a skin-specific microarray. Promoter sequence analysis of differentially expressed genes allowed us to reverse-engineer a regulatory network. Likewise, by testing two-loci models derived from all pair-wise comparisons across piebald-associated SNP, we generated an epistatic network. At the intersection of both networks, we identified thirteen genes with insulin-like growth factor binding protein 7 (IGFBP7), platelet-derived growth factor alpha (PDGFRA) and the tetraspanin platelet activator CD9 at the kernel of the intersection. Further, we report a number of differentially expressed genes in regions containing highly associated SNP including ATRN, DOCK7, FGFR1OP, GLI3, SILV and TBX15. The application of network theory facilitated co-analysis of genetic variation with gene expression, recapitulated aspects of the known molecular biology of skin pigmentation and provided insights into the transcription regulation and epistatic interactions involved in piebald Merino sheep.

A whole-genome association study for pig reproductive traits

A whole-genome association study was performed for reproductive traits in commercial sows using the PorcineSNP60 BeadChip and Bayesian statistical methods. The traits included total number born (TNB), number born alive (NBA), number of stillborn (SB), number of mummified foetuses at birth (MUM) and gestation length (GL) in each of the first three parities. We report the associations of informative QTL and the genes within the QTL for each reproductive trait in different parities. These results provide evidence of gene effects having temporal impacts on reproductive traits in different parities. Many QTL identified in this study are new for pig reproductive traits. Around 48% of total genes located in the identified QTL regions were predicted to be involved in placental functions. The genomic regions containing genes important for foetal developmental (e.g. MEF2C) and uterine functions (e.g. PLSCR4) were associated with TNB and NBA in the first two parities. Similarly, QTL in other foetal developmental (e.g. HNRNPD and AHR) and placental (e.g. RELL1 and CD96) genes were associated with SB and MUM in different parities. The QTL with genes related to utero-placental blood flow (e.g. VEGFA) and hematopoiesis (e.g. MAFB) were associated with GL differences among sows in this population. Pathway analyses using genes within QTL identified some modest underlying biological pathways, which are interesting candidates (e.g. the nucleotide metabolism pathway for SB) for pig reproductive traits in different parities. Further validation studies on large populations are warranted to improve our understanding of the complex genetic architecture for pig reproductive traits.

Generation of genome-scale gene-associated SNPs in catfish for the construction of a high-density SNP array

Background

Single nucleotide polymorphisms (SNPs) have become the marker of choice for genome-wide association studies. In order to provide the best genome coverage for the analysis of performance and production traits, a large number of relatively evenly distributed SNPs are needed. Gene-associated SNPs may fulfill these requirements of large numbers and genome wide distribution. In addition, gene-associated SNPs could themselves be causative SNPs for traits. The objective of this project was to identify large numbers of gene-associated SNPs using high-throughput next generation sequencing.

Results

Transcriptome sequencing was conducted for channel catfish and blue catfish using Illumina next generation sequencing technology. Approximately 220 million reads (15.6 Gb) for channel catfish and 280 million reads (19.6 Gb) for blue catfish were obtained by sequencing gene transcripts derived from various tissues of multiple individuals from a diverse genetic background. A total of over 35 billion base pairs of expressed short read sequences were generated. Over two million putative SNPs were identified from channel catfish and almost 2.5 million putative SNPs were identified from blue catfish. Of these putative SNPs, a set of filtered SNPs were identified including 342,104 intra-specific SNPs for channel catfish, 366,269 intra-specific SNPs for blue catfish, and 420,727 inter-specific SNPs between channel catfish and blue catfish. These filtered SNPs are distributed within 16,562 unique genes in channel catfish and 17,423 unique genes in blue catfish.

Conclusions

For aquaculture species, transcriptome analysis of pooled RNA samples from multiple individuals using Illumina sequencing technology is both technically efficient and cost-effective for generating expressed sequences. Such an approach is most effective when coupled to existing EST resources generated using traditional sequencing approaches because the reference ESTs facilitate effective assembly of the expressed short reads. When multiple individuals with different genetic backgrounds are used, RNA-Seq is very effective for the identification of SNPs. The SNPs identified in this report will provide a much needed resource for genetic studies in catfish and will contribute to the development of a high-density SNP array. Validation and testing of these SNPs using SNP arrays will form the material basis for genome association studies and whole genome-based selection in catfish.