The Relationship between Gene Network Structure and Expression Variation among Individuals and Species

Progress and opportunities through use of genomics in animal production

The rearing of farmed animals is a vital component of global food production systems, but its impact on the environment, human health, animal welfare, and biodiversity is being increasingly challenged. Developments in genetic and genomic technologies have had a key role in improving the productivity of farmed animals for decades. Advances in genome sequencing, annotation, and editing offer a means not only to continue that trend, but also, when combined with advanced data collection, analytics, cloud computing, appropriate infrastructure, and regulation, to take precision livestock farming (PLF) and conservation to an advanced level. Such an approach could generate substantial additional benefits in terms of reducing use of resources, health treatments, and environmental impact, while also improving animal health and welfare.

Patterns of selection across gene regulatory networks

Gene regulatory networks (GRNs) are the core engine of organismal development. If we would like to understand the origin and diversification of phenotypes, it is necessary to consider the structure of GRNs in order to reconstruct the links between genetic mutations and phenotypic change. Much of the progress in evolutionary developmental biology, however, has occurred without a nuanced consideration of the evolution of functional relationships between genes, especially in the context of their broader network interactions. Characterizing and comparing GRNs across traits and species in a more detailed way will allow us to determine how network position influences what genes drive adaptive evolution. In this perspective paper, we consider the architecture of developmental GRNs and how positive selection strength may vary across a GRN. We then propose several testable models for these patterns of selection and experimental approaches to test these models.

Rethinking positional information and digit identity: The role of late interdigit signaling

Seminal work from John Fallon's lab has illuminated how digit identity determination involves ongoing late regulation and occurs progressively during phalanx formation. Complementary genetic analyses in mice and several papers in this special issue have begun to flesh out how interdigit signaling accomplishes this, but major questions remain unaddressed, including how uncommitted progenitors from which phalanges arise are maintained, and what factors set limits on digit extension and phalanx number, particularly in mammals. This review summarizes what has been learned in the two decades since control of digit identity by late interdigit signals was first identified and what remains poorly understood, and will hopefully spark renewed interest in a process that is critical to evolutionary limb adaptations but nevertheless remains enigmatic.

Postnatal development in a marsupial model, the fat-tailed dunnart (Sminthopsis crassicaudata; Dasyuromorphia: Dasyuridae)

Marsupials exhibit unique biological features that provide fascinating insights into many aspects of mammalian development. These include their distinctive mode of reproduction, altricial stage at birth, and the associated heterochrony that is required for their crawl to the pouch and teat attachment. Marsupials are also an invaluable resource for mammalian comparative biology, forming a distinct lineage from the extant placental and egg-laying monotreme mammals. Despite their unique biology, marsupial resources are lagging behind those available for placentals. The fat-tailed dunnart (Sminthopsis crassicaudata) is a laboratory based marsupial model, with simple and robust husbandry requirements and a short reproductive cycle making it amenable to experimental manipulations. Here we present a detailed staging series for the fat-tailed dunnart, focusing on their accelerated development of the forelimbs and jaws. This study provides the first skeletal developmental series on S. crassicaudata and provides a fundamental resource for future studies exploring mammalian diversification, development and evolution.

Non-model systems in mammalian forelimb evo-devo

Mammal forelimbs are highly diverse, ranging from the elongated wing of a bat to the stout limb of the mole. The mammal forelimb has been a long-standing system for the study of early developmental patterning, proportional variation, shape change, and the reduction of elements. However, most of this work has been performed in mice, which neglects the wide variation present across mammal forelimbs. This review emphasizes the critical role of non-model systems in limb evo-devo and highlights new emerging models and their potential. We discuss the role of gene networks in limb evolution, and touch on functional analyses that lay the groundwork for further developmental studies. Mammal limb evo-devo is a rich field, and here we aim to synthesize the findings of key recent works and the questions to which they lead.

Gene Regulatory and Expression Differences between Mouse and Pig Limb Buds Provide Insights into the Evolutionary Emergence of Artiodactyl Traits

