Positional Candidate Genes Research Articles

Molecular basis for a new bovine model of Niemann-Pick type C disease.

Niemann-Pick type C disease is a lysosomal storage disease affecting primarily the nervous system that results in premature death. Here we present the first report and investigation of Niemann-Pick type C disease in Australian Angus/Angus-cross calves. After a preliminary diagnosis of Niemann-Pick type C, samples from two affected calves and two obligate carriers were analysed using single nucleotide polymorphism genotyping and homozygosity mapping, and NPC1 was considered as a positional candidate gene. A likely causal missense variant on chromosome 24 in the NPC1 gene (NM_174758.2:c.2969C>G) was identified by Sanger sequencing of cDNA. SIFT analysis, protein alignment and protein modelling predicted the variant to be deleterious to protein function. Segregation of the variant with disease was confirmed in two additional affected calves and two obligate carrier dams. Genotyping of 403 animals from the original herd identified an estimated allele frequency of 3.5%. The Niemann-Pick type C phenotype was additionally confirmed via biochemical analysis of Lysotracker Green, cholesterol, sphingosine and glycosphingolipids in fibroblast cell cultures originating from two affected calves. The identification of a novel missense variant for Niemann-Pick type C disease in Angus/Angus-cross cattle will enable improved breeding and management of this disease in at-risk populations. The results from this study offer a unique opportunity to further the knowledge of human Niemann-Pick type C disease through the potential availability of a bovine model of disease.

Molecular basis of a new ovine model for human 3M syndrome-2

BackgroundBrachygnathia, cardiomegaly and renal hypoplasia syndrome (BCRHS, OMIA 001595–9940) is a previously reported recessively inherited disorder in Australian Poll Merino/Merino sheep. Affected lambs are stillborn with various congenital defects as reflected in the name of the disease, as well as short stature, a short and broad cranium, a small thoracic cavity, thin ribs and brachysternum. The BCRHS phenotype shows similarity to certain human short stature syndromes, in particular the human 3M syndrome-2. Here we report the identification of a likely disease-causing variant and propose an ovine model for human 3M syndrome-2.ResultsEight positional candidate genes were identified among the 39 genes in the approximately 1 Mb interval to which the disease was mapped previously. Obscurin like cytoskeletal adaptor 1 (OBSL1) was selected as a strong positional candidate gene based on gene function and the resulting phenotypes observed in humans with mutations in this gene. Whole genome sequencing of an affected lamb (BCRHS3) identified a likely causal variant ENSOARG00000020239:g.220472248delC within OBSL1. Sanger sequencing of seven affected, six obligate carrier, two phenotypically unaffected animals from the original flock and one unrelated control animal validated the variant. A genotyping assay was developed to genotype 583 animals from the original flock, giving an estimated allele frequency of 5%.ConclusionsThe identification of a likely disease-causing variant resulting in a frameshift (p.(Val573Trpfs*119)) in the OBSL1 protein has enabled improved breeding management of the implicated flock. The opportunity for an ovine model for human 3M syndrome and ensuing therapeutic research is promising given the availability of carrier ram semen for BCRHS.

Single-step genome-wide association study for social genetic effects and direct genetic effects on growth in Landrace pigs

In livestock social interactions, social genetic effects (SGE) represent associations between phenotype of one individual and genotype of another. Such associations occur when the trait of interest is affected by transmissible phenotypes of social partners. The aim of this study was to estimate SGE and direct genetic effects (DGE, genetic effects of an individual on its own phenotype) on average daily gain (ADG) in Landrace pigs, and to conduct single-step genome-wide association study using SGE and DGE as dependent variables to identify quantitative trait loci (QTLs) and their positional candidate genes. A total of 1,041 Landrace pigs were genotyped using the Porcine SNP 60K BeadChip. Estimates of the two effects were obtained using an extended animal model. The SGE contributed 16% of the total heritable variation of ADG. The total heritability estimated by the extended animal model including both SGE and DGE was 0.52. The single-step genome-wide association study identified a total of 23 QTL windows for the SGE on ADG distributed across three chromosomes (i.e., SSC1, SSC2, and SSC6). Positional candidate genes within these QTL regions included PRDM13, MAP3K7, CNR1, HTR1E, IL4, IL5, IL13, KIF3A, EFHD2, SLC38A7, mTOR, CNOT1, PLCB2, GABRR1, and GABRR2, which have biological roles in neuropsychiatric processes. The results of biological pathway and gene network analyses also support the association of the neuropsychiatric processes with SGE on ADG in pigs. Additionally, a total of 11 QTL windows for DGE on ADG in SSC2, 3, 6, 9, 10, 12, 14, 16, and 17 were detected with positional candidate genes such as ARL15. We found a putative pleotropic QTL for both SGE and DGE on ADG on SSC6. Our results in this study provide important insights that can help facilitate a better understanding of the molecular basis of SGE for socially affected traits.

