Identifying risk variants for embryo aneuploidy using ultra-low coverage whole-genome sequencing from preimplantation genetic testing.

Aneuploidy frequently arises during human meiosis and is the primary cause of early miscarriage and in vitro fertilization (IVF) failure. Individuals undergoing IVF exhibit significant variability in aneuploidy rates, although the exact genetic causes of the variability in aneuploid egg production remain unclear. Preimplantation genetic testing for aneuploidy (PGT-A) using next-generation sequencing is a standard test for identifying and selecting IVF-derived euploid embryos. The wealth of embryo aneuploidy data and ultra-low coverage whole-genome sequencing (ulc-WGS) data from PGT-A have the potential to discover variants in parental genomes that are associated with aneuploidy risk in their embryos. Using ulc-WGS data from ∼10,000 PGT-A biopsies, we imputed genotype likelihoods of genetic variants in embryo genomes. We then used the imputed variants and embryo aneuploidy calls to perform a genome-wide association study of aneuploidy incidence. Finally, we carried out functional evaluation of the identified candidate gene in a mouse oocyte system. We identified one locus on chromosome 3 that is significantly associated with meiotic aneuploidy risk. One candidate gene, CCDC66, encompassed by this locus, is involved in chromosome segregation during meiosis. Using mouse oocytes, we showed that CCDC66 regulates meiotic progression and chromosome segregation fidelity, especially in older mice. Our work extended the research utility of PGT-A ulc-WGS data by allowing robust association testing and improved the understanding of the genetic contribution to maternal meiotic aneuploidy risk. Importantly, we introduce a generalizable method that has potential to be leveraged for similar association studies that use ulc-WGS data.

A Massive Proteogenomic Screen Identifies Thousands of Novel Peptides From the Human "Dark" Proteome.

