The objective of this study was to estimate economic returns and costs associated with 4 scenarios of genetic evaluation that combine genotypes, phenotypes, and pedigree information from a vertically integrated purebred (PB) and commercial (CM) beef cattle system. Inference was to a genetic evaluation for a production system producing Charolais terminal sires for 10,000 CM cows. The first genetic evaluation scenario, denoted PB_A, modeled a genetic evaluation in which pedigree information and phenotypes are available for PB seedstock animals. Scenario PB_H contained the same information as PB_A with the addition of 25K density (GeneSeek Genomic Profiler LD) single nucleotide polymorphism (SNP) genotypes from PB animals. Scenario PBCM_A contained pedigree records and phenotypes from PB and CM cattle. Scenario PBCM_H contained phenotypes, pedigree, and genotypes from the PB and CM animals. Estimates of prediction error variance, (co)variance, and selection index parameters were used to estimate accuracy of selection candidates (rTI) and genetic gain resulting from selection on an economic index in US dollars (ΔG). Annual costs and incomes were used to determine the 30-yr cumulative net present value (CNPV) per CM calf resulting from selection in these genetic evaluation scenarios. Adding genotypes and CM production phenotypes to genetic evaluation increased the rTI of selection candidates and ΔG across all 4 scenarios. Scenario PBCM_H produced the highest annual ΔG in the PB herd at US$11.91 per head. Including CM phenotypes and parentage testing in the genetic evaluation increased the time to breakeven from 12 yr in PB_A to 19 years in PBCM_A after accounting for the cost of that information. Adding CM phenotypes and genotypes increased the breakeven time from 12 yr in PB_H to 18 yr in PBCM_H. Scenario PB_H produced the highest 30-yr CNPV per slaughtered CM calf at US$371.16. These results using field data indicate that economically relevant rTI and ΔG can be realized by adding 25K SNP genotypes and CM phenotypes to genetic evaluation, but the additional cost of that data significantly delays the economic return to the enterprise.

Interferon β (IFN-β) signaling activates the transcription factor complex ISGF3 to induce gene expression programs critical for antiviral defense and host immune responses. It has also been observed that IFN-β activates a second transcription factor complex, γ-activated factor (GAF), but the significance of this coordinated activation is unclear. We report that in murine lung epithelial cells (MLE12) high doses of IFN-β indeed activate both ISGF3 and GAF, which bind to distinct genomic locations defined by their respective DNA sequence motifs. In contrast, low doses of IFN-β preferentially activate ISGF3 but not GAF. Surprisingly, in MLE12 cells GAF binding does not induce nearby gene expression even when strongly bound to the promoter. Yet expression of interferon stimulated genes is enhanced when GAF and ISGF3 are both active compared to ISGF3 alone. We propose that GAF may function as a dose-sensitive amplifier of ISG expression to enhance antiviral immunity and establish pro-inflammatory states.

Loss of heterochromatin has been implicated as a cause of pre-mature aging and age-associated decline in organ functions in mammals; however, the specific cell types and gene loci affected by this type of epigenetic change have remained unclear. To address this knowledge gap, we probed chromatin accessibility at single-cell resolution in the brains, hearts, skeletal muscles, and bone marrows from young, middle-aged, and old mice, and assessed age-associated changes at 353,126 candidate cis-regulatory elements (cCREs) across 32 major cell types. Unexpectedly, we detected increased chromatin accessibility within specific heterochromatin domains in old mouse excitatory neurons. The gain of chromatin accessibility at these genomic loci was accompanied by the cell-type-specific loss of heterochromatin and activation of LINE1 elements. Immunostaining further confirmed the loss of the heterochromatin mark H3K9me3 in the excitatory neurons but not in inhibitory neurons or glial cells. Our results reveal the cell-type-specific changes in chromatin landscapes in old mice and shed light on the scope of heterochromatin loss in mammalian aging.

Structural and transcriptional changes during early brain maturation follow fixed developmental programs defined by genetics. However, whether this is true for functional network activity remains unknown, primarily due to experimental inaccessibility of the initial stages of the living human brain. Here, we developed human cortical organoids that dynamically change cellular populations during maturation and exhibited consistent increases in electrical activity over the span of several months. The spontaneous network formation displayed periodic and regular oscillatory events that were dependent on glutamatergic and GABAergic signaling. The oscillatory activity transitioned to more spatiotemporally irregular patterns, and synchronous network events resembled features similar to those observed in preterm human electroencephalography. These results show that the development of structured network activity in a human neocortex model may follow stable genetic programming. Our approach provides opportunities for investigating and manipulating the role of network activity in the developing human cortex.

Misregulated gene expression in human hearts can result in cardiovascular diseases that are leading causes of mortality worldwide. However, the limited information on the genomic location of candidate cis-regulatory elements (cCREs) such as enhancers and promoters in distinct cardiac cell types has restricted the understanding of these diseases. Here, we defined >287,000 cCREs in the four chambers of the human heart at single-cell resolution, which revealed cCREs and candidate transcription factors associated with cardiac cell types in a region-dependent manner and during heart failure. We further found cardiovascular disease-associated genetic variants enriched within these cCREs including 38 candidate causal atrial fibrillation variants localized to cardiomyocyte cCREs. Additional functional studies revealed that two of these variants affect a cCRE controlling KCNH2/HERG expression and action potential repolarization. Overall, this atlas of human cardiac cCREs provides the foundation for illuminating cell type-specific gene regulation in human hearts during health and disease.

