A sister of NANOG regulates genes expressed in pre-implantation human development

The identification and classification of candidate genes during the zygotic genome activation in the mammals

Zygotic genome activation (ZGA) is a critical event in early embryonic development, and thousands of genes are involved in this delicate and sophisticated biological process. To date, however, only a handful of these genes have revealed their core functions in this special process, and therefore the roles of other genes still remain unclear. In the present study, we used previously published transcriptome profiling to identify potential key genes (candidate genes) in minor ZGA and major ZGA in both human and mouse specimens, and further identified the conserved genes across species. Our results showed that 887 and 760 genes, respectively, were thought to be specific to human and mouse in major ZGA, and the other 135 genes were considered to be orthologous genes. Moreover, the conserved genes were most enriched in rRNA processing in the nucleus and cytosol, ribonucleoprotein complex biogenesis, ribonucleoprotein complex assembly and ribosome large subunit biogenesis. The findings of this first comprehensive identification and characterization of candidate genes in minor and major ZGA provide relevant insights for future studies on ZGA.

NANOGP1, a tandem duplicate of NANOG, exhibits partial functional conservation in human naïve pluripotent stem cells

Gene duplication events can drive evolution by providing genetic material for new gene functions, and they create opportunities for diverse developmental strategies to emerge between species. To study the contribution of duplicated genes to human early development, we examined the evolution and function of NANOGP1, a tandem duplicate of the transcription factor NANOG. We found that NANOGP1 and NANOG have overlapping but distinct expression profiles, with high NANOGP1 expression restricted to early epiblast cells and naïve-state pluripotent stem cells. Sequence analysis and epitope-tagging revealed that NANOGP1 is protein coding with an intact homeobox domain. The duplication that created NANOGP1 occurred earlier in primate evolution than previously thought and has been retained only in great apes, whereas Old World monkeys have disabled the gene in different ways, including homeodomain point mutations. NANOGP1 is a strong inducer of naïve pluripotency; however, unlike NANOG, it is not required to maintain the undifferentiated status of human naïve pluripotent cells. By retaining expression, sequence and partial functional conservation with its ancestral copy, NANOGP1 exemplifies how gene duplication and subfunctionalisation can contribute to transcription factor activity in human pluripotency and development.

A Random Forest-Based Genome-Wide Scan Reveals Fertility-Related Candidate Genes and Potential Inter-Chromosomal Epistatic Regions Associated With Age at First Calving in Nellore Cattle

This study aimed to perform a genome-wide association analysis (GWAS) using the Random Forest (RF) approach for scanning candidate genes for age at first calving (AFC) in Nellore cattle. Additionally, potential epistatic effects were investigated using linear mixed models with pairwise interactions between all markers with high importance scores within the tree ensemble non-linear structure. Data from Nellore cattle were used, including records of animals born between 1984 and 2015 and raised in commercial herds located in different regions of Brazil. The estimated breeding values (EBV) were computed and used as the response variable in the genomic analyses. After quality control, the remaining number of animals and SNPs considered were 3,174 and 360,130, respectively. Five independent RF analyses were carried out, considering different initialization seeds. The importance score of each SNP was averaged across the independent RF analyses to rank the markers according to their predictive relevance. A total of 117 SNPs associated with AFC were identified, which spanned 10 autosomes (2, 3, 5, 10, 11, 17, 18, 21, 24, and 25). In total, 23 non-overlapping genomic regions embedded 262 candidate genes for AFC. Enrichment analysis and previous evidence in the literature revealed that many candidate genes annotated close to the lead SNPs have key roles in fertility, including embryo pre-implantation and development, embryonic viability, male germinal cell maturation, and pheromone recognition. Furthermore, some genomic regions previously associated with fertility and growth traits in Nellore cattle were also detected in the present study, reinforcing the effectiveness of RF for pre-screening candidate regions associated with complex traits. Complementary analyses revealed that many SNPs top-ranked in the RF-based GWAS did not present a strong marginal linear effect but are potentially involved in epistatic hotspots between genomic regions in different autosomes, remarkably in the BTAs 3, 5, 11, and 21. The reported results are expected to enhance the understanding of genetic mechanisms involved in the biological regulation of AFC in this cattle breed.

Causal and Candidate Gene Variants in a Large Cohort of Women With Primary Ovarian Insufficiency

CONTEXT

A genetic etiology likely accounts for the majority of unexplained primary ovarian insufficiency (POI).

OBJECTIVE

We hypothesized that heterozygous rare variants and variants in enhanced categories are associated with POI.

DESIGN

The study was an observational study.

SETTING

Subjects were recruited at academic institutions.

PATIENTS

Subjects from Boston (n = 98), the National Institutes of Health and Washington University (n = 98), Pittsburgh (n = 20), Italy (n = 43), and France (n = 32) were diagnosed with POI (amenorrhea with an elevated follicle-stimulating hormone level). Controls were recruited for health in old age or were from the 1000 Genomes Project (total n = 233).

INTERVENTION