A genome-wide association study for clinical mastitis in the dual-purpose German Black Pied cattle breed

The dual-purpose German Black Pied Cattle (DSN) has become an endangered breed of approximately 2,550 registered cows in Germany. The breed is genetically related to Holstein-Friesian cattle because the old DSN breed contributed to the selection of the modern Holstein dairy cow. In dairy farms, breeders aim to improve animal health and well-being by reducing the number of mastitis cases, which would also reduce milk losses and treatment costs. On the genomic level, no markers associated with clinical mastitis have been reported in DSN. Therefore, we performed a genome-wide association study on 1,062 DSN cows using a univariate linear mixed model that included a relatedness matrix to correct for population stratification. Although the statistical power was limited by the small population size, 3 markers were significantly associated, and 2 additional markers showed a suggestive association with clinical mastitis. Those markers accounted for 1 to 3% of the variance of clinical mastitis in the examined DSN population. One marker was found in the intragenic region of NEURL1 on BTA26, and the other 4 markers in intergenic regions on BTA3, BTA6, and BTA9. Further analyses identified 23 positional candidate genes. Among them is BMPR1B, which has been previously associated with clinical mastitis in other dairy cattle breeds. The markers presented here can be used for selection for mastitis-resistant animals in the endangered DSN population, and can broadly contribute to a better understanding of mastitis determinants in dairy cattle breeds.

Characterization of Celiac Disease-Related Epitopes and Gluten Fractions, and Identification of Associated Loci in Durum Wheat

While durum wheat is a major food source in Mediterranean countries, storage (i.e., gluten) proteins are however responsible for celiac disease (CD), a serious autoimmune disease that occurs in genetically predisposed subjects. Different gluten epitopes—defined as “immunogenic” (IP) and “toxic” (TP) peptides—are involved in the pathology and their content in wheat grain depends on environmental and genetic factors. Detection of IP and TP is not trivial, and no work has been conducted so far to identify the genomic regions associated with their accumulation in wheat. In the present study, a genome-wide association study was performed on a durum wheat collection to identify marker–trait associations (MTAs) between 5730 high quality SNPs and the accumulation of CD-related peptides and gluten protein composition measured in two consecutive cropping seasons (2015/2016 and 2016/2017). High-molecular-weight glutenin subunits (HMW-GS) were more stable between the two years, and differences in total gluten proteins were mainly due to low-molecular-weight glutenin subunits (LMW-GS) and accumulation of gliadins. In the first instance, association tests were conducted on yellow pigment content (YP), a highly inheritable trait with a well-known genetic basis, and several significant MTAs were found corresponding to loci already known for being related to YP. These findings showed that MTAs found for the rest of the measured traits were reliable. In total, 28 significant MTAs were found for gluten composition, while 14 were found to be associated with IP and TP. Noteworthy, neither significant (−log10p > 4.7) nor suggestive (−log10p > 3.3) MTAs for the accumulation of CD-triggering epitopes were found on Gli-A1/Glu-A3 and Gli-B1/Glu-B3 loci, thus suggesting regulatory rather than structural gene effect. A PBF transcription factor on chromosome 5B, known to be involved in the regulation of the expression of CD-related peptides, was identified among the positional candidate genes in the LD-decay range around significant SNPs. Results obtained in the present study provide useful insights and resources for the long-term objective of selecting low-toxic durum wheat varieties while maintaining satisfactory gluten quality.

Genome‐wide association analysis of resistance to Pythium ultimum in common bean (Phaseolus vulgaris)