Respiratory failure associated with COVID-19 has placed focus on the lungs. Here, we present single-nucleus accessible chromatin profiles of 90,980 nuclei and matched single-nucleus transcriptomes of 46,500 nuclei in non-diseased lungs from donors of ~30 weeks gestation,~3 years and ~30 years. We mapped candidate cis-regulatory elements (cCREs) and linked them to putative target genes. We identified distal cCREs with age-increased activity linked to SARS-CoV-2 host entry gene TMPRSS2 in alveolar type 2 cells, which had immune regulatory signatures and harbored variants associated with respiratory traits. At the 3p21.31 COVID-19 risk locus, a candidate variant overlapped a distal cCRE linked to SLC6A20, a gene expressed in alveolar cells and with known functional association with the SARS-CoV-2 receptor ACE2. Our findings provide insight into regulatory logic underlying genes implicated in COVID-19 in individual lung cell types across age. More broadly, these datasets will facilitate interpretation of risk loci for lung diseases.

The fetal period is a critical time for brain development, characterized by neurogenesis, neural migration, and synaptogenesis 1-3 . Microglia, the tissue resident macrophages of the brain, are observed as early as the fourth week of gestation 4 and are thought to engage in a variety of processes essential for brain development and homeostasis 5-11 . Conversely, microglia phenotypes are highly regulated by the brain environment 12-14 . Mechanisms by which human brain development influences the maturation of microglia and microglia potential contribution to neurodevelopmental disorders remain poorly understood. Here, we performed transcriptomic analysis of human fetal and postnatal microglia and corresponding cortical tissue to define age-specific brain environmental factors that may drive microglia phenotypes. Comparative analysis of open chromatin profiles using bulk and single-cell methods in conjunction with a new computational approach that integrates epigenomic and single-cell RNA-seq data allowed decoding of cellular heterogeneity with inference of subtype- and development stage-specific transcriptional regulators. Interrogation of in vivo and in vitro iPSC-derived microglia models provides evidence for roles of putative instructive signals and downstream gene regulatory networks which establish human-specific fetal and postnatal microglia gene expression programs and potentially contribute to neurodevelopmental disorders.

Microglia phenotypes are highly regulated by the brain environment, but the transcriptional networks that specify the maturation of human microglia are poorly understood. Here, we characterized stage-specific transcriptomes and epigenetic landscapes of fetal and postnatal human microglia and acquired corresponding data in induced pluripotent stem cell (iPSC)-derived microglia, in cerebral organoids, and following engraftment into humanized mice. Parallel development of computational approaches that considered transcription factor (TF) co-occurrence and enhancer activity allowed prediction of shared and state-specific gene regulatory networks associated with fetal and postnatal microglia. Additionally, many features of the human fetal-to-postnatal transition were recapitulated in a time-dependent manner following the engraftment of iPSC cells into humanized mice. These data and accompanying computational approaches will facilitate further efforts to elucidate mechanisms by which human microglia acquire stage- and disease-specific phenotypes.

The evolutionarily conserved splicing regulator neuro-oncological ventral antigen 1 (NOVA1) plays a key role in neural development and function. NOVA1 also includes a protein-coding difference between the modern human genome and Neanderthal and Denisovan genomes. To investigate the functional importance of an amino acid change in humans, we reintroduced the archaic allele into human induced pluripotent cells using genome editing and then followed their neural development through cortical organoids. This modification promoted slower development and higher surface complexity in cortical organoids with the archaic version of NOVA1 Moreover, levels of synaptic markers and synaptic protein coassociations correlated with altered electrophysiological properties in organoids expressing the archaic variant. Our results suggest that the human-specific substitution in NOVA1, which is exclusive to modern humans since divergence from Neanderthals, may have had functional consequences for our species' evolution.

Angiotensin-converting enzyme 2 (ACE2) and accessory proteases (TMPRSS2 and CTSL) are needed for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cellular entry, and their expression may shed light on viral tropism and impact across the body. We assessed the cell-type-specific expression of ACE2, TMPRSS2 and CTSL across 107 single-cell RNA-sequencing studies from different tissues. ACE2, TMPRSS2 and CTSL are coexpressed in specific subsets of respiratory epithelial cells in the nasal passages, airways and alveoli, and in cells from other organs associated with coronavirus disease 2019 (COVID-19) transmission or pathology. We performed a meta-analysis of 31 lung single-cell RNA-sequencing studies with 1,320,896 cells from 377 nasal, airway and lung parenchyma samples from 228 individuals. This revealed cell-type-specific associations of age, sex and smoking with expression levels of ACE2, TMPRSS2 and CTSL. Expression of entry factors increased with age and in males, including in airway secretory cells and alveolar type 2 cells. Expression programs shared by ACE2+TMPRSS2+ cells in nasal, lung and gut tissues included genes that may mediate viral entry, key immune functions and epithelial-macrophage cross-talk, such as genes involved in the interleukin-6, interleukin-1, tumor necrosis factor and complement pathways. Cell-type-specific expression patterns may contribute to the pathogenesis of COVID-19, and our work highlights putative molecular pathways for therapeutic intervention.