Search

Background Feed accounts for about 70% of the total cost of poultry meat production. Residual feed intake (RFI) has become the preferred measure of feed efficiency because it is phenotypically independent of growth rate and body weight. In this study, our aim was to estimate genetic parameters and identify quantitative trait loci (QTL) for feed efficiency in 3314 purebred broilers using a genome-wide association study. Broilers were genotyped using a custom 55 K single nucleotide polymorphism (SNP) array. Results Estimates of genomic heritability for seven growth and feed efficiency traits, including body weight at 28 days of age (BW28), BW42, average daily feed intake (ADFI), RFI, and RFI adjusted for weight of abdominal fat (RFIa), ranged from 0.12 to 0.26. Eleven genome-wide significant SNPs and 15 suggestively significant SNPs were detected, of which 19 clustered around two genomic regions. A region on chromosome 16 (2.34–2.66 Mb) was associated with both BW28 and BW42, and the most significant SNP in this region, AX_101003762, accounted for 7.6% of the genetic variance of BW28. The other region, on chromosome 1 (91.27–92.43 Mb) was associated with RFI and ADFI, and contains the NSUN3 and EPHA6 as candidate genes. The most significant SNP in this region, AX_172588157, accounted for 4.4% of the genetic variance of RFI. In addition, a genomic region containing the gene AGK on chromosome 1 was found to be associated with RFIa. The NSUN3 and AGK genes were found to be differentially expressed in breast muscle, thigh muscle, and abdominal fat between male broilers with high and low RFI. Conclusions We identified QTL regions for BW28 and BW42 (spanning 0.32 Mb) and RFI (spanning 1.16 Mb). The NSUN3, EPHA6, and AGK were identified as the most likely candidate genes for these QTL. These genes are involved in mitochondrial function and behavioral regulation. These results contribute to the identification of candidate genes and variants for growth and feed efficiency in poultry.

Background The pond snail Lymnaea stagnalis (L. stagnalis) has been widely used as a model organism in neurobiology, ecotoxicology, and parasitology due to the relative simplicity of its central nervous system (CNS). However, its usefulness is restricted by a limited availability of transcriptome data. While sequence information for the L. stagnalis CNS transcripts has been obtained from EST libraries and a de novo RNA-seq assembly, the quality of these assemblies is limited by a combination of low coverage of EST libraries, the fragmented nature of de novo assemblies, and lack of reference genome. Results In this study, taking advantage of the recent availability of a preliminary L. stagnalis genome, we generated an RNA-seq library from the adult L. stagnalis CNS, using a combination of genome-guided and de novo assembly programs to identify 17,832 protein-coding L. stagnalis transcripts. We combined our library with existing resources to produce a transcript set with greater sequence length, completeness, and diversity than previously available ones. Using our assembly and functional domain analysis, we profiled L. stagnalis CNS transcripts encoding ion channels and ionotropic receptors, which are key proteins for CNS function, and compared their sequences to other vertebrate and invertebrate model organisms. Interestingly, L. stagnalis transcripts encoding numerous putative Ca2+ channels showed the most sequence similarity to those of Mus musculus, Danio rerio, Xenopus tropicalis, Drosophila melanogaster, and Caenorhabditis elegans, suggesting that many calcium channel-related signaling pathways may be evolutionarily conserved. Conclusions Our study provides the most thorough characterization to date of the L. stagnalis transcriptome and provides insights into differences between vertebrates and invertebrates in CNS transcript diversity, according to function and protein class. Furthermore, this study provides a complete characterization of the ion channels of Lymnaea stagnalis, opening new avenues for future research on fundamental neurobiological processes in this model system.

Background Susceptibility of European sea bass (Dicentrarchus labrax L.) to viral nervous necrosis (VNN) is well-known. Interest towards selective breeding as a tool to enhance genetic resistance in this species has increased sharply due to the major threat represented by VNN for farmed sea bass and limitations concerning specific therapeutical measures. A sea bass experimental population (N = 650) was challenged with nervous necrosis virus (NNV) to investigate genetic variation in VNN mortality. In addition, relationships of this trait with serum cortisol concentration after stress exposure, antibody titer against NNV antigens, and body weight at a fixed age were studied to identify potential indicator traits of VNN resistance. Results The estimate of heritability for VNN mortality was moderate and ranged from 0.15 (HPD95%, 95% highest posterior density interval: 0.02, 0.31) to 0.23 (HPD95%: 0.06, 0.47). Heritability estimates for cortisol concentration, antibody titer, and body weight were 0.19 (HPD95%: 0.07, 0.34), 0.36 (HPD95%: 0.16, 0.59) and 0.57 (HPD95%: 0.33, 0.84), respectively. Phenotypic relationships between traits were trivial and not statistically significant, except for the estimated correlation between antibody titer and body weight (0.24). Genetic correlations of mortality with body weight or antibody titer (− 0.39) exhibited a 0.89 probability of being negative. A negligible genetic correlation between mortality and cortisol concentration was detected. Antibody titer was estimated to be positively correlated with body weight (0.49). Conclusions Antibody titer against NNV offers the opportunity to use indirect selection to enhance resistance, while the use of cortisol concentration as an indicator trait in breeding programs for VNN resistance is questionable. The estimate of heritability for VNN mortality indicates the feasibility of selective breeding to enhance resistance to NNV and raises attention to the development of genomic prediction tools to simplify testing procedures for selection candidates.

The integration of a viral genome into the host genome has a major impact on the trajectory of the infected cell. Integration location and variation within the associated viral genome can influence both clonal expansion and persistence of infected cells. Methods based on short-read sequencing can identify viral insertion sites, but the sequence of the viral genomes within remains unobserved. We develop PCIP-seq, a method that leverages long reads to identify insertion sites and sequence their associated viral genome. We apply the technique to exogenous retroviruses HTLV-1, BLV, and HIV-1, endogenous retroviruses, and human papillomavirus., SCOPUS: ar.j, info:eu-repo/semantics/published

ERANET grant BBSRC BB/M028291/1, Swedish Research Council Formas 2014-01840, ARN (Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria) 291815, European Cooperation in Science and Technology (COST) CA15224, UK Research & Innovation (UKRI) Biotechnology and Biological Sciences Research Council (BBSRC) BB/P013759/1

BackgroundProduct quality and production efficiency of Atlantic salmon are, to a large extent, influenced by the deposition and depletion of lipid reserves. Fillet lipid content is a heritable trait and is unfavourably correlated with growth, thus genetic management of fillet lipid content is needed for sustained genetic progress in these two traits. The laboratory-based reference method for recording fillet lipid content is highly accurate and precise but, at the same time, expensive, time-consuming, and destructive. Here, we test the use of rapid and cheaper vibrational spectroscopy methods, namely near-infrared (NIR) and Raman spectroscopy both as individual phenotypes and phenotypic predictors of lipid content in Atlantic salmon.ResultsRemarkably, 827 of the 1500 individual Raman variables (i.e. Raman shifts) of the Raman spectrum were significantly heritable (heritability (h2) ranging from 0.15 to 0.65). Similarly, 407 of the 2696 NIR spectral landscape variables (i.e. wavelengths) were significantly heritable (h2= 0.27–0.40). Both Raman and NIR spectral landscapes had significantly heritable regions, which are also informative in spectroscopic predictions of lipid content. Partial least square predicted lipid content using Raman and NIR spectra were highly concordant and highly genetically correlated with the lipid content values ($${r}_{\text{g}}$$rg= 0.91–0.98) obtained with the reference method using Lin’s concordance correlation coefficient (CCC = 0.63–0.90), and were significantly heritable ($${h}^{2}$$h2= 0.52–0.67).ConclusionsBoth NIR and Raman spectral landscapes show substantial additive genetic variation and are highly genetically correlated with the reference method. These findings lay down the foundation for rapid spectroscopic measurement of lipid content in salmonid breeding programmes.

Background Single nucleotide polymorphism (SNP) array technology has been increasingly used to generate large quantities of SNP data for use in genetic studies. As new arrays are developed to take advantage of new technology and of improved probe design using new genome sequence and panel data, a need to integrate data from different arrays and array platforms has arisen. This study was undertaken in view of our need for an integrated high-quality dataset of Illumina Infinium® 20 K and Affymetrix Axiom® 480 K SNP array data in apple (Malus × domestica). In this study, we qualify and quantify the compatibility of SNP calling, defined as SNP calls that are both accurate and concordant, across both arrays by two approaches. First, the concordance of SNP calls was evaluated using a set of 417 duplicate individuals genotyped on both arrays starting from a set of 10,295 robust SNPs on the Infinium array. Next, the accuracy of the SNP calls was evaluated on additional germplasm (n = 3141) from both arrays using Mendelian inconsistent and consistent errors across thousands of pedigree links. While performing this work, we took the opportunity to evaluate reasons for probe failure and observed discordant SNP calls. Results Concordance among the duplicate individuals was on average of 97.1% across 10,295 SNPs. Of these SNPs, 35% had discordant call(s) that were further curated, leading to a final set of 8412 (81.7%) SNPs that were deemed compatible. Compatibility was highly influenced by the presence of alternate probe binding locations and secondary polymorphisms. The impact of the latter was highly influenced by their number and proximity to the 3′ end of the probe. Conclusions The Infinium and Axiom SNP array data were mostly compatible. However, data integration required intense data filtering and curation. This work resulted in a workflow and information that may be of use in other data integration efforts. Such an in-depth analysis of array concordance and accuracy as ours has not been previously described in the literature and will be useful in future work on SNP array data integration and interpretation, and in probe/platform development.

