Sequenced-based GWAS for linear classification traits in Belgian Blue beef cattle reveals new coding variants in genes regulating body size in mammals.

BACKGROUND: Cohorts of individuals that have been genotyped and phenotyped for genomic selection programs offer the opportunity to better understand genetic variation associated with complex traits. Here, we performed an association study for traits related to body size and muscular development in intensively selected beef cattle. We leveraged multiple trait information to refine and interpret the significant associations. RESULTS: After a multiple-step genotype imputation to the sequence-level for 14,762 Belgian Blue beef (BBB) cows, we performed a genome-wide association study (GWAS) for 11 traits related to muscular development and body size. The 37 identified genome-wide significant quantitative trait loci (QTL) could be condensed in 11 unique QTL regions based on their position. Evidence for pleiotropic effects was found in most of these regions (e.g., correlated association signals, overlap between credible sets (CS) of candidate variants). Thus, we applied a multiple-trait approach to combine information from different traits to refine the CS. In several QTL regions, we identified strong candidate genes known to be related to growth and height in other species such as LCORL-NCAPG or CCND2. For some of these genes, relevant candidate variants were identified in the CS, including three new missense variants in EZH2, PAPPA2 and ADAM12, possibly two additional coding variants in LCORL, and candidate regulatory variants linked to CCND2 and ARMC12. Strikingly, four other QTL regions associated with dimension or muscular development traits were related to five (recessive) deleterious coding variants previously identified. CONCLUSIONS: Our study further supports that a set of common genes controls body size across mammalian species. In particular, we added new genes to the list of those associated with height in both humans and cattle. We also identified new strong candidate causal variants in some of these genes, strengthening the evidence of their causality. Several breed-specific recessive deleterious variants were identified in our QTL regions, probably as a result of the extreme selection for muscular development in BBB cattle.

An organism-wide ATAC-seq peak catalog for the bovine and its use to identify regulatory variants.

We report the generation of an organism-wide catalog of 976,813 cis-acting regulatory elements for the bovine detected by the assay for transposase accessible chromatin using sequencing (ATAC-seq). We regroup these regulatory elements in 16 components by nonnegative matrix factorization. Correlation between the genome-wide density of peaks and transcription start sites, correlation between peak accessibility and expression of neighboring genes, and enrichment in transcription factor binding motifs support their regulatory potential. Using a previously established catalog of 12,736,643 variants, we show that the proportion of single-nucleotide polymorphisms mapping to ATAC-seq peaks is higher than expected and that this is owing to an approximately 1.3-fold higher mutation rate within peaks. Their site frequency spectrum indicates that variants in ATAC-seq peaks are subject to purifying selection. We generate eQTL data sets for liver and blood and show that variants that drive eQTL fall into liver- and blood-specific ATAC-seq peaks more often than expected by chance. We combine ATAC-seq and eQTL data to estimate that the proportion of regulatory variants mapping to ATAC-seq peaks is approximately one in three and that the proportion of variants mapping to ATAC-seq peaks that are regulatory is approximately one in 25. We discuss the implication of these findings on the utility of ATAC-seq information to improve the accuracy of genomic selection.

High-resolution structural variants catalogue in a large-scale whole genome sequenced bovine family cohort data.

