Comprehensive transcriptomic analysis of long non-coding RNAs in bovine ovarian follicles and early embryos.

Long non-coding RNAs (lncRNAs) have been the subject of numerous studies over the past decade. First thought to come from aberrant transcriptional events, lncRNAs are now considered a crucial component of the genome with roles in multiple cellular functions. However, the functional annotation and characterization of bovine lncRNAs during early development remain limited. In this comprehensive analysis, we review lncRNAs expression in bovine ovarian follicles and early embryos, based on a unique database comprising 468 microarray hybridizations from a single platform designed to target 7,724 lncRNA transcripts, of which 5,272 are intergenic (lincRNA), 958 are intronic, and 1,524 are antisense (lncNAT). Compared to translated mRNA, lncRNAs have been shown to be more tissue-specific and expressed in low copy numbers. This analysis revealed that protein-coding genes and lncRNAs are both expressed more in oocytes. Differences between the oocyte and the 2-cell embryo are also more apparent in terms of lncRNAs than mRNAs. Co-expression network analysis using WGCNA generated 25 modules with differing proportions of lncRNAs. The modules exhibiting a higher proportion of lncRNAs were found to be associated with fewer annotated mRNAs and housekeeping functions. Functional annotation of co-expressed mRNAs allowed attribution of lncRNAs to a wide array of key cellular events such as meiosis, translation initiation, immune response, and mitochondrial related functions. We thus provide evidence that lncRNAs play diverse physiological roles that are tissue-specific and associated with key cellular functions alongside mRNAs in bovine ovarian follicles and early embryos. This contributes to add lncRNAs as active molecules in the complex regulatory networks driving folliculogenesis, oogenesis and early embryogenesis all of which are necessary for reproductive success.

Cytogenetic screening of a Canadian swine breeding nucleus using a newly developed karyotyping method named oligo-banding.

BACKGROUND: The frequency of chromosomal rearrangements in Canadian breeding boars has been estimated at 0.91 to 1.64%. These abnormalities are widely recognized as a potential cause of subfertility in livestock production. Since artificial insemination is practiced in almost all intensive pig production systems, the use of elite boars carrying cytogenetic defects that have an impact on fertility can lead to major economic losses. To avoid keeping subfertile boars in artificial insemination centres and spreading chromosomal defects within populations, cytogenetic screening of boars is crucial. Different techniques are used for this purpose, but several issues are frequently encountered, i.e. environmental factors can influence the quality of results, the lack of genomic information outputted by these techniques, and the need for prior cytogenetic skills. The aim of this study was to develop a new pig karyotyping method based on fluorescent banding patterns. RESULTS: The use of 207,847 specific oligonucleotides generated 96 fluorescent bands that are distributed across the 18 autosomes and the sex chromosomes. Tested alongside conventional G-banding, this oligo-banding method allowed us to identify four chromosomal translocations and a rare unbalanced chromosomal rearrangement that was not detected by conventional banding. In addition, this method allowed us to investigate chromosomal imbalance in spermatozoa. CONCLUSIONS: The use of oligo-banding was found to be appropriate for detecting chromosomal aberrations in a Canadian pig nucleus and its convenient design and use make it an interesting tool for livestock karyotyping and cytogenetic studies.

DNA methylation dynamic in male rat germ cells during gametogenesis.

In mammals, a near complete resetting of DNA methylation (DNAme) is observed during germline establishment. This wave of epigenetic reprogramming is sensitive to the environment, which could impair the establishment of an optimal state of the gamete epigenome, hence proper embryo development. Yet, we lack a comprehensive understanding of DNAme dynamics during spermatogenesis, especially in rats, the model of choice for toxicological studies. Using a combination of cell sorting and DNA methyl-seq capture, we generated a stage-specific mapping of DNAme in nine populations of differentiating germ cells from perinatal life to spermiogenesis. DNAme was found to reach its lowest level at gestational day 18, the last demethylated coding regions being associated with negative regulation of cell movement. The following de novo DNAme displayed three different kinetics with common and distinct genomic enrichments, suggesting a non-random process. DNAme variations were also detected at key steps of chromatin remodeling during spermiogenesis, revealing potential sensitivity. These methylome datasets for coding sequences during normal spermatogenesis in rat provide an essential reference for studying epigenetic-related effects of disease or environmental factors on the male germline.

Characterization of the genetic pool of the Canadienne dairy cattle breed.