Although the human gene annotation has been continuously improved over the past 2 decades, numerous studies demonstrated the existence of a "dark proteome", consisting of proteins that were critical for biological processes but not included in widely used gene catalogs. The Genotype-Tissue Expression project generated more than 15,000 RNA-seq datasets from multiple tissues, which modeled 30 million transcripts in the human genome. To provide a resource of high-confidence novel proteins from the dark proteome, we screened 50,000 mass spectrometry runs from over 900 projects to identify proteins translated from the Genotype-Tissue Expression transcript model with proteomic support. We also integrated 3.8 million common genetic variants from the gnomAD database to improve peptide identification. As a result, we identified 170,529 novel peptides with proteomic evidence, of which 6048 passed the strictest standard we defined and were supported by PepQuery. We provided a user-friendly website (https://ncorf.genes.fun/) for researchers to check the evidence of novel peptides from their studies. The findings will improve our understanding of coding genes and facilitate genomic data interpretation in biomedical research.

Genomic Evidence for the Complex Evolutionary History of Macaques (Genus Macaca).

The genus Macaca is widely distributed, occupies a variety of habitats, shows diverse phenotypic characteristics, and is one of the best-studied genera of nonhuman primates. Here, we reported five re-sequencing Macaca genomes, including one M. cyclopis, one M. fuscata, one M. thibetana, one M. silenus, and one M. sylvanus. Together with published genomes of other macaque species, we combined 20 genome sequences of 10 macaque species to investigate the gene introgression and genetic differences among the species. The network analysis of the SNV-fragment trees indicates a reticular phylogeny of macaque species. Combining the results from various analytical methods, we identified extensive ancient introgression events among macaque species. The multiple introgression signals between different species groups were also observed, such as between fascicularis group species and silenus group species. However, gene flow signals between fascicularis and sinica group were not as strong as those between fascicularis group and silenus group. On the other hand, the unidirect gene flow in M. arctoides probably occurred between the progenitor of M. arctoides and the common ancestor of fascicularis group. Our study also shows that the genetic backgrounds and genetic diversity of different macaques vary dramatically among species, even among populations of the same species. In conclusion, using whole genome sequences and multiple methods, we have studied the evolutionary history of the genus Macaca and provided evidence for extensive introgression among the species.

Common genetic risk factors in ASD and ADHD co-occurring families.

Autism spectrum disorder (ASD) and attention-deficit/hyperactivity disorder (ADHD) are two major neurodevelopmental disorders that frequently co-occur. However, the genetic mechanism of the co-occurrence remains unclear. The New Jersey Language and Autism Genetics Study (NJLAGS) collected more than 100 families with at least one member affected by ASD. NJLAGS families show a high prevalence of ADHD and provide a good opportunity to study shared genetic risk factors for ASD and ADHD. The linkage study of the NJLAGS families revealed regions on chromosomes 12 and 17 that are significantly associated with ADHD. Using whole-genome sequencing data on 272 samples from 73 NJLAGS families, we identified potential risk genes for ASD and ADHD. Within the linkage regions, we identified 36 genes that are associated with ADHD using a pedigree-based gene prioritization approach. KDM6B (Lysine Demethylase 6B) is the highest-ranking gene, which is a known risk gene for neurodevelopmental disorders, including ASD and ADHD. At the whole-genome level, we identified 207 candidate genes from the analysis of both small variants and structure variants, including both known and novel genes. Using enrichment and protein-protein interaction network analyses, we identified gene ontology terms and pathways enriched for ASD and ADHD candidate genes, such as cilia function and cation channel activity. Candidate genes and pathways identified in our study improve the understanding of the genetic etiology of ASD and ADHD and will lead to new diagnostic or therapeutic interventions for ASD and ADHD in the future.

Mathematical modeling of human oocyte aneuploidy.

Aneuploidy is the leading contributor to pregnancy loss, congenital anomalies, and in vitro fertilization (IVF) failure in humans. Although most aneuploid conceptions are thought to originate from meiotic division errors in the female germline, quantitative studies that link the observed phenotypes to underlying error mechanisms are lacking. In this study, we developed a mathematical modeling framework to quantify the contribution of different mechanisms of erroneous chromosome segregation to the production of aneuploid eggs. Our model considers the probabilities of all possible chromosome gain/loss outcomes that arise from meiotic errors, such as nondisjunction (NDJ) in meiosis I and meiosis II, and premature separation of sister chromatids (PSSC) and reverse segregation (RS) in meiosis I. To understand the contributions of different meiotic errors, we fit our model to aneuploidy data from 11,157 blastocyst-stage embryos. Our best-fitting model captures several known features of female meiosis, for instance, the maternal age effect on PSSC. More importantly, our model reveals previously undescribed patterns, including an increased frequency of meiosis II errors among eggs affected by errors in meiosis I. This observation suggests that the occurrence of NDJ in meiosis II is associated with the ploidy status of an egg. We further demonstrate that the model can be used to identify IVF patients who produce an extreme number of aneuploid embryos. The dynamic nature of our mathematical model makes it a powerful tool both for understanding the relative contributions of mechanisms of chromosome missegregation in human female meiosis and for predicting the outcomes of assisted reproduction.

Structural Variations Contribute to the Genetic Etiology of Autism Spectrum Disorder and Language Impairments.

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by restrictive interests and/or repetitive behaviors and deficits in social interaction and communication. ASD is a multifactorial disease with a complex polygenic genetic architecture. Its genetic contributing factors are not yet fully understood, especially large structural variations (SVs). In this study, we aimed to assess the contribution of SVs, including copy number variants (CNVs), insertions, deletions, duplications, and mobile element insertions, to ASD and related language impairments in the New Jersey Language and Autism Genetics Study (NJLAGS) cohort. Within the cohort, ~77% of the families contain SVs that followed expected segregation or de novo patterns and passed our filtering criteria. These SVs affected 344 brain-expressed genes and can potentially contribute to the genetic etiology of the disorders. Gene Ontology and protein-protein interaction network analysis suggested several clusters of genes in different functional categories, such as neuronal development and histone modification machinery. Genes and biological processes identified in this study contribute to the understanding of ASD and related neurodevelopment disorders.

Predicting embryonic aneuploidy rate in IVF patients using whole-exome sequencing.

Infertility is a major reproductive health issue that affects about 12% of women of reproductive age in the United States. Aneuploidy in eggs accounts for a significant proportion of early miscarriage and in vitro fertilization failure. Recent studies have shown that genetic variants in several genes affect chromosome segregation fidelity and predispose women to a higher incidence of egg aneuploidy. However, the exact genetic causes of aneuploid egg production remain unclear, making it difficult to diagnose infertility based on individual genetic variants in mother's genome. In this study, we evaluated machine learning-based classifiers for predicting the embryonic aneuploidy risk in female IVF patients using whole-exome sequencing data. Using two exome datasets, we obtained an area under the receiver operating curve of 0.77 and 0.68, respectively. High precision could be traded off for high specificity in classifying patients by selecting different prediction score cutoffs. For example, a strict prediction score cutoff of 0.7 identified 29% of patients as high-risk with 94% precision. In addition, we identified MCM5, FGGY, and DDX60L as potential aneuploidy risk genes that contribute the most to the predictive power of the model. These candidate genes and their molecular interaction partners are enriched for meiotic-related gene ontology categories and pathways, such as microtubule organizing center and DNA recombination. In summary, we demonstrate that sequencing data can be mined to predict patients' aneuploidy risk thus improving clinical diagnosis. The candidate genes and pathways we identified are promising targets for future aneuploidy studies.

Pedigree-based estimation of human mobile element retrotransposition rates.

Germline mutation rates in humans have been estimated for a variety of mutation types, including single-nucleotide and large structural variants. Here, we directly measure the germline retrotransposition rate for the three active retrotransposon elements: L1, Alu, and SVA. We used three tools for calling mobile element insertions (MEIs) (MELT, RUFUS, and TranSurVeyor) on blood-derived whole-genome sequence (WGS) data from 599 CEPH individuals, comprising 33 three-generation pedigrees. We identified 26 de novo MEIs in 437 births. The retrotransposition rate estimates for Alu elements, one in 40 births, is roughly half the rate estimated using phylogenetic analyses, a difference in magnitude similar to that observed for single-nucleotide variants. The L1 retrotransposition rate is one in 63 births and is within range of previous estimates (1:20-1:200 births). The SVA retrotransposition rate, one in 63 births, is much higher than the previous estimate of one in 900 births. Our large, three-generation pedigrees allowed us to assess parent-of-origin effects and the timing of insertion events in either gametogenesis or early embryonic development. We find a statistically significant paternal bias in Alu retrotransposition. Our study represents the first in-depth analysis of the rate and dynamics of human retrotransposition from WGS data in three-generation human pedigrees.

The lineage-specific transcription factor CDX2 navigates dynamic chromatin to control distinct stages of intestine development.

Lineage-restricted transcription factors, such as the intestine-specifying factor CDX2, often have dual requirements across developmental time. Embryonic loss of CDX2 triggers homeotic transformation of intestinal fate, whereas adult-onset loss compromises crucial physiological functions but preserves intestinal identity. It is unclear how such diverse requirements are executed across the developmental continuum. Using primary and engineered human tissues, mouse genetics, and a multi-omics approach, we demonstrate that divergent CDX2 loss-of-function phenotypes in embryonic versus adult intestines correspond to divergent CDX2 chromatin-binding profiles in embryonic versus adult stages. CDX2 binds and activates distinct target genes in developing versus adult mouse and human intestinal cells. We find that temporal shifts in chromatin accessibility correspond to these context-specific CDX2 activities. Thus, CDX2 is not sufficient to activate a mature intestinal program; rather, CDX2 responds to its environment, targeting stage-specific genes to contribute to either intestinal patterning or mature intestinal function. This study provides insights into the mechanisms through which lineage-specific regulatory factors achieve divergent functions over developmental time.

PrecisionProDB: improving the proteomics performance for precision medicine.

SUMMARY: As the next-generation sequencing technology becomes broadly applied, genomics and transcriptomics are becoming more commonly used in both research and clinical settings. However, proteomics is still an obstacle to be conquered. For most peptide search programs in proteomics, a standard reference protein database is used. Because of the thousands of coding DNA variants in each individual, a standard reference database does not provide perfect match for many proteins/peptides of an individual. A personalized reference database can improve the detection power and accuracy for individual proteomics data. To connect genomics and proteomics, we designed a Python package PrecisionProDB that is specialized for generating a personized protein database for proteomics applications. PrecisionProDB supports multiple popular file formats and reference databases, and can generate a personized database in minutes. To demonstrate the application of PrecisionProDB, we generated human population-specific reference protein databases with PrecisionProDB, which improves the number of identified peptides by 0.34% on average. In addition, by incorporating cell line-specific variants into the protein database, we demonstrated a 0.71% improvement for peptide identification in the Jurkat cell line. With PrecisionProDB and these datasets, researchers and clinicians can improve their peptide search performance by adopting the more representative protein database or adding population and individual-specific proteins to the search database with minimum increase of efforts. AVAILABILITYAND IMPLEMENTATION: PrecisionProDB and pre-calculated protein databases are freely available at https://github.com/ATPs/PrecisionProDB and https://github.com/ATPs/PrecisionProDB_references. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Whole-exome sequencing identifies genes associated with Tourette's disorder in multiplex families.

Tourette's Disorder (TD) is a neurodevelopmental disorder (NDD) that affects about 0.7% of the population and is one of the most heritable NDDs. Nevertheless, because of its polygenic nature and genetic heterogeneity, the genetic etiology of TD is not well understood. In this study, we combined the segregation information in 13 TD multiplex families with high-throughput sequencing and genotyping to identify genes associated with TD. Using whole-exome sequencing and genotyping array data, we identified both small and large genetic variants within the individuals. We then combined multiple types of evidence to prioritize candidate genes for TD, including variant segregation pattern, variant function prediction, candidate gene expression, protein-protein interaction network, candidate genes from previous studies, etc. From the 13 families, 71 strong candidate genes were identified, including both known genes for NDDs and novel genes, such as HtrA Serine Peptidase 3 (HTRA3), Cadherin-Related Family Member 1 (CDHR1), and Zinc Finger DHHC-Type Palmitoyltransferase 17 (ZDHHC17). The candidate genes are enriched in several Gene Ontology categories, such as dynein complex and synaptic membrane. Candidate genes and pathways identified in this study provide biological insight into TD etiology and potential targets for future studies.

Analysis of DNA variants in miRNAs and miRNA 3'UTR binding sites in female infertility patients.

Early human embryogenesis relies on maternal gene products accumulated during oocyte growth and maturation, until around day-3 post-fertilization when human zygotic genome activation occurs. The maternal-to-zygotic transition (MZT) is a tightly coordinated process of selective maternal transcript clearance and new zygotic transcript production. If MZT is disrupted, it will lead to developmental arrest and pregnancy loss. It is well established that microRNA (miRNA) mutations disrupt regulation of their target transcripts. We hypothesize that some cases of embryonic arrest and pregnancy loss could be explained by the mutations in the maternal genome that affect miRNA-target transcript pairs. To this end, we examined mutations within miRNAs or miRNA binding sites in the 3' untranslated regions (3'UTR) of target transcripts. Using whole-exome sequencing data from 178 women undergoing in vitro fertilization (IVF) procedures, we identified 1197 variants in miRNA genes, including 93 single nucleotide variants (SNVs) and 19 small insertions/deletions (INDELs) within the seed region of 100 miRNAs. Eight miRNA seed-region variants were significantly enriched among our patients when compared to a normal population. Within predicted 3'UTR miRNA binding sites, we identified 7393 SNVs and 1488 INDELs. Between our patients and a normal population, 52 SNVs and 30 INDELs showed significant association in the single-variant testing, whereas 51 genes showed significant association in the gene-burden analysis for genes that are expressed in preimplantation embryos. Interestingly, we found that many genes with disrupted 3'UTR miRNA binding sites follow gene expression patterns resembling MZT. In addition, some of these variants showed dramatic allele frequency difference between the patient and the normal group, offering potential utility as biomarkers for screening patients prior to IVF procedures.