BACKGROUND: Structural variants (SVs) are chromosomal segments that differ between genomes, such as deletions, duplications, insertions, inversions and translocations. The genomics revolution enabled the discovery of sub-microscopic SVs via array and whole-genome sequencing (WGS) data, paving the way to unravel the functional impact of SVs. Recent human expression QTL mapping studies demonstrated that SVs play a disproportionally large role in altering gene expression, underlining the importance of including SVs in genetic analyses. Therefore, this study aimed to generate and explore a high-quality bovine SV catalogue exploiting a unique cattle family cohort data (total 266 samples, forming 127 trios). RESULTS: We curated 13,731 SVs segregating in the population, consisting of 12,201 deletions, 1,509 duplications, and 21 multi-allelic CNVs (> 50-bp). Of these, we validated a subset of copy number variants (CNVs) utilising a direct genotyping approach in an independent cohort, indicating that at least 62% of the CNVs are true variants, segregating in the population. Among gene-disrupting SVs, we prioritised two likely high impact duplications, encompassing ORM1 and POPDC3 genes, respectively. Liver expression QTL mapping results revealed that these duplications are likely causing altered gene expression, confirming the functional importance of SVs. Although most of the accurately genotyped CNVs are tagged by single nucleotide polymorphisms (SNPs) ascertained in WGS data, most CNVs were not captured by individual SNPs obtained from a 50K genotyping array. CONCLUSION: We generated a high-quality SV catalogue exploiting unique whole genome sequenced bovine family cohort data. Two high impact duplications upregulating the ORM1 and POPDC3 are putative candidates for postpartum feed intake and hoof health traits, thus warranting further investigation. Generally, CNVs were in low LD with SNPs on the 50K array. Hence, it remains crucial to incorporate CNVs via means other than tagging SNPs, such as investigation of tagging haplotypes, direct imputation of CNVs, or direct genotyping as done in the current study. The SV catalogue and the custom genotyping array generated in the current study will serve as valuable resources accelerating utilisation of full spectrum of genetic variants in bovine genomes.

Comparing the transcriptome of developing native and iPSC-derived mouse retinae by single cell RNA sequencing.

We report the generation and analysis of single-cell RNA-Seq data (> 38,000 cells) from mouse native retinae and induced pluripotent stem cell (iPSC)-derived retinal organoids at four matched stages of development spanning the emergence of the major retinal cell types. We combine information from temporal sampling, visualization of 3D UMAP manifolds, pseudo-time and RNA velocity analyses, to show that iPSC-derived 3D retinal organoids broadly recapitulate the native developmental trajectories. However, we observe relaxation of spatial and temporal transcriptome control, premature emergence and dominance of photoreceptor precursor cells, and susceptibility of dynamically regulated pathways and transcription factors to culture conditions in retinal organoids. We demonstrate that genes causing human retinopathies are enriched in cell-type specifying genes and identify a subset of disease-causing genes with expression profiles that are highly conserved between human retinae and murine retinal organoids. This study provides a resource to the community that will be useful to assess and further improve protocols for ex vivo recapitulation and study of retinal development.

A 12 kb multi-allelic copy number variation encompassing a GC gene enhancer is associated with mastitis resistance in dairy cattle.

Clinical mastitis (CM) is an inflammatory disease occurring in the mammary glands of lactating cows. CM is under genetic control, and a prominent CM resistance QTL located on chromosome 6 was reported in various dairy cattle breeds. Nevertheless, the biological mechanism underpinning this QTL has been lacking. Herein, we mapped, fine-mapped, and discovered the putative causal variant underlying this CM resistance QTL in the Dutch dairy cattle population. We identified a ~12 kb multi-allelic copy number variant (CNV), that is in perfect linkage disequilibrium with a lead SNP, as a promising candidate variant. By implementing a fine-mapping and through expression QTL mapping, we showed that the group-specific component gene (GC), a gene encoding a vitamin D binding protein, is an excellent candidate causal gene for the QTL. The multiplicated alleles are associated with increased GC expression and low CM resistance. Ample evidence from functional genomics data supports the presence of an enhancer within this CNV, which would exert cis-regulatory effect on GC. We observed that strong positive selection swept the region near the CNV, and haplotypes associated with the multiplicated allele were strongly selected for. Moreover, the multiplicated allele showed pleiotropic effects for increased milk yield and reduced fertility, hinting that a shared underlying biology for these effects may revolve around the vitamin D pathway. These findings together suggest a putative causal variant of a CM resistance QTL, where a cis-regulatory element located within a CNV can alter gene expression and affect multiple economically important traits.

Global transcriptomic profiles of circulating leucocytes in early lactation cows with clinical or subclinical mastitis.

