Malignant features of minipig melanomas prior to spontaneous regression.

In MeLiM minipigs, melanomas develop around birth, can metastasize, and have histopathologic characteristics similar to humans. Interestingly, MeLiM melanomas eventually regress. This favorable outcome raises the question of their malignancy, which we investigated. We clinically followed tens of tumors from onset to first signs of regression. Transcriptome analysis revealed an enrichment of all cancer hallmarks in melanomas, although no activating or suppressing somatic mutation were found in common driver genes. Analysis of tumor cell genomes revealed high mutation rates without UV signature. Canonical proliferative, survival and angiogenic pathways were detected in MeLiM tumor cells all along progression stages. Functionally, we show that MeLiM melanoma cells are capable to grow in immunocompromised mice, with serial passages and for a longer time than in MeLiM pigs. Pigs set in place an immune response during progression with dense infiltration by myeloid cells while melanoma cells are deficient in B2M expression. To conclude, our data on MeLiM melanomas reveal several malignancy characteristics. The combination of these features with the successful spontaneous regression of these tumors make it an outstanding model to study an efficient anti-tumor immune response.

Multi-omics data integration for the identification of biomarkers for bull fertility.

Bull fertility is an important economic trait, and the use of subfertile semen for artificial insemination decreases the global efficiency of the breeding sector. Although the analysis of semen functional parameters can help to identify infertile bulls, no tools are currently available to enable precise predictions and prevent the commercialization of subfertile semen. Because male fertility is a multifactorial phenotype that is dependent on genetic, epigenetic, physiological and environmental factors, we hypothesized that an integrative analysis might help to refine our knowledge and understanding of bull fertility. We combined -omics data (genotypes, sperm DNA methylation at CpGs and sperm small non-coding RNAs) and semen parameters measured on a large cohort of 98 Montbéliarde bulls with contrasting fertility levels. Multiple Factor Analysis was conducted to study the links between the datasets and fertility. Four methodologies were then considered to identify the features linked to bull fertility variation: Logistic Lasso, Random Forest, Gradient Boosting and Neural Networks. Finally, the features selected by these methods were annotated in terms of genes, to conduct functional enrichment analyses. The less relevant features in -omics data were filtered out, and MFA was run on the remaining 12,006 features, including the 11 semen parameters and a balanced proportion of each type of-omics data. The results showed that unlike the semen parameters studied the-omics datasets were related to fertility. Biomarkers related to bull fertility were selected using the four methodologies mentioned above. The most contributory CpGs, SNPs and miRNAs targeted genes were all found to be involved in development. Interestingly, fragments derived from ribosomal RNAs were overrepresented among the selected features, suggesting roles in male fertility. These markers could be used in the future to identify subfertile bulls in order to increase the global efficiency of the breeding sector.

Bioinformatic Mining and Structure-Activity Profiling of Baeyer-Villiger Monooxygenases from Mycobacterium tuberculosis.

Mycobacterium tuberculosis is the etiological agent of tuberculosis (TB), one of the deadliest infectious diseases. The alarming health context coupled with the emergence of resistant M. tuberculosis strains highlights the urgent need to expand the range of anti-TB antibiotics. A subset of anti-TB drugs in use are prodrugs that require bioactivation by a class of M. tuberculosis enzymes called Baeyer-Villiger monooxygenases (BVMOs), which remain understudied. To examine the prevalence and the molecular function of BVMOs in mycobacteria, we applied a comprehensive bioinformatic analysis that identified six BVMOs in M. tuberculosis, including Rv3083 (MymA), Rv3854c (EthA), Rv0565c, and Rv0892, which were selected for further characterization. Homology modeling and substrate docking analysis, performed on this subset, suggested that Rv0892 is closer to the cyclohexanone BVMO, while Rv0565c and EthA are structurally and functionally similar to MymA, which is by far the most prominent type I BVMO enzyme. Thanks to an unprecedented purification and assay optimization, biochemical studies confirmed that all four BVMOs display BV-oxygenation activity. We also showed that MymA displays a distinctive substrate preference that we further investigated by kinetic parameter determination and that correlates with in silico modeling. We provide insights into distribution of BVMOs and the structural basis of their substrate profiling, and we discuss their possible redundancy in M. tuberculosis, raising questions about their versatility in prodrug activation and their role in physiology and infection. IMPORTANCE Tuberculosis (TB), caused by Mycobacterium tuberculosis, is one of the leading causes of death worldwide. The rise in drug resistance highlights the urgent need for innovation in anti-TB drug development. Many anti-TB drugs require bioactivation by Baeyer-Villiger monooxygenases (BVMOs). Despite their emerging importance, BVMO structural and functional features remain enigmatic. We applied a comprehensive bioinformatic analysis and confirmed the presence of six BVMOs in M. tuberculosis, including MymA, EthA, and Rv0565c-activators of the second-line prodrug ethionamide-and the novel BVMO Rv0892. Combining in silico characterization with in vitro validation, we outlined their structural framework and substrate preference. Markedly, MymA displayed an enhanced capacity and a distinct selectivity profile toward ligands, in agreement with its catalytic site topology. These features ground the molecular basis for structure-function comprehension of the specificity in these enzymes and expand the repertoire of BVMOs with selective and/or overlapping activity for application in the context of improving anti-TB therapy.