Meiosis interrupted: the genetics of female infertility via meiotic failure.

Idiopathic or 'unexplained' infertility represents as many as 30% of infertility cases worldwide. Conception, implantation, and term delivery of developmentally healthy infants require chromosomally normal (euploid) eggs and sperm. The crux of euploid egg production is error-free meiosis. Pathologic genetic variants dysregulate meiotic processes that occur during prophase I, meiotic resumption, chromosome segregation, and in cell cycle regulation. This dysregulation can result in chromosomally abnormal (aneuploid) eggs. In turn, egg aneuploidy leads to a broad range of clinical infertility phenotypes, including primary ovarian insufficiency and early menopause, egg fertilization failure and embryonic developmental arrest, or recurrent pregnancy loss. Therefore, maternal genetic variants are emerging as infertility biomarkers, which could allow informed reproductive decision-making. Here, we select and deeply examine human genetic variants that likely cause dysregulation of critical meiotic processes in 14 female infertility-associated genes: SYCP3, SYCE1, TRIP13, PSMC3IP, DMC1, MCM8, MCM9, STAG3, PATL2, TUBB8, CEP120, AURKB, AURKC, andWEE2. We discuss the function of each gene in meiosis, explore genotype-phenotype relationships, and delineate the frequencies of infertility-associated variants.

