RESUMEN
Most deleterious variants are recessive and segregate at relatively low frequency. Therefore, high sample sizes are required to identify these variants. In this study we report a large-scale sequence based genome-wide association study (GWAS) in pigs, with a total of 120,000 Large White and 80,000 Synthetic breed animals imputed to sequence using a reference population of approximately 1,100 whole genome sequenced pigs. We imputed over 20 million variants with high accuracies (R2>0.9) even for low frequency variants (1-5% minor allele frequency). This sequence-based analysis revealed a total of 14 additive and 9 non-additive significant quantitative trait loci (QTLs) for growth rate and backfat thickness. With the non-additive (recessive) model, we identified a deleterious missense SNP in the CDHR2 gene reducing growth rate and backfat in homozygous Large White animals. For the Synthetic breed, we revealed a QTL on chromosome 15 with a frameshift variant in the OBSL1 gene. This QTL has a major impact on both growth rate and backfat, resembling human 3M-syndrome 2 which is related to the same gene. With the additive model, we confirmed known QTLs on chromosomes 1 and 5 for both breeds, including variants in the MC4R and CCND2 genes. On chromosome 1, we disentangled a complex QTL region with multiple variants affecting both traits, harboring 4 independent QTLs in the span of 5 Mb. Together we present a large scale sequence-based association study that provides a key resource to scan for novel variants at high resolution for breeding and to further reduce the frequency of deleterious alleles at an early stage in the breeding program.
Asunto(s)
Estudio de Asociación del Genoma Completo , Polimorfismo de Nucleótido Simple , Humanos , Animales , Porcinos/genética , Polimorfismo de Nucleótido Simple/genética , Sitios de Carácter Cuantitativo/genética , Fenotipo , Frecuencia de los Genes , Genotipo , Proteínas del Citoesqueleto/genéticaRESUMEN
Climate change is a threat to sustainable livestock production and livelihoods in the tropics. It has adverse impacts on feed and water availability, disease prevalence, production, environmental temperature, and biodiversity. Unravelling the drivers of local adaptation and understanding the underlying genetic variation in random mating indigenous livestock populations informs the design of genetic improvement programmes that aim to increase productivity and resilience. In the present study, we combined environmental, genomic, and phenotypic information of Ethiopian indigenous chickens to investigate their environmental adaptability. Through a hybrid sampling strategy, we captured wide biological and ecological variabilities across the country. Our environmental dataset comprised mean values of 34 climatic, vegetation and soil variables collected over a thirty-year period for 260 geolocations. Our biological dataset included whole genome sequences and quantitative measurements (on eight traits) from 513 individuals, representing 26 chicken populations spread along 4 elevational gradients (6-7 populations per gradient). We performed signatures of selection analyses ([Formula: see text] and XP-EHH) to detect footprints of natural selection, and redundancy analyses (RDA) to determine genotype-environment and genotype-phenotype-associations. RDA identified 1909 outlier SNPs linked with six environmental predictors, which have the highest contributions as ecological drivers of adaptive phenotypic variation. The same method detected 2430 outlier SNPs that are associated with five traits. A large overlap has been observed between signatures of selection identified by[Formula: see text]and XP-EHH showing that both methods target similar selective sweep regions. Average genetic differences measured by [Formula: see text] are low between gradients, but XP-EHH signals are the strongest between agroecologies. Genes in the calcium signalling pathway, those associated with the hypoxia-inducible factor (HIF) transcription factors, and sports performance (GALNTL6) are under selection in high-altitude populations. Our study underscores the relevance of landscape genomics as a powerful interdisciplinary approach to dissect adaptive phenotypic and genetic variation in random mating indigenous livestock populations.
Asunto(s)
Pollos , Genómica , Humanos , Animales , Pollos/genética , Genómica/métodos , Genotipo , Genoma , Selección Genética , Polimorfismo de Nucleótido Simple , Variación GenéticaRESUMEN
BACKGROUND: Integration of high throughput DNA genotyping and RNA-sequencing data enables the discovery of genomic regions that regulate gene expression, known as expression quantitative trait loci (eQTL). In pigs, efforts to date have been mainly focused on purebred lines for traits with commercial relevance as such growth and meat quality. However, little is known on genetic variants and mechanisms associated with the robustness of an animal, thus its overall health status. Here, the liver, lung, spleen, and muscle transcriptomes of 100 three-way crossbred female finishers were studied, with the aim of identifying novel eQTL regulatory regions and transcription factors (TFs) associated with regulation of porcine metabolism and health-related traits. RESULTS: An expression genome-wide association study with 535,896 genotypes and the expression of 12,680 genes in liver, 13,310 genes in lung, 12,650 genes in spleen, and 12,595 genes in muscle resulted in 4,293, 10,630, 4,533, and 6,871 eQTL regions for each of these tissues, respectively. Although only a small fraction of the eQTLs were annotated as cis-eQTLs, these presented a higher number of polymorphisms per region and significantly stronger associations with their target gene compared to trans-eQTLs. Between 20 and 115 eQTL hotspots were identified across the four tissues. Interestingly, these were all enriched for immune-related biological processes. In spleen, two TFs were identified: ERF and ZNF45, with key roles in regulation of gene expression. CONCLUSIONS: This study provides a comprehensive analysis with more than 26,000 eQTL regions identified that are now publicly available. The genomic regions and their variants were mostly associated with tissue-specific regulatory roles. However, some shared regions provide new insights into the complex regulation of genes and their interactions that are involved with important traits related to metabolism and immunity.
Asunto(s)
Estudio de Asociación del Genoma Completo , Sitios de Carácter Cuantitativo , Animales , Porcinos/genética , Polimorfismo de Nucleótido Simple , Femenino , Factores de Transcripción/genética , Factores de Transcripción/metabolismo , Hígado/metabolismo , Especificidad de Órganos/genética , Bazo/metabolismo , Transcriptoma , Regulación de la Expresión Génica , Pulmón/metabolismo , Pulmón/inmunología , GenotipoRESUMEN
In livestock breeding, single nucleotide polymorphism arrays have become a cornerstone of modern livestock breeding. SNP arrays facilitate the identification of genetic markers linked to economically important traits and provide a powerful tool for predicting breeding values. However, conventional breeding programs often overlook additional genomic features contained in the SNP array data that can provide valuable insights into the genetic diversity, copy number variation, inbreeding levels and potential challenges in breeding lines. Here we present genomeprofile, a tool using SNP array-based genomic data, offering a comprehensive profile of breeding animals including the identification of copy number variants and runs of homozygosity, and screening for aneuploidy. By integrating these features into the breeding landscape, genomeprofile enables a more comprehensive picture of genomic variation, ultimately enhancing precision breeding strategies. To illustrate the practicality and efficacy of genomeprofile, we applied the tool to a dataset of four pig breeding lines. The genomeprofile tool is a user-friendly tool that processes genotype data in finalreport or plink ped format efficiently into useful output. The output contains copy number variations, runs of homozygosity, selection signatures, aneuploidy and inbreeding per individual and across populations. This allows breeding companies and researchers to identify unique individuals or regions in the genome of interest based on routinely collected data.
RESUMEN
In general, the relationship between the predicted functional consequences of missense mutations mapping to genes known to be involved in human diseases and the severity of disease manifestations is weak. In this study, we tested in pigs whether missense single nucleotide polymorphisms (SNPs), predicted to have consequences on the function of genes related to lipid metabolism are associated with lipid phenotypes. Association analysis demonstrated that nine out of 72 nominally associated SNPs were classified as "highly" or "very highly consistent" in silico-predicted functional mutations and did not show association with lipid traits expected to be affected by inactivation of the corresponding gene. Although the lack of endophenotypes and the limited sample size of certain genotypic classes might have limited to some extent the reach of the current study, our data indicate that present-day bioinformatic tools have a modest ability to predict the impact of missense mutations on complex phenotypes.
Asunto(s)
Metabolismo de los Lípidos , Mutación Missense , Porcinos , Humanos , Animales , Metabolismo de los Lípidos/genética , Fenotipo , Genotipo , Lípidos , Polimorfismo de Nucleótido SimpleRESUMEN
In the past, small population sizes and unequal ancestor contributions have resulted in high inbreeding rates (ΔF) in the Friesian horse. Two decades ago, the studbook implemented a mating quota and started publishing individual kinships and reduced ΔF below 1% per generation. However, since then, the breeding population size has decreased and this raises the question whether current breeding strategies are sufficient to keep ΔF below desired rates. The aim of this study was to (1) reflect on past inbreeding trends and their main determinants, using pedigree analysis and (2) evaluate the effectiveness of the current and additional breeding strategies using stochastic simulations. We estimated the current ΔF (2013-2022) at 0.72% per generation. While the total contribution of the top 10 sires to the number of offspring per year has decreased from 75% in 1980 to 35% in 2022, this was mainly due to an increased number of approved studbook sires, and not due to more equalized contributions among sires. Of the simulated breeding strategies, selecting only breeding stallions with a below average mean kinship (i.e., "mean kinship selection") was most effective to decrease ΔF (from 0.66% to 0.33%). Increasing the number of breeding sires only had an effect when also a mating quota was applied. However, its effect remained limited. For example, a ~1.5 fold increase, combined with a mating quota of 80 offspring per sire per year, reduced ΔF from 0.55% to 0.51%. When increasing the number of breeding mares, a practically unfeasible large increase was needed for a meaningful reduction in ΔF (e.g. twice as many mares were needed to reduce ΔF from 0.66% to 0.56%). Stratified mating quotas, a novel approach in which we assigned each sire a mating quota (of 60, 80, 100 or 120 offspring per year) based on its mean kinship to recently born foals, resulted in a lower ΔF (0.43%) than a general mating quota of 90 offspring per sire per year (0.55%). Overall, while the current ΔF is below 1%, we recommend to implement additional strategies to further reduce ΔF below 0.5% in the Friesian horse population. For this breed and similar populations, we recommend to focus on breeding strategies based on kinship levels to effectively reduce ΔF.
RESUMEN
Lethal recessive alleles cause pre- or postnatal death in homozygous affected individuals, reducing fertility. Especially in small size domestic and wild populations, those alleles might be exposed by inbreeding, caused by matings between related parents that inherited the same recessive lethal allele from a common ancestor. In this study we report five relatively common (up to 13.4% carrier frequency) recessive lethal haplotypes in two commercial pig populations. The lethal haplotypes have a large effect on carrier-by-carrier matings, decreasing litter sizes by 15.1 to 21.6%. The causal mutations are of different type including two splice-site variants (affecting POLR1B and TADA2A genes), one frameshift (URB1), and one missense (PNKP) variant, resulting in a complete loss-of-function of these essential genes. The recessive lethal alleles affect up to 2.9% of the litters within a single population and are responsible for the death of 0.52% of the total population of embryos. Moreover, we provide compelling evidence that the identified embryonic lethal alleles contribute to the observed heterosis effect for fertility (i.e. larger litters in crossbred offspring). Together, this work marks specific recessive lethal variation describing its functional consequences at the molecular, phenotypic, and population level, providing a unique model to better understand fertility and heterosis in livestock.
Asunto(s)
Genes Letales , Mutación con Pérdida de Función , Sus scrofa/embriología , Sus scrofa/genética , Secuencia de Aminoácidos , Animales , Femenino , Fertilidad/genética , Genes Recesivos , Flujo Genético , Genética de Población , Haplotipos , Vigor Híbrido/genética , Hibridación Genética/genética , Tamaño de la Camada/genética , Masculino , Embarazo , ARN Polimerasa I/genética , Análisis de Secuencia de ARN , Secuenciación Completa del GenomaRESUMEN
Structural variations (SVs) are an important source of phenotypic diversity in cattle. Here, 72 whole genome sequences representing taurine and zebu cattle were used to identify SVs. Applying multiple approaches, 16,738 SVs were identified. A comparison against the Database of Genomic Variants archives revealed that 1575 SVs were novel in our data. A novel duplication covering the entire GALNT15 gene, was observed only in N'Dama. A duplication, which was previously reported only in zebu and associated with navel length, was also observed in N'Dama. Investigation of a novel deletion located upstream of CAST13 gene and identified only in Italian cattle and zebu, revealed its introgressed origin in the former. Overall, our data highlights how the SVs distribution in cattle is also shaped by forces such as demographical differences and gene flow. The cattle SVs of this study and its meta-data can be visualized on an interactive genome browser at https://tinyurl.com/svCowArs.
Asunto(s)
Enfermedades de los Bovinos , Animales , Bovinos/genética , Enfermedades de los Bovinos/genética , Genoma , GenómicaRESUMEN
The genotype-phenotype link is a major research topic in the life sciences but remains highly complex to disentangle. Part of the complexity arises from the number of genes contributing to the observed phenotype. Despite the vast increase of molecular data, pinpointing the causal variant underlying a phenotype of interest is still challenging. In this study, we present an approach to map causal variation and molecular pathways underlying important phenotypes in pigs. We prioritize variation by utilizing and integrating predicted variant impact scores (pCADD), functional genomic information, and associated phenotypes in other mammalian species. We demonstrate the efficacy of our approach by reporting known and novel causal variants, of which many affect non-coding sequences. Our approach allows the disentangling of the biology behind important phenotypes by accelerating the discovery of novel causal variants and molecular mechanisms affecting important phenotypes in pigs. This information on molecular mechanisms could be applicable in other mammalian species, including humans.
Asunto(s)
Variación Genética , Genómica , Animales , Genotipo , Mamíferos , Fenotipo , Porcinos/genéticaRESUMEN
Understanding the genetic basis of similar phenotypes shared between lineages is a long-lasting research interest. Even though animal evolution offers many examples of parallelism, for many phenotypes little is known about the underlying genes and mutations. We here use a combination of whole-genome sequencing, expression analyses, and comparative genomics to study the parallel genetic origin of ptilopody (Pti) in chicken. Ptilopody (or foot feathering) is a polygenic trait that can be observed in domesticated and wild avian species and is characterized by the partial or complete development of feathers on the ankle and feet. In domesticated birds, ptilopody is easily selected to fixation, though extensive variation in the type and level of feather development is often observed. By means of a genome-wide association analysis, we identified two genomic regions associated with ptilopody. At one of the loci, we identified a 17-kb deletion affecting PITX1 expression, a gene known to encode a transcription regulator of hindlimb identity and development. Similarly to pigeon, at the second loci, we observed ectopic expression of TBX5, a gene involved in forelimb identity and a key determinant of foot feather development. We also observed that the trait evolved only once as foot-feathered birds share the same haplotype upstream TBX5. Our findings indicate that in chicken and pigeon ptilopody is determined by the same set of genes that affect similar molecular pathways. Our study confirms that ptilopody has evolved through parallel evolution in chicken and pigeon.
Asunto(s)
Evolución Biológica , Pollos/genética , Plumas/crecimiento & desarrollo , Factores de Transcripción Paired Box/genética , Proteínas de Dominio T Box/genética , Animales , Pollos/crecimiento & desarrollo , Pollos/metabolismo , Columbidae/genética , Pie , Haplotipos , Herencia Multifactorial , Factores de Transcripción Paired Box/metabolismo , Proteínas de Dominio T Box/metabolismo , Secuenciación Completa del GenomaRESUMEN
Livestock populations can be used to study recessive defects caused by deleterious alleles. The frequency of deleterious alleles including recessive lethal alleles can stay at high or moderate frequency within a population, especially if recessive lethal alleles exhibit an advantage for favourable traits in heterozygotes. In this study, we report such a recessive lethal deletion of 212kb (del) within the BBS9 gene in a breeding population of pigs. The deletion produces a truncated BBS9 protein expected to cause a complete loss-of-function, and we find a reduction of approximately 20% on the total number of piglets born from carrier by carrier matings. Homozygous del/del animals die mid- to late-gestation, as observed from high increase in numbers of mummified piglets resulting from carrier-by-carrier crosses. The moderate 10.8% carrier frequency (5.4% allele frequency) in this pig population suggests an advantage on a favourable trait in heterozygotes. Indeed, heterozygous carriers exhibit increased growth rate, an important selection trait in pig breeding. Increased growth and appetite together with a lower birth weight for carriers of the BBS9 null allele in pigs is analogous to the phenotype described in human and mouse for (naturally occurring) BBS9 null-mutants. We show that fetal death, however, is induced by reduced expression of the downstream BMPER gene, an essential gene for normal foetal development. In conclusion, this study describes a lethal 212kb deletion with pleiotropic effects on two different genes, one resulting in fetal death in homozygous state (BMPER), and the other increasing growth (BBS9) in heterozygous state. We provide strong evidence for balancing selection resulting in an unexpected high frequency of a lethal allele in the population. This study shows that the large amounts of genomic and phenotypic data routinely generated in modern commercial breeding programs deliver a powerful tool to monitor and control lethal alleles much more efficiently.
Asunto(s)
Regulación del Desarrollo de la Expresión Génica , Frecuencia de los Genes , Genes Letales/fisiología , Endogamia , Sus scrofa/genética , Animales , Conjuntos de Datos como Asunto , Femenino , Fertilidad/genética , Genes Recesivos/fisiología , Técnicas de Genotipaje , Heterocigoto , Homocigoto , Masculino , Modelos Animales , Sus scrofa/crecimiento & desarrolloRESUMEN
BACKGROUND: A balanced constitutional reciprocal translocation (RT) is a mutual exchange of terminal segments of two non-homologous chromosomes without any loss or gain of DNA in germline cells. Carriers of balanced RTs are viable individuals with no apparent phenotypical consequences. These animals produce, however, unbalanced gametes and show therefore reduced fertility and offspring with congenital abnormalities. This cytogenetic abnormality is usually detected using chromosome staining techniques. The aim of this study was to test the possibilities of using paired end short read sequencing for detection of balanced RTs in boars and investigate their breakpoints and junctions. RESULTS: Balanced RTs were recovered in a blinded analysis, using structural variant calling software DELLY, in 6 of the 7 carriers with 30 fold short read paired end sequencing. In 15 non-carriers we did not detect any RTs. Reducing the coverage to 20 fold, 15 fold and 10 fold showed that at least 20 fold coverage is required to obtain good results. One RT was not detected using the blind screening, however, a highly likely RT was discovered after unblinding. This RT was located in a repetitive region, showing the limitations of short read sequence data. The detailed analysis of the breakpoints and junctions suggested three junctions showing microhomology, three junctions with blunt-end ligation, and three micro-insertions at the breakpoint junctions. The RTs detected also showed to disrupt genes. CONCLUSIONS: We conclude that paired end short read sequence data can be used to detect and characterize balanced reciprocal translocations, if sequencing depth is at least 20 fold coverage. However, translocations in repetitive areas may require large fragments or even long read sequence data.
Asunto(s)
Aberraciones Cromosómicas , Translocación Genética , Animales , ADN , Heterocigoto , Masculino , Porcinos/genéticaRESUMEN
Genome sequences nowadays play a central role in molecular biology and bioinformatics. These sequences are shared with the scientific community through sequence databases. The sequence repositories of the International Nucleotide Sequence Database Collaboration (INSDC, comprising GenBank, ENA and DDBJ) are the largest in the world. Preparing an annotated sequence in such a way that it will be accepted by the database is challenging because many validation criteria apply. In our opinion, it is an undesirable situation that researchers who want to submit their sequence need either a lot of experience or help from partners to get the job done. To save valuable time and money, we list a number of recommendations for people who want to submit an annotated genome to a sequence database, as well as for tool developers, who could help to ease the process.
Asunto(s)
Genoma , Bases de Datos de Ácidos Nucleicos , Humanos , Internet , National Library of Medicine (U.S.) , Estados UnidosRESUMEN
The rate of directional environmental change may have profound consequences for evolutionary dynamics and outcomes. Yet, most evolution experiments impose a sudden large change in the environment, after which the environment is kept constant. We previously cultured replicate Saccharomyces cerevisiae populations for 500 generations in the presence of either gradually increasing or constant high concentrations of the heavy metals cadmium, nickel, and zinc. Here, we investigate how each of these treatments affected genomic evolution. Whole-genome sequencing of evolved clones revealed that adaptation occurred via a combination of SNPs, small indels, and whole-genome duplications and other large-scale structural changes. In contrast to some theoretical predictions, gradual and abrupt environmental change caused similar numbers of genomic changes. For cadmium, which is toxic already at comparatively low concentrations, mutations in the same genes were used for adaptation to both gradual and abrupt increase in concentration. Conversely, for nickel and zinc, which are toxic at high concentrations only, mutations in different genes were used for adaptation depending on the rate of change. Moreover, evolution was more repeatable following a sudden change in the environment, particularly for nickel and zinc. Our results show that the rate of environmental change and the nature of the selection pressure are important drivers of evolutionary dynamics and outcomes, which has implications for a better understanding of societal problems such as climate change and pollution.
Asunto(s)
Adaptación Fisiológica/genética , Saccharomyces cerevisiae/genética , Selección Genética/genética , Aclimatación , Adaptación Biológica , Evolución Biológica , Evolución Molecular Dirigida , Ambiente , Genoma , Genómica , Mutación , Níquel/metabolismo , Zinc/metabolismoRESUMEN
BACKGROUND: Sex-linked slow (SF) and fast (FF) feathering rates at hatch have been widely used in poultry breeding for autosexing at hatch. In chicken, the sex-linked K (SF) and k+ (FF) alleles are responsible for the feathering rate phenotype. Allele K is dominant and a partial duplication of the prolactin receptor gene has been identified as the causal mutation. Interestingly, some domesticated turkey lines exhibit similar slow- and fast-feathering phenotypes, but the underlying genetic components and causal mutation have never been investigated. In this study, our aim was to investigate the molecular basis of feathering rate at hatch in domestic turkey. RESULTS: We performed a sequence-based case-control association study and detected a genomic region on chromosome Z, which is statistically associated with rate of feathering at hatch in turkey. We identified a 5-bp frameshift deletion in the prolactin receptor (PRLR) gene that is responsible for slow feathering at hatch. All female cases (SF turkeys) were hemizygous for this deletion, while 188 controls (FF turkeys) were hemizygous or homozygous for the reference allele. This frameshift mutation introduces a premature stop codon and six novel amino acids (AA), which results in a truncated PRLR protein that lacks 98 C-terminal AA. CONCLUSIONS: We present the causal mutation for feathering rate in turkey that causes a partial C-terminal loss of the prolactin receptor, and this truncated PRLR protein is strikingly similar to the protein encoded by the slow feathering K allele in chicken.
Asunto(s)
Pollos/genética , Mutación del Sistema de Lectura , Receptores de Prolactina/genética , Análisis de Secuencia de ADN/veterinaria , Pavos/genética , Alelos , Secuencia de Aminoácidos , Animales , Pollos/metabolismo , Duplicación Cromosómica , Plumas , Femenino , Estudios de Asociación Genética/veterinaria , Hemicigoto , Masculino , Fenotipo , Receptores de Prolactina/metabolismo , Pavos/metabolismoRESUMEN
BACKGROUND: Deleterious genetic variation can increase in frequency as a result of mutations, genetic drift, and genetic hitchhiking. Although individual effects are often small, the cumulative effect of deleterious genetic variation can impact population fitness substantially. In this study, we examined the genome of commercial purebred chicken lines for deleterious and functional variations, combining genotype and whole-genome sequence data. RESULTS: We analysed over 22,000 animals that were genotyped on a 60 K SNP chip from four purebred lines (two white egg and two brown egg layer lines) and two crossbred lines. We identified 79 haplotypes that showed a significant deficit in homozygous carriers. This deficit was assumed to stem from haplotypes that potentially harbour lethal recessive variations. To identify potentially deleterious mutations, a catalogue of over 10 million variants was derived from 250 whole-genome sequenced animals from three purebred white-egg layer lines. Out of 4219 putative deleterious variants, 152 mutations were identified that likely induce embryonic lethality in the homozygous state. Inferred deleterious variation showed evidence of purifying selection and deleterious alleles were generally overrepresented in regions of low recombination. Finally, we found evidence that mutations, which were inferred to be evolutionally intolerant, likely have positive effects in commercial chicken populations. CONCLUSIONS: We present a comprehensive genomic perspective on deleterious and functional genetic variation in egg layer breeding lines, which are under intensive selection and characterized by a small effective population size. We show that deleterious variation is subject to purifying selection and that there is a positive relationship between recombination rate and purging efficiency. In addition, multiple putative functional coding variants were discovered in selective sweep regions, which are likely under positive selection. Together, this study provides a unique molecular perspective on functional and deleterious variation in commercial egg-laying chickens, which can enhance current genomic breeding practices to lower the frequency of undesirable variants in the population.
Asunto(s)
Pollos/genética , Polimorfismo de Nucleótido Simple , Eliminación de Secuencia , Secuenciación Completa del Genoma/veterinaria , Animales , Animales Domésticos , Cruzamiento , Variación Genética , Genotipo , Haplotipos , Recombinación Genética , Selección GenéticaRESUMEN
BACKGROUND: Lethal recessive variation can cause prenatal death of homozygous offspring. Although usually present at low-frequency in populations, the impact on individual fitness can be substantial. Until recently, the presence of recessive embryonic lethal variation could only be measured indirectly through reduced fertility. In this study, we estimate the presence of genetic loci associated with both early and late termination of development during gestation in pigs from the wealth of genome data routinely generated by a commercial breeding company. RESULTS: We examined three commercial pig (Sus scrofa) populations for potentially deleterious genetic variation based on 80 K SNP-chip genotypes, and estimate the effects on reproductive traits. 24,000 pigs from three populations were analyzed for missing or depletion of homozygous haplotypes. We identified 145 haplotypes (ranging from 0.5-4 Mb in size) in the genome with complete absence or depletion of homozygous animals. Thirty-five haplotypes show a negative effect on at least one of the analysed reproductive traits (total number born, number of stillborn, and number of mummified piglets). One variant in particular appeared to result in relative late termination of development of fetuses, responsible for a significant fraction of observed stillborn piglets ('mummies'), as they die mid-gestation. Moreover, we identified the BMPER gene as a likely candidate underlying this phenomenon. CONCLUSIONS: Our study shows that although lethal recessive variation is present, the frequency of these alleles is invariably low in these highly managed populations. Nevertheless, due to cumulative effects of deleterious variants, large numbers of affected offspring are produced. Furthermore, our study demonstrates the use of a large-scale commercial genetic experiment to systematically screen for 'natural knockouts' that can increase understanding of gene function.
Asunto(s)
Genes Recesivos/genética , Polimorfismo de Nucleótido Simple , Sus scrofa/genética , Animales , Haplotipos , Homocigoto , Encuestas y CuestionariosRESUMEN
BACKGROUND: Folsomia candida is a model in soil biology, belonging to the family of Isotomidae, subclass Collembola. It reproduces parthenogenetically in the presence of Wolbachia, and exhibits remarkable physiological adaptations to stress. To better understand these features and adaptations to life in the soil, we studied its genome in the context of its parthenogenetic lifestyle. RESULTS: We applied Pacific Bioscience sequencing and assembly to generate a reference genome for F. candida of 221.7 Mbp, comprising only 162 scaffolds. The complete genome of its endosymbiont Wolbachia, was also assembled and turned out to be the largest strain identified so far. Substantial gene family expansions and lineage-specific gene clusters were linked to stress response. A large number of genes (809) were acquired by horizontal gene transfer. A substantial fraction of these genes are involved in lignocellulose degradation. Also, the presence of genes involved in antibiotic biosynthesis was confirmed. Intra-genomic rearrangements of collinear gene clusters were observed, of which 11 were organized as palindromes. The Hox gene cluster of F. candida showed major rearrangements compared to arthropod consensus cluster, resulting in a disorganized cluster. CONCLUSIONS: The expansion of stress response gene families suggests that stress defense was important to facilitate colonization of soils. The large number of HGT genes related to lignocellulose degradation could be beneficial to unlock carbohydrate sources in soil, especially those contained in decaying plant and fungal organic matter. Intra- as well as inter-scaffold duplications of gene clusters may be a consequence of its parthenogenetic lifestyle. This high quality genome will be instrumental for evolutionary biologists investigating deep phylogenetic lineages among arthropods and will provide the basis for a more mechanistic understanding in soil ecology and ecotoxicology.
Asunto(s)
Artrópodos/genética , Artrópodos/fisiología , Genómica , Suelo , Animales , Antibacterianos/biosíntesis , Artrópodos/metabolismo , Reordenamiento Génico , Transferencia de Gen Horizontal , Familia de Multigenes/genética , FilogeniaRESUMEN
Grain legume improvement is currently impeded by a lack of genomic resources. The paucity of genome information for faba bean can be attributed to the intrinsic difficulties of assembling/annotating its giant (~13 Gb) genome. In order to address this challenge, RNA-sequencing analysis was performed on faba bean (cv. Wizard) leaves. Read alignment to the faba bean reference transcriptome identified 16 300 high quality unigenes. In addition, Illumina paired-end sequencing was used to establish a baseline for genomic information assembly. Genomic reads were assembled de novo into contigs with a size range of 50-5000 bp. Over 85% of sequences did not align to known genes, of which ~10% could be aligned to known repetitive genetic elements. Over 26 000 of the reference transcriptome unigenes could be aligned to DNA-sequencing (DNA-seq) reads with high confidence. Moreover, this comparison identified 56 668 potential splice points in all identified unigenes. Sequence length data were extended at 461 putative loci through alignment of DNA-seq contigs to full-length, publicly available linkage marker sequences. Reads also yielded coverages of 3466× and 650× for the chloroplast and mitochondrial genomes, respectively. Inter- and intraspecies organelle genome comparisons established core legume organelle gene sets, and revealed polymorphic regions of faba bean organelle genomes.
Asunto(s)
Productos Agrícolas/genética , Genoma Mitocondrial , Genoma de Planta , Genómica/métodos , Fitomejoramiento/métodos , Análisis de Secuencia de ARN/métodos , Vicia faba/genéticaRESUMEN
BACKGROUND: Studies on vertebrate DNA methylomes have revealed a regulatory role of tissue specific DNA methylation in relation to gene expression. However, it is not well known how tissue-specific methylation varies between different functional and structural components of genes and genomes. Using whole-genome bisulfite sequencing data we here describe both CpG and non-CpG methylation profiles of whole blood and brain tissue in relation to gene features, CpG-islands (CGIs), transposable elements (TE), and their functional roles in an ecological model species, the great tit (Parus major). RESULTS: We show that hypomethylation at the transcription start site (TSS) is enriched in genes with functional classes that relate directly to processes specific to each tissue type. We find that 6877 (~21 %) of the CGIs are differentially methylated between blood and brain, of which 1186 and 2055 are annotated to promoter and intragenic regions, respectively. We observe that CGI methylation in promoter regions is more conserved between tissues compared to CGI methylation in intra and inter-genic regions. Differentially methylated CGIs in promoter and intragenic regions are overrepresented in genomic loci linked to development, suggesting a distinct role for CGI methylation in regulating expression during development. Additionally, we find significant non-CpG methylation in brain but not in blood with a strong preference for methylation at CpA dinucleotide sites. Finally, CpG hypermethylation of TEs is significantly stronger in brain compared to blood, but does not correlate with TE activity. Surprisingly, TEs showed significant hypomethylation in non-CpG contexts which was negatively correlated with TE expression. CONCLUSION: The discovery that TSS methylation levels are directly linked to functional classes related to each tissue provides new insights in the regulatory role of DNA-methylation patterns. The dominant sequence motifs for brain non-CpG methylation, similar to those found in mammals, suggests that a conserved non-CpG regulatory mechanism was already present in the amniote ancestor. The negative correlation between brain non-CpG methylation and TE activity (not found for CpG methylation) suggests that non-CpG is the dominant regulatory form of methylation in TE silencing.