International audience; While virtually all animals show certain abilities for regeneration after an injury, these abilities vary greatly among metazoans. Porifera (Sponges) is basal metazoans characterized by a wide variety of different regenerative processes, including whole-body regeneration (WBR). Considering phylogenetic position and unique body organization, sponges are highly promising models, as they can shed light on the origin and early evolution of regeneration in general and WBR in particular. The present review summarizes available data on the morphogenetic and cellular mechanisms accompanying different types of WBR in sponges. Sponges show a high diversity of WBR, which principally could be divided into (1) WBR from a body fragment and (2) WBR by aggregation of dissociated cells. Sponges belonging to different phylogenetic clades and even to different species and/or differing in the anatomical structure undergo different morphogeneses after similar operations. A common characteristic feature of WBR in sponges is the instability of the main body axis: a change of the organism polarity is described during all types of WBR. The cellular mechanisms of WBR are different across sponge classes, while cell dedifferentiations and transdifferentiations are involved in regeneration processes in all sponges. Data considering molecular regulation of WBR in sponges are extremely scarce. However, the possibility to achieve various types of WBR ensured by common morphogenetic and cellular basis in a single species makes sponges highly accessible for future comprehensive physiological, biochemical, and molecular studies of regeneration processes.

BackgroundAedes aegyptimosquito, the principal global vector of arboviral diseases, lays eggs and undergoes larval and pupal development to become adult mosquitoes in fresh water (FW). It has recently been observed to develop in coastal brackish water (BW) habitats of up to 50% sea water, and such salinity tolerance shown to be an inheritable trait. Genomics of salinity tolerance inAe. aegyptihas not been previously studied, but it is of fundamental biological interest and important for controlling arboviral diseases in the context of rising sea levels increasing coastal ground water salinity.ResultsBW- and FW-Ae. aegyptiwere compared by RNA-seq analysis on the gut, anal papillae and rest of the carcass in fourth instar larvae (L4), proteomics of cuticles shed when L4 metamorphose into pupae, and transmission electron microscopy of cuticles in L4 and adults. Genes for specific cuticle proteins, signalling proteins, moulting hormone-related proteins, membrane transporters, enzymes involved in cuticle metabolism, and cytochrome P450 showed different mRNA levels in BW and FW L4 tissues. The salinity-tolerantAe. aegyptiwere also characterized by altered L4 cuticle proteomics and changes in cuticle ultrastructure of L4 and adults.ConclusionsThe findings provide new information on molecular and ultrastructural changes associated with salinity adaptation in FW mosquitoes. Changes in cuticles of larvae and adults of salinity-tolerantAe. aegyptiare expected to reduce the efficacy of insecticides used for controlling arboviral diseases. Expansion of coastal BW habitats and their neglect for control measures facilitates the spread of salinity-tolerantAe. aegyptiand genes for salinity tolerance. The transmission of arboviral diseases can therefore be amplified in multiple ways by salinity-tolerantAe. aegyptiand requires appropriate mitigating measures. The findings inAe. aegypti have attendant implications for the development of salinity tolerance in other fresh water mosquito vectors and the diseases they transmit.

International audience; In mammals, milk is essential for the growth, development, and health. Milk quantity and quality are dependent on mammary development, strongly influenced by nutrition. This review provides an overview of the data on nutritional regulations of mammary development and gene expression involved in milk component synthesis. Mammary development is described related to rodents, rabbits, and pigs, common models in mammary biology. Molecular mechanisms of the nutritional regulation of milk synthesis are reported in ruminants regarding the importance of ruminant milk in human health. The effects of dietary quantitative and qualitative alterations are described considering the dietary composition and in regard to the periods of nutritional susceptibly. During lactation, the effects of lipid supplementation and feed restriction or deprivation are discussed regarding gene expression involved in milk biosynthesis, in ruminants. Moreover, nutrigenomic studies underline the role of the mammary structure and the potential influence of microRNAs. Knowledge from three lactating and three dairy livestock species contribute to understanding the variety of phenotypes reported in this review and highlight (1) the importance of critical physiological stages, such as puberty gestation and early lactation and (2) the relative importance of the various nutrients besides the total energetic value and their interaction.

BackgroundBarley is one of the founder crops of Neolithic agriculture and is among the most-grown cereals today. The only trait that universally differentiates the cultivated and wild subspecies is ‘non-brittleness’ of the rachis (the stem of the inflorescence), which facilitates harvesting of the crop. Other phenotypic differences appear to result from facultative or regional selective pressures. The population structure resulting from these regional events has been interpreted as evidence for multiple domestications or a mosaic ancestry involving genetic interaction between multiple wild or proto-domesticated lineages. However, each of the three mutations that confer non-brittleness originated in the western Fertile Crescent, arguing against multiregional origins for the crop.ResultsWe examined exome data for 310 wild, cultivated and hybrid/feral barley accessions and showed that cultivated barley is structured into six genetically-defined groups that display admixture, resulting at least in part from two or more significant passages of gene flow with distinct wild populations. The six groups are descended from a single founding population that emerged in the western Fertile Crescent. Only a few loci were universally targeted by selection, the identity of these suggesting that changes in seedling emergence and pathogen resistance could represent crucial domestication switches. Subsequent selection operated on a regional basis and strongly contributed to differentiation of the genetic groups.ConclusionsIdentification of genetically-defined groups provides clarity to our understanding of the population history of cultivated barley. Inference of population splits and mixtures together with analysis of selection sweeps indicate descent from a single founding population, which emerged in the western Fertile Crescent. This founding population underwent relatively little genetic selection, those changes that did occur affecting traits involved in seedling emergence and pathogen resistance, indicating that these phenotypes should be considered as ‘domestication traits’. During its expansion out of the western Fertile Crescent, the crop underwent regional episodes of gene flow and selection, giving rise to a modern genetic signature that has been interpreted as evidence for multiple domestications, but which we show can be rationalized with a single origin.

International audience; This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BYCyanobacteria are widely-diverse, environmentally crucial photosynthetic prokaryotes of great interests for basic and applied science. Work to date has focused mostly on the three non-nitrogen fixing unicellular species Synechocystis PCC 6803, Synechococcus PCC 7942, and Synechococcus PCC 7002, which have been selected for their genetic and physiological interests summarized in this review. Extensive “omics” data sets have been generated, and genome-scale models (GSM) have been developed for the rational engineering of these cyanobacteria for biotechnological purposes. We presently discuss what should be done to improve our understanding of the genotype-phenotype relationships of these models and generate robust and predictive models of their metabolism. Furthermore, we also emphasize that because Synechocystis PCC 6803, Synechococcus PCC 7942, and Synechococcus PCC 7002 represent only a limited part of the wide biodiversity of cyanobacteria, other species distantly related to these three models, should be studied. Finally, we highlight the need to strengthen the communication between academic researchers, who know well cyanobacteria and can engineer them for biotechnological purposes, but have a limited access to large photobioreactors, and industrial partners who attempt to use natural or engineered cyanobacteria to produce interesting chemicals at reasonable costs, but may lack knowledge on cyanobacterial physiology and metabolism. View Full-TextKeywords: Synechocystis PCC 6803; Synechococcus PCC 7942; Synechococcus PCC 7002; biodiversity; genotype-phenotype relationships; conjugation; transformation; RSF1010 plasmid; neutral cloning sites

Ornithine transcarbamylase deficiency (OTCD) is an X-linked liver disorder caused by partial or total loss of OTC enzyme activity. It is characterized by elevated plasma ammonia, leading to neurological impairments, coma, and death in the most severe cases. OTCD is managed by combining dietary restrictions, essential amino acids, and ammonia scavengers. However, to date, liver transplantation provides the best therapeutic outcome. AAV-mediated gene-replacement therapy represents a promising curative strategy. Here, we generated an AAV2/8 vector expressing a codon-optimized human OTC cDNA by the α1-AAT liver-specific promoter. Unlike standard codon-optimization approaches, we performed multiple codon-optimization rounds via common algorithms and ortholog sequence analysis that significantly improved mRNA translatability and therapeutic efficacy. AAV8-hOTC-CO (codon optimized) vector injection into adult OTCSpf-Ash mice (5.0E11 vg/kg) mediated long-term complete correction of the phenotype. Adeno-Associated viral (AAV) vector treatment restored the physiological ammonia detoxification liver function, as indicated by urinary orotic acid normalization and by conferring full protection against an ammonia challenge. Removal of liver-specific transcription factor binding sites from the AAV backbone did not affect gene expression levels, with a potential improvement in safety. These results demonstrate that AAV8-hOTC-CO gene transfer is safe and results in sustained correction of OTCD in mice, supporting the translation of this approach to the clinic., Graphical Abstract, This study presents an efficient AAV-mediated liver-specific gene therapy for OTCD, with improved vector safety features, resulting in steady hOTC liver expression and complete phenotype normalization of Spf-Ash mice, even after an ammonia challenge, at relatively low vector doses. These results support the potential translation of this strategy to the clinic.