Dynamics of cattle sperm sncRNAs during maturation, from testis to ejaculated sperm.

BACKGROUND: During epididymal transit, spermatozoa go through several functional maturation steps, resulting from interactions with epididymal secretomes specific to each region. In particular, the sperm membrane is under constant remodeling, with sequential attachment and shedding of various molecules provided by the epididymal lumen fluid and epididymosomes, which also deliver sncRNA cargo to sperm. As a result, the payload of sperm sncRNAs changes during the transit from the epididymis caput to the cauda. This work was designed to study the dynamics of cattle sperm sncRNAs from spermatogenesis to final maturation. RESULTS: Comprehensive catalogues of sperm sncRNAs were obtained from testicular parenchyma, epididymal caput, corpus and cauda, as well as ejaculated semen from three Holstein bulls. The primary cattle sncRNA sperm content is markedly remodeled as sperm mature along the epididymis. Expression of piRNAs, which are abundant in testis parenchyma, decreases dramatically at epididymis. Conversely, sperm progressively acquires miRNAs, rsRNAs, and tsRNAs along epididymis, with regional specificities. For instance, miRNAs and tsRNAs are enriched in epididymis cauda and ejaculated sperm, while rsRNA expression peaks at epididymis corpus. In addition, epididymis corpus contains mainly 20 nt long piRNAs, instead of 30 nt in all other locations. Beyond the bulk differences in abundance of sncRNAs classes, K-means clustering was performed to study their spatiotemporal expression profile, highlighting differences in specific sncRNAs and providing insights into their putative biological role at each maturation stage. For instance, Gene Ontology analyses using miRNA targets highlighted enriched processes such as cell cycle regulation, response to stress and ubiquitination processes in testicular parenchyma, protein metabolism in epididymal sperm, and embryonic morphogenesis in ejaculated sperm. CONCLUSIONS: Our findings confirm that the sperm sncRNAome does not simply reflect a legacy of spermatogenesis. Instead, sperm sncRNA expression shows a remarkable level of plasticity resulting probably from the combination of multiple factors such as loss of the cytoplasmic droplet, interaction with epididymosomes, and more surprisingly, the putative in situ production and/or modification of sncRNAs by sperm. Given the suggested role of sncRNA in epigenetic trans-generational inheritance, our detailed spatiotemporal analysis may pave the way for a study of sperm sncRNAs role in embryo development.

A comprehensive overview of bull sperm-borne small non-coding RNAs and their diversity across breeds.

BACKGROUND: Mature sperm carry thousands of RNAs, including mRNAs, lncRNAs, tRNAs, rRNAs and sncRNAs, though their functional significance is still a matter of debate. Growing evidence suggests that sperm RNAs, especially sncRNAs, are selectively retained during spermiogenesis or specifically transferred during epididymis maturation, and are thus delivered to the oocyte at fertilization, providing resources for embryo development. However , a deep characterization of the sncRNA content of bull sperm and its expression profile across breeds is currently lacking. To fill this gap, we optimized a guanidinium-Trizol total RNA extraction protocol to prepare high-quality RNA from frozen bull sperm collected from 40 representative bulls from six breeds. Deep sequencing was performed (40 M single 50-bp reads per sample) to establish a comprehensive repertoire of cattle sperm sncRNA. RESULTS: Our study showed that it comprises mostly piRNAs (26%), rRNA fragments (25%), miRNAs (20%) and tRNA fragments (tsRNA, 14%). We identified 5p-halves as the predominant tsRNA subgroup in bull sperm, originating mostly from Gly and Glu isoacceptors. Our study also increased by ~ 50% the sperm repertoire of known miRNAs and identified 2022 predicted miRNAs. About 20% of sperm miRNAs were located within genomic clusters, expanding the list of known polycistronic pri-miRNA clusters and defining several networks of co-expressed miRNAs. Strikingly, our study highlighted the great diversity of isomiRs, resulting mainly from deletions and non-templated additions (A and U) at the 3p end. Substitutions within miRNA sequence accounted for 40% of isomiRs, with G>A, U>C and C>U substitutions being the most frequent variations. In addition, many sncRNAs were found to be differentially expressed across breeds. CONCLUSIONS: Our study provides a comprehensive overview of cattle sperm sncRNA, and these findings will pave the way for future work on the role of sncRNAs in embryo development and their relevance as biomarkers of semen fertility.