Bovine mastitis, an inflammatory disease of the mammary gland, is classified as subclinical or clinical. Circulating neutrophils are recruited to the udder to combat infection. We compared the transcriptomic profiles in circulating leukocytes between healthy cows and those with naturally occurring subclinical or clinical mastitis. Holstein Friesian dairy cows from six farms in EU countries were recruited. Based on milk somatic cell count and clinical records, cows were classified as healthy (n = 147), subclinically (n = 45) or clinically mastitic (n = 22). Circulating leukocyte RNA was sequenced with Illumina NextSeq single end reads (30 M). Differentially expressed genes (DEGs) between the groups were identified using CLC Genomics Workbench V21, followed by GO enrichment analysis. Both subclinical and clinical mastitis caused significant changes in the leukocyte transcriptome, with more intensive changes attributed to clinical mastitis. We detected 769 DEGs between clinical and healthy groups, 258 DEGs between subclinical and healthy groups and 193 DEGs between clinical and subclinical groups. Most DEGs were associated with cell killing and immune processes. Many upregulated DEGs in clinical mastitis encoded antimicrobial peptides (AZU1, BCL3, CAMP, CATHL1, CATHL2, CATHL4,CATHL5, CATHL6, CCL1, CXCL2, CXCL13, DEFB1, DEFB10, DEFB4A, DEFB7, LCN2, PGLYRP1, PRTN3, PTX3, S100A8, S100A9, S100A12, SLC11A1, TF and LTF) which were not upregulated in subclinical mastitis. The use of transcriptomic profiles has identified a much greater up-regulation of genes encoding antimicrobial peptides in circulating leukocytes of cows with naturally occurring clinical compared with subclinical mastitis. These could play a key role in combatting disease organisms.

Elevated Fibroblast Growth Factor Signaling Is Critical for the Pathogenesis of the Dwarfism in Evc2/Limbin Mutant Mice.

Ellis-van Creveld (EvC) syndrome is a skeletal dysplasia, characterized by short limbs, postaxial polydactyly, and dental abnormalities. EvC syndrome is also categorized as a ciliopathy because of ciliary localization of proteins encoded by the two causative genes, EVC and EVC2 (aka LIMBIN). While recent studies demonstrated important roles for EVC/EVC2 in Hedgehog signaling, there is still little known about the pathophysiological mechanisms underlying the skeletal dysplasia features of EvC patients, and in particular why limb development is affected, but not other aspects of organogenesis that also require Hedgehog signaling. In this report, we comprehensively analyze limb skeletogenesis in Evc2 mutant mice and in cell and tissue cultures derived from these mice. Both in vivo and in vitro data demonstrate elevated Fibroblast Growth Factor (FGF) signaling in Evc2 mutant growth plates, in addition to compromised but not abrogated Hedgehog-PTHrP feedback loop. Elevation of FGF signaling, mainly due to increased Fgf18 expression upon inactivation of Evc2 in the perichondrium, critically contributes to the pathogenesis of limb dwarfism. The limb dwarfism phenotype is partially rescued by inactivation of one allele of Fgf18 in the Evc2 mutant mice. Taken together, our data uncover a novel pathogenic mechanism to understand limb dwarfism in patients with Ellis-van Creveld syndrome.

Expression of Evc2 in craniofacial tissues and craniofacial bone defects in Evc2 knockout mouse.

OBJECTIVE: Our objectives were to determine the expression of EVC2 in craniofacial tissues and investigate the effect of Evc2 deficiency on craniofacial bones using Evc2 knockout (KO) mouse model. DESIGN: Evc2 KO mice were generated by introducing a premature stop codon followed by the Internal Ribosomal Entry Site fused to ß-galactosidase (LacZ). Samples from wild-type (WT), heterozygous (Het) and homozygous Evc2 KO mice were prepared. LacZ staining and immunohistochemistry (IHC) with anti-ß-galactosidase, anti-EVC2 and anti-SOX9 antibodies were performed. The craniofacial bones were stained with alcian blue and alizarin red. RESULTS: The LacZ activity in KO was mainly observed in the anterior parts of viscerocranium. The Evc2-expressing cells were identified in many cartilageous regions by IHC with anti-ß-galactosidase antibody in KO and Het embryos. The endogenous EVC2 protein was observed in these areas in WT embryos. Double labeling with anti-SOX9 antibody showed that these cells were mainly chondrocytes. At adult stages, the expression of EVC2 was found in chondrocytes of nasal bones and spheno-occipital synchondrosis, and osteocytes and endothelial-like cells of the premaxilla and mandible. The skeletal double staining demonstrated that craniofacial bones, where the expression of EVC2 was observed, in KO had the morphological defects as compared to WT. CONCLUSION: To our knowledge, our study was the first to identify the types of Evc2-expressing cells in craniofacial tissues. Consistent with the expression pattern, abnormal craniofacial bone morphology was found in the Evc2 KO mice, suggesting that EVC2 may be important during craniofacial growth and development.