Cystic fibrosis (CF) is a chronic genetic disease that mainly affects the respiratory and gastrointestinal systems. No curative treatments are available, but the follow-up in specialized centers has greatly improved the patient life expectancy. Robust biomarkers are required to monitor the disease, guide treatments, stratify patients, and provide outcome measures in clinical trials. In the present study, we outline a strategy to select putative DNA methylation biomarkers of lung disease severity in cystic fibrosis patients. In the discovery step, we selected seven potential biomarkers using a genome-wide DNA methylation dataset that we generated in nasal epithelial samples from the MethylCF cohort. In the replication step, we assessed the same biomarkers using sputum cell samples from the MethylBiomark cohort. Of interest, DNA methylation at the cg11702988 site (ATP11A gene) positively correlated with lung function and BMI, and negatively correlated with lung disease severity, P. aeruginosa chronic infection, and the number of exacerbations. These results were replicated in prospective sputum samples collected at four time points within an 18-month period and longitudinally. To conclude, (i) we identified a DNA methylation biomarker that correlates with CF severity, (ii) we provided a method to easily assess this biomarker, and (iii) we carried out the first longitudinal analysis of DNA methylation in CF patients. This new epigenetic biomarker could be used to stratify CF patients in clinical trials.

International audience; The last decade has seen a significant increase of cell therapy protocols using effector T cells (Teffs) in particular, but also, more recently, non-engineered and expanded polyclonal regulatory T cells (Tregs) to control pathological immune responses such as cancer, autoimmune diseases, or transplantation rejection. However, limitations, such as stability, migration, and specificity of the cell products, have been seen. Thus, genetic engineering of these cell subsets is expected to provide the next generation of T cell therapy products. Lentiviral vectors are commonly used to modify Teffs; however, Tregs are more sensitive to mechanical stress and require specific culture conditions. Also, there is a lack of reproducible and efficient protocols to expand and genetically modify Tregs without affecting their growth and function. Due to smaller number of cells and poorer viability upon culture in vitro, mouse Tregs are more difficult to transduce and amplify in vitro than human Tregs. Here we propose a step-by-step protocol to produce both human and mouse genetically modified CD8+ and CD4+ Tregs in sufficient amounts to assess their therapeutic efficacy in humanized immunocompromised mouse models and murine models of disease and to establish pre-clinical proofs of concept. We report, for the first time, an efficient and reproducible method to isolate Tregs from human blood or mouse spleen, transduce with a lentiviral vector, and culture, in parallel, CD8+ and CD4+ Tregs while preserving their function. Beyond chimeric antigen receptor (CAR)-Treg cell therapy, this protocol will promote the development of potential new engineered T cell therapies to treat autoimmune diseases and transplantation rejection.

Hemizygous pathogenic variants in CACNA1F lead to defective signal transmission from retinal photoreceptors to bipolar cells and cause incomplete congenital stationary night blindness in humans. Although the primary defect is at the terminal end of first-order neurons (photoreceptors), there is limited knowledge of higher-order neuronal changes (inner retinal) in this disorder. This study aimed to investigate inner retinal changes in CACNA1F-retinopathy by analyzing macular ganglion cell layer-inner plexiform layer (GCL-IPL) thickness and optic disc pallor in 22 subjects with molecularly confirmed CACNA1F-retinopathy. Detailed ocular phenotypic data including distance and color vision, refraction and electroretinogram (ERG) were collected. Distance vision was universally reduced (mean: 0.42 LogMAR), six had abnormal color vision and myopia was common (n = 15, mean: −6.32 diopters). Mean GCL-IPL thickness was significantly lower in patients (55.00 µm) compared to age-matched controls (n = 87, 84.57 µm, p &lt, &lt, 0.001). The GCL-IPL thickness correlated with scotopic standard (p = 0.04) and bright-flash (p = 0.014) ERG b/a ratios and photopic b-wave amplitudes (p = 0.05). Twenty-one patients had some degree of disc pallor (bilateral in 19). Fifteen putative disease-causing, including five novel variants were identified. This study establishes macular inner retinal thinning and optic atrophy as characteristic features of CACNA1F-retinopathy, which are independent of myopia and could impact potential future treatment strategies.

Seed development needs the coordination of multiple molecular mechanisms to promote correct tissue development, seed filling and the acquisition of germination capacity, desiccation tolerance, longevity and dormancy. Heat stress can negatively impact these processes and upon the increasing of global mean temperatures, global food security is threatened. Here, we explored the impact of heat stress on seed physiology, morphology, gene expression and methylation on three stages of seed development. Notably, Arabidopsis Col-0 plants under heat stress presented a decrease in germination capacity and also a decrease in longevity. We observed that upon mild stress, gene expression and DNA methylation were moderately affected. Nevertheless, upon severe heat stress during seed development, gene expression was intensively modified, promoting heat stress response mechanisms, including the activation of ABA pathway. By analyzing candidate epigenetic marks using mutants’ physiological assays, we observed that the lack of DNA demethylation by ROS1 gene impaired seed germination by affecting germination-related genes expression. On the other hand, we also observed that upon severe stress, a large proportion of differentially methylated regions (DMRs) were located in promoters and gene sequences of germination-related genes. To conclude, our results indicate that DNA (de)methylation could be a key regulatory process to ensure proper seed germination of seeds produced under heat stress.Graphic Abstract

Populations in North Africa (NA) are characterized by a high rate of consanguinity. Consequently, the proportion of founder mutations might be higher than expected and could be a major cause for the high prevalence of recessive genetic disorders like Fanconi anemia (FA). We report clinical, cytogenetic, and molecular characterization ofFANCAin 29 North African FA patients from Tunisia, Libya, and Algeria. Cytogenetic tests revealed high rates of spontaneous chromosome breakages for all patients except two of them.FANCAmolecular analysis was performed using three different molecular approaches which allowed us to identify causal mutations as homozygous or compound heterozygous forms. It included a nonsense mutation (c.2749C > T; p.Arg917Ter), one reported missense mutation (c.1304G > A; p.Arg435His), a novel missense variant (c.1258G > A; p.Asp409Glu), and theFANCAmost common reported mutation (c.3788_3790delTCT; p.Phe1263del). Furthermore, three founder mutations were identified in 86.7% of the 22 Tunisian patients: (1) a deletion of exon 15, in 36.4% patients (8/22); (2), a deletion of exons 4 and 5 in 23% (5/22) and (3) an intronic mutation c.2222 + 166G > A, in 27.3% (6/22). Despite the relatively small number of patients studied, our results depict the mutational landscape of FA among NA populations and it should be taken into consideration for appropriate genetic counseling.

Inherited metabolic disorders (IMDs) can present with psychiatric signs that vary widely from one disease to another. This picture is further complicated by the fact that these features occur at very different illness time points, which may further delay appropriate diagnosis and treatment. In this case series of 62 children and adolescents suffering from IMDs, we clustered psychiatric signs (on the basis of the fifth edition of the Diagnostic and Statistical Manual for Mental Disorders classification) as well as impaired cognitive domains (on the basis of the Research Domain Criteriamatrix) according to their mean age of onset (5.7 ± 4 years). We observed consistent patterns of occurrence across disorders. Externalizing symptoms, sleep problems, and cross‐domain self‐regulation deficits were found to precede the IMD diagnosis. Repetitive thoughts and behaviors as well as emotional dysregulation were found to occur around the disease onset. Finally, late‐onset features included dissociative or eating disorders, together with impaired emotion knowledge. Clinicians should specifically look for the co‐occurrence of age‐specific atypical signs, such as treatment resistance or worsening with psychotropic medication in the earliest stages and symptom fluctuation, confusion, catatonia, or isolated visual hallucinations. We believe that the combined characterizations of psychiatric signs and impaired neurocognitive domains may enable the earliest detection of IMDs and the appropriate care of these particular manifestations. Key Points Psychiatric signs are common in inherited metabolic disorders (IMDs) and may occur in the same age‐range as other clinical manifestations. Three clusters of psychiatric signs and two clusters of neurocognitive domains can be defined according to their mean age of onset. Warning signs to be used in liaison psychiatry should include age‐specific cognitive impairments