BACKGROUND: Canadienne cattle are the oldest breed of dairy cattle in North America. The Canadienne breed originates from cattle that were brought to America by the mid-seventeenth century French settlers. The herd book was established in 1886 and the current breed characteristics include dark coat color, small size compared to the modern Holstein breed, and overall rusticity shaped by the harsh environmental conditions that were prevalent during the settlement of North America. The Canadienne breed is an invaluable genetic resource due to its high resilience, longevity and fertility. However, it is heavily threatened with a current herd limited to an estimated 1200 registered animals, of which less than 300 are fullblood. To date, no effort has been made to document the genetic pool of this heritage breed in order to preserve it. RESULTS: In this project, we used genomic data, which allow a precise description of the genetic makeup of a population, to provide valuable information on the genetic diversity of this heritage breed and suggest management options for its long-term viability. Using a panel that includes 640,000 single nucleotide polymorphisms (SNPs), we genotyped 190 animals grouped into six purity ranges. Unsupervised clustering analyses revealed three genetically distinct groups among those with the higher levels of purity. The observed heterozygosity was higher than expected even in the 100% purebreds. Comparison with Holstein genotypes showed significantly shorter runs of homozygosity for the Canadienne breed, which was unexpected due to the high inbreeding value calculated from pedigree data. CONCLUSIONS: Overall, our data indicate that the fullblood gene pool of the Canadienne breed is more diversified than expected and that bloodline management could promote breed sustainability. In its current state, the Canadienne is not a dead-end breed but remains highly vulnerable due to its small population size.

Epitranscriptome marks detection and localization of RNA modifying proteins in mammalian ovarian follicles.

BACKGROUND: Most of the resources that support the early development of the embryo are stored in the oocyte. Clearing of maternal resources and activation of the embryonic genome to produce its own mRNA transcripts marks the maternal-to-embryo transition. Dependence on stored mRNA can last from a few hours to several days, depending on animal species. The mechanisms regulating stabilization and recruitment of stored maternal transcripts have not yet been described in full detail but are known to involve reversible polyadenylation and modulation of 3'UTR-mediated elements. RNA epigenetic modifications, new players in this field, have an important role in RNA regulation and stabilization. RESULTS: The objectives of this study were first to determine if some of post-transcriptional methylation of stored mRNA is greater in oocytes than in somatic cells. We found that m6A, known to be the most prevalent and involved in various aspects of RNA metabolism and physiological functions, is particularly abundant in porcine oocyte mRNA compared to liver used as a somatic tissue reference. The second objective was to compare the epitranscriptome machinery, such as methyltransferases ("writers"), binding proteins ("readers") and demethylases ("erasers") catalyzing the different process, in follicles and oocytes of different mammalian species by immunofluorescence and confocal microscopy. The expression and localization patterns of these proteins differ between mice, pigs and cows ovaries and oocytes. m5C-associated proteins were generally less abundant. In contrast, m6A-associated proteins were expressed strongly during the early and late stages of folliculogenesis. Transzonal projections were found to contain more granules bearing the m5C mark in mice but both m5C and m6A methylation marks in association with mature oocytes of pigs and cows. Eraser proteins showed the greatest interspecies diversity in terms of distribution in the germinal tissues. CONCLUSIONS: So far, few studies have looked at the oocyte and ovarian epitranscriptomic profile. Our findings indicate that a hitherto unrecognized species-specific layer of transcript regulation occurs at the RNA level and might be consequential during the oocyte transcriptional silencing period.

Observation of Genetic Gain with Instrumental Insemination of Honeybee Queens.

Controlling mating in the honeybee (Apis mellifera) is part of one of the greatest challenges for the beekeeping industry's genetic selection programs due to specific characteristics of their reproduction. Several techniques for supervising honeybee mating with relative effective control have been developed over the years to allow honeybee selection. As part of this project, we compared the genetic gains for several colony performance traits, obtained using the BLUP-animal method, according to the selection pressure applied in controlled reproduction (directed fertilization versus instrumental insemination). Our results show similar genetic gains for hygienic behavior and honey production between colonies whether queens were fertilized naturally or via instrumental insemination, as well as similar or lower genetic gains for colonies with queens inseminated for spring development. In addition, we noticed greater fragility in queens following insemination. These findings show that instrumental insemination is an effective tool for reproductive control in genetic selection and for estimating breeding values more precisely. However, this technique does not result in queens of superior genetic quality for commercial purposes.

Thirty-year survey of bibliometrics used in the research literature of pain: Analysis, evolution, and pitfalls.