Ellis Van Creveld2 is Required for Postnatal Craniofacial Bone Development.

Ellis-van Creveld (EvC) syndrome is a genetic disorder with mutations in either EVC or EVC2 gene. Previous case studies reported that EvC patients underwent orthodontic treatment, suggesting the presence of craniofacial bone phenotypes. To investigate whether a mutation in EVC2 gene causes a craniofacial bone phenotype, Evc2 knockout (KO) mice were generated and cephalometric analysis was performed. The heads of wild type (WT), heterozygous (Het) and homozygous Evc2 KO mice (1-, 3-, and 6-week-old) were prepared and cephalometric analysis based on the selected reference points on lateral X-ray radiographs was performed. The linear and angular bone measurements were then calculated, compared between WT, Het and KO and statistically analyzed at each time point. Our data showed that length of craniofacial bones in KO was significantly lowered by ∼20% to that of WT and Het, the growth of certain bones, including nasal bone, palatal length, and premaxilla was more affected in KO, and the reduction in these bone length was more significantly enhanced at later postnatal time points (3 and 6 weeks) than early time point (1 week). Furthermore, bone-to-bone relationship to cranial base and cranial vault in KO was remarkably changed, i.e. cranial vault and nasal bone were depressed and premaxilla and mandible were developed in a more ventral direction. Our study was the first to show the cause-effect relationship between Evc2 deficiency and craniofacial defects in EvC syndrome, demonstrating that Evc2 is required for craniofacial bone development and its deficiency leads to specific facial bone growth defect. Anat Rec, 299:1110-1120, 2016. © 2016 Wiley Periodicals, Inc.

Ectopic Expression of Retrotransposon-Derived PEG11/RTL1 Contributes to the Callipyge Muscular Hypertrophy.

The callipyge phenotype is an ovine muscular hypertrophy characterized by polar overdominance: only heterozygous +Mat/CLPGPat animals receiving the CLPG mutation from their father express the phenotype. +Mat/CLPGPat animals are characterized by postnatal, ectopic expression of Delta-like 1 homologue (DLK1) and Paternally expressed gene 11/Retrotransposon-like 1 (PEG11/RTL1) proteins in skeletal muscle. We showed previously in transgenic mice that ectopic expression of DLK1 alone induces a muscular hypertrophy, hence demonstrating a role for DLK1 in determining the callipyge hypertrophy. We herein describe newly generated transgenic mice that ectopically express PEG11 in skeletal muscle, and show that they also exhibit a muscular hypertrophy phenotype. Our data suggest that both DLK1 and PEG11 act together in causing the muscular hypertrophy of callipyge sheep.

Generation of Evc2/Limbin global and conditional KO mice and its roles during mineralized tissue formation.

Ellis-van Creveld (EvC) syndrome (OMIM 225500) is an autosomal recessive disease characterized with chondrodysplastic dwarfism in association with abnormalities in oral cavity. Ciliary proteins EVC and EVC2 have been identified as causative genes and they play an important role on Hedgehog signal transduction. We have also identified a causative gene LIMBIN for bovine chondrodysplastic dwarfism (bcd) that is later identified as the bovine ortholog of EVC2. Here, we report generation of conventional and conditional mutant Evc2/Limbin alleles that mimics mutations found in EvC patients and bcd cattle. Resulted homozygous mice showed no ciliary localization of EVC2 and EVC and displayed reduced Hedgehog signaling activity in association with skeletal and oral defects similar to the EvC patients. Cartilage-specific disruption of Evc2/Limbin resulted in similar but milder skeletal defects, whereas osteoblast-specific disruption did not cause overt changes in skeletal system. Neural crest-specific disruption of Evc2/Limbin resulted in defective incisor growth similar to that seen in conventional knockouts; however, differentiation of amelobolasts was relatively normal in the conditional knockouts. These results showcased functions of EVC2/LIMBIN during formation of mineralized tissues. Availability of the conditional allele for this gene should facilitate further detailed analyses of the role of EVC2/LIMBIN in pathogenesis of EvC syndrome. genesis 53:612-626, 2015. © 2015 Wiley Periodicals, Inc.