Digit loss/reductions are evolutionary adaptations in cursorial mammals such as pigs. To gain mechanistic insight into these processes, we performed a comparative molecular analysis of limb development in mouse and pig embryos, which revealed a loss of anterior-posterior polarity during distal progression of pig limb bud development. These alterations in pig limb buds are paralleled by changes in the mesenchymal response to Sonic hedgehog (SHH) signaling, which is altered upstream of the reduction and loss of Fgf8 expression in the ectoderm that overlaps the reduced and vestigial digit rudiments of the pig handplate, respectively. Furthermore, genome-wide open chromatin profiling using equivalent developmental stages of mouse and pig limb buds reveals the functional divergence of about one-third of the regulatory genome. This study uncovers widespread alterations in the regulatory landscapes of genes essential for limb development that likely contributed to the morphological diversion of artiodactyl limbs from the pentadactyl archetype of tetrapods.

X chromosome-dependent disruption of placental regulatory networks in hybrid dwarf hamsters

Embryonic development in mammals is highly sensitive to changes in gene expression within the placenta. The placenta is also highly enriched for genes showing parent-of-origin or imprinted expression, which is predicted to evolve rapidly in response to parental conflict. However, little is known about the evolution of placental gene expression, or if divergence of placental gene expression plays an important role in mammalian speciation. We used crosses between two species of dwarf hamsters (Phodopus sungorus and Phodopus campbelli) to examine the genetic and regulatory underpinnings of severe placental overgrowth in their hybrids. Using quantitative genetic mapping and mitochondrial substitution lines, we show that overgrowth of hybrid placentas was primarily caused by genetic differences on the maternally inherited P. sungorus X chromosome. Mitochondrial interactions did not contribute to abnormal hybrid placental development, and there was only weak correspondence between placental disruption and embryonic growth. Genome-wide analyses of placental transcriptomes from the parental species and first- and second-generation hybrids revealed a central group of co-expressed X-linked and autosomal genes that were highly enriched for maternally biased expression. Expression of this gene network was strongly correlated with placental size and showed widespread misexpression dependent on epistatic interactions with X-linked hybrid incompatibilities. Collectively, our results indicate that the X chromosome is likely to play a prominent role in the evolution of placental gene expression and the accumulation of hybrid developmental barriers between mammalian species.

Development of the Proximal-Anterior Skeletal Elements in the Mouse Hindlimb Is Regulated by a Transcriptional and Signaling Network Controlled by Sall4

The vertebrate limb serves as an experimental paradigm to study mechanisms that regulate development of the stereotypical skeletal elements. In this study, we simultaneously inactivated Sall4 using Hoxb6Cre and Plzf in mouse embryos, and found that their combined function regulates development of the proximal-anterior skeletal elements in hindlimbs. The Sall4; Plzf double knockout exhibits severe defects in the femur, tibia, and anterior digits, distinct defects compared to other allelic series of Sall4; Plzf We found that Sall4 regulates Plzf expression prior to hindlimb outgrowth. Further expression analysis indicated that Hox10 genes and GLI3 are severely downregulated in the Sall4; Plzf double knockout hindlimb bud. In contrast, PLZF expression is reduced but detectable in Sall4; Gli3 double knockout limb buds, and SALL4 is expressed in the Plzf; Gli3 double knockout limb buds. These results indicate that Plzf, Gli3, and Hox10 genes downstream of Sall4, regulate femur and tibia development. In the autopod, we show that Sall4 negatively regulates Hedgehog signaling, which allows for development of the most anterior digit. Collectively, our study illustrates genetic systems that regulate development of the proximal-anterior skeletal elements in hindlimbs.

Genotype Components as Predictors of Phenotype in Model Gene Regulatory Networks