During the last decades, the emergence of Bibliometrics and the progress in Pain research have led to a proliferation of bibliometric studies on the medical and scientific literature of pain (B/P). This study charts the evolution of the B/P literature published during the last 30 years. Using various searching techniques, 189 B/P studies published from 1993 to August 2022 were collected for analysis-half were published since 2018. Most of the selected B/P publications use classic bibliometric analysis of Pain in toto, while some focus on specific types of Pain with Headache/Migraine, Low Back Pain, Chronic Pain, and Cancer Pain dominating. Each study is characterized by the origin (geographical, economical, institutional, …) and the medical/scientific context over a specified time span to provide a detailed landscape of the Pain research literature. Some B/P studies have been developed to pinpoint difficulties in appropriately identifying the Pain literature or to highlight some general publishing pitfalls. Having observed that most of the recent B/P studies have integrated newly emergent software visualization tools (SVTs), we found an increase of anomalies and suggest that readers exercise caution when interpreting results in the B/P literature details.

Chromosome-level assembly of the Rangifer tarandus genome and validation of cervid and bovid evolution insights.

BACKGROUND: Genome assembly into chromosomes facilitates several analyses including cytogenetics, genomics and phylogenetics. Despite rapid development in bioinformatics, however, assembly beyond scaffolds remains challenging, especially in species without closely related well-assembled and available reference genomes. So far, four draft genomes of Rangifer tarandus (caribou or reindeer, a circumpolar distributed cervid species) have been published, but none with chromosome-level assembly. This emblematic northern species is of high interest in ecological studies and conservation since most populations are declining. RESULTS: We have designed specific probes based on Oligopaint FISH technology to upgrade the latest published reindeer and caribou chromosome-level genomes. Using this oligonucleotide-based method, we found six mis-assembled scaffolds and physically mapped 68 of the largest scaffolds representing 78% of the most recent R. tarandus genome assembly. Combining physical mapping and comparative genomics, it was possible to document chromosomal evolution among Cervidae and closely related bovids. CONCLUSIONS: Our results provide validation for the current chromosome-level genome assembly as well as resources to use chromosome banding in studies of Rangifer tarandus.

Genomic health is dependent on long-term population demographic history.

Current genetic methods of population assessment in conservation biology have been challenged by genome-scale analyses due to their quantitatively novel insights. These analyses include assessments of runs-of-homozygosity (ROH), genomic evolutionary rate profiling (GERP), and mutational load. Here, we aim to elucidate the relationships between these measures using three divergent ungulates: white-tailed deer, caribou, and mountain goat. The white-tailed deer is currently expanding, while caribou are in the midst of a significant decline. Mountain goats remain stable, having suffered a large historical bottleneck. We assessed genome-wide signatures of inbreeding using the inbreeding coefficient F and %ROH (FROH ) and identified evolutionarily constrained regions with GERP. Mutational load was estimated by identifying mutations in highly constrained elements (CEs) and sorting intolerant from tolerant (SIFT) mutations. Our results showed that F and FROH are higher in mountain goats than in caribou and white-tailed deer. Given the extended bottleneck and low Ne of the mountain goat, this supports the idea that the genome-wide effects of demographic change take time to accrue. Similarly, we found that mountain goats possess more highly constrained CEs and the lowest dN/dS values, both of which are indicative of greater purifying selection; this is also reflected by fewer mutations in CEs and deleterious mutations identified by SIFT. In contrast, white-tailed deer presented the highest mutational load with both metrics, in addition to dN/dS, while caribou were intermediate. Our results demonstrate that extended bottlenecks may lead to reduced diversity and increased FROH in ungulates, but not necessarily an increase in mutational load, probably due to the purging of deleterious alleles in small populations.

Hedgehog signaling regulates Wolffian duct development through the primary cilium†.

Primary cilia play pivotal roles in embryonic patterning and organogenesis through transduction of the Hedgehog signaling pathway (Hh). Although mutations in Hh morphogens impair the development of the gonads and trigger male infertility, the contribution of Hh and primary cilia in the development of male reproductive ductules, including the epididymis, remains unknown. From a Pax2Cre; IFT88fl/fl knock-out mouse model, we found that primary cilia deletion is associated with imbalanced Hh signaling and morphometric changes in the Wolffian duct (WD), the embryonic precursor of the epididymis. Similar effects were observed following pharmacological blockade of primary cilia formation and Hh modulation on WD organotypic cultures. The expression of genes involved in extracellular matrix, mesenchymal-epithelial transition, canonical Hh and WD development was significantly altered after treatments. Altogether, we identified the primary cilia-dependent Hh signaling as a master regulator of genes involved in WD development. This provides new insights regarding the etiology of sexual differentiation and male infertility issues.