Experimental evaluation does not reveal a direct effect of microRNA from the callipyge locus on DLK1 expression.

BACKGROUND: Polar overdominance at the ovine callipyge (CLPG) locus involves the post-transcriptional trans-inhibition of DLK1 in skeletal muscle of CLPG/CLPG sheep. The abundant maternally expressed microRNAs (miRNAs) mapping to the imprinted DLK1-GTL2 domain are prime candidate mediators of this trans-effect. RESULTS: We have tested the affinity of 121 miRNAs processed from this locus for DLK1 by co-transfecting COS1 cells with a vector expressing the full-length ovine DLK1 with corresponding mimic miRNAs. None of the tested miRNAs was able to down regulate DLK1 to the extent observed in vivo. CONCLUSIONS: This suggests that other factors, with or without these miRNAs, are involved in mediating the observed trans-effect.

Deep sequencing reveals abundant noncanonical retroviral microRNAs in B-cell leukemia/lymphoma.

Viral tumor models have significantly contributed to our understanding of oncogenic mechanisms. How transforming delta-retroviruses induce malignancy, however, remains poorly understood, especially as viral mRNA/protein are tightly silenced in tumors. Here, using deep sequencing of broad windows of small RNA sizes in the bovine leukemia virus ovine model of leukemia/lymphoma, we provide in vivo evidence of the production of noncanonical RNA polymerase III (Pol III)-transcribed viral microRNAs in leukemic B cells in the complete absence of Pol II 5'-LTR-driven transcriptional activity. Processed from a cluster of five independent self-sufficient transcriptional units located in a proviral region dispensable for in vivo infectivity, bovine leukemia virus microRNAs represent ∼40% of all microRNAs in both experimental and natural malignancy. They are subject to strong purifying selection and associate with Argonautes, consistent with a critical function in silencing of important cellular and/or viral targets. Bovine leukemia virus microRNAs are strongly expressed in preleukemic and malignant cells in which structural and regulatory gene expression is repressed, suggesting a key role in tumor onset and progression. Understanding how Pol III-dependent microRNAs subvert cellular and viral pathways will contribute to deciphering the intricate perturbations that underlie malignant transformation.

Variants modulating the expression of a chromosome domain encompassing PLAG1 influence bovine stature.

We report mapping of a quantitative trait locus (QTL) with a major effect on bovine stature to a ∼780-kb interval using a Hidden Markov Model-based approach that simultaneously exploits linkage and linkage disequilibrium. We re-sequenced the interval in six sires with known QTL genotype and identified 13 clustered candidate quantitative trait nucleotides (QTNs) out of >9,572 discovered variants. We eliminated five candidate QTNs by studying the phenotypic effect of a recombinant haplotype identified in a breed diversity panel. We show that the QTL influences fetal expression of seven of the nine genes mapping to the ∼780-kb interval. We further show that two of the eight candidate QTNs, mapping to the PLAG1-CHCHD7 intergenic region, influence bidirectional promoter strength and affect binding of nuclear factors. By performing expression QTL analyses, we identified a splice site variant in CHCHD7 and exploited this naturally occurring null allele to exclude CHCHD7 as single causative gene.

Demonstrating polymorphic miRNA-mediated gene regulation in vivo: application to the g+6223G->A mutation of Texel sheep.

We herein describe the development of a biochemical method to evaluate the effect of single nucleotide polymorphisms (SNPs) in target genes on their regulation by microRNAs in vivo. The method is based on the detection of allelic imbalance in RNAs coimmunoprecipitated with AGO proteins from tissues of heterozygous individuals. We characterize the performances of our approach using a model system in a cell culture, and then apply it successfully to prove that the 3'UTR g+6223G-->A mutation operates by promoting RISC-dependent down-regulation of myostatin (MSTN) in skeletal muscle of Texel sheep.

Methylation analysis of the imprinted DLK1-GTL2 domain supports the random parental origin of the IGH-involving del(14q) in B-cell malignancies.

