Evolution of beak morphology in the Ground Tit revealed by comparative transcriptomics

Genetic parameter estimation and molecular foundation of chicken beak shape

The bird beak is mainly functioned as feeding and attacking, and its shape has extremely important significance for survival and reproduction. In chickens, since beak shape could lead to some disadvantages including pecking and waste of feed, it is important to understand the inheritance of chicken beak shape. In the present study, we firstly established 4 indicators to describe the chicken beak shapes, including upper beak length (UL), lower beak length (LL), distance between upper and lower beak tips (DB) and upper beak curvature (BC). And then, we measured the 4 beak shape indicators as well as some production traits including body weight (BW), shank length (SL), egg weight (EW), eggshell strength (ES) of a layer breed, Rhode Island Red (RIR), in order to estimate genetic parameters of chicken beak shape. The heritabilities of UL and LL were 0.41 and 0.37, and the heritabilities of DB and BC were 0.22 and 0.21, indicating that beak shape was a highly or mediumly heritable. There were significant positive genetic and phenotypic correlations among UL, LL, and DB. And UL was positively correlated with body weight (BW18) and shank length (SL18) at 18 weeks of age in genetics, and DB was positively correlated with BC in terms of genetics and phenotype. We also found that layers of chicken cages played a role on beak shape, which could be attributed to the difference of lightness in different cage layers. By a genome-wide association study (GWAS) for the chicken UL, we identified 9 significant candidate genes associated with UL in RIR. For the variants with low minor allele frequencies (MAF <0.01) and outside of high linkage disequilibrium (LD) regions, we also conducted rare variant association studies (RVA) and GWAS to find the association between genotype and phenotype. We also analyzed transcriptomic data from multiple tissues of chicken embryos and revealed that all of the 9 genes were highly expressed in beak of chicken embryos, indicating their potential function for beak development. Our results provided the genetic foundation of chicken beak shape, which could help chicken breeding on beak related traits.

LSD1 knockdown confers protection against osteoclast formation by reducing histone 3 lysine 9 monomethylation and dimethylation in ITGB3 promoter

ITGB3, an osteoclast marker, is involved in osteoclast formation. Nevertheless, its related mechanism remains poorly characterized. Herein, this study examines the mechanisms affecting osteoclast formation with the involvement of ITGB3. Osteoclast formation was induced with macrophage colony-stimulating factor (M-CSF) and receptor activator of nuclear factor-kappa B ligand (RANKL), followed by measurement of the mRNA and protein expression of ITGB3 and LSD1. After gain- and loss-of-function assays, cell viability and the expression of osteoclast marker genes (NFATc1, ACP5, and CTSK) were assessed, and osteoclast formation was evaluated with TRAP staining. ChIP assays were used to examine histone 3 lysine 9 (H3K9) monomethylation (H3K9me1) and H3K9 dimethylation (H3K9me2) modifications and LSD1 protein enrichment in the ITGB3 promoter. During osteoclast formation, ITGB3 and LSD1 were gradually augmented. Knockdown of LSD1 or ITGB3 curbed cell viability, the expression of osteoclast marker genes, and osteoclast formation. Moreover, overexpression of ITGB3 nullified the suppressive impact of LSD1 knockdown on osteoclast formation. Mechanistically, LSD1 promoted ITGB3 expression by reducing H3K9 levels in the ITGB3 promoter. LSD1 enhanced ITGB3 expression by decreasing H3K9me1 and H3K9me2 levels in ITGB3 promoter to boost osteoclast formation.

Divergence in interspecific and intersubspecific gene expression between two closely related horseshoe bats (Rhinolophus)

Closely related species have been used as representative systems to investigate the genetic mechanisms involved in the early stages of species differentiation. Previous studies have indicated that variation in gene expression might be a sensitive indicator of initial species divergence, although the role of expression divergence, and especially that associated with phenotypic variation remained relatively undefined. For three organs (cochlea, brain, and liver) from two closely related bat species (Rhinolophus siamensis and R. episcopus), the interspecific and intersubspecific gene expression profiles were compared using transcriptomics in this study. Striking organ specificity of expression was observed, and expression profiles exhibited similarities between cochlea and brain tissues. Numerous differentially expressed genes (DEGs) were identified for each organ in the interspecific comparison (cochlea/brain/liver: 1,069/647/692) and intersubspecific comparison (608/528/368). Functional enrichment analysis indicated vital variation in expression related to the immune system, ion activities, neuronal function, and multisensory system regulation in both comparisons. DEGs relevant to the variation in echolocation calls (RF) were found, and some of them were involved in the pivotal patterns of expression variation. The regulation of immune, ion channel, neural activity, and sophisticated sensory functions at the expression level might be key mechanisms in the early species divergence of bats, and the expression variation related to acoustical signal could have played a crucial part. This study expands our knowledge of gene expression and patterns of variation for three key organs to echolocation at both the interspecific and intersubspecific levels. Further, the framework described here provides insight into the genetic basis of phenotypic variation during the incipient stage of species differentiation.