AbstractPythium root rot (PRR) caused by Pythium spp. is an important root rot disease affecting common bean productivity. The objective of this study was to conduct a genome‐wide association analysis of resistance to PRR in the common bean of Andean gene pool using single nucleotide polymorphism (SNP) markers. About 260 genotypes of the Andean diversity panel (ADP) were evaluated under screen house conditions using Pythium ultimum isolate MS61 in Uganda. Sixteen significant signals for resistance to PRR were detected on chromosomes Pv01, Pv02, Pv04, Pv05 and Pv09 using 260K GBS‐based and 6K Beadchip SNPs. Common significant signals were detected on Pv02 and Pv09 for PRR. Positional candidate genes associated with significant SNPs on Pv02 were Phvul.002G119700, 16.97 kb near marker S02_25507837 (25.50 Mb), encoding Subtilase family protein, and Phvul.002G278400 near marker ss715645959 (44.79 Mb) encoding Defensin‐like (DEFL) protein involved in plant defence responses. Based on the relatively high heritability estimates observed for PRR in this study, significant SNP markers associated with genomic regions for resistance to PRR could be validated for marker‐assisted breeding in Andean beans.

Genetic Basis of Response of Ghanaian Local Chickens to Infection With a Lentogenic Newcastle Disease Virus.

Newcastle disease (ND) is a global threat to domestic poultry, especially in rural areas of Africa and Asia, where the loss of entire backyard local chicken flocks often threatens household food security and income. To investigate the genetics of Ghanaian local chicken ecotypes to Newcastle disease virus (NDV), in this study, three popular Ghanaian chicken ecotypes (regional populations) were challenged with a lentogenic NDV strain at 28 days of age. This study was conducted in parallel with a similar study that used three popular Tanzanian local chicken ecotypes and after two companion studies in the United States, using Hy-line Brown commercial laying birds. In addition to growth rate, NDV response traits were measured following infection, including anti-NDV antibody levels [pre-infection and 10 days post-infection (dpi)], and viral load (2 and 6 dpi). Genetic parameters were estimated, and two genome-wide association study analysis methods were used on data from 1,440 Ghanaian chickens that were genotyped on a chicken 600K Single Nucleotide Polymorphism (SNP) chip. Both Ghana and Tanzania NDV challenge studies revealed moderate to high (0.18 – 0.55) estimates of heritability for all traits, except viral clearance where the heritability estimate was not different from zero for the Tanzanian ecotypes. For the Ghana study, 12 quantitative trait loci (QTL) for growth and/or response to NDV from single-SNP analyses and 20 genomic regions that explained more than 1% of genetic variance using the Bayes B method were identified. Seven of these windows were also identified as having at least one significant SNP in the single SNP analyses for growth rate, anti-NDV antibody levels, and viral load at 2 and 6 dpi. An important gene for growth during stress, CHORDC1 associated with post-infection growth rate was identified as a positional candidate gene, as well as other immune related genes, including VAV2, IL12B, DUSP1, and IL17B. The QTL identified in the Ghana study did not overlap with those identified in the Tanzania study. However, both studies revealed QTL with genes vital for growth and immune response during NDV challenge. The Tanzania parallel study revealed an overlapping QTL on chromosome 24 for viral load at 6 dpi with the US NDV study in which birds were challenged with NDV under heat stress. This QTL region includes genes related to immune response, including TIRAP, ETS1, and KIRREL3. The moderate to high estimates of heritability and the identified QTL suggest that host response to NDV of local African chicken ecotypes can be improved through selective breeding to enhance increased NDV resistance and vaccine efficacy.

Deciphering the mode of action and position of genetic variants impacting on egg number in broiler breeders

BackgroundAim of the present study was first to identify genetic variants associated with egg number (EN) in female broilers, second to describe the mode of their gene action (additive and/or dominant) and third to provide a list with implicated candidate genes for the trait. A number of 2586 female broilers genotyped with the high density (~ 600 k) SNP array and with records on EN (mean = 132.4 eggs, SD = 29.8 eggs) were used. Data were analyzed with application of additive and dominant multi-locus mixed models.ResultsA number of 7 additive, 4 dominant and 6 additive plus dominant marker-trait significant associations were detected. A total number of 57 positional candidate genes were detected within 50 kb downstream and upstream flanking regions of the 17 significant markers. Functional enrichment analysis pinpointed two genes (BHLHE40 and CRTC1) to be involved in the ‘entrainment of circadian clock by photoperiod’ biological process. Gene prioritization analysis of the positional candidate genes identified 10 top ranked genes (GDF15, BHLHE40, JUND, GDF3, COMP, ITPR1, ELF3, ELL, CRLF1 and IFI30). Seven prioritized genes (GDF15, BHLHE40, JUND, GDF3, COMP, ELF3, CRTC1) have documented functional relevance to reproduction, while two more prioritized genes (ITPR1 and ELL) are reported to be related to egg quality in chickens.ConclusionsPresent results have shown that detailed exploration of phenotype-marker associations can disclose the mode of action of genetic variants and help in identifying causative genes associated with reproductive traits in the species.