Models of gene regulatory networks (GRN) have proven useful for understanding many aspects of the highly complex behavior of biological control networks. Randomly generated non-Boolean networks were used in experimental simulations to generate data on dynamic phenotypes as a function of several genotypic parameters. We found that predictive relationships between some phenotypes and quantitative genotypic parameters such as number of network genes, interaction density, and initial condition could be derived depending on the strength of the topological (positional) genotype on specific phenotypes. We quantitated the strength of the topological genotype effect (TGE) on a number of phenotypes in multi-gene networks. For phenotypes with a low influence of topological genotype, derived and empirical relationships using quantitative genotype parameters were accurate in phenotypic outcomes. We found a number of dynamic network properties, including oscillation behaviors, that were largely dependent on genotype topology, and for which no such general quantitative relationships were determinable. It remains to be determined if these results are applicable to biological gene regulatory networks.

Gene expression variability: the other dimension in transcriptome analysis

Transcriptome sequencing is a powerful technique to study molecular changes that underlie the differences in physiological conditions and disease progression. A typical question that is posed in such studies is finding genes with significant changes between sample groups. In this respect expression variability is regarded as a nuisance factor that is primarily of technical origin and complicates the data analysis. However, it is becoming apparent that the biological variation in gene expression might be an important molecular phenotype that can affect physiological parameters. In this review we explore the recent literature on technical and biological variability in gene expression, sources of expression variability, (epi-)genetic hallmarks, and evolutionary constraints in genes with robust and variable gene expression. We provide an overview of recent findings on effects of external cues, such as diet and aging, on expression variability and on other biological phenomena that can be linked to it. We discuss metrics and tools that were developed for quantification of expression variability and highlight the importance of future studies in this direction. To assist the adoption of expression variability analysis, we also provide a detailed description and computer code, which can easily be utilized by other researchers. We also provide a reanalysis of recently published data to highlight the value of the analysis method.

Gene Network Analysis of Interstitial Macrophages After Treatment with Induced Pluripotent Stem Cells Secretome (iPSC-cm) in the Bleomycin Injured Rat Lung

Idiopathic pulmonary fibrosis (IPF) is a complex disease involving various cell types. Macrophages are essential in maintenance of physiological homeostasis, wound repair and fibrosis in the lung. Macrophages play a crucial role in repair and remodeling by altering their phenotype and secretory pattern in response to injury. The secretome of induced pluripotent stem cells (iPSC-cm) attenuates injury and fibrosis in bleomycin injured rat lungs. In the current study, we evaluate the effect of iPSC-cm on gene expression and phenotype of interstitial macrophage in bleomycin injured rat lungs in vivo. iPSC-cm was intratracheally instilled 7 days after bleomycin induced lung injury and assessed 7 days later and single cell isolation was performed. Macrophages were FACS sorted and microarray analysis was performed. We characterized changes in the rat lung interstitial macrophages using transcriptional profiling. iPSC-cm reduced the total collagen content of the lung and reduced different macrophage populations. Gene set enrichment analysis revealed involvement of three essential pathways (a) immune modulation, (b) branching morphogenesis and (c) canonical Wnt signaling. This study demonstrates that iPSC-cm reduces fibrosis in bleomycin injured rat lung by partially altering the macrophages and regulating their gene expression.

The Combinatorial Creature: Cortical Phenotypes within and across Lifetimes

The neocortex is one of the most distinctive structures of the mammalian brain, yet also one of the most varied in terms of both size and organization. Multiple processes have contributed to this variability, including evolutionary mechanisms (i.e., alterations in gene sequence) that alter the size, organization, and connections of neocortex, and activity dependent mechanisms that can also modify these same features. Thus, changes to the neocortex can occur over different time-scales, including within a single generation. This combination of genetic and activity dependent mechanisms that create a given cortical phenotype allows the mammalian neocortex to rapidly and flexibly adjust to different body and environmental contexts, and in humans permits culture to impact brain construction.

Regulatory integration of Hox factor activity with T-box factors in limb development