Extensive Adaptive Variation in Gene Expression within and between Closely Related Horseshoe Bats (Chiroptera, Rhinolophus) Revealed by Three Organs

In the process of species differentiation and adaption, the relative influence of natural selection on gene expression variation often remains unclear (especially its impact on phenotypic divergence). In this study, we used differentially expressed genes from brain, cochlea, and liver samples collected from two species of bats to determine the gene expression variation forced by natural selection when comparing at the interspecific (Rhinolophus siamensis and R. episcopus episcopus) and the intraspecific (R. e. episcopus and R. episcopus spp.) levels. In both cases, gene expression variation was extensively adaptive (>66.0%) and mainly governed by directional selection, followed by stabilizing selection, and finally balancing selection. The expression variation related to acoustic signals (resting frequency, RF) and body size (forearm length, FA) was also widely governed by natural selection (>69.1%). Different functional patterns of RF- or FA-related adaptive expression variation were found between the two comparisons, which manifested as abundant immune-related regulations between subspecies (indicating a relationship between immune response and phenotypic adaption). Our study verifies the extensive adaptive expression variation between both species and subspecies and provides insight into the effects of natural selection on species differentiation and adaptation as well as phenotypic divergence at the expression level.

microRNA-25-3p suppresses osteogenic differentiation of BMSCs in patients with osteoporosis by targeting ITGB3

This study was conducted to investigate the impact of the microRNA (miR)-25-3p/ITGB3 axis on the osteogenic differentiation of bone marrow mesenchymal stem cells (BMSCs) from patients with osteoporosis (OP). BMSCs isolated from the bone marrow of healthy controls and OP patients were identified by flow cytometry, in which ITGB3 mRNA and miR-25-3p expression was detected by RT-qPCR and ITGB3, Runx2, OPN, ALP, and OSX protein expression by western blot. The binding between ITGB3 and miR-25-3p was assessed by dual-luciferase reporter gene and Ago2-RIP assays. BMSC osteogenic differentiation was observed by alizarin red staining and ALP activity. The differentiation of BMSCs to adipocytes and chondrocytes was measured by oil red O staining and alcian blue staining, respectively. BMSCs were successfully isolated from the bone marrow of healthy controls (normal-BMSCs) and OP patients (OP-BMSCs). ITGB3, Runx2, OPN, ALP, and OSX expression was poorer and miR-25-3p expression was higher in OP-BMSCs than in normal-BMSCs. Mechanistically, ITGB3 was negatively targeted by miR-25-3p. After osteogenic, adipogenic, and chondrogenic differentiation of BMSCs were successfully induced, adipogenic differentiation increased and osteogenic and chondrogenic differentiation decreased in OP-BMSCs compared with normal-BMSCs. Overexpression of ITGB3 facilitated mineralized nodule formation and elevated ALP activity and Runx2, OPN, and ALP expression in OP-BMSCs. miR-25-3p upregulation diminished mineralized nodule formation, ALP activity, and Runx2, OPN, and ALP expression in OP-BMSCs and normal-BMSCs, which was annulled by additional ITGB3 overexpression. miR-25-3p targets ITGB3, thereby suppressing osteogenic differentiation of BMSCs from OP patients.

Identification of Key Candidate Genes for Beak Length Phenotype by Whole-Genome Resequencing in Geese

The domestic goose is an important economic animal in agriculture and its beak, a trait with high heritability, plays an important role in promoting food intake and defending against attacks. In this study, we sequenced 772 420-day-old Xingguo gray geese (XGG) using a low-depth (~1 ×) whole-genome resequencing strategy. We detected 12,490,912 single nucleotide polymorphisms (SNPs) using the standard GATK and imputed with STITCH. We then performed a genome-wide association study on the beak length trait in XGG. The results indicated that 57 SNPs reached genome-wide significance levels for the beak length trait and were assigned to seven genes, including TAPT1, DHX15, CCDC149, LGI2, SEPSECS, ANAPC4, and Slc34a2. The different genotypes of the most significant SNP (top SNP), which was located upstream of LGI2 and explained 7.24% of the phenotypic variation in beak length, showed significant differences in beak length. Priority-based significance analysis concluded that CCDC149, LGI2, and SEPSECS genes in the most significant quantitative trait locus interval were the most plausible positional and functional candidate genes for beak length development in the XGG population. These findings not only enhance our understanding of the genetic mechanism of the beak length phenotype in geese, but also lay the foundation for further studies to facilitate the genetic selection of traits in geese.