Abstract Background In fish, minimally invasive blood sampling is widely used to monitor physiological stress with blood plasma biomarkers. As fish blood cells are nucleated, they might be a source a potential new markers derived from ‘omics technologies. We modified the epiGBS (epiGenotyping By Sequencing) technique to explore changes in genome-wide cytosine methylation in the red blood cells (RBCs) of challenged European sea bass (Dicentrarchus labrax), a species widely studied in both natural and farmed environments. Results We retrieved 501,108,033 sequencing reads after trimming, with a mean mapping efficiency of 73.0% (unique best hits). Minor changes in RBC methylome appeared to manifest after the challenge test and a family-effect was detected. Only fifty-seven differentially methylated cytosines (DMCs) close to 51 distinct genes distributed on 17 of 24 linkage groups (LGs) were detected between RBCs of pre- and post-challenge individuals. Thirty-seven of these genes were previously reported as differentially expressed in the brain of zebrafish, most of them involved in stress coping differences. While further investigation remains necessary, few DMC-related genes associated to the Brain Derived Neurotrophic Factor, a protein that favors stress adaptation and fear memory, appear relevant to integrate a centrally produced stress response in RBCs. Conclusion Our modified epiGBS protocol was powerful to analyze patterns of cytosine methylation in RBCs of D. labrax and to evaluate the impact of a challenge using minimally invasive blood samples. This study is the first approximation to identify epigenetic biomarkers of exposure to stress in fish.

International audience; The mammary gland undergoes important anatomical and physiological changes from embryogenesis through puberty, pregnancy, lactation and involution. These steps are under the control of a complex network of molecular factors, in which epigenetic mechanisms play a role that is increasingly well described. Recently, studies investigating epigenetic modifications and their impacts on gene expression in the mammary gland have been performed at different physiological stages and in different mammary cell types. This has led to the establishment of a role for epigenetic marks in milk component biosynthesis. This review aims to summarize the available knowledge regarding the involvement of the four main molecular mechanisms in epigenetics: DNA methylation, histone modifications, polycomb protein activity and non-coding RNA functions.

Accurately measuring the phenotype at the individual level is critical to the success of selective breeding programs. Feed efficiency is a key sustainability trait and is typically approached through feed conversion ratio (FCR). This requires measurements of body weight gain (BWG) and feed intake (FI), the latter being technically challenging in fish. We assessed two of the principal methods for measuring feed intake in fish over consecutive days: (1) group rearing 10 fish per group and video recording the meals and (2) rearing fish individually on a restricted ration. Juvenile Nile tilapia (Oreochromis niloticus) from the Genetically Improved Farmed Tilapia (GIFT) strain and the Cirad strain were entered into the study (128 GIFT and 109 Cirad). The GIFT strain were reared over three consecutive periods of 7 days each under different feeding, recording, and rearing scenarios (i) in groups fed an optimal ration (g100) or (ii) fed a 50% restricted ration (g50) both with video records of all meals and (iii) reared in isolation and fed a 50% restrictive ration. The Cirad strain were tested similarly but only for scenarios (i) and (iii). All fish were fed twice daily with a calculated ration. Correlations showed the same trends for the GIFT and the Cirad strains. For the GIFT strain, correlations were positive and significant for BWG and FI measured in scenarios (i) and (ii), 0.49 and 0.63, respectively, and FI measured in scenarios (i) and (iii) (0.50) but not for BWG measured in scenarios (i) and (iii) (0.29, NS). The phenotypic correlation estimated for FCR between scenarios (i) and (ii) with fish fed an optimal or a 50% restricted ration was low and not significant (0.22). Feed Conversion Ratio for GIFT fish reared in groups or in isolation and fed with a restricted ration [scenarios (ii) and (iii)] were not significantly correlated either. Social interactions between fish, potentially impacting their efficiency, may explain the results. Therefore, selective breeding programs seeking to improve feed efficiency will need to carefully plan the feeding rate and the rearing system used to estimate FCR in order to optimize selection for the targeted production system.

Background Variation in locomotor capacity among animals often reflects adaptations to different environments. Despite evidence that physical performance is heritable, the molecular basis of locomotor performance and performance trade-offs remains poorly understood. In this study we identify the genes, signaling pathways, and regulatory processes possibly responsible for the trade-off between burst performance and endurance observed in Xenopus allofraseri, using a transcriptomic approach. Results We obtained a total of about 121 million paired-end reads from Illumina RNA sequencing and analyzed 218,541 transcripts obtained from a de novo assembly. We identified 109 transcripts with a significant differential expression between endurant and burst performant individuals (FDR ≤ 0.05 and logFC ≥2), and blast searches resulted in 103 protein-coding genes. We found major differences between endurant and burst-performant individuals in the expression of genes involved in the polymerization and ATPase activity of actin filaments, cellular trafficking, proteoglycans and extracellular proteins secreted, lipid metabolism, mitochondrial activity and regulators of signaling cascades. Remarkably, we revealed transcript isoforms of key genes with functions in metabolism, apoptosis, nuclear export and as a transcriptional corepressor, expressed in either burst-performant or endurant individuals. Lastly, we find two up-regulated transcripts in burst-performant individuals that correspond to the expression of myosin-binding protein C fast-type (mybpc2). This suggests the presence of mybpc2 homoeologs and may have been favored by selection to permit fast and powerful locomotion. Conclusion These results suggest that the differential expression of genes belonging to the pathways of calcium signaling, endoplasmic reticulum stress responses and striated muscle contraction, in addition to the use of alternative splicing and effectors of cellular activity underlie locomotor performance trade-offs. Ultimately, our transcriptomic analysis offers new perspectives for future analyses of the role of single nucleotide variants, homoeology and alternative splicing in the evolution of locomotor performance trade-offs. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07517-1.

Background Streptococcus pneumoniae serotype 1 remains a leading cause of invasive pneumococcal diseases, even in countries with PCV-10/PCV-13 vaccine implementation. The main objective of this study, which is part of the Pneumococcal African Genome project (PAGe), was to determine the phylogenetic relationships of serotype 1 isolates recovered from children patients in Casablanca (Morocco), compared to these from other African countries; and to investigate the contribution of accessory genes and recombination events to the genetic diversity of this serotype. Results The genome average size of the six-pneumococcus serotype 1 from Casablanca was 2,227,119 bp, and the average content of coding sequences was 2113, ranging from 2041 to 2161. Pangenome analysis of the 80 genomes used in this study revealed 1685 core genes and 1805 accessory genes. The phylogenetic tree based on core genes and the hierarchical bayesian clustering analysis revealed five sublineages with a phylogeographic structure by country. The Moroccan strains cluster in two different lineages, the five invasive strains clusters altogether in a divergent clade distantly related to the non-invasive strain, that cluster with all the serotype 1 genomes from Africa. Conclusions The whole genome sequencing provides increased resolution analysis of the highly virulent serotype 1 in Casablanca, Morocco. Our results are concordant with previous works, showing that the phylogeography of S. pneumoniae serotype 1 is structured by country, and despite the small size (six isolates) of the Moroccan sample, our analysis shows the genetic cohesion of the Moroccan invasive isolates.

Recently, a disorder caused by the heterozygous de novo c.1267C>T (p.R423*) substitution in SLC37A4 has been described. This causes mislocalization of the glucose-6-phosphate transporter to the Golgi leading to a congenital disorder of glycosylation type II (SLC37A4-CDG). Only one patient has been reported showing liver disease that improved with age and mild dysmorphism. Here we report the second patient with a type II CDG caused by the same heterozygous de novo c.1267C>T (p.R423*) mutation thereby confirming the pathogenicity of this variant and expanding the clinical picture with type 1 diabetes, severe scoliosis, and membranoproliferative glomerulonephritis. Additional clinical and biochemical data provide further insight into the mechanism and prognosis of SLC37A4-CDG. keywords: CDG, G6PT1, glycogen storage disease, glycosylation, hepatopathy, SLC37A4 ispartof: JIMD Reports vol:58 issue:1 pages:122-128 ispartof: location:United States status: published