Association Analysis between SPP1, POFUT1 and PRLR Gene Variation and Milk Yield, Composition and Coagulation Traits in Sarda Sheep.

Simple SummaryThe purpose of our research was to analyze the association between the three candidate genes secreted phosphoprotein 1 (SPP1), protein O-fucosyltransferase 1 (POFUT1) and prolactin receptor (PRLR) with milk production, quality and coagulation properties in 380 Sarda dairy sheep. Results revealed an association between SPP1 and somatic cells count, in line with the function of this gene and with its genomic position. We revealed an association of POFUT1 variation with milk coagulation properties, and PRLR with quality. This information can be useful for future breeding schemes in sheep.Many studies focus on the identification of genomic regions that undergo selective processes, where evidence of selection is revealed and positional candidate genes are identified. The aim of the research was to evaluate the association between positional candidate genes, namely secreted phosphoprotein 1 (SPP1, sheep chromosome Ovis aries OAR6, 36.651–36.658 Mb), protein O-fucosyltransferase 1 (POFUT1, OAR13, 61.006–61.027 Mb) and prolactin receptor (PRLR, OAR16, 38.969–39.028 Mb) with milk yield, composition and coagulation traits. Eight single nucleotide polymorphisms (SNPs) mapping to the three genes were genotyped in 380 Sarda dairy sheep. Statistical analysis revealed an association between SNP rs161844011 at SPP1 (chromosome position Oar_v3 OAR6:36651870, gene region exon 7) and somatic cell score, while POFUT1 SNP rs424501869 (OAR13:61007495, intron 1) was associated with curd firmness both 45 and 60 min after rennet addition (p = 0.015 and p = 0.007, respectively). SNP rs400874750 at PRLR gene (OAR16:39004070, intron 2) had a significant association with lactose content (p = 0.020), somatic cell score (p = 0.038), rennet coagulation time (p = 0.018) and curd firming time (p = 0.047). The outcome of this research confirmed predictions based on genomic studies, producing new information regarding the SPP1, POFUT1 and PRLR genes, which may be useful for future breeding schemes.

Identification of Loci and Pathways Associated with Heifer Conception Rate in U.S. Holsteins.

Heifer conception rate (HCR) is defined as the percentage of inseminated heifers that become pregnant at each service. The genome-wide association analyses in this study focused on identifying the loci associated with Holstein heifer (n = 2013) conception rate at first service (HCR1) and the number of times bred (TBRD) to achieve a pregnancy. There were 348 unique loci associated (p < 5 × 10−8) with HCR1 and 615 unique loci associated (p < 5 × 10−8) with TBRD. The two phenotypes shared 302 loci, and 56 loci were validated in independent cattle populations. There were 52 transcription factor binding sites (TFBS) and 552 positional candidate genes identified in the HCR1- and TBRD-associated loci. The positional candidate genes and the TFBS associated with HCR1 and TBRD were used in the ingenuity pathway analysis (IPA). In the IPA, 11 pathways, 207 master regulators and 11 upstream regulators were associated (p < 1.23 × 10−5) with HCR1 and TBRD. The validated loci associated with both HCR1 and TBRD make good candidates for genomic selection and further investigations to elucidate the mechanisms associated with subfertility and infertility.

The morbid cutaneous anatomy of the human genome revealed by a bioinformatic approach

Computational approaches have been developed to prioritize candidate genes in disease gene identification. They are based on different pieces of evidences associating each gene with the given disease. In this study, 648 genes underlying genodermatoses have been compared to 1808 genes involved in other genetic diseases using a bioinformatic approach. These genes were studied at the structural, evolutionary and functional levels. Results show that genes underlying genodermatoses present longer CDS and have more exons. Significant differences were observed in nucleotide motif and amino-acid compositions. Evolutionary conservation analysis revealed that genodermatoses genes have less paralogs, more orthologs in Mouse and Dog and are less conserved. Functional analysis revealed that genodermatosis genes seem to be involved in immune system and skin layers. The Bayesian network model returned a rate of good classification of around 80%. This computational approach could help investigators working in the field of dermatology by prioritizing positional candidate genes for mutation screening.