Determining the variation in premaxillary and dentary bone morphology that may underlie beak shape between two pure layer lines

Beak treatment is an effective method of reducing the damage inflicted by severe feather pecking (SFP) but there is significant pressure to eliminate these treatments and rely solely on alternative strategies. Substantial variation in beak shape exists within non-beak treated layer flocks and beak shape appears to be heritable. There is the potential to use this pre-existing variation and genetically select for hens whose beak shapes are less apt to cause damage during SFP. To do this, we must first understand the range of phenotypes that exist for both the external beak shape and the bones that provide its structure. The objective of this study was to determine the variation in premaxillary (within the top beak) and dentary (within the bottom beak) bone morphology that exists in 2 non-beak treated pure White Leghorn layer lines using geometric morphometrics to analyze radiographs. Lateral head radiographs were taken of 825 hens and the premaxillary and dentary bones were landmarked. Landmark coordinates were standardized by Procrustes superimposition and the covariation was analyzed by principal components analysis and multivariate regression using Geomorph (an R package). Three principal components (PCs) explained 85% of total premaxillary bone shape variation and showed that the shape ranged from long and narrow with pointed bone tips to short and wide with more curved tips. Two PCs explained 81% of total dentary bone shape variation. PC1 described the dentary bone length and width and PC2 explained the angle between the bone tip and its articular process. For both bones, shape was significantly associated with bone size and differed significantly between the two lines. Bone size accounted for 42% of the total shape variation for both bones. Together, the results showed a range of phenotypic variation in premaxillary and dentary bone shape, which in turn may influence beak shape. These bone phenotypes will guide further quantitative genetic and behavioral analyses that will help identify which beaks shapes cause the least damage when birds engage in SFP.

Comparative Genomics Reveals Evolution of a Beak Morphology Locus in a High-Altitude Songbird

The Ground Tit (Pseudopodoces humilis) has lived on the Qinghai-Tibet Plateau for ∼5.7 My and has the highest altitudinal distribution among all parids. This species has evolved an elongated beak in response to long-term selection imposed by ground-foraging and cavity-nesting habits, yet the genetic basis for beak elongation remains unknown. Here, we perform genome-wide analyses across 14 parid species and identify 25 highly divergent genomic regions that are significantly associated with beak length, finding seven candidate genes involved in bone morphogenesis and remolding. Neutrality tests indicate that a model allowing for a selective sweep in the highly conserved COL27A1 gene best explains variation in beak length. We also identify two nonsynonymous fixed mutations in the collagen domain that are predicted to be functionally deleterious yet may have facilitated beak elongation. Our study provides evidence of adaptive alleles in COL27A1 with major effects on beak elongation of Ps. humilis.

The Genomic Landscape of Divergence Across the Speciation Continuum in Island-Colonising Silvereyes (Zosterops lateralis)

Inferring the evolutionary dynamics at play during the process of speciation by analyzing the genomic landscape of divergence is a major pursuit in population genomics. However, empirical assessments of genomic landscapes under varying evolutionary scenarios that are known a priori are few, thereby limiting our ability to achieve this goal. Here we combine RAD-sequencing and individual-based simulations to evaluate the genomic landscape of divergence in the silvereye (Zosterops lateralis). Using pairwise comparisons that differ in divergence timeframe and the presence or absence of gene flow, we document how genomic patterns accumulate along the speciation continuum. In contrast to previous predictions, our results provide limited support for the idea that divergence accumulates around loci under divergent selection or that genomic islands widen with time. While a small number of genomic islands were found in populations diverging with and without gene flow, in few cases were SNPs putatively under selection tightly associated with genomic islands. The transition from localized to genome-wide levels of divergence was captured using individual-based simulations that considered only neutral processes. Our results challenge the ubiquity of existing verbal models that explain the accumulation of genomic differences across the speciation continuum and instead support the idea that divergence both within and outside of genomic islands is important during the speciation process.

The role of core and variable Gene Regulatory Network modules in tooth development and evolution