Roles of the cumulus-oocyte transzonal network and the Fragile X protein family in oocyte competence.

In brief: RNA granules travel through the cumulus cell network of transzonal projections which is associated with oocyte developmental competence, and RNA packaging involves RNA-binding proteins of the Fragile X protein family. Abstract: The determinants of oocyte developmental competence have puzzled scientists for decades. It is known that follicular conditions can nurture the production of a high-quality oocyte, but the underlying mechanisms remain unknown. Somatic cumulus cells most proximal to the oocyte are known to have cellular extensions that reach across the zona pellucida and contact with the oocyte plasma membrane. Herein, it was found that transzonal projections (TZPs) network quality is associated with developmental competence. Knowing that ribonucleoparticles are abundant within TZPs, the distribution of RNA-binding proteins was studied. The Fragile X-related proteins (FXR1P and FXR2P) and two partnering protein families, namely cytoplasmic FMRP-interacting protein and nuclear FMRP-interacting protein, exhibited distinctive patterns consistent with roles in regulating mRNA packaging, transport, and translation. The expression of green fluorescent protein (GFP)-FMRP fusion protein in cumulus cells showed active granule formation and their transport and transfer through filipodia connecting with neighboring cells. Near the projections' ends was found the cytoskeletal anchoring protein Filamin A and active protein synthesis sites. This study highlights key proteins involved in delivering mRNA to the oocyte. Thus, cumulus cells appear to indeed support the development of high-quality oocytes via the transzonal network.

Demographic history and conservation genomics of caribou (Rangifer tarandus) in Québec.

The loss of genetic diversity is a challenge many species are facing, with genomics being a potential tool to inform and prioritize decision-making. Most caribou (Rangifer tarandus) populations have experienced significant recent declines throughout Québec, Canada, and are considered of concern, threatened or endangered. Here, we calculated the ancestral and contemporary patterns of genomic diversity of five representative caribou populations and applied a comparative population genomics framework to assess the interplay between demographic events and genomic diversity. We first calculated a caribou specific mutation rate, µ, by extracting orthologous genes from related ungulates and estimating the rate of synonymous mutations. Whole genome re-sequencing was then completed on 67 caribou: from these data we calculated nucleotide diversity, θ π and estimated the coalescent or ancestral effective population size (N e), which ranged from 12,030 to 15,513. When compared to the census size, N C, the endangered Gaspésie Mountain caribou population had the highest ancestral N e:N C ratio which is consistent with recent work suggesting high ancestral N e:N C is of conservation concern. In contrast, values of contemporary N e, estimated from linkage-disequilibrium, ranged from 11 to 162, with Gaspésie having among the highest contemporary N e:N C ratio. Importantly, classic conservation genetics theory would predict this population to be of less concern based on this ratio. Interestingly, F varied only slightly between populations, and despite evidence of bottlenecks across the province, runs of homozygosity were not abundant in the genome. Tajima's D estimates mirrored the demographic models and current conservation status. Our study highlights how genomic patterns are nuanced and potentially misleading if viewed only through a contemporary lens; we argue a holistic conservation genomics view should integrate ancestral N e and Tajima's D into management decisions.

Design and validation of a 63K genome-wide SNP-genotyping platform for caribou/reindeer (Rangifer tarandus).

BACKGROUND: Development of large single nucleotide polymorphism (SNP) arrays can make genomic data promptly available for conservation problematic. Medium and high-density panels can be designed with sufficient coverage to offer a genome-wide perspective and the generated genotypes can be used to assess different genetic metrics related to population structure, relatedness, or inbreeding. SNP genotyping could also permit sexing samples with unknown associated metadata as it is often the case when using non-invasive sampling methods favored for endangered species. Genome sequencing of wild species provides the necessary information to design such SNP arrays. We report here the development of a SNP-array for endangered Rangifer tarandus using a multi-platform sequencing approach from animals found in diverse populations representing the entire circumpolar distribution of the species. RESULTS: From a very large comprehensive catalog of SNPs detected over the entire sample set (N = 894), a total of 63,336 SNPs were selected. SNP selection accounted for SNPs evenly distributed across the entire genome (~ every 50Kb) with known minor alleles across populations world-wide. In addition, a subset of SNPs was selected to represent rare and local alleles found in Eastern Canada which could be used for ecotype and population assignments - information urgently needed for conservation planning. In addition, heterozygosity from SNPs located in the X-chromosome and genotyping call-rate of SNPs located into the SRY gene of the Y-chromosome yielded an accurate and robust sexing assessment. All SNPs were validated using a high-throughput SNP-genotyping chip. CONCLUSION: This design is now integrated into the first genome-wide commercially available genotyping platform for Rangifer tarandus. This platform would pave the way to future genomic investigation of populations for this endangered species, including estimation of genetic diversity parameters, population assignments, as well as animal sexing from genetic SNP data for non-invasive samples.