MYO, a Candidate Gene for Haploid Induction in Maize Causes Male Sterility.

Doubled haploid technology is highly successful in maize breeding programs and is contingent on the ability of maize inducers to efficiently produce haploids. Knowledge of the genes involved in haploid induction is important for not only developing better maize inducers, but also to create inducers in other crops. The main quantitative trait loci involved in maize haploid induction are qhir1 and qhir8. The gene underlying qhir1 has been discovered and validated by independent research groups. Prior to initiation of this study, the gene associated with qhir8 had yet to be recognized. Therefore, this research focused on characterizing positional candidate genes underlying qhir8. Pursuing this goal, a strong candidate for qhir8, GRMZM2G435294 (MYO), was silenced by RNAi. Analysis of crosses with these heterozygous RNAi-transgenic lines for haploid induction rate revealed that the silencing of MYO significantly enhanced haploid induction rate by an average of 0.6% in the presence of qhir1. Recently, GRMZM2G465053 (ZmDMP) was identified by map-based gene isolation and shown to be responsible for qhir8. While our results suggest that MYO may contribute to haploid induction rate, results were inconsistent and only showing minor increases in haploid induction rate compared to ZmDMP. Instead, reciprocal crosses clearly revealed that the silencing of MYO causes male sterility.

Genome-wide association study of milk production traits in a crossbred dairy sheep population using three statistical models.

Milk production is one of the most important characteristics of dairy sheep, and the identification of genes affecting milk production traits is critical to understanding the genetics and improve milk production in future generations. Three statistical techniques, namely GWAS, ridge-regression BLUP and BayesC , were used to identify SNPs in significant association with three milk production traits (milk yield, fat yield and protein yield) in a crossbred dairy sheep population. The results suggested that chromosomes 1, 3, 4, 5, 7 and 11 were likely to harbor genes important to milk production because these chromosomes had the greatest top-100-SNP variance contributions on the three milk production traits. The GWAS analysis identified between 74 and 288 genome-wide significant SNP (P<0.05) whereas the BayesCπ model revealed between six and 63 SNPs, each with >95% posterior probability of inclusion as having a non-zero association effect on at least one of the three milk production traits. Positional candidate genes for milk production in sheep were searched, based on the sheep genomic assembly OAR version 3.1, such as those which map position coincided with or was located within 0.1 Mbp of a genome-wide suggestive or significant SNP. These identified SNPs and candidate genes supported some previous findings and also added new information about genetic markers for genetic improvement of lactation in dairy sheep, but keeping in mind that the majority of these positional candidate genes are not necessarily true causative loci for these traits and future validations are thus necessary.

A comprehensive analysis of AHRR gene as a candidate for cleft lip with or without cleft palate

Cleft lip and palate (CL/P) is among the most common congenital malformations and affects 1 in 700 newborns. CL/P is caused by genetic and environmental factors (maternal smoking, alcohol or drug use and others). Many genes and loci were associated with cleft lip/palate but the amount of heterogeneity justifies identifying new causal genes and variants. AHRR (Aryl-Hydrocarbon Receptor Repressor) gene has recently been related to CL/P however, few functional studies analyze the genotypephenotype interaction of AHRR with CL/P. Several studies associate the molecular pathway of AHRR to CL/P which indicates this gene as a functional candidate in CL/P etiology. Methods: Systematic Literature Review was performed using PUBMED database with the keywords cleft lip, cleft palate, orofacial cleft, AHRR and synonyms. SLR resulted in 37 included articles. Results: AHRR is a positional and functional candidate gene for CL/P. In silico analysis detected interactions with other genes previously associated to CL/P like ARNT and CYP1A1. AHRR protein regulates cellular toxicity through TCDD mediated AHR pathway. Exposure to TCDD in animal embryos is AHR mediated and lead to cleft palate due to palate fusion failure and post fusion rupture. AHRR regulates cellular growth and differentiation, fundamental to lip and palatogenesis. AHRR decreases carcinogenesis and recently a higher tumor risk has been described in CL/P patients and families. AHRR is also a smoking biomarker due to changed methylation sites found in smokers DNA although folate intake may partially revert these methylation alterations. This corroborates the role of maternal smoking and lack of folate supplementation as risk factors for CL/P. Conclusion: This research identified the importance of AHRR in dioxin response and demonstrated an example of genetic and environmental interaction, indispensable in the development of many complex diseases.