Leukemias/lymphomas with IGH-involving del(14q)(1) commonly lose the DLK1-GTL2 imprinted domain that comprises several paternally and maternally expressed genes, including a cluster of microRNAs. Given that deletion of this region could lead to inactivation of a monoallelically expressed tumor suppressor gene, our study aimed at determination of the parental origin of del(14q/IGH). The designed allele-specific methylation study of the DLK1/GTL2 intergenic differentially methylated region allowed us to determine the parental origin of del(14q/IGH) in 9/20 analyzed cases. In six cases del(14q/IGH) was of the paternal origin and in three cases of the maternal origin. These findings argue against the concept that a TSG/anti-oncomir located in the imprinted region is systematically inactivated by a targeted deletion of its functional allele.

A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep.

Texel sheep are renowned for their exceptional meatiness. To identify the genes underlying this economically important feature, we performed a whole-genome scan in a Romanov x Texel F2 population. We mapped a quantitative trait locus with a major effect on muscle mass to chromosome 2 and subsequently fine-mapped it to a chromosome interval encompassing the myostatin (GDF8) gene. We herein demonstrate that the GDF8 allele of Texel sheep is characterized by a G to A transition in the 3' UTR that creates a target site for mir1 and mir206, microRNAs (miRNAs) that are highly expressed in skeletal muscle. This causes translational inhibition of the myostatin gene and hence contributes to the muscular hypertrophy of Texel sheep. Analysis of SNP databases for humans and mice demonstrates that mutations creating or destroying putative miRNA target sites are abundant and might be important effectors of phenotypic variation.

The callipyge mutation enhances bidirectional long-range DLK1-GTL2 intergenic transcription in cis.

The callipyge mutation (CLPG) is an A to G transition that affects a muscle-specific long-range control element located in the middle of the 90-kb DLK1-GTL2 intergenic (IG) region. It causes ectopic expression of a 327-kb cluster of imprinted genes in skeletal muscle, resulting in the callipyge muscular hypertrophy and its non-Mendelian inheritance pattern known as polar overdominance. We herein demonstrate that the CLPG mutation alters the muscular epigenotype of the DLK1-GTL2 IG region in cis, including hypomethylation, acquisition of novel DNase-I hypersentivite sites, and, most strikingly, strongly enhanced bidirectional, long-range IG transcription. The callipyge phenotype thus emerges as a unique model to study the functional significance of IG transcription, which recently has proven to be a widespread, yet elusive, feature of the mammalian genome.

Optimized beta-galactosidase staining method for simultaneous detection of endogenous gene expression in early mouse embryos.

beta-Galactosidase (beta-gal) is one of the popular reporters for detecting the expression of endogenous or exogenous genes. Here we report 6-chloro-3-indoxyl-beta-D-galactopyranoside (S-gal) is more sensitive for beta-gal activity than 5-bromo-4-chloro-3-indolyl-beta-D-galactoside (X-gal), particularly during the early developmental stages of mouse embryos. Further, we successfully combined beta-gal staining with S-gal and in situ hybridization using DIG-labeled probes in both whole and sections of early stage embryos due to the sensitivity and color compatibility of S-gal.

A comprehensive genetic map of the cattle genome based on 3802 microsatellites.

A microsatellite-based high-density genetic map facilitates for fine mapping of hereditary traits of interest, characterization of meiosis, and providing a foundation for physical map construction. Here, we developed a comprehensive genetic map on the basis of >880,000 genotypes across the USDA MARC cattle reference families with a potential genetic resolution of 0.8 cM at the 95% confidence level ( approximately 800 kb in the bovine genome). We incorporated 2325 microsatellites into the second-generation genetic map by linkage analysis based on sex-averaged two-point LOD scores (>3.0), of which 2293 were fine-mapped by multipoint linkage analysis. The new 3160-cM map comprised of 29 sex-averaged autosomal linkage groups and a sex-specific X-chromosome linkage group includes 3960 markers with 2389 positions, resulting in an average interval size of 1.4 cM. More than half (51%) of the total length of the map is covered with intervals of 2.0 cM or less, and the largest gap is a 10.2-cM interval on the X-linkage group. The new map should accelerate fine mapping and positional cloning of genes for genetic diseases and economically important traits in cattle, as well as related livestock species, such as sheep and goat.