Mammalian cumulus-oocyte complex communication: a dialog through long and short distance messaging.

Communications are crucial to ovarian follicle development and to ovulation, and while both folliculogenesis and oogenesis are distinct processes, they share highly interdependent signaling pathways. Signals from distant organs such as the brain must be processed and compartments within the follicle have to be synchronized. The hypothalamic-pituitary-gonadal (HPG) axis relies on long-distance signalling analogous to wireless communication by which data is disseminated in the environment and cells equipped with the appropriate receptors receive and interpret the messages. In contrast, direct cell-to-cell transfer of molecules is a very targeted, short distance messaging system. Numerous signalling pathways have been identified and proven to be essential for the production of a developmentally competent egg. The development of the cumulus-oocyte complex relies largely on short distance communications or direct transfer type via extensions of corona radiata cells through the zona pellucida. The type of information transmitted through these transzonal projections is still largely uncharacterized. This review provides an overview of current understanding of the mechanisms by which the gamete receives and transmits information within the follicle. Moreover, it highlights the fact that in addition to the well-known systemic long-distance based communications from the HPG axis, these mechanisms acting more locally should also be considered as important targets for controlling/optimizing oocyte quality.

Corn or wheat-based diet to manipulate lipid content in early embryos of Jersey cows.

Numerous field observations and previous reports have described a breed effect on embryonic survival to cryopreservation. Embryos collected from Jersey cows show a poor survival rate compared to Holstein embryos due to their rich intracellular lipid content. Less well documented in the scientific literature, but also reported from field observations, it was proposed that survival rates can be improved by removing corn from the feedstuff. This project was aimed to test whether corn replacement would modify intracellular lipid content in Jersey cow embryos. Non-lactating purebred Jersey cows were fed with an isoenergetic diet supplemented with either corn or wheat. Embryos were collected in vivo following ovarian stimulation and insemination. Results show only marginal impact on blood lipids. Overall, embryos were morphologically indistinguishable. Some lipids known to be associated with lower cryotolerance were significantly impacted by the treatment namely triacylglycerol, sphingomyelin and phosphatidylcholine. Corn supplementation significantly diminished mitochondrial activity in blastomeres. These treatment effects suggest an indirect mechanism potentially impacting embryo quality through mitochondrial dysfunction rather than direct lipid uptake.

Alcohol induced increases in sperm Histone H3 lysine 4 trimethylation correlate with increased placental CTCF occupancy and altered developmental programming.

Using a mouse model, studies by our group reveal that paternal preconception alcohol intake affects offspring fetal-placental growth, with long-lasting consequences on adult metabolism. Here, we tested the hypothesis that chronic preconception male alcohol exposure impacts histone enrichment in sperm and that these changes are associated with altered developmental programming in the placenta. Using chromatin immunoprecipitation, we find alcohol-induced increases in sperm histone H3 lysine 4 trimethylation (H3K4me3) that map to promoters and presumptive enhancer regions enriched in genes driving neurogenesis and craniofacial development. Given the colocalization of H3K4me3 with the chromatin binding factor CTCF across both sperm and embryos, we next examined CTCF localization in the placenta. We find global changes in CTCF binding within placentae derived from the male offspring of alcohol-exposed sires. Furthermore, altered CTCF localization correlates with dysregulated gene expression across multiple gene clusters; however, these transcriptional changes only occur in male offspring. Finally, we identified a correlation between genomic regions exhibiting alcohol-induced increases in sperm H3K4me3 and increased CTCF binding in male placentae. Collectively, our analysis demonstrates that the chromatin landscape of sperm is sensitive to chronic alcohol exposure and that a subset of these affected regions exhibits increased placental CTCF enrichment.