Mutations in ASH1L confer susceptibility to Tourette syndrome.

Tourette syndrome (TS) is a childhood-onset neuropsychiatric disorder characterized by repetitive motor movements and vocal tics. The clinical manifestations of TS are complex and often overlap with other neuropsychiatric disorders. TS is highly heritable; however, the underlying genetic basis and molecular and neuronal mechanisms of TS remain largely unknown. We performed whole-exome sequencing of a hundred trios (probands and their parents) with detailed records of their clinical presentations and identified a risk gene, ASH1L, that was both de novo mutated and associated with TS based on a transmission disequilibrium test. As a replication, we performed follow-up targeted sequencing of ASH1L in additional 524 unrelated TS samples and replicated the association (P value = 0.001). The point mutations in ASH1L cause defects in its enzymatic activity. Therefore, we established a transgenic mouse line and performed an array of anatomical, behavioral, and functional assays to investigate ASH1L function. The Ash1l+/- mice manifested tic-like behaviors and compulsive behaviors that could be rescued by the tic-relieving drug haloperidol. We also found that Ash1l disruption leads to hyper-activation and elevated dopamine-releasing events in the dorsal striatum, all of which could explain the neural mechanisms for the behavioral abnormalities in mice. Taken together, our results provide compelling evidence that ASH1L is a TS risk gene.