Identification of Susceptibility Loci for Spontaneous Coronary Artery Dissection

Spontaneous coronary artery dissection (SCAD), an idiopathic disorder that predominantly affects young to middle-aged women, has emerged as an important cause of acute coronary syndrome, myocardial infarction, and sudden cardiac death. To identify common single-nucleotide variants (SNVs) associated with SCAD susceptibility. This single-center genome-wide association study examined approximately 5 million genotyped and imputed SNVs and subsequent SNV-targeted replication analysis results in individuals enrolled in the Mayo Clinic SCAD registry from August 30, 2011, to August 2, 2018. Data analysis was performed from June 21, 2017, to December 30, 2019. Genetic loci and positional candidate genes associated with SCAD. This study included 484 white women with SCAD (mean [SD] age, 46.6 [9.2] years) and 1477 white female controls in the discovery cohort (mean [SD] age, 64.0 [14.5] years) and 183 white women with SCAD (mean [SD] age, 47.1 [9.9] years) and 340 white female controls in the replication cohort (mean [SD] age, 51.0 [15.3] years). Associations with SCAD risk reached genome-wide significance at 3 loci (1q21.3 [OR, 1.78; 95% CI, 1.51-2.09; P = 2.63 × 10-12], 6p24.1 [OR, 1.77; 95% CI, 1.51-2.09; P = 7.09 × 10-12], and 12q13.3 [OR, 1.67; 95% CI, 1.42-1.97; P = 3.62 × 10-10]), and 7 loci had evidence suggestive of an association (1q24.2 [OR, 2.10; 95% CI, 1.58-2.79; P = 2.88 × 10-7], 3q22.3 [OR, 1.47; 95% CI, 1.26-1.71; P = 6.65 × 10-7], 4q34.3 [OR, 1.84; 95% CI, 1.44-2.35; P = 9.80 × 10-7], 8q24.3 [OR, 2.57; 95% CI, 1.76-3.75; P = 9.65 × 10-7], 15q21.1 [OR, 1.75; 95% CI, 1.40-2.18; P = 7.23 × 10-7], 16q24.1 [OR, 1.91; 95% CI, 1.49-2.44; P = 2.56 × 10-7], and 21q22.11 [OR, 2.11; 95% CI, 1.59-2.82; P = 3.12 × 10-7]) after adjusting for the top 5 principal components. Associations were validated for 5 of the 10 risk alleles in the replication cohort. In a meta-analysis of the discovery and replication cohorts, associations for the 5 SNVs were significant, with relatively large effect sizes (1q21.3 [OR, 1.77; 95% CI, 1.54-2.03; P = 3.26 × 10-16], 6p24.1 [OR, 1.71; 95% CI, 1.49-1.97; P = 4.59 × 10-14], 12q13.3 [OR, 1.69; 95% CI, 1.47-1.94; P = 1.42 × 10-13], 15q21.1 [OR, 1.79; 95% CI, 1.48-2.17; P = 2.12 × 10-9], and 21q22.11 [OR, 2.18; 95% CI, 1.70-2.81; P = 1.09 × 10-9]). Each index SNV was within or near a gene highly expressed in arterial tissue and previously linked to SCAD (PHACTR1) and/or other vascular disorders (LRP1, LINC00310, and FBN1). This study revealed 5 replicated risk loci and positional candidate genes for SCAD, most of which are associated with extracoronary arteriopathies. Moreover, the alternate alleles of 3 SNVs have been previously associated with atherosclerotic coronary artery disease, further implicating allelic susceptibility to coronary artery atherosclerosis vs dissection.

Genomic Scan for Selection Signature Reveals Fat Deposition in Chinese Indigenous Sheep with Extreme Tail Types.