We performed whole exome sequencing (WES), and data were analyzed using a rare variant scoring method and a Bayes factor-based framework for identifying genes harboring pathogenic variants. We performed functional studies on identified genes that were not previously implicated in POI in a D. melanogaster model.

MAIN OUTCOME

Genes with rare pathogenic variants and gene sets with increased burden of deleterious variants were identified.

RESULTS

Candidate heterozygous variants were identified in known genes and genes with functional evidence. Gene sets with increased burden of deleterious alleles included the categories transcription and translation, DNA damage and repair, meiosis and cell division. Variants were found in novel genes from the enhanced categories. Functional evidence supported 7 new risk genes for POI (USP36, VCP, WDR33, PIWIL3, NPM2, LLGL1, and BOD1L1).

CONCLUSIONS

Candidate causative variants were identified through WES in women with POI. Aggregating clinical data and genetic risk with a categorical approach may expand the genetic architecture of heterozygous rare gene variants causing risk for POI.

Genome activation in equine in vitro-produced embryos

Embryonic genome activation is a critical event in embryo development, in which the transcriptional program of the embryo is initiated. The timing and regulation of this process are species-specific. In vitro embryo production is becoming an important clinical and research tool in the horse; however, very little is known about genome activation in this species. The objective of this work was to identify the timing of genome activation, and the transcriptional networks involved, in in vitro-produced horse embryos. RNA-Seq was performed on oocytes and embryos at eight stages of development (MII, zygote, 2-cell, 4-cell, 8-cell, 16-cell, morula, blastocyst; n = 6 per stage, 2 from each of 3 mares). Transcription of seven genes was initiated at the 2-cell stage. The first substantial increase in gene expression occurred at the 4-cell stage (minor activation), followed by massive gene upregulation and downregulation at the 8-cell stage (major activation). An increase in intronic nucleotides, indicative of transcription initiation, was also observed at the 4-cell stage. Co-expression network analyses identified groups of genes that appeared to be regulated by common mechanisms. Investigation of hub genes and binding motifs enriched in the promoters of co-expressed genes implicated several transcription factors. This work represents, to the best of our knowledge, the first genomic evaluation of embryonic genome activation in horse embryos.

Allele-specific expression analysis reveals conserved and unique features of preimplantation development in equine ICSI embryos

Embryonic genome activation and dosage compensation are major genetic events in early development. Combined analysis of single embryo RNA-seq data and parental genome sequencing was used to evaluate parental contributions to early development and investigate X-chromosome dynamics. In addition, we evaluated dimorphism in gene expression between male and female embryos. Evaluation of parent-specific gene expression revealed a minor increase in paternal expression at the 4-cell stage that increased at the 8-cell stage. We also detected eight genes with allelic expression bias that may have an important role in early development, notably NANOGNB. The main actor in X-chromosome inactivation, XIST, was significantly upregulated at the 8-cell, morula, and blastocyst stages in female embryos, with high expression at the latter. Sexual dimorphism in gene expression was identified at all stages, with strong representation of the X-chromosome in females from the 16-cell to the blastocyst stage. Female embryos showed biparental X-chromosome expression at all stages after the 4-cell stage, demonstrating the absence of imprinted X-inactivation at the embryo level. The analysis of gene dosage showed incomplete dosage compensation (0.5 < X:A < 1) in MII oocytes and embryos up to the 4-cell stage, an increase of the X:A ratio at the 16-cell and morula stages after genome activation, and a decrease of the X:A ratio at the blastocyst stage, which might be associated with the beginning of X-chromosome inactivation. This study represents the first critical analysis of parent- and sex-specific gene expression in early equine embryos produced in vitro.

Dynamic Molecular Evolution of Mammalian Homeobox Genes: Duplication, Loss, Divergence and Gene Conversion Sculpt PRD Class Repertoires

The majority of homeobox genes are highly conserved across animals, but the eutherian-specific ETCHbox genes, embryonically expressed and highly divergent duplicates of CRX, are a notable exception. Here we compare the ETCHbox genes of 34 mammalian species, uncovering dynamic patterns of gene loss and tandem duplication, including the presence of a large tandem array of LEUTX loci in the genome of the European rabbit (Oryctolagus cuniculus). Despite extensive gene gain and loss, all sampled species possess at least two ETCHbox genes, suggesting their collective role is indispensable. We find evidence for positive selection and show that TPRX1 and TPRX2 have been the subject of repeated gene conversion across the Boreoeutheria, homogenising their sequences and preventing divergence, especially in the homeobox region. Together, these results are consistent with a model where mammalian ETCHbox genes are dynamic in evolution due to functional overlap, yet have collective indispensable roles.

Evidence for novel epigenetic marks within plants