Origins and mechanisms leading to aneuploidy in human eggs.

The gain or loss of a chromosome-or aneuploidy-acts as one of the major triggers for infertility and pregnancy loss in humans. These chromosomal abnormalities affect more than 40% of eggs in women at both ends of the age spectrum, that is, young girls as well as women of advancing maternal age. Recent studies in human oocytes and embryos using genomics, cytogenetics, and in silico modeling all provide new insight into the rates and potential genetic and cellular factors associated with aneuploidy at varying stages of development. Here, we review recent studies that are shedding light on potential molecular mechanisms of chromosome missegregation in oocytes and embryos across the entire female reproductive life span.

SIRT7 mediates L1 elements transcriptional repression and their association with the nuclear lamina.

Long interspersed elements-1 (LINE-1, L1) are retrotransposons that hold the capacity of self-propagation in the genome with potential mutagenic outcomes. How somatic cells restrict L1 activity and how this process becomes dysfunctional during aging and in cancer cells is poorly understood. L1s are enriched at lamin-associated domains, heterochromatic regions of the nuclear periphery. Whether this association is necessary for their repression has been elusive. Here we show that the sirtuin family member SIRT7 participates in the epigenetic transcriptional repression of L1 genome-wide in both mouse and human cells. SIRT7 depletion leads to increased L1 expression and retrotransposition. Mechanistically, we identify a novel interplay between SIRT7 and Lamin A/C in L1 repression. Our results demonstrate that SIRT7-mediated H3K18 deacetylation regulates L1 expression and promotes L1 association with elements of the nuclear lamina. The failure of such activity might contribute to the observed genome instability and compromised viability in SIRT7 knockout mice. Overall, our results reveal a novel function of SIRT7 on chromatin organization by mediating the anchoring of L1 to the nuclear envelope, and a new functional link of the nuclear lamina with transcriptional repression.