Simple SummaryAccording to the tail types, sheep can be briefly classified into three groups (fat-tailed, fat-rumped, and thin-tailed sheep). In this study, we used these three typical breeds from Chinese indigenous sheep breeds to perform a genome scan for selective sweeps using Ovine Infinium HD SNP BeadChip genotype data. Results showed that 25 genomic regions exhibited selection signals and harbored 73 positional candidate genes. These genes were documented not only to be associated with tail fat formation, but also be related to reproduction, body conformation, and appearance. Our findings contributed to understanding genetic basis of fat deposition in sheep tail and provide a reference for developing new sheep breeds with an ideal tail type.It is a unique feature that fat can be deposited in sheep tails and rumps. To elucidate the genetic mechanism underlying this trait, we collected 120 individuals from three Chinese indigenous sheep breeds with extreme tail types, namely large fat-tailed sheep (n = 40), Altay sheep (n = 40), and Tibetan sheep (n = 40), and genotyped them using the Ovine Infinium HD SNP BeadChip. Then genomic scan for selection signatures was performed using the hapFLK. In total, we identified 25 genomic regions exhibiting evidence of having been under selection. Bioinformatic analysis of the genomic regions showed that selection signatures related to multiple candidate genes had a demonstrated role in phenotypic variation. Nine genes have documented association with sheep tail types, including WDR92, TBX12, WARS2, BMP2, VEGFA, PDGFD, HOXA10, ALX4, and ETAA1. Moreover, a number of genes were of particular interest, including RXFP2 associated with the presence/absence and morphology of horns; MITF involved in coat color; LIN52 and SYNDIG1L related to the number of teats; MSRB3 gene associated with ear sizes; LTBP2 considered as a positional candidate genes for number of ribs; JAZF1 regulating lipid metabolism; PGRMC2, SPAG17, TSHR, GTF2A1, and LARP1B implicated with reproductive traits. Our findings provide insights into fat tail formation and a reference for carrying out molecular breeding and conservation in sheep.

Use of gene expression profile to identify potentially relevant transcripts to myofibrillar fragmentation index trait.

The myofibrillar fragmentation index (MFI) is an indicative trait for meat tenderness. Longissimus thoracis muscle samples from the 20 most extreme bulls (out of 80 bulls set) for MFI (high (n = 10) and low (n = 10) groups) trait were used to perform transcriptomic analysis, using RNA Sequencing (RNA-Seq). An average of 24.616 genes was expressed in the Nellore muscle transcriptome analysis. A total of 96 genes were differentially expressed (p value ≤ 0.001) between the two groups of divergent bulls for MFI. The HEBP2 and BDH1 genes were overexpressed in animals with high MFI. The MYBPH and MYL6, myosin encoders, were identified. The differentially expressed genes were related to increase mitochondria efficiency, especially in cells under oxidative stress conditions, and these also were related to zinc and calcium binding, membrane transport, and muscle constituent proteins, such as actin and myosin. Most of those genes were involved in metabolic pathways of oxidation-reduction, transport of lactate in the plasma membrane, and muscle contraction. This is the first study applying MFI phenotypes in transcriptomic studies to identify and understand differentially expressed genes for beef tenderness. These results suggest that differences detected in gene expression between high and low MFI animals are related to reactive mechanisms and structural components of oxidative fibers under the condition of cellular stress. Some genes may be selected as positional candidate genes to beef tenderness, MYL6, MYBPH, TRIM63, TRIM55, TRIOBP, and CHRNG genes. The use of MFI phenotypes could enhance results of meat tenderness studies.

Loci associated with conception rate in crossbred beef heifers.

The inability of beef cattle to maintain full term pregnancies has become an economic concern for the beef industry. Herd management and nutritional improvements have alleviated environmental impacts on embryonic and fetal loss, yet additional gains can be made through genomic selection. The objectives of this study were to identify loci and gene-sets in crossbred beef heifers associated with the number of services required to become pregnant (TBRD) and heifer conception rate at first service (HCR1). Heifers (n = 709) from a commercial beef operation underwent one round of artificial insemination, before exposure to bulls for natural service for 50 days. Pregnancy and time of conception was determined by ultrasound 35 days after the breeding season. Heifers were genotyped using the GeneSeek (Lincoln, NE) Bovine GGP50K BeadChip prior to genome-wide association analyses (GWAA) conducted using an EIGENSTRAT-like model to identify loci associated (P < 1 × 10−5) with TBRD and HCR1. One locus was associated (P = 8.97 × 10−6) with TBRD on BTA19 and included the positional candidate gene ASIC2, which is differentially expressed in the endometrium of fertility classified heifers, and the positional candidate gene, SPACA3. Gene-set enrichment analyses using SNP (GSEA-SNP) data, was performed and identified one gene-set, oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen as enriched (NES = 3.15) with TBRD and contained nine leading edge genes that contributed to the enrichment of the gene set. The enriched gene-set is involved in catalyzing oxidation-reduction reactions, which have been associated with oxidative stressors impacting pregnancy success. No loci were associated nor gene-sets enriched with HCR1. Identification of loci, positional candidate genes, gene-sets and leading edge genes enriched for fertility facilitate genomic selection that allows producers to select for reproductively superior cattle, reduce costs associated with infertility, and increase percent calf crop.