In tetrapods, Tbx4, Tbx5 and Hox cluster genes are crucial for forelimb and hindlimb development and mutations in these genes are responsible for congenital limb defects. The molecular basis of their integrated mechanisms of action in the context of limb development remains poorly understood. We studied Tbx4 and Hoxc10 owing to their overlapping loss-of-function phenotypes and colocalized expression in mouse hindlimb buds. We report an extensive overlap between Tbx4 and Hoxc10 genome occupancy and their putative target genes. Tbx4 and Hoxc10 interact directly with each other, have the ability to bind to a previously unrecognized T-box-Hox composite DNA motif and show synergistic activity when acting on reporter genes. Pitx1, the master regulator for hindlimb specification, also shows extensive genomic colocalization with Tbx4 and Hoxc10. Genome occupancy by Tbx4 in hindlimb buds is similar to Tbx5 occupancy in forelimbs. By contrast, another Hox factor, Hoxd13, also interacts with Tbx4/Tbx5 but antagonizes Tbx4/Tbx5-dependent transcriptional activity. Collectively, the modulation of Tbx-dependent activity by Hox factors acting on common DNA targets may integrate different developmental processes for the balanced formation of proportionate limbs.

Epigenetic and Transcriptional Variability Shape Phenotypic Plasticity

Epigenetic and transcriptional variability contribute to the vast diversity of cellular and organismal phenotypes and are key in human health and disease. In this review, we describe different types, sources, and determinants of epigenetic and transcriptional variability, enabling cells and organisms to adapt and evolve to a changing environment. We highlight the latest research and hypotheses on how chromatin structure and the epigenome influence gene expression variability. Further, we provide an overview of challenges in the analysis of biological variability. An improved understanding of the molecular mechanisms underlying epigenetic and transcriptional variability, at both the intra- and inter-individual level, provides great opportunity for disease prevention, better therapeutic approaches, and personalized medicine.

Timing the developmental origins of mammalian limb diversity

Mammals have highly diverse limbs that have contributed to their occupation of almost every niche. Researchers have long been investigating the development of these diverse limbs, with the goals of identifying developmental processes and potential biases that shape mammalian limb diversity. To date, researchers have used techniques ranging from the genomic to the anatomic to investigate the developmental processes shaping the limb morphology of mammals from five orders (Marsupialia, Chiroptera, Rodentia, Cetartiodactyla, and Perissodactyla). Results of these studies suggest that the differential expression of genes controlling diverse cellular processes underlies mammalian limb diversity. Results also suggest that the earliest development of the limb tends to be conserved among mammalian species, while later limb development tends to be more variable. This research has established the mammalian limb as a model system for evolutionary developmental biology, and set the stage for more in-depth, cross-disciplinary research into the genetic controls, tissue-level cellular behaviors, and selective pressures that have driven the developmental evolution of mammalian limbs. Ideally, these studies will be performed in a diverse suite of mammalian species within a comparative, phylogenetic framework.

Transcriptomic insights into the genetic basis of mammalian limb diversity

Background

From bat wings to whale flippers, limb diversification has been crucial to the evolutionary success of mammals. We performed the first transcriptome-wide study of limb development in multiple species to explore the hypothesis that mammalian limb diversification has proceeded through the differential expression of conserved shared genes, rather than by major changes to limb patterning. Specifically, we investigated the manner in which the expression of shared genes has evolved within and among mammalian species.

Results

We assembled and compared transcriptomes of bat, mouse, opossum, and pig fore- and hind limbs at the ridge, bud, and paddle stages of development. Results suggest that gene expression patterns exhibit larger variation among species during later than earlier stages of limb development, while within species results are more mixed. Consistent with the former, results also suggest that genes expressed at later developmental stages tend to have a younger evolutionary age than genes expressed at earlier stages. A suite of key limb-patterning genes was identified as being differentially expressed among the homologous limbs of all species. However, only a small subset of shared genes is differentially expressed in the fore- and hind limbs of all examined species. Similarly, a small subset of shared genes is differentially expressed within the fore- and hind limb of a single species and among the forelimbs of different species.

Conclusions

Taken together, results of this study do not support the existence of a phylotypic period of limb development ending at chondrogenesis, but do support the hypothesis that the hierarchical nature of development translates into increasing variation among species as development progresses.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0902-6) contains supplementary material, which is available to authorized users.

A new developmental mechanism for the separation of the mammalian middle ear ossicles from the jaw