Variation in patterns of gene expression can result from modifications in the genome that occur without a change in the sequence of the DNA; such modifications include methylation of cytosine to generate 5-methylcytosine (5mC) resulting in the generation of heritable epimutation and novel epialleles. This type of non-sequence variation is called epigenetics. The enzymes responsible for generation of such DNA modifications in mammals are named DNA methyltransferases (DNMT) including DNMT1, DNMT2 and DNMT3. The later stages of oxidations to these modifications are catalyzed by Ten Eleven Translocation (TET) proteins, which contain catalytic domains belonging to the 2-oxoglutarate dependent dioxygenase family. In various mammalian cells/tissues including embryonic stem cells, cancer cells and brain tissues, it has been confirmed that these proteins are able to induce the stepwise oxidization of 5-methyl cytosine to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and finally 5-carboxylcytosine (5caC). Each stage from initial methylation until the end of the DNA demethylation process is considered as a specific epigenetic mark that may regulate gene expression. This review discusses controversial evidence for the presence of such oxidative products, particularly 5hmC, in various plant species. Whereas some reports suggest no evidence for enzymatic DNA demethylation, other reports suggest that the presence of oxidative products is followed by the active demethylation and indicate the contribution of possible TET-like proteins in the regulation of gene expression in plants. The review also summarizes the results obtained by expressing the human TET conserved catalytic domain in transgenic plants.

Expression dynamics of repetitive DNA in early human embryonic development

BACKGROUND

The last decade witnessed a number of genome-wide studies on human pre-implantation, which mostly focused on genes and provided only limited information on repeats, excluding the satellites. Considering the fact that repeats constitute a large portion of our genome with reported links to human physiology and disease, a thorough understanding of their spatiotemporal regulation during human embryogenesis will give invaluable clues on chromatin dynamics across time and space. Therefore, we performed a detailed expression analysis of all repetitive DNA elements including the satellites across stages of human pre-implantation and embryonic stem cells.

RESULTS

We uncovered stage-specific expressions of more than a thousand repeat elements whose expressions fluctuated with a mild global decrease at the blastocyst stage. Most satellites were highly expressed at the 4-cell level and expressions of ACRO1 and D20S16 specifically peaked at this point. Whereas all members of the SVA elements were highly upregulated at 8-cell and morula stages, other transposons and small RNA repeats exhibited a high level of variation among their specific subtypes. Our repeat enrichment analysis in gene promoters coupled with expression correlations highlighted potential links between repeat expressions and nearby genes, emphasising mostly 8-cell and morula specific genes together with SVA_D, LTR5_Hs and LTR70 transposons. The DNA methylation analysis further complemented the understanding on the mechanistic aspects of the repeatome's regulation per se and revealed critical stages where DNA methylation levels are negatively correlating with repeat expression.

CONCLUSIONS

Taken together, our study shows that specific expression patterns are not exclusive to genes and long non-coding RNAs but the repeatome also exhibits an intriguingly dynamic pattern at the global scale. Repeats identified in this study; particularly satellites, which were historically associated with heterochromatin, and those with potential links to nearby gene expression provide valuable insights into the understanding of key events in genomic regulation and warrant further research in epigenetics, genomics and developmental biology.

Effect of BMP15 and/or AMH during in vitro maturation of oocytes from involuntarily culled dairy cows

The high metabolic activity to which the dairy cattle are exposed to maintain milk production altered steroid metabolism that affects reproductive physiology and reduce oocyte competence. Our aims were (a) to characterize the competence of immature oocytes collected from dairy cattle based on the expression of genes in cumulus cells (CCs) and (b) to improve oocyte competence to support preimplantation embryo development by the supplementation of maturation medium with bone morphogenetic protein 15 (BMP15) and/or anti-mullerian hormone (AMH). Oocyte donors were identified at the moment of ovary collection and grouped by involuntarily culled dairy cows (Holstein breed) or beef cattle. The embryo development speed to blastocyst of the cull dairy cattle versus beef cattle (control group) was lower. Besides, <10% of oocytes (with CC biopsies) derived from dairy cattle were able to develop to the blastocyst stage. In addition, a higher level of expression and a positive correlation were observed in the expression of most of the genes evaluated (LUM, KRT18, KRT8, CLIC3, BMPR1B, and SLC38A3) in the cumulus-oocyte complexes that produced blastocysts versus those which did not develop correctly (arrested development). Further, use of BMP15 in the maturation of oocytes from dairy cattle seems to increase competence, modulating the expression of OCT4, SOX2, CDX2, GATA6, and TP1 in resulting blastocysts.

Novel and divergent genes in the evolution of placental mammals

Analysis of genome sequences within a phylogenetic context can give insight into the mode and tempo of gene and protein evolution, including inference of gene ages. This can reveal whether new genes arose on particular evolutionary lineages and were recruited for new functional roles. Here, we apply MCL clustering with all-versus-all reciprocal BLASTP to identify and phylogenetically date 'Homology Groups' among vertebrate proteins. Homology Groups include new genes and highly divergent duplicate genes. Focusing on the origin of the placental mammals within the Eutheria, we identify 357 novel Homology Groups that arose on the stem lineage of Placentalia, 87 of which are deduced to play core roles in mammalian biology as judged by extensive retention in evolution. We find the human homologues of novel eutherian genes are enriched for expression in preimplantation embryo, brain, and testes, and enriched for functions in keratinization, reproductive development, and the immune system.