Characterization of the primate TRIM gene family reveals the recent evolution in primates.

The tripartite motif (TRIM) gene family encodes diverse distinct proteins that play important roles in many biological processes. However, the molecular evolution and phylogenetic relationships of TRIM genes in primates are still elusive. We performed a genomic approach to identify and characterize TRIM genes in human and other six primate genomes. In total, 537 putative functional TRIM genes were identified and TRIM members varied among primates. A neighbor joining (NJ) tree based on the protein sequences of 82 human TRIM genes indicates seven TRIM groups, which is consistent with the results based on the architectural motifs. Many TRIM gene duplication events were identified, indicating a recent expansion of TRIM family in primate lineages. Interestingly, the chimpanzee genome shows the greatest TRIM gene expansion among the primates; however, its congeneric species, bonobo, has the least number of TRIM genes and no duplication event. Moreover, we identified a ~ 200 kb deletion on chromosome 11 of bonobos that results in a loss of cluster3 TRIM genes. The loss of TRIM genes might have occurred within the last 2 mys. Analysis of positive selection recovered 9 previously reported and 21 newly identified positively selected TRIM genes. In particular, most positive selected sites are located in the B30.2 domains. Our results have provided new insight into the evolution of primate TRIM genes and will broaden our understanding on the functions of the TRIM family.

Exome sequencing links CEP120 mutation to maternally derived aneuploid conception risk.