Gene expression dynamics is a key component of polyploid evolution, varying in nature, intensity, and temporal scales, most particularly in allopolyploids, where two or more sub-genomes from differentiated parental species and different repeat contents are merged. Here, we investigated transcriptome evolution at different evolutionary time scales among tetraploid, hexaploid, and neododecaploidSpartinaspecies (Poaceae, Chloridoideae) that successively diverged in the last 6–10 my, at the origin of differential phenotypic and ecological traits. Of particular interest are the recent (19th century) hybridizations between the two hexaploidsSpartina alterniflora(2n= 6x= 62) andS. maritima(2n= 6x= 60) that resulted in two sterile F1 hybrids:Spartina×townsendii(2n= 6x= 62) in England andSpartina×neyrautii(2n= 6x= 62) in France. Whole genome duplication ofS.×townsendiigave rise to the invasive neo-allododecaploid speciesSpartina anglica(2n= 12x= 124). New transcriptome assemblies and annotations for tetraploids and the enrichment of previously published reference transcriptomes for hexaploids and the allododecaploid allowed identifying 42,423 clusters of orthologs and distinguishing 21 transcribed transposable element (TE) lineages across the seven investigatedSpartinaspecies. In 4xand 6xmesopolyploids, gene and TE expression changes were consistent with phylogenetic relationships and divergence, revealing weak expression differences in the tetraploid sister speciesSpartina bakeriandSpartina versicolor(S. alternifloraandS. maritimathat diverged 2–4 mya. Differentially expressed genes were involved in glycolysis, post-transcriptional protein modifications, epidermis development, biosynthesis of carotenoids. Most detected TE lineages (exceptSINEelements) were found more expressed in hexaploids than in tetraploids, in line with their abundance in the corresponding genomes. Comparatively, an astonishing (52%) expression repatterning and deviation from parental additivity were observed following recent reticulate evolution (involving the F1 hybrids and the neo-allododecaploidS. anglica), with various patterns of biased homoeologous gene expression, including genes involved in epigenetic regulation. Downregulation of TEs was observed in both hybrids and accentuated in the neo-allopolyploid. Our results reinforce the view that allopolyploidy represents springboards to new regulatory patterns, offering to worldwide invasive species, such asS. anglica, the opportunity to colonize stressful and fluctuating environments on saltmarshes.

Today, it is common knowledge that environmental factors can change the color of many animals. Studies have shown that the molecular mechanisms underlying such modifications could involve epigenetic factors. Since 2013, the pearl oyster Pinctada margaritifera var. cumingii has become a biological model for questions on color expression and variation in Mollusca. A previous study reported color plasticity in response to water depth variation, specifically a general darkening of the nacre color at greater depth. However, the molecular mechanisms behind this plasticity are still unknown. In this paper, we investigate the possible implication of epigenetic factors controlling shell color variation through a depth variation experiment associated with a DNA methylation study performed at the whole genome level with a constant genetic background. Our results revealed six genes presenting differentially methylated CpGs in response to the environmental change, among which four are linked to pigmentation processes or regulations (GART, ABCC1, MAPKAP1, and GRL101), especially those leading to darker phenotypes. Interestingly, the genes perlucin and MGAT1, both involved in the biomineralization process (deposition of aragonite and calcite crystals), also showed differential methylation, suggesting that a possible difference in the physical/spatial organization of the crystals could cause darkening (iridescence or transparency modification of the biomineral). These findings are of great interest for the pearl production industry, since wholly black pearls and their opposite, the palest pearls, command a higher value on several markets. They also open the route of epigenetic improvement as a new means for pearl production improvement.

Background In recent years, the breeding of honeybees has gained significant scientific interest, and numerous theoretical and practical improvements have been made regarding the collection and processing of their performance data. It is now known that the selection of high-quality drone material is crucial for mid to long-term breeding success. However, there has been no conclusive mathematical theory to explain these findings. Methods We derived mathematical formulas to describe the response to selection of a breeding population and an unselected passive population of honeybees that benefits indirectly from genetic improvement in the breeding population via migration. This was done under the assumption of either controlled or uncontrolled mating of queens in the breeding population. Results Our model equations confirm what has been observed in simulation studies. In particular, we have proven that the breeding population and the passive population will show parallel genetic gain after some years and we were able to assess the responses to selection for different breeding strategies. Thus, we confirmed the crucial importance of controlled mating for successful honeybee breeding. When compared with data from simulation studies, the derived formulas showed high coefficients of determination \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$> 0.95$$\end{document}>0.95 in cases where many passive queens had dams from the breeding population. For self-sufficient passive populations, the coefficients of determination were lower (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 0.8$$\end{document}∼0.8) if the breeding population was under controlled mating. This can be explained by the limited simulated time-frame and lower convergence rates. Conclusion The presented theoretical derivations allow extrapolation of honeybee-specific simulation results for breeding programs to a wide range of population parameters. Furthermore, they provide general insights into the genetic dynamics of interdependent populations, not only for honeybees but also in a broader context.

International audience; Nuclear Receptors (NRs) are a superfamily of transcription factors specific to metazoans that have the unique ability to directly translate the message of a signaling molecule into a transcriptional response. In vertebrates, NRs are pivotal players in countless processes of both embryonic and adult physiology, with embryonic development being one of the most dynamic periods of NR activity. Accumulating evidence suggests that NR signaling is also a major regulator of development in marine invertebrates, although ligands and transactivation dynamics are not necessarily conserved with respect to vertebrates. The explosion of genome sequencing projects and the interpretation of the resulting data in a phylogenetic context allowed significant progress toward an understanding of NR superfamily evolution, both in terms of molecular activities and developmental functions. In this context, marine invertebrates have been crucial for characterizing the ancestral states of NR-ligand interactions, further strengthening the importance of these organisms in the field of evolutionary developmental biology.

International audience; The shell color of the Mollusca has attracted naturalists and collectors for hundreds of years, while the molecular pathways regulating pigment production and the pigments themselves remain poorly described. In this study, our aim was to identify the main pigments and their molecular pathways in the pearl oyster Pinctada margaritifera—the species displaying the broadest range of colors. Three inner shell colors were investigated—red, yellow, and green. To maximize phenotypic homogeneity, a controlled population approach combined with common garden conditioning was used. Comparative analysis of transcriptomes (RNA-seq) of P. margaritifera with different shell colors revealed the central role of the heme pathway, which is involved in the production of red (uroporphyrin and derivates), yellow (bilirubin), and green (biliverdin and cobalamin forms) pigments. In addition, the Raper–Mason, and purine metabolism pathways were shown to produce yellow pigments (pheomelanin and xanthine) and the black pigment eumelanin. The presence of these pigments in pigmented shell was validated by Raman spectroscopy. This method also highlighted that all the identified pathways and pigments are expressed ubiquitously and that the dominant color of the shell is due to the preferential expression of one pathway compared with another. These pathways could likely be extrapolated to many other organisms presenting broad chromatic variation

Background Knowledge about potential functional relationships among traits of interest offers a unique opportunity to understand causal mechanisms and to optimize breeding goals, management practices, and prediction accuracy. In this study, we inferred the phenotypic causal networks among five traits in a turkey population and assessed the effect of the use of such causal structures on the accuracy of predictions of breeding values. Methods Phenotypic data on feed conversion ratio, residual feed intake, body weight, breast meat yield, and walking score in addition to genotype data from a commercial breeding population were used. Causal links between the traits were detected using the inductive causation algorithm based on the joint distribution of genetic effects obtained from a standard Bayesian multiple trait model. Then, a structural equation model was implemented to infer the magnitude of causal structure coefficients among the phenotypes. Accuracies of predictions of breeding values derived using pedigree- and blending-based multiple trait models were compared to those obtained with the pedigree- and blending-based structural equation models. Results In contrast to the two unconditioned traits (i.e., feed conversion ratio and breast meat yield) in the causal structures, the three conditioned traits (i.e., residual feed intake, body weight, and walking score) showed noticeable changes in estimates of genetic and residual variances between the structural equation model and the multiple trait model. The analysis revealed interesting functional associations and indirect genetic effects. For example, the structural coefficient for the path from body weight to walking score indicated that a 1-unit genetic improvement in body weight is expected to result in a 0.27-unit decline in walking score. Both structural equation models outperformed their counterpart multiple trait models for the conditioned traits. Applying the causal structures led to an increase in accuracy of estimated breeding values of approximately 7, 6, and 20% for residual feed intake, body weight, and walking score, respectively, and different rankings of selection candidates for the conditioned traits. Conclusions Our results suggest that structural equation models can improve genetic selection decisions and increase the prediction accuracy of breeding values of selection candidates. The identified causal relationships between the studied traits should be carefully considered in future turkey breeding programs.

MPE and SERO were funded by MINECO grant AGL2016-78709-R and from the EU through the BFU2016-77236-P (MINECO/AEI/FEDER, EU) and the “Centro de Excelencia Severo Ochoa 2016-2019″ award SEV-2015-0533. LCRA is funded by “Don Carlos Antonio López” Graduate program (BECAL) from Paraguay. YRC was funded by the European Union H2020 Research and Innovation programme under Marie Skłodowska-Curie grant (P-Sphere) agreement No 6655919. This work was partially funded by the Animal Genetics Division of INRAE as part of the CHACUBA project. Y. Rodriguez-Valera was funded by the France Embassy in Havana (Cuba) and INRAE.

Cas13a, an effector of type VI CRISPR-Cas systems, is an RNA guided RNase with multiplexing and therapeutic potential. This study employs the Leptotrichia shahii (Lsh) Cas13a and a repeat-based CRISPR RNA (crRNA) to track and eliminate toxic RNA aggregates in myotonic dystrophy type 1 (DM1) – a neuromuscular disease caused by CTG expansion in the DMPK gene. We demonstrate that LshCas13a cleaves CUG repeat RNA in biochemical assays and reduces toxic RNA load in patient-derived myoblasts. As a result, LshCas13a reverses the characteristic adult-to-embryonic missplicing events in several key genes that contribute to DM1 phenotype. The deactivated LshCas13a can further be repurposed to track RNA-rich organelles within cells. Our data highlights the reprogrammability of LshCas13a and the possible use of Cas13a to target expanded repeat sequences in microsatellite expansion diseases.

Background Carbapenem-resistant Enterobacteriaceae are considered by WHO as “critical” priority pathogens for which novel antibiotics are urgently needed. The dissemination of carbapenemase-producing Escherichia coli (CP-Ec) in the community is a major public health concern. However, the global molecular epidemiology of CP-Ec isolates remains largely unknown as well as factors contributing to the acquisition of carbapenemase genes. Methods We first analyzed the whole-genome sequence and the evolution of the E. coli sequence type (ST) 410 and its disseminated clade expressing the carbapenemase OXA-181. We reconstructed the phylogeny of 19 E. coli ST enriched in CP-Ec and corresponding to a total of 2026 non-redundant isolates. Using the EpiCs software, we determined the significance of the association between specific mutations and the acquisition of a carbapenemase gene and the most probable order of events. The impact of the identified mutations was assessed experimentally by genetic manipulations and phenotypic testing. Results In 13 of the studied STs, acquisition of carbapenemase genes occurred in multidrug-resistant lineages characterized by a combination of mutations in ftsI encoding the penicillin-binding protein 3 and in the porin genes ompC and ompF. Mutated ftsI genes and a specific ompC allele related to that from ST38 inducing reduced susceptibility to diverse β-lactams spread across the species by recombination. We showed that these mutations precede in most cases the acquisition of a carbapenemase gene. The ompC allele from ST38 might have contributed to the selection of CP-Ec disseminated lineages within this ST. On the other hand, in the pandemic ST131 lineage, CP-Ec were not associated with mutations in ompC or ftsI and show no signs of dissemination. Conclusions Lineages of CP-Ec have started to disseminate globally. However, their selection is a multistep process involving mutations, recombination, acquisition of antibiotic resistance genes, and selection by β-lactams from diverse families. This process did not yet occur in the high-risk lineage ST131.

Background Thousands of years of natural and artificial selection have resulted in indigenous cattle breeds that are well-adapted to the environmental challenges of their local habitat and thereby are considered as valuable genetic resources. Understanding the genetic background of such adaptation processes can help us design effective breeding objectives to preserve local breeds and improve commercial cattle. To identify regions under putative selection, GGP HD 150 K single nucleotide polymorphism (SNP) arrays were used to genotype 106 individuals representing five Swedish breeds i.e. native to different regions and covering areas with a subarctic cold climate in the north and mountainous west, to those with a continental climate in the more densely populated south regions. Results Five statistics were incorporated within a framework, known as de-correlated composite of multiple signals (DCMS) to detect signatures of selection. The obtained p-values were adjusted for multiple testing (FDR DCAF8, PPP1R12A, SLC16A3, UCP2, UCP3, TIGAR), cold acclimation (AQP3, AQP7, HSPB8), body size and stature (PLAG1, KCNA6, NDUFA9, AKAP3, C5H12orf4, RAD51AP1, FGF6, TIGAR, CCND2, CSMD3), resistance to disease and bacterial infection (CHI3L2, GBP6, PPFIBP1, REP15, CYP4F2, TIGD2, PYURF, SLC10A2, FCHSD2, ARHGEF17, RELT, PRDM2, KDM5B), reproduction (PPP1R12A, ZFP36L2, CSPP1), milk yield and components (NPC1L1, NUDCD3, ACSS1, FCHSD2), growth and feed efficiency (TMEM68, TGS1, LYN, XKR4, FOXA2, GBP2, GBP5, FGD6), and polled phenotype (URB1, EVA1C). Conclusions We identified genomic regions that may provide background knowledge to understand the mechanisms that are involved in economic traits and adaptation to cold climate in cattle. Incorporating p-values of different statistics in a single DCMS framework may help select and prioritize candidate genes for further analyses.

Background Positional weight matrix (PWM) is a de facto standard model to describe transcription factor (TF) DNA binding specificities. PWMs inferred from in vivo or in vitro data are stored in many databases and used in a plethora of biological applications. This calls for comprehensive benchmarking of public PWM models with large experimental reference sets. Results Here we report results from all-against-all benchmarking of PWM models for DNA binding sites of human TFs on a large compilation of in vitro (HT-SELEX, PBM) and in vivo (ChIP-seq) binding data. We observe that the best performing PWM for a given TF often belongs to another TF, usually from the same family. Occasionally, binding specificity is correlated with the structural class of the DNA binding domain, indicated by good cross-family performance measures. Benchmarking-based selection of family-representative motifs is more effective than motif clustering-based approaches. Overall, there is good agreement between in vitro and in vivo performance measures. However, for some in vivo experiments, the best performing PWM is assigned to an unrelated TF, indicating a binding mode involving protein-protein cooperativity. Conclusions In an all-against-all setting, we compute more than 18 million performance measure values for different PWM-experiment combinations and offer these results as a public resource to the research community. The benchmarking protocols are provided via a web interface and as docker images. The methods and results from this study may help others make better use of public TF specificity models, as well as public TF binding data sets.

© The Author(s). 2020. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data., Background: How intestinal epithelial cells interact with the microbiota and how this is regulated at the gene expression level are critical questions. Smarcad1 is a conserved chromatin remodeling factor with a poorly understood tissue function. As this factor is highly expressed in the stem and proliferative zones of the intestinal epithelium, we explore its role in this tissue. Results: Specific deletion of Smarcad1 in the mouse intestinal epithelium leads to colitis resistance and substantial changes in gene expression, including a striking increase of expression of several genes linked to innate immunity. Absence of Smarcad1 leads to changes in chromatin accessibility and significant changes in histone H3K9me3 over many sites, including genes that are differentially regulated upon Smarcad1 deletion. We identify candidate members of the gut microbiome that elicit a Smarcad1-dependent colitis response, including members of the poorly understood TM7 phylum. Conclusions: Our study sheds light onto the role of the chromatin remodeling machinery in intestinal epithelial cells in the colitis response and shows how a highly conserved chromatin remodeling factor has a distinct role in anti-microbial defense. This work highlights the importance of the intestinal epithelium in the colitis response and the potential of microbial species as pharmacological and probiotic targets in the context of inflammatory diseases., This work was supported by BBSRC Institute Strategic Programme Grants (BBS/E/B/000C0404, BBS/E/B/000C0405), a MRC project grant to PVW and MV (MR/N009398/1), a BBSRC-Fundação de Amparo à Pesquisa do Estado de São Paulo/FAPESP-Brazil Pump-priming award (BB/N013565/1,FAPESP#2015/50379-1), and a Newton Advanced Fellowship from the Royal Society to MARV and PVW (NAF\R1\180116).

International audience; Background : X-chromosomal loci present different inheritance patterns compared to autosomal loci and must be modeled accordingly. Sexual chromosomes are not systematically considered in whole-genome relationship matrices although rules based on genealogical or marker information have been derived. Loci on the X-chromosome could have a significant contribution to the additive genetic variance, in particular for some traits such as those related to reproduction. Thus, accounting for the X-chromosome relationship matrix might be informative to better understand the architecture of complex traits (e.g., by estimating the variance associated to this chromosome) and to improve their genomic prediction. For such applications, previous studies have shown the benefits of combining information from genotyped and ungenotyped individuals.Results: In this paper, we start by presenting rules to compute a genomic relationship matrix (GRM) for the X-chromosome (GX) without making any assumption on dosage compensation, and based on coding of gene content with 0/1 for males and 0/1/2 for females. This coding adjusts naturally to previously derived pedigree-based relationships (S) for the X-chromosome. When needed, we propose to accommodate and estimate dosage compensation and genetic heterogeneity across sexes via multiple trait models. Using a Holstein dairy cattle dataset, including males and females, we then empirically illustrate that realized relationships (GX) matches expectations (S). However, GX presents high deviations from S. GX has also a lower dimensionality compared to the autosomal GRM. In particular, individuals are frequently identical along the entire chromosome. Finally, we confirm that the heritability of gene content for markers on the X-chromosome that are estimated by using S is 1, further demonstrating that S and GX can be combined. For the pseudo-autosomal region, we demonstrate that the expected relationships vary according to position because of the sex-gradient. We end by presenting the rules to construct the 'H matrix’ by combining both relationship matrices.Conclusions: This work shows theoretically and empirically that a pedigree-based relationship matrix built with rules specifically developed for the X-chromosome (S) matches the realized GRM for the X-chromosome. Therefore, applications that combine expected relationships and genotypes for markers on the X-chromosome should use S and GX.

BackgroundPolyploidy is widespread in animals and especially in plants. Different kinds of ploidies exist, for example, hexaploidy in wheat, octaploidy in strawberries, and diploidy, triploidy, tetraploidy, and pseudo-tetraploidy (partly tetraploid) in fish. Triploid offspring from diploid parents occur frequently in the wild in Atlantic salmon (Salmo salar) and, as with triploidy in general, the triploid individuals are sterile. Induced triploidy in Atlantic salmon is common practice to produce sterile fish. In Norwegian aquaculture, production of sterile triploid fish is an attempt by government and industry to limit genetic introgression between wild and farmed fish. However, triploid fish may have traits and properties that differ from those of diploids. Investigating the genetics behind traits in triploids has proved challenging because genotype calling of genetic markers in triploids is not supported by standard software. Our aim was to develop a method that can be used for genotype calling of genetic markers in triploid individuals.ResultsAllele signals were produced for 381 triploid Atlantic salmon offspring using a 56 K Thermo Fisher GeneTitan genotyping platform. Genotypes were successfully called by applying finite normal mixture models to the (transformed) allele signals. Subsets of markers were filtered by quality control statistics for use with downstream analyses. The quality of the called genotypes was sufficient to allow for assignment of diploid parents to the triploid offspring and to discriminate between maternal and paternal parents from autosomal inheritance patterns. In addition, as the maternal inheritance in triploid offspring is identical to gynogenetic inheritance, the maternal recombination pattern for each chromosome could be mapped by using a similar approach as that used in gene-centromere mapping.ConclusionsWe show that calling of dense marker genotypes for triploid individuals is feasible. The resulting genotypes can be used in parentage assignment of triploid offspring to diploid parents, to discriminate between maternal and paternal parents using autosomal inheritance patterns, and to map the maternal recombination pattern using an approach similar to gene-centromere mapping. Genotyping of triploid individuals is important both for selective breeding programs and unravelling the underlying genetics of phenotypes recorded in triploids. In principle, the developed method can be used for genotype calling of other polyploid organisms.

Drosophila has been studied as a biological model for many years and many discoveries in biology rely on this species. Research on transposable elements (TEs) is not an exception. Drosophila has contributed significantly to our knowledge on the mechanisms of transposition and their regulation, but above all, it was one of the first organisms on which genetic and genomic studies of populations were done. In this review article, in a very broad way, we will approach the TEs of Drosophila with a historical hindsight as well as recent discoveries in the field.

Background Cattle international genetic evaluations allow the comparison of estimated breeding values (EBV) across different environments, i.e. countries. For international evaluations, across-country genetic correlations (rg) need to be estimated. However, lack of convergence of the estimated parameters and high standard errors of the rg are often experienced for beef cattle populations due to limited across-country genetic connections. Furthermore, using all available genetic connections to estimate rg is prohibitive due to computational constraints, thus sub-setting the data is necessary. Our objective was to investigate and compare the impact of strategies of data sub-setting on estimated across-country rg and their computational requirements. Methods Phenotype and pedigree information for age-adjusted weaning weight was available for ten European countries and 3,128,338 Limousin beef cattle males and females. Using a Monte Carlo based expectation–maximization restricted maximum likelihood (MC EM REML) methodology, we estimated across-country rg by using a multi-trait animal model where countries are modelled as different correlated traits. Values of rg were estimated using the full data and four different sub-setting strategies that aimed at selecting the most connected herds from the largest population. Results Using all available data, direct and maternal rg (standard errors in parentheses) were on average equal to 0.79 (0.14) and 0.71 (0.19), respectively. Direct-maternal within-country and between-country rg were on average equal to − 0.12 (0.09) and 0.00 (0.14), respectively. Data sub-setting scenarios gave similar results: on average, estimated rg were smaller compared to using all data for direct (0.02) and maternal (0.05) genetic effects. The largest differences were obtained for the direct-maternal within-country and between-country rg, which were, on average 0.13 and 0.12 smaller compared to values obtained by using all data. Standard errors always increased when reducing the data, by 0.02 to 0.06, on average. The proposed sub-setting strategies reduced the required computing time up to 22% compared to using all data. Conclusions Estimating all 120 across-country rg that are required for beef cattle international evaluations, using a multi-trait MC EM REML approach, is feasible but involves long computing time. We propose four strategies to reduce computational requirements while keeping a multi-trait estimation approach. In all scenarios with data sub-setting, the estimated rg were consistently smaller (mainly for direct-maternal rg) and had larger standard errors.

Background Bisphenol S (BPS) is a common bisphenol A (BPA) substitute, since BPA is virtually banned worldwide. However, BPS and BPA have both endocrine disrupting properties. Their effects appear mostly in adulthood following perinatal exposures. The objective of the present study was to investigate the impact of perinatal and chronic exposure to BPS at the low dose of 1.5 μg/kg body weight/day on the transcriptome and methylome of the liver in 23 weeks-old C57BL6/J male mice. Results This multi-omic study highlights a major impact of BPS on gene expression (374 significant deregulated genes) and Gene Set Enrichment Analysis show an enrichment focused on several biological pathways related to metabolic liver regulation. BPS exposure also induces a hypomethylation in 58.5% of the differentially methylated regions (DMR). Systematic connections were not found between gene expression and methylation profile excepted for 18 genes, including 4 genes involved in lipid metabolism pathways (Fasn, Hmgcr, Elovl6, Lpin1), which were downregulated and featured differentially methylated CpGs in their exons or introns. Conclusions This descriptive study shows an impact of BPS on biological pathways mainly related to an integrative disruption of metabolism (energy metabolism, detoxification, protein and steroid metabolism) and, like most high-throughput studies, contributes to the identification of potential exposure biomarkers. Graphical abstract Supplementary Information The online version contains supplementary material available at 10.1186/s12864-020-07294-3.

Background Native Swedish sheep breeds are part of the North European short-tailed sheep group; characterized in part by their genetic uniqueness. Our objective was to study the population structure of native Swedish sheep. Five breeds were genotyped using the 600 K SNP array. Dalapäls and Klövsjö sheep are from the middle of Sweden; Gotland and Gute sheep from Gotland, an island in the Baltic Sea; and Fjällnäs sheep from northern Sweden. We studied population structure by: principal component analysis (PCA), cluster-based analysis of admixture, and an estimated population tree. Results The analyses of the five Swedish breeds revealed that these breeds are five distinct breeds, while Gute and Gotland are more closely related to each other as seen in all analyses. All breeds had long branch lengths in the population tree indicating they’ve been subjected to drift. We repeated our analyses using 39 K SNP and including 50 K SNP genotypes from other European and southwestern Asian breeds from the Sheep HapMap project and 600 K SNP genotypes from a dataset of French sheep. Results arranged breeds into five groups: south-west Asia, south-west Europe, central Europe, north Europe and north European short-tailed sheep. Within this last group, Norwegian and Icelandic breeds, Finn and Romanov sheep, Scottish breeds, and Gute and Gotland sheep were more closely related while the remaining Swedish breeds and Ouessant sheep were distinct from all breeds and had longer branches in the population tree. Conclusions We showed population structure of five Swedish breeds and their structure within European and southwestern Asian breeds. Swedish breeds are unique, distinct breeds that have been subjected to drift but group with other north European short-tailed sheep.

BACKGROUND: We describe the latest improvements to the long-range phasing (LRP) and haplotype library imputation (HLI) algorithms for successful phasing of both datasets with one million individuals and datasets genotyped using different sets of single nucleotide polymorphisms (SNPs). Previous publicly available implementations of the LRP algorithm implemented in AlphaPhase could not phase large datasets due to the computational cost of defining surrogate parents by exhaustive all-against-all searches. Furthermore, the AlphaPhase implementations of LRP and HLI were not designed to deal with large amounts of missing data that are inherent when using multiple SNP arrays.METHODS: We developed methods that avoid the need for all-against-all searches by performing LRP on subsets of individuals and then concatenating the results. We also extended LRP and HLI algorithms to enable the use of different sets of markers, including missing values, when determining surrogate parents and identifying haplotypes. We implemented and tested these extensions in an updated version of AlphaPhase, and compared its performance to the software package Eagle2.RESULTS: A simulated dataset with one million individuals genotyped with the same 6711 SNPs for a single chromosome took less than a day to phase, compared to more than seven days for Eagle2. The percentage of correctly phased alleles at heterozygous loci was 90.2 and 99.9% for AlphaPhase and Eagle2, respectively. A larger dataset with one million individuals genotyped with 49,579 SNPs for a single chromosome took AlphaPhase 23 days to phase, with 89.9% of alleles at heterozygous loci phased correctly. The phasing accuracy was generally lower for datasets with different sets of markers than with one set of markers. For a simulated dataset with three sets of markers, 1.5% of alleles at heterozygous positions were phased incorrectly, compared to 0.4% with one set of markers.CONCLUSIONS: The improved LRP and HLI algorithms enable AlphaPhase to quickly and accurately phase very large and heterogeneous datasets. AlphaPhase is an order of magnitude faster than the other tested packages, although Eagle2 showed a higher level of phasing accuracy. The speed gain will make phasing achievable for very large genomic datasets in livestock, enabling more powerful breeding and genetics research and application.

BackgroundAPOBEC-driven mutagenesis and functional positive selection of mutated genes may synergistically drive the higher frequency of some hotspot driver mutations compared to other mutations within the same gene, as we reported forFGFR3S249C. Only a few APOBEC-associated driver hotspot mutations have been identified in bladder cancer (BCa). Here, we systematically looked for and characterised APOBEC-associated hotspots in BCa.MethodsWe analysed 602 published exome-sequenced BCas, for part of which gene expression data were also available. APOBEC-associated hotspots were identified by motif-mapping, mutation signature fitting and APOBEC-mediated mutagenesis comparison. Joint analysis of DNA hairpin stability and gene expression was performed to predict driver or passenger hotspots. Aryl hydrocarbon receptor (AhR) activity was calculated based on its target genes expression. Effects of AhR knockout/inhibition on BCa cell viability were analysed.ResultsWe established a panel of 44 APOBEC-associated hotspot mutations in BCa, which accounted for about half of the hotspot mutations. Fourteen of them overlapped with the hotspots found in other cancer types with high APOBEC activity. They mostly occurred in the DNA lagging-strand templates and the loop of DNA hairpins. APOBEC-associated hotspots presented systematically a higher prevalence than the other mutations within each APOBEC-target gene, independently of their functional impact. A combined analysis of DNA loop stability and gene expression allowed to distinguish known passenger from known driver hotspot mutations in BCa, including loss-of-function mutations affecting tumour suppressor genes, and to predict new candidate drivers, such asAHRQ383H. We further characterisedAHRQ383H as an activating driver mutation associated with high AhR activity in luminal tumours. High AhR activity was also found in tumours presenting amplifications ofAHRand its co-receptorARNT. We finally showed that BCa cells presenting those different genetic alterations were sensitive to AhR inhibition.ConclusionsOur study identified novel potential drivers within APOBEC-associated hotspot mutations in BCa reinforcing the importance of APOBEC mutagenesis in BCa. It could allow a better understanding of BCa biology and aetiology and have clinical implications such as AhR as a potential therapeutic target. Our results also challenge the dogma that all hotspot mutations are drivers and mostly gain-of-function mutations affecting oncogenes.

International audience; Prostate cancer (PCa) is a major concern in public health, with many genetically distinct subsets. Genomic alterations in PCa are extraordinarily complex, and both germline and somatic mutations are of great importance in the development of this tumor. The aim of this review is to provide an overview of genetic changes that can occur in the development of PCa and their role in potential therapeutic approaches. Various pathways and mechanisms proposed to play major roles in PCa are described in detail to provide an overview of current knowledge.

Background Production and health traits are central in cattle breeding. Advances in next-generation sequencing technologies and genotype imputation have increased the resolution of gene mapping based on genome-wide association studies (GWAS). Thus, numerous candidate genes that affect milk yield, milk composition, and mastitis resistance in dairy cattle are reported in the literature. Effect-bearing variants often affect multiple traits. Because the detection of overlapping quantitative trait loci (QTL) regions from single-trait GWAS is too inaccurate and subjective, multi-trait analysis is a better approach to detect pleiotropic effects of variants in candidate genes. However, large sample sizes are required to achieve sufficient power. Multi-trait meta-analysis is one approach to deal with this problem. Thus, we performed two multi-trait meta-analyses, one for three milk production traits (milk yield, protein yield and fat yield), and one for milk yield and mastitis resistance. Results For highly correlated traits, the power to detect pleiotropy was increased by multi-trait meta-analysis compared with the subjective assessment of overlapping of single-trait QTL confidence intervals. Pleiotropic effects of lead single nucleotide polymorphisms (SNPs) that were detected from the multi-trait meta-analysis were confirmed by bivariate association analysis. The previously reported pleiotropic effects of variants within the DGAT1 and MGST1 genes on three milk production traits, and pleiotropic effects of variants in GHR on milk yield and fat yield were confirmed. Furthermore, our results suggested that variants in KCTD16, KCNK18 and ENSBTAG00000023629 had pleiotropic effects on milk production traits. For milk yield and mastitis resistance, we identified possible pleiotropic effects of variants in two genes, GC and DGAT1. Conclusions Multi-trait meta-analysis improves our ability to detect pleiotropic interactions between milk production traits and identifies variants with pleiotropic effects on milk production traits and mastitis resistance. In particular, this should contribute to better understand the biological mechanisms that underlie the unfavorable genetic correlation between milk yield and mastitis.

Background In this study, we assessed the accuracy of genomic prediction for carcass weight (CWT), marbling score (MS), eye muscle area (EMA) and back fat thickness (BFT) in Hanwoo cattle when using genomic best linear unbiased prediction (GBLUP), weighted GBLUP (wGBLUP), and a BayesR model. For these models, we investigated the potential gain from using pre-selected single nucleotide polymorphisms (SNPs) from a genome-wide association study (GWAS) on imputed sequence data and from gene expression information. We used data on 13,717 animals with carcass phenotypes and imputed sequence genotypes that were split in an independent GWAS discovery set of varying size and a remaining set for validation of prediction. Expression data were used from a Hanwoo gene expression experiment based on 45 animals. Results Using a larger number of animals in the reference set increased the accuracy of genomic prediction whereas a larger independent GWAS discovery dataset improved identification of predictive SNPs. Using pre-selected SNPs from GWAS in GBLUP improved accuracy of prediction by 0.02 for EMA and up to 0.05 for BFT, CWT, and MS, compared to a 50 k standard SNP array that gave accuracies of 0.50, 0.47, 0.58, and 0.47, respectively. Accuracy of prediction of BFT and CWT increased when BayesR was applied with the 50 k SNP array (0.02 and 0.03, respectively) and was further improved by combining the 50 k array with the top-SNPs (0.06 and 0.04, respectively). By contrast, using BayesR resulted in limited improvement for EMA and MS. wGBLUP did not improve accuracy but increased prediction bias. Based on the RNA-seq experiment, we identified informative expression quantitative trait loci, which, when used in GBLUP, improved the accuracy of prediction slightly, i.e. between 0.01 and 0.02. SNPs that were located in genes, the expression of which was associated with differences in trait phenotype, did not contribute to a higher prediction accuracy. Conclusions Our results show that, in Hanwoo beef cattle, when SNPs are pre-selected from GWAS on imputed sequence data, the accuracy of prediction improves only slightly whereas the contribution of SNPs that are selected based on gene expression is not significant. The benefit of statistical models to prioritize selected SNPs for estimating genomic breeding values is trait-specific and depends on the genetic architecture of each trait.

Background Tilapias (Family Cichlidae) are the second most important group of aquaculture species in the world. They have been the subject of much research on sex determination due to problems caused by early maturation in culture and their complex sex-determining systems. Different sex-determining loci (linkage group 1, 20 and 23) have been detected in various tilapia stocks. The ‘genetically improved farmed tilapia’ (GIFT) stock, founded from multiple Nile tilapia (Oreochromis niloticus) populations, with some likely to have been introgressed with O. mossambicus, is a key resource for tilapia aquaculture. The sex-determining mechanism in the GIFT stock was unknown, but potentially complicated due to its multiple origins. Results A bulk segregant analysis (BSA) version of double-digest restriction-site associated DNA sequencing (BSA-ddRADseq) was developed and used to detect and position sex-linked single nucleotide polymorphism (SNP) markers in 19 families from the GIFT strain breeding nucleus and two Stirling families as controls (a single XY locus had been previously mapped to LG1 in the latter). About 1500 SNPs per family were detected across the genome. Phenotypic sex in Stirling families showed strong association with LG1, whereas only SNPs located in LG23 showed clear association with sex in the majority of the GIFT families. No other genomic regions linked to sex determination were apparent. This region was validated using a series of LG23-specific DNA markers (five SNPs with highest association to sex from this study, the LG23 sex-associated microsatellite UNH898 and ARO172, and the recently isolated amhy marker for individual fish (n = 284). Conclusions Perhaps surprisingly given its multiple origins, sex determination in the GIFT strain breeding nucleus was associated only with a locus in LG23. BSA-ddRADseq allowed cost-effective analysis of multiple families, strengthening this conclusion. This technique has potential to be applied to other complex traits. The sex-linked SNP markers identified will be useful for potential marker-assisted selection (MAS) to control sex-ratio in GIFT tilapia to suppress unwanted reproduction during growout.

Background Increased contamination of European and Asian wheat and barley crops with “emerging” mycotoxins such as enniatins or beauvericin, produced by Fusarium avenaceum and Fusarium tricinctum, suggest that these phylogenetically close species could be involved in future food-safety crises. Results The mitochondrial genomes of F. tricinctum strain INRA104 and F. avenaceum strain FaLH27 have been annotated. A comparative analysis was carried out then extended to a set of 25 wild strains. Results show that they constitute two distinct species, easily distinguished by their mitochondrial sequences. The mitochondrial genetic variability is mainly located within the intergenic regions. Marks of variations show they have evolved (i) by Single Nucleotide Polymorphisms (SNPs), (ii) by length variations mediated by insertion/deletion sequences (Indels), and (iii) by length mutations generated by DNA sliding events occurring in mononucleotide (A)n or (T)n microsatellite type sequences arranged in a peculiar palindromic organization. The optionality of these palindromes between both species argues for their mobility. The presence of Indels and SNPs in palindrome neighbouring regions suggests their involvement in these observed variations. Moreover, the intraspecific and interspecific variations in the presence/absence of group I introns suggest a high mobility, resulting from several events of gain and loss during short evolution periods. Phylogenetic analyses of intron orthologous sequences suggest that most introns could have originated from lateral transfers from phylogenetically close or distant species belonging to various Ascomycota genera and even to the Basidiomycota fungal division. Conclusions Mitochondrial genome evolution between F. tricinctum and F. avenaceum is mostly driven by two types of mobile genetic elements, implicated in genome polymorphism. The first one is represented by group I introns. Indeed, both genomes harbour optional (inter- or intra-specifically) group I introns, all carrying putatively functional hegs, arguing for a high mobility of these introns during short evolution periods. The gain events were shown to involve, for most of them, lateral transfers between phylogenetically distant species. This study has also revealed a new type of mobile genetic element constituted by a palindromic arrangement of (A) n and (T) n microsatellite sequences whose presence was related to occurrence of SNPs and Indels in the neighbouring regions.