Multiple mammalian lineages independently evolved a definitive mammalian middle ear (DMME) through breakdown of Meckel's cartilage (MC). However, the cellular and molecular drivers of this evolutionary transition remain unknown for most mammal groups. Here, we identify such drivers in the living marsupial opossum Monodelphis domestica, whose MC transformation during development anatomically mirrors the evolutionary transformation observed in fossils. Specifically, we link increases in cellular apoptosis and TGF-BR2 signalling to MC breakdown in opossums. We demonstrate that a simple change in TGF-β signalling is sufficient to inhibit MC breakdown during opossum development, indicating that changes in TGF-β signalling might be key during mammalian evolution. Furthermore, the apoptosis that we observe during opossum MC breakdown does not seemingly occur in mouse, consistent with homoplastic DMME evolution in the marsupial and placental lineages.

Limb development: a paradigm of gene regulation

The limb is a commonly used model system for developmental biology. Given the need for precise control of complex signalling pathways to achieve proper patterning, the limb is also becoming a model system for gene regulation studies. Recent developments in genomic technologies have enabled the genome-wide identification of regulatory elements that control limb development, yielding insights into the determination of limb morphology and forelimb versus hindlimb identity. The modulation of regulatory interactions - for example, through the modification of regulatory sequences or chromatin architecture - can lead to morphological evolution, acquired regeneration capacity or limb malformations in diverse species, including humans.

A new mammalian model system for thalidomide teratogenesis: Monodelphis domestica

From 1957 to 1962, thalidomide caused birth defects in >10,000 children. While the drug was pulled from the market, thalidomide is currently prescribed to treat conditions including leprosy. As a result, a new generation of babies with thalidomide defects is being born in the developing world. This represents a serious problem, as the mechanisms by which thalidomide disrupts development remain unresolved. This lack of resolution is due, in part, to the absence of an appropriate mammalian model for thalidomide teratogenesis. We test the hypothesis that opossum (Monodelphis domestica) is well suited to model human thalidomide defects. Results suggest that opossum embryos exposed to thalidomide display a range of phenotypes (e.g., heart, craniofacial, limb defects) and penetrance similar to humans. Furthermore, all opossums with thalidomide defects exhibit vascular disruptions. Results therefore support the hypotheses that opossums make a good mammalian model for thalidomide teratogenesis, and that thalidomide can severely disrupt angiogenesis in mammals.

Use of a Conditional Ubr5 Mutant Allele to Investigate the Role of an N-End Rule Ubiquitin-Protein Ligase in Hedgehog Signalling and Embryonic Limb Development

Hedgehog (Hh) signalling is a potent regulator of cell fate and function. While much is known about the events within a Hh-stimulated cell, far less is known about the regulation of Hh-ligand production. Drosophila Hyperplastic Discs (Hyd), a ubiquitin-protein ligase, represents one of the few non-transcription factors that independently regulates both hh mRNA expression and pathway activity. Using a murine embryonic stem cell system, we revealed that shRNAi of the mammalian homologue of hyd, Ubr5, effectively prevented retinoic-acid-induced Sonic hedgehog (Shh) expression. We next investigated the UBR5:Hh signalling relationship in vivo by generating and validating a mouse bearing a conditional Ubr5 loss-of-function allele. Conditionally deleting Ubr5 in the early embryonic limb-bud mesenchyme resulted in a transient decrease in Indian hedgehog ligand expression and decreased Hh pathway activity, around E13.5. Although Ubr5-deficient limbs and digits were, on average, shorter than control limbs, the effects were not statistically significant. Hence, while loss of UBR5 perturbed Hedgehog signalling in the developing limb, there were no obvious morphological defects. In summary, we report the first conditional Ubr5 mutant mouse and provide evidence for a role for UBR5 in influencing Hh signalling, but are uncertain to whether the effects on Hedgehog signaling were direct (cell autonomous) or indirect (non-cell-autonomous). Elaboration of the cellular/molecular mechanism(s) involved may help our understanding on diseases and developmental disorders associated with aberrant Hh signalling.