STUDY QUESTION: What are the genetic factors that increase the risk of aneuploid egg production? SUMMARY ANSWER: A non-synonymous variant rs2303720 within centrosomal protein 120 (CEP120) disrupts female meiosis in vitro in mouse. WHAT IS KNOWN ALREADY: The production of aneuploid eggs, with an advanced maternal age as an established contributing factor, is the major cause of IVF failure, early miscarriage and developmental anomalies. The identity of maternal genetic variants contributing to egg aneuploidy irrespective of age is missing. STUDY DESIGN, SIZE, DURATION: Patients undergoing fertility treatment (n = 166) were deidentified and selected for whole-exome sequencing. PARTICIPANTS/MATERIALS, SETTING, METHODS: Patients self-identified their ethnic groups and their ages ranged from 22 to 49 years old. The study was performed using genomes from White, non-Hispanic patients divided into controls (97) and cases (69) according to the number of aneuploid blastocysts derived during each IVF procedure. Following a gene prioritization strategy, a mouse oocyte system was used to validate the functional significance of the discovered associated genetic variants. MAIN RESULTS AND THE ROLE OF CHANCE: Patients producing a high proportion of aneuploid blastocysts (considered aneuploid if they missed any of the 40 chromatids or had extra copies) were found to carry a higher mutational burden in genes functioning in cytoskeleton and microtubule pathways. Validation of the functional significance of a non-synonymous variant rs2303720 within Cep120 on mouse oocyte meiotic maturation revealed that ectopic expression of CEP120:p.Arg947His caused decreased spindle microtubule nucleation efficiency and increased incidence of aneuploidy. LIMITATIONS, REASONS FOR CAUTION: Functional validation was performed using the mouse oocyte system. Because spindle building pathways differ between mouse and human oocytes, the defects we observed upon ectopic expression of the Cep120 variant may alter mouse oocyte meiosis differently than human oocyte meiosis. Further studies using knock-in 'humanized' mouse models and in human oocytes will be needed to translate our findings to human system. Possible functional differences of the variant between ethnic groups also need to be investigated. WIDER IMPLICATIONS OF THE FINDINGS: Variants in centrosomal genes appear to be important contributors to the risk of maternal aneuploidy. Functional validation of these variants will eventually allow prescreening to select patients that have better chances to benefit from preimplantation genetic testing. STUDY FUNDING/COMPETING INTEREST(S): This study was funded through R01-HD091331 to K.S. and J.X. and EMD Serono Grant for Fertility Innovation to N.R.T. N.R.T. is a shareholder and an employee of Genomic Prediction. TRIAL REGISTRATION NUMBER: N/A.

Worldwide patterns of genomic variation and admixture in gray wolves.

The gray wolf (Canis lupus) is a widely distributed top predator and ancestor of the domestic dog. To address questions about wolf relationships to each other and dogs, we assembled and analyzed a data set of 34 canine genomes. The divergence between New and Old World wolves is the earliest branching event and is followed by the divergence of Old World wolves and dogs, confirming that the dog was domesticated in the Old World. However, no single wolf population is more closely related to dogs, supporting the hypothesis that dogs were derived from an extinct wolf population. All extant wolves have a surprisingly recent common ancestry and experienced a dramatic population decline beginning at least ∼30 thousand years ago (kya). We suggest this crisis was related to the colonization of Eurasia by modern human hunter-gatherers, who competed with wolves for limited prey but also domesticated them, leading to a compensatory population expansion of dogs. We found extensive admixture between dogs and wolves, with up to 25% of Eurasian wolf genomes showing signs of dog ancestry. Dogs have influenced the recent history of wolves through admixture and vice versa, potentially enhancing adaptation. Simple scenarios of dog domestication are confounded by admixture, and studies that do not take admixture into account with specific demographic models are problematic.

A YY1-dependent increase in aerobic metabolism is indispensable for intestinal organogenesis.

During late gestation, villi extend into the intestinal lumen to dramatically increase the surface area of the intestinal epithelium, preparing the gut for the neonatal diet. Incomplete development of the intestine is the most common gastrointestinal complication in neonates, but the causes are unclear. We provide evidence in mice that Yin Yang 1 (Yy1) is crucial for intestinal villus development. YY1 loss in the developing endoderm had no apparent consequences until late gestation, after which the intestine differentiated poorly and exhibited severely stunted villi. Transcriptome analysis revealed that YY1 is required for mitochondrial gene expression, and ultrastructural analysis confirmed compromised mitochondrial integrity in the mutant intestine. We found increased oxidative phosphorylation gene expression at the onset of villus elongation, suggesting that aerobic respiration might function as a regulator of villus growth. Mitochondrial inhibitors blocked villus growth in a fashion similar to Yy1 loss, thus further linking oxidative phosphorylation with late-gestation intestinal development. Interestingly, we find that necrotizing enterocolitis patients also exhibit decreased expression of oxidative phosphorylation genes. Our study highlights the still unappreciated role of metabolic regulation during organogenesis, and suggests that it might contribute to neonatal gastrointestinal disorders.