Among the developmental processes that have been proposed to influence the direction of evolution, the modular organization of developmental gene regulatory networks (GRNs) has shown particular promise. In theory, GRNs have core modules comprised of essential, conserved circuits of genes, and sub-modules of downstream, secondary circuits of genes that are more susceptible to variation. While this idea has received considerable interest as of late, the field of evo-devo lacks the experimental systems needed to rigorously evaluate this hypothesis. Here, we introduce an experimental system, the vertebrate tooth, that has great potential as a model for testing this hypothesis. Tooth development and its associated GRN have been well studied and modeled in both model and non-model organisms. We propose that the existence of modules within the tooth GRN explains both the conservation of developmental mechanisms and the extraordinary diversity of teeth among vertebrates. Based on experimental data, we hypothesize that there is a conserved core module of genes that is absolutely necessary to ensure tooth or cusp initiation and development. In regard to tooth shape variation between species, we suggest that more relaxed sub-modules activated at later steps of tooth development, e.g., during the morphogenesis of the tooth and its cusps, control the different axes of tooth morphological variation.

Rapid phenotypic evolution with shallow genomic differentiation during early stages of high elevation adaptation in Eurasian Tree Sparrows

Known as the ‘third polar region’, the Qinghai-Tibet Plateau represents one of the harshest highland environments in the world and yet a number of organisms thrive there. Previous studies of birds, animals and humans have focused on well-differentiated populations in later stages of phenotypic divergence. The adaptive processes during the initial phase of highland adaptation remain poorly understood. We studied a human commensal, the Eurasian Tree Sparrow, which has followed human agriculture to the Qinghai-Tibet Plateau. Despite strong phenotypic differentiation at multiple levels, in particular in muscle-related phenotypes, highland and lowland populations show shallow genomic divergence and the colonization event occurred within the past few thousand years. In a one-month acclimation experiment investigating phenotypic plasticity, we exposed adult lowland tree sparrows to a hypoxic environment and did not observe muscle changes. Through population genetic analyses, we identified a signature of polygenic adaptation, whereby shifts in allele frequencies are spread across multiple loci, many of which are associated with muscle-related processes. Our results reveal a case of positive selection in which polygenic adaptation appears to drive rapid phenotypic evolution, shedding light on early stages of adaptive evolution to a novel environment.

Comparative Transcriptome Analysis of Unusual Localized Skin Laxity in Sika Deer (Cervus nippon)

Cutis laxa represents a heterogeneous group of rare, inherited, or acquired connective tissue disorders with the common feature of loose and redundant skin with decreased elasticity. The skin of affected deer showed abnormal collagen fiber morphology. To identify the differentially expressed genes of the unusual localized skin laxity in sika deer, we performed transcriptome analysis in the affected and control sika deer. The transcriptome analysis showed 700 genes with significant differential expression in the affected skin as compared with normal skin. Pathway analysis revealed an enrichment of genes involved in tumor necrosis factor signaling, the extracellular matrix-receptor interaction, platelet activation, and Huntington's disease. A gene network was constructed, and the hub nodes such as PTGS2, THBS1, COL1A1, FOS, and NOS3 were found through PPI network analysis, which may contributed to the unusual localized skin laxity in sika deer. Abnormal expression patterns of genes during the development of the affected sika deer were successfully uncovered in the present study, which provides a reference for revealing the related mechanism underlying cutis laxa in sika deer and human beings.

Neural crest and the origin of species-specific pattern

For well over half of the 150 years since the discovery of the neural crest, the special ability of these cells to function as a source of species-specific pattern has been clearly recognized. Initially, this observation arose in association with chimeric transplant experiments among differentially pigmented amphibians, where the neural crest origin for melanocytes had been duly noted. Shortly thereafter, the role of cranial neural crest cells in transmitting species-specific information on size and shape to the pharyngeal arch skeleton as well as in regulating the timing of its differentiation became readily apparent. Since then, what has emerged is a deeper understanding of how the neural crest accomplishes such a presumably difficult mission, and this includes a more complete picture of the molecular and cellular programs whereby neural crest shapes the face of each species. This review covers studies on a broad range of vertebrates and describes neural-crest-mediated mechanisms that endow the craniofacial complex with species-specific pattern. A major focus is on experiments in quail and duck embryos that reveal a hierarchy of cell-autonomous and non-autonomous signaling interactions through which neural crest generates species-specific pattern in the craniofacial integument, skeleton, and musculature. By controlling size and shape throughout the development of these systems, the neural crest underlies the structural and functional integration of the craniofacial complex during evolution.