Multi-species annotation of transcriptome and chromatin structure in domesticated animals.

BACKGROUND: Comparative genomics studies are central in identifying the coding and non-coding elements associated with complex traits, and the functional annotation of genomes is a critical step to decipher the genotype-to-phenotype relationships in livestock animals. As part of the Functional Annotation of Animal Genomes (FAANG) action, the FR-AgENCODE project aimed to create reference functional maps of domesticated animals by profiling the landscape of transcription (RNA-seq), chromatin accessibility (ATAC-seq) and conformation (Hi-C) in species representing ruminants (cattle, goat), monogastrics (pig) and birds (chicken), using three target samples related to metabolism (liver) and immunity (CD4+ and CD8+ T cells). RESULTS: RNA-seq assays considerably extended the available catalog of annotated transcripts and identified differentially expressed genes with unknown function, including new syntenic lncRNAs. ATAC-seq highlighted an enrichment for transcription factor binding sites in differentially accessible regions of the chromatin. Comparative analyses revealed a core set of conserved regulatory regions across species. Topologically associating domains (TADs) and epigenetic A/B compartments annotated from Hi-C data were consistent with RNA-seq and ATAC-seq data. Multi-species comparisons showed that conserved TAD boundaries had stronger insulation properties than species-specific ones and that the genomic distribution of orthologous genes in A/B compartments was significantly conserved across species. CONCLUSIONS: We report the first multi-species and multi-assay genome annotation results obtained by a FAANG project. Beyond the generation of reference annotations and the confirmation of previous findings on model animals, the integrative analysis of data from multiple assays and species sheds a new light on the multi-scale selective pressure shaping genome organization from birds to mammals. Overall, these results emphasize the value of FAANG for research on domesticated animals and reinforces the importance of future meta-analyses of the reference datasets being generated by this community on different species.

A missense mutation in TYRP1 causes the chocolate plumage color in chicken and alters melanosome structure.

The chocolate plumage color in chickens is due to a sex-linked recessive mutation, choc, which dilutes eumelanin pigmentation. Because TYRP1 is sex-linked in chickens, and TYRP1 mutations determine brown coat color in mammals, TYRP1 appeared as the obvious candidate gene for the choc mutation. By combining gene mapping with gene capture, a complete association was identified between the chocolate phenotype and a missense mutation leading to a His214Asn change in the ZnA zinc-binding domain of the protein. A diagnostic test confirmed complete association by screening 428 non-chocolate chickens of various origins. This is the first TYRP1 mutation described in the chicken. Electron microscopy analysis showed that melanosomes were more numerous in feather follicles of chocolate chickens but exhibited an abnormal structure characterized by a granular content and an irregular shape. A similar altered morphology was observed on melanosomes of another TYRP1 mutant in birds, the roux mutation of the quail.

Inferring the evolution of the major histocompatibility complex of wild pigs and peccaries using hybridisation DNA capture-based sequencing.

The major histocompatibility complex (MHC) is a key genomic model region for understanding the evolution of gene families and the co-evolution between host and pathogen. To date, MHC studies have mostly focused on species from major vertebrate lineages. The evolution of MHC classical (Ia) and non-classical (Ib) genes in pigs has attracted interest because of their antigen presentation roles as part of the adaptive immune system. The pig family Suidae comprises over 18 extant species (mostly wild), but only the domestic pig has been extensively sequenced and annotated. To address this, we used a DNA-capture approach, with probes designed from the domestic pig genome, to generate MHC data for 11 wild species of pigs and their closest living family, Tayassuidae. The approach showed good efficiency for wild pigs (~80% reads mapped, ~87× coverage), compared to tayassuids (~12% reads mapped, ~4× coverage). We retrieved 145 MHC loci across both families. Phylogenetic analyses show that the class Ia and Ib genes underwent multiple duplications and diversifications before suids and tayassuids diverged from their common ancestor. The histocompatibility genes mostly form orthologous groups and there is genetic differentiation for most of these genes between Eurasian and sub-Saharan African wild pigs. Tests of selection showed that the peptide-binding region of class Ib genes was under positive selection. These findings contribute to better understanding of the evolutionary history of the MHC, specifically, the class I genes, and provide useful data for investigating the immune response of wild populations against pathogens.