Association of Melanophilin (MLPH) gene polymorphism with coat colour in Rex rabbits

Rex rabbit, with multiple phenotypes and colourful fur, is an interesting model for assessing the effect of coat colour gene mutations on characteristic pigmentation phenotype. Based on previous study, the &lt;em&gt;melanophilin&lt;/em&gt; (&lt;em&gt;MLPH&lt;/em&gt;) gene is a positional candidate gene related coat colour dilution. The fur colours are a lighter shade, e.g. grey instead of black. We sequenced 1689 base pairs of the &lt;em&gt;MLPH&lt;/em&gt; gene in Chinchilla and black Rex rabbit. A total of 13 polymorphisms were identified, including seven missense mutations. The rabbit &lt;em&gt;MLPH&lt;/em&gt; gene has a very high GC content and the protein shows 64.87% identity to the orthologous human protein (lack of homologous amino acids encoded by human MLPH exon 9). Hardy-Weinberg test showed that, except for the g.606C&amp;gt;A single nucleotid polymorphism (SNP), all other SNPs were in Hardy-Weinberg equilibrium. Haplotype analysis revealed that the seven missense mutation SNPs of two strains of Rex rabbits formed 10 haplotypes, but there were only seven major types of haplotypes (haplotype frequency &lt;em&gt;P&lt;/em&gt;&amp;gt;0.05). The major haplotypes of the Chinchilla and black Rex rabbits were H1/H2/H3/H4/H5 and H1/H2/H3/H6/H8, respectively. The special haplotypes of Chinchilla Rex rabbit (H4, H5, H7) were consistently associated with the Chinchilla phenotype. This study provides evidence that different coat colour formation may be caused by one or more mutations within &lt;em&gt;MLPH&lt;/em&gt; gene in several Rex rabbit strains. The data on polymorphisms that are associated with the Chinchilla phenotype facilitate the breeding of rabbits with defined coat colours.

Systems Genetics of Optic Nerve Axon Necrosis During Glaucoma.

In this study, we identify genomic regions that modulate the number of necrotic axons in optic nerves of a family of mice, some of which have severe glaucoma, and define a set of high priority positional candidate genes that modulate retinal ganglion cell (RGC) axonal degeneration. A large cohort of the BXD family were aged to greater than 13 months of age. Optic nerves from 74 strains and the DBA/2J (D2) parent were harvested, sectioned, and stained with p-phenylenediamine. Numbers of necrotic axons per optic nerve cross-section were counted from 1 to 10 replicates per genotype. Strain means and standard errors were uploaded into GeneNetwork 2 for mapping and systems genetics analyses (Trait 18614). The number of necrotic axons per nerve ranged from only a few hundred to more than 4,000. Using conventional interval mapping as well as linear mixed model mapping, we identified a single locus on chromosome 12 between 109 and 112.5 Mb with a likelihood ratio statistic (LRS) of ~18.5 (p genome-wide ~0.1). Axon necrosis is not linked to locations of major known glaucoma genes in this family, including Gpnmb, Tyrp1, Cdh11, Pou6f2, and Cacna2d1. This indicates that although these genes contribute to pigmentary dispersion or elevated IOP, none directly modulates axon necrosis. Of 156 positional candidates, eight genes—CDC42 binding protein kinase beta (Cdc42bpb); eukaryotic translation initiation factor 5 (Eif5); BCL2-associated athanogene 5 (Bag5); apoptogenic 1, mitochondrial (Apopt1); kinesin light chain 1 (Klc1); X-ray repair cross complementing 3 (Xrcc3); protein phosphatase 1, regulatory subunit 13B (Ppp1r13b); and transmembrane protein 179 (Tmem179)—passed stringent criteria and are high priority candidates. Several candidates are linked to mitochondria and/or axons, strengthening their plausible role as modulators of ON necrosis. Additional studies are required to validate and/or eliminate plausible candidates. Surprisingly, IOP and ON necrosis are inversely correlated across the BXD family in mice >13 months of age and these two traits share few genes among their top ocular and retinal correlates. These data suggest that the two traits are independently modulated or that a more complex and multifaceted approach is required to reveal their association.