Sunflower oil supplementation affects the expression of miR-20a-5p and miR-142-5p in the lactating bovine mammary gland.

Oil supplementation in dairy cattle diets is used to modulate milk fat composition, as well as the expression of mammary lipogenic genes, whose regulation remains unclear. MiRNAs are small non-coding RNA considered as crucial regulators of gene expression, offering clues to explain the mechanism underlying gene nutriregulation. The present study was designed to identify miRNAs whose expression in the cow mammary gland is modulated by sunflower oil supplementation. MiRNomes were obtained using RNAseq technology from the mammary gland of lactating cows receiving a low forage diet, supplemented or not with 4% sunflower oil. Among the 272 miRNAs characterized, eight were selected for RT-qPCR validations, showing the significant down-regulation of miR-142-5p and miR-20a-5p by sunflower supplementation. These two miRNAs are predicted to target genes whose expression was reported as differentially expressed by sunflower supplementation. Among their putative targets, ELOVL6 gene involved in lipid metabolism has been studied. However, a first analysis did not show its significant down-regulation, in response to the over-expression of miR-142-5p, of miR-20a-5p, or both, in a bovine mammary epithelial cell line. However, a clearer understanding of the miRNA expression by lipid supplementation would help to decipher the regulation of lactating cow mammary gland in response to nutrition.

Draft Anaplasma phagocytophilum Genome Sequences from Five Cows, Two Horses, and One Roe Deer Collected in Europe.

Anaplasma phagocytophilum is a zoonotic tick-borne intracellular bacterium responsible for granulocytic anaplasmosis. As it is difficult to isolate and cultivate, only 20 A. phagocytophilum genomes have been sequenced to date. Here, we present eight A. phagocytophilum genome sequences obtained using alternative approaches based on sequence capture technology.

Food Deprivation Affects the miRNome in the Lactating Goat Mammary Gland.

BACKGROUND: Nutrition affects milk composition thus influencing its nutritional properties. Nutrition also modifies the expression of mammary genes, whose regulation is not fully understood. MicroRNAs (miRNA) are small non coding RNA which are important post-transcriptional regulators of gene expression by targeting messenger RNAs. Our goal was to characterize miRNA whose expression is regulated by nutrition in the lactating goat mammary gland, which may provide clues to deciphering regulations of the biosynthesis and secretion of milk components. METHODOLOGY/PRINCIPAL FINDINGS: Using high-throughput sequencing technology, miRNomes of the lactating mammary gland were established from lactating goats fed ad libitum or deprived of food for 48 h affecting milk production and composition. High throughput miRNA sequencing revealed 30 miRNA with an expression potentially modulated by food deprivation; 16 were down-regulated and 14 were up-regulated. Diana-microT predictive tools suggested a potential role for several nutriregulated miRNA in lipid metabolism. Among the putative targets, 19 were previously identified as differently expressed genes (DEG). The functions of these 19 DEG revealed, notably, their involvement in tissue remodelling. CONCLUSION/SIGNIFICANCE: In conclusion, this study offers the first evidence of nutriregulated miRNA in the ruminant mammary gland. Characterization of these 30 miRNA could contribute to a clearer understanding of gene regulation in the mammary gland in response to nutrition.

Annotation of the goat genome using next generation sequencing of microRNA expressed by the lactating mammary gland: comparison of three approaches.

BACKGROUND: MicroRNAs (miRNA) are small endogenous non-coding RNA involved in the post-transcriptional regulation of specific mRNA targets. The first whole goat genome sequence became available in 2013, with few annotations. Our goal was to establish a list of the miRNA expressed in the mammary gland of lactating goats, thus enabling implementation of the goat miRNA repertoire and considerably enriching annotation of the goat genome. RESULTS: Here, we performed high throughput RNA sequencing on 10 lactating goat mammary glands. The bioinformatic detection of miRNA was carried out using miRDeep2 software. Three different methods were used to predict, quantify and annotate the sequenced reads. The first was a de novo approach based on the prediction of miRNA from the goat genome only. The second approach used bovine miRNA as an external reference whereas the last one used recently available goat miRNA. The three methods enabled the prediction and annotation of hundreds of miRNA, more than 95% were commonly identified. Using bovine miRNA, 1,178 distinct miRNA were detected, together with the annotation of 88 miRNA for which corresponding precursors could not be retrieved in the goat genome, and which were not detected using the de novo approach or with the use of goat miRNA. Each chromosomal coordinate of the precursors determined here were generated and depicted on a reference localisation map. Forty six goat miRNA clusters were also reported. The study revealed 263 precursors located in goat protein-coding genes, amongst which the location of 43 precursors was conserved between human, mouse and bovine, revealing potential new gene regulations in the goat mammary gland. Using the publicly available cattle QTL database, and cow precursors conserved in the goat and expressed in lactating mammary gland, 114 precursors were located within known QTL regions for milk production and composition. CONCLUSIONS: The results reported here represent the first major identification study on miRNA expressed in the goat mammary gland at peak lactation. The elements generated by this study will now be used as references to decipher the regulation of miRNA expression in the goat mammary gland and to clarify their involvement in the lactation process.

Comparative genome analyses reveal distinct structure in the saltwater crocodile MHC.

The major histocompatibility complex (MHC) is a dynamic genome region with an essential role in the adaptive immunity of vertebrates, especially antigen presentation. The MHC is generally divided into subregions (classes I, II and III) containing genes of similar function across species, but with different gene number and organisation. Crocodylia (crocodilians) are widely distributed and represent an evolutionary distinct group among higher vertebrates, but the genomic organisation of MHC within this lineage has been largely unexplored. Here, we studied the MHC region of the saltwater crocodile (Crocodylus porosus) and compared it with that of other taxa. We characterised genomic clusters encompassing MHC class I and class II genes in the saltwater crocodile based on sequencing of bacterial artificial chromosomes. Six gene clusters spanning ∼452 kb were identified to contain nine MHC class I genes, six MHC class II genes, three TAP genes, and a TRIM gene. These MHC class I and class II genes were in separate scaffold regions and were greater in length (2-6 times longer) than their counterparts in well-studied fowl B loci, suggesting that the compaction of avian MHC occurred after the crocodilian-avian split. Comparative analyses between the saltwater crocodile MHC and that from the alligator and gharial showed large syntenic areas (>80% identity) with similar gene order. Comparisons with other vertebrates showed that the saltwater crocodile had MHC class I genes located along with TAP, consistent with birds studied. Linkage between MHC class I and TRIM39 observed in the saltwater crocodile resembled MHC in eutherians compared, but absent in avian MHC, suggesting that the saltwater crocodile MHC appears to have gene organisation intermediate between these two lineages. These observations suggest that the structure of the saltwater crocodile MHC, and other crocodilians, can help determine the MHC that was present in the ancestors of archosaurs.

Comparative genomics of first available bovine Anaplasma phagocytophilum genome obtained with targeted sequence capture.

BACKGROUND: Anaplasma phagocytophilum is a zoonotic and obligate intracellular bacterium transmitted by ticks. In domestic ruminants, it is the causative agent of tick-borne fever, which causes significant economic losses in Europe. As A. phagocytophilum is difficult to isolate and cultivate, only nine genome sequences have been published to date, none of which originate from a bovine strain.Our goals were to; 1/ develop a sequencing methodology which efficiently circumvents the difficulties associated with A. phagocytophilum isolation and culture; 2/ describe the first genome of a bovine strain; and 3/ compare it with available genomes, in order to both explore key genomic features at the species level, and to identify candidate genes that could be specific to bovine strains. RESULTS: DNA was extracted from a bovine blood sample infected by A. phagocytophilum. Following a whole genome capture approach, A. phagocytophilum DNA was enriched 197-fold in the sample and then sequenced using Illumina technology. In total, 58.9% of obtained reads corresponded to the A. phagocytophilum genome, covering 85.3% of the HZ genome. Then by performing comparisons with nine previously-sequenced A. phagocytophilum genomes, we determined the core genome of these ten strains. Following analysis, 1281 coding DNA sequences, including 1001 complete sequences, were detected in the A. phagocytophilum bovine genome, of which four appeared to be unique to the bovine isolate. These four coding DNA sequences coded for "hypothetical proteins of unknown function" and require further analysis. We also identified nine proteins common to both European domestic ruminants tested. CONCLUSION: Using a whole genome capture approach, we have sequenced the first A. phagocytophilum genome isolated from a cow. To the best of our knowledge, this is the first time that this method has been used to selectively enrich pathogenic bacterial DNA from samples also containing host DNA. The four proteins unique to the A. phagocytophilum bovine genome could be involved in host tropism, therefore their functions need to be explored.

Characterisation and comparison of lactating mouse and bovine mammary gland miRNomes.

BACKGROUND: The mammary gland is a dynamic organ that undergoes important physiological changes during reproductive cycles. Until now, data regarding the characterisation of miRNA in the mammary gland have been scarce and mainly focused on their abnormal expression in breast cancer. Our goal was to characterise the microRNA (miRNA) involved in mechanisms regulating the mammary function, with particular focus on the lactation stage. METHODOLOGY/PRINCIPAL FINDINGS: Using high-throughput sequencing technology, the exhaustive repertoires of miRNA expressed (miRNome) in mouse and bovine mammary glands during established lactation were identified, characterized and compared. Furthermore, in order to obtain more information on miRNA loading in the RNA-induced silencing complex (RISC), the miRNome was compared with that obtained from RNA associated with the AGO2 protein (AGO2-miRNome) in mouse lactating mammary gland. This study enabled the identification of 164 and 167 miRNA in mouse and bovine, respectively. Among the 30 miRNA most highly expressed in each species, 24 were common to both species and six of them were preferentially highly expressed in lactating than non-lactating mammary gland. The potential functional roles of these 24 miRNA were deduced using DIANA-miRPath software, based on miRNA/mRNA interactions. Moreover, seven putative novel miRNA were identified. Using DAVID analysis, it was concluded that the predicted targets of two of these putative novel miRNA are involved in mammary gland morphogenesis. CONCLUSION/SIGNIFICANCE: Our study provides an overview of the characteristics of lactating mouse and bovine mammary gland miRNA expression profiles. Moreover, species-conserved miRNA involved in this fundamental biological function were identified. These miRNomes will now be used as references for further studies during which the impact of animal breeding on the miRNA expression will be analysed.

Prion protein and Shadoo are involved in overlapping embryonic pathways and trophoblastic development.

The potential requirement of either the Prion or Shadoo protein for early mouse embryogenesis was recently suggested. However, the current data did not allow to precise the developmental process that was affected in the absence of both proteins and that led to the observed early lethal phenotype. In the present study, using various Prnp transgenic mouse lines and lentiviral vectors expressing shRNAs that target the Shadoo-encoding mRNA, we further demonstrate the specific requirement of at least one of these two PrP-related proteins at early developmental stages. Histological analysis reveals developmental defect of the ectoplacental cone and important hemorrhage surrounding the Prnp-knockout-Sprn-knockdown E7.5 embryos. By restricting the RNA interference to the trophoblastic cell lineages, the observed lethal phenotype could be attributed to the sole role of these proteins in this trophectoderm-derived compartment. RNAseq analysis performed on early embryos of various Prnp and Sprn genotypes indicated that the simultaneous down-regulation of these two proteins affects cell-adhesion and inflammatory pathways as well as the expression of ectoplacental-specific genes. Overall, our data provide biological clues in favor of a crucial and complementary embryonic role of the prion protein family in Eutherians and emphasizes the need to further evaluate its implication in normal and pathological human placenta biology.

Genes under positive selection in a model plant pathogenic fungus, Botrytis.

The rapid evolution of particular genes is essential for the adaptation of pathogens to new hosts and new environments. Powerful methods have been developed for detecting targets of selection in the genome. Here we used divergence data to compare genes among four closely related fungal pathogens adapted to different hosts to elucidate the functions putatively involved in adaptive processes. For this goal, ESTs were sequenced in the specialist fungal pathogens Botrytis tulipae and Botrytis ficariarum, and compared with genome sequences of Botrytis cinerea and Sclerotinia sclerotiorum, responsible for diseases on over 200 plant species. A maximum likelihood-based analysis of 642 predicted orthologs detected 21 genes showing footprints of positive selection. These results were validated by resequencing nine of these genes in additional Botrytis species, showing they have also been rapidly evolving in other related species. Twenty of the 21 genes had not previously been identified as pathogenicity factors in B. cinerea, but some had functions related to plant-fungus interactions. The putative functions were involved in respiratory and energy metabolism, protein and RNA metabolism, signal transduction or virulence, similarly to what was detected in previous studies using the same approach in other pathogens. Mutants of B. cinerea were generated for four of these genes as a first attempt to elucidate their functions.

Analysis of porcine MHC using microarrays.

The major histocompatibility complex (MHC) in Mammals is one of the most gene dense regions of the genome and contains the polymorphic histocompatibility gene families known to be involved in pathogen response and control of auto-immunity. The MHC is a complex genetic system that provides an interesting model system to study genome expression regulation and genetic diversity at the megabase scale. The pig MHC or SLA (Swine Leucocyte Antigen) complex spans 2.4 megabases and 151 loci have been annotated. We will review key results from previous RNA expression studies using microarrays containing probes specific to annotated loci within SLA and in addition present novel data obtained using high-density tiling arrays encompassing the whole SLA complex. We have focused on transcriptome modifications of porcine peripheral blood mononuclear cells stimulated with a mixture of phorbol myristate acetate and ionomycin known to activate B and T cell proliferation. Our results show that numerous loci mapping to the SLA complex are affected by the treatment. A general decreased level of expression for class I and II genes and an up-regulation of genes involved in peptide processing and transport were observed. Tiling array-based experiments contributed to refined gene annotations as presented for one SLA class I gene referred to as SLA-11. In conclusion, high-density tiling arrays can serve as an excellent tool to draw comprehensive transcription maps, and improve genome annotations for the SLA complex. We are currently studying their relevance to characterize SLA genetic diversity in combination with high throughput next generation sequencing.

Transcriptomic analysis brings new insight into the biological role of the prion protein during mouse embryogenesis.

The biological function of the Prion protein remains largely unknown but recent data revealed its implication in early zebrafish and mammalian embryogenesis. To gain further insight into its biological function, comparative transcriptomic analysis between FVB/N and FVB/N Prnp knockout mice was performed at early embryonic stages. RNAseq analysis revealed the differential expression of 73 and 263 genes at E6.5 and E7.5, respectively. The related metabolic pathways identified in this analysis partially overlap with those described in PrP1 and PrP2 knockdown zebrafish embryos and prion-infected mammalian brains and emphasize a potentially important role for the PrP family genes in early developmental processes.

Identification of differentially expressed genes in chickens differing in muscle glycogen content and meat quality.

BACKGROUND: The processing ability of poultry meat is highly related to its ultimate pH, the latter being mainly determined by the amount of glycogen in the muscle at death. The genetic determinism of glycogen and related meat quality traits has been established in the chicken but the molecular mechanisms involved in variations in these traits remain to be fully described. In this study, Chicken Genome Arrays (20 K) were used to compare muscle gene expression profiles of chickens from Fat (F) and Lean (L) lines that exhibited high and low muscle glycogen content, respectively, and of individuals exhibiting extremely high (G+) or low (G-) muscle glycogen content originating from the F2 cross between the Fat and Lean lines. Real-time RT-PCR was subsequently performed to validate the differential expression of genes either selected from the microarray analysis or whose function in regulating glycogen metabolism was well known. RESULTS: Among the genes found to be expressed in chicken P. major muscle, 197 and 254 transcripts appeared to be differentially expressed on microarrays for the F vs. L and the G+ vs. G- comparisons, respectively. Some involved particularly in lipid and carbohydrate metabolism were selected for further validation studies by real-time RT-PCR. We confirmed that, as in mammals, the down-regulation of CEBPB and RGS2 coincides with a decrease in peripheral adiposity in the chicken, but these genes are also suggested to affect muscle glycogen turnover through their role in the cAMP-dependent signalling pathway. Several other genes were suggested to have roles in the regulation of glycogen storage in chicken muscle. PDK4 may act as a glycogen sensor in muscle, UGDH may compete for glycogen synthesis by using UDP-glucose for glucoronidation, and PRKAB1, PRKAG2, and PHKD may impact on glycogen turnover in muscle, through AMP-activated signalling pathways. CONCLUSIONS: This study is the first stage in the understanding of molecular mechanisms underlying variations in poultry meat quality. Large scale analyses are now required to validate the role of the genes identified and ultimately to find molecular markers that can be used for selection or to optimize rearing practices.