A pangenome reference of 36 Chinese populations.

Human genomics is witnessing an ongoing paradigm shift from a single reference sequence to a pangenome form, but populations of Asian ancestry are underrepresented. Here we present data from the first phase of the Chinese Pangenome Consortium, including a collection of 116 high-quality and haplotype-phased de novo assemblies based on 58 core samples representing 36 minority Chinese ethnic groups. With an average 30.65× high-fidelity long-read sequence coverage, an average contiguity N50 of more than 35.63 megabases and an average total size of 3.01 gigabases, the CPC core assemblies add 189 million base pairs of euchromatic polymorphic sequences and 1,367 protein-coding gene duplications to GRCh38. We identified 15.9 million small variants and 78,072 structural variants, of which 5.9 million small variants and 34,223 structural variants were not reported in a recently released pangenome reference1. The Chinese Pangenome Consortium data demonstrate a remarkable increase in the discovery of novel and missing sequences when individuals are included from underrepresented minority ethnic groups. The missing reference sequences were enriched with archaic-derived alleles and genes that confer essential functions related to keratinization, response to ultraviolet radiation, DNA repair, immunological responses and lifespan, implying great potential for shedding new light on human evolution and recovering missing heritability in complex disease mapping.

Genomes, fossils, and the concurrent rise of modern birds and flowering plants in the Late Cretaceous.

The phylogeny and divergence timing of the Neoavian radiation remain controversial despite recent progress. We analyzed the genomes of 124 species across all Neoavian orders, using data from 25,460 loci spanning four DNA classes, including 5,756 coding sequences, 12,449 conserved nonexonic elements, 4,871 introns, and 2,384 intergenic segments. We conducted a comprehensive sensitivity analysis to account for the heterogeneity across different DNA classes, leading to an optimal tree of Neoaves with high resolution. This phylogeny features a novel Neoavian dichotomy comprising two monophyletic clades: a previously recognized Telluraves (land birds) and a newly circumscribed Aquaterraves (waterbirds and relatives). Molecular dating analyses with 20 fossil calibrations indicate that the diversification of modern birds began in the Late Cretaceous and underwent a constant and steady radiation across the KPg boundary, concurrent with the rise of angiosperms as well as other major Cenozoic animal groups including placental and multituberculate mammals. The KPg catastrophe had a limited impact on avian evolution compared to the Paleocene-Eocene Thermal Maximum, which triggered a rapid diversification of seabirds. Our findings suggest that the evolution of modern birds followed a slow process of gradualism rather than a rapid process of punctuated equilibrium, with limited interruption by the KPg catastrophe. This study places bird evolution into a new context within vertebrates, with ramifications for the evolution of the Earth's biota.

Phylogenetic inference of inter-population transmission rates for infectious diseases.

Estimating transmission rates is a challenging yet essential aspect of comprehending and controlling the spread of infectious diseases. Various methods exist for estimating transmission rates, each with distinct assumptions, data needs, and constraints. This study introduces a novel phylogenetic approach called transRate, which integrates genetic information with traditional epidemiological approaches to estimate inter-population transmission rates. The phylogenetic method is statistically consistent as the sample size (i.e. the number of pathogen genomes) approaches infinity under the multi-population susceptible-infected-recovered model. Simulation analyses indicate that transRate can accurately estimate the transmission rate with a sample size of 200 ~ 400 pathogen genomes. Using transRate, we analyzed 40,028 high-quality sequences of SARS-CoV-2 in human hosts during the early pandemic. Our analysis uncovered significant transmission between populations even before widespread travel restrictions were implemented. The development of transRate provides valuable insights for scientists and public health officials to enhance their understanding of the pandemic's progression and aiding in preparedness for future viral outbreaks. As public databases for genomic sequences continue to expand, transRate is increasingly vital for tracking and mitigating the spread of infectious diseases.

Historic samples reveal loss of wild genotype through domestic chicken introgression during the Anthropocene.

Human activities have precipitated a rise in the levels of introgressive gene flow among animals. The investigation of conspecific populations at different time points may shed light on the magnitude of human-mediated introgression. We used the red junglefowl Gallus gallus, the wild ancestral form of the chicken, as our study system. As wild junglefowl and domestic chickens readily admix, conservationists fear that domestic introgression into junglefowl may compromise their wild genotype. By contrasting the whole genomes of 51 chickens with 63 junglefowl from across their natural range, we found evidence of a loss of the wild genotype across the Anthropocene. When comparing against the genomes of junglefowl from approximately a century ago using rigorous ancient-DNA protocols, we discovered that levels of domestic introgression are not equal among and within modern wild populations, with the percentage of domestic ancestry around 20-50%. We identified a number of domestication markers in which chickens are deeply differentiated from historic junglefowl regardless of breed and/or geographic provenance, with eight genes under selection. The latter are involved in pathways dealing with development, reproduction and vision. The wild genotype is an allelic reservoir that holds most of the genetic diversity of G. gallus, a species which is immensely important to human society. Our study provides fundamental genomic infrastructure to assist in efforts to prevent a further loss of the wild genotype through introgression of domestic alleles.

Gene flow and an anomaly zone complicate phylogenomic inference in a rapidly radiated avian family (Prunellidae).

BACKGROUND: Resolving the phylogeny of rapidly radiating lineages presents a challenge when building the Tree of Life. An Old World avian family Prunellidae (Accentors) comprises twelve species that rapidly diversified at the Pliocene-Pleistocene boundary. RESULTS: Here we investigate the phylogenetic relationships of all species of Prunellidae using a chromosome-level de novo assembly of Prunella strophiata and 36 high-coverage resequenced genomes. We use homologous alignments of thousands of exonic and intronic loci to build the coalescent and concatenated phylogenies and recover four different species trees. Topology tests show a large degree of gene tree-species tree discordance but only 40-54% of intronic gene trees and 36-75% of exonic genic trees can be explained by incomplete lineage sorting and gene tree estimation errors. Estimated branch lengths for three successive internal branches in the inferred species trees suggest the existence of an empirical anomaly zone. The most common topology recovered for species in this anomaly zone was not similar to any coalescent or concatenated inference phylogenies, suggesting presence of anomalous gene trees. However, this interpretation is complicated by the presence of gene flow because extensive introgression was detected among these species. When exploring tree topology distributions, introgression, and regional variation in recombination rate, we find that many autosomal regions contain signatures of introgression and thus may mislead phylogenetic inference. Conversely, the phylogenetic signal is concentrated to regions with low-recombination rate, such as the Z chromosome, which are also more resistant to interspecific introgression. CONCLUSIONS: Collectively, our results suggest that phylogenomic inference should consider the underlying genomic architecture to maximize the consistency of phylogenomic signal.

Troponin T1 in tumorigenesis and immune modulation: Insights into multiple cancers and kidney renal clear cell carcinoma.

Troponin T1 (TNNT1) plays a crucial role in muscle contraction but its role in cancer, particularly in kidney renal clear cell carcinoma (KIRC), is not well-understood. This study explores the expression, clinical significance and biological functions of TNNT1 in various cancers, with an emphasis on its involvement in KIRC. We analysed TNNT1 expression in cancers using databases like TCGA and GTEx, assessing its prognostic value, mutation patterns, methylation status and functional implications. The study also examined TNNT1's effect on the tumour microenvironment and drug sensitivity in KIRC, complemented by in vitro TNNT1 knockdown experiments in KIRC cells. TNNT1 is overexpressed in several cancers and linked to adverse outcomes, showing frequent upregulation mutations and abnormal methylation. Functionally, TNNT1 connects to muscle and cancer pathways, affects immune infiltration and drug responses, and its overexpression in KIRC is associated with advanced disease and reduced survival. Knocking down TNNT1 curbed KIRC cell growth. TNNT1's aberrant expression plays a significant role in tumorigenesis and immune modulation, highlighting its value as a prognostic biomarker and a potential therapeutic target in KIRC and other cancers. Further studies are essential to understand TNNT1's oncogenic mechanisms in KIRC.

Fusobacterium nucleatum-derived small extracellular vesicles facilitate tumor growth and metastasis via TLR4 in breast cancer.

BACKGROUND: The contributive role of the microbiome in tumor progression has been reported in multiple studies, such as the Fusobacterium nucleatum (F. nucleatum) in breast cancer (BC). This study aimed to explore the role of F. nucleatum-derived small extracellular vesicles (Fn-EVs) in BC and preliminarily uncover the mechanism. METHODS: Ten normal and 20 cancerous breast tissues were harvested to investigate the gDNA expression of F. nucleatum and its relation with the clinical characteristics of BC patients. After isolating Fn-EVs by ultracentrifugation from F. nucleatum (ATCC 25,586), both MDA-MB-231 and MCF-7 cells were treated with PBS, Fn, or Fn-EVs, followed by being subjected to CCK-8, Edu staining, wound healing, and Transwell assays to detect their cell viability, proliferation, migration, and invasion. TLR4 expression in BC cells with diverse treatments was assessed by western blot. In vivo experiments were performed to verify its role in tumor growth and liver metastasis. RESULTS: The F. nucleatum gDNA levels of breast tissues in BC patients were significantly higher than those in normal subjects, and positively associated with tumor size and metastasis. Fn-EVs administration significantly enhanced the cell viability, proliferation, migration, and invasion of BC cells, while knocking down TLR4 in BC cells could block these effects. Furthermore, in vivo study verified the contributive role of Fn-EVs in tumor growth and metastasis of BC, which might rely on its regulation of TLR4. CONCLUSIONS: Collectively, our results suggest that F. nucleatum plays an important role in BC tumor growth and metastasis by regulating TLR4 through Fn-EVs. Thus, a better understanding of this process may aid in the development of novel therapeutic agents.

The combinational nano-immunotherapy of ferumoxytol and poly(I:C) inhibits melanoma via boosting anti-angiogenic immunity.

Angiogenesis plays a key role in the progression and metastasis of melanoma, and the pro-angiogenic effect of macrophages is one major reason for the failure of current anti-angiogenic therapies. Here, a nano-immunotherapy combining ferumoxytol and poly(I:C) (ferumoxytol/poly(I:C)) has been developed to boost the anti-angiogenic activities of macrophages to inhibit melanoma. Our findings demonstrated that ferumoxytol/poly(I:C) was a highly efficacious anti-tumor therapy with limited toxicity. Both in vivo and in vitro experiments indicated that this combination was successful in impeding angiogenesis. Ferumoxytol/poly(I:C) was demonstrated to reduce the viability of endothelial cells, thus hindering tube formation. Particularly, ferumoxytol/poly(I:C) was able to polarize macrophages to the M1 phenotype and decrease the expression of vascular endothelial growth factor, which in turn amplified the anti-angiogenic properties of ferumoxytol/poly(I:C). This combination of ferumoxytol/poly(I:C) nano-immunotherapy enriches the anti-angiogenic therapeutic nature of ferumoxytol and will shed new light on the treatment of melanoma.

Optimization scheme of machine learning model for genetic division between northern Han, southern Han, Korean and Japanese.

Han Chinese, Korean and Japanese are the main populations of East Asia, and Han Chinese presents a gradient admixture from north to south. There are differences among the East Asian populations in genetic structure. To achieve fine-scale genetic classification of southern (S-) and northern (N-) Han Chinese, Korean and Japanese individuals in this study, we collected and analyzed 1185 ancestry informative SNPs (AISNPs) from previous literature reports and our laboratory findings. First, two machine learning algorithms, softmax and randomForest, were used to build genetic classification models. Then, phylogenetic tree, STRUCTURE and principal component analysis were used to evaluate the performance of classification for different AISNP panels. The 234-AISNP panel achieved a fine-scale differentiation among the target populations in four classification schemes. The accuracy of the softmax model was 92%, which realized the accurate classification of the S-Han, N-Han, Korean and Japanese individuals. The two machine learning models tested in this study provided important references for the high-resolution discrimination of close-range populations and will be useful tools to optimize marker panels for developing forensic DNA ancestry inference systems.

Genetic expression and mutational profile analysis in different pathologic stages of hepatocellular carcinoma patients.

BACKGROUND: The clinical pathologic stages (stage I, II, III-IV) of hepatocellular carcinoma (HCC) are closely linked to the clinical prognosis of patients. This study aims at investigating the gene expression and mutational profile in different clinical pathologic stages of HCC. METHODS: Based on the TCGA-LIHC cohort, we utilized a series of analytical approaches, such as statistical analysis, random forest, decision tree, principal component analysis (PCA), to identify the differential gene expression and mutational profiles. The expression patterns of several targeting genes were also verified by analyzing the Chinese HLivH060PG02 HCC cohort, several GEO datasets, HPA database, and diethylnitrosamine-induced HCC mouse model. RESULTS: We identified a series of targeting genes with copy number variation, which is statistically associated with gene expression. Non-synonymous mutations mainly existed in some genes (e.g.,TTN, TP53, CTNNB1). Nevertheless, no association between gene mutation frequency and pathologic stage distribution was detected. The random forest and decision tree modeling analysis data showed a group of genes related to different HCC pathologic stages, including GAS2L3 and SEMA3F. Additionally, our PCA data indicated several genes associated with different pathologic stages, including SNRPA and SNRPD2. Compared with adjacent normal tissues, we observed a highly expressed level of GAS2L3, SNRPA, and SNRPD2 (P = 0.002) genes in HCC tissues of our HLivH060PG02 cohort. We also detected the high expression pattern of GAS2L3, SEMA3F, SNRPA, and SNRPD2 in the datasets of GSE102079, GSE76427, GSE64041, GSE121248, GSE84005, and the qPCR assay using diethylnitrosamine-induced HCC mouse model. Moreover, SEMA3F and SNRPD2 protein were highly stained in the HCC tissues of the HPA database. The high expression level of these four genes was associated with the poor survival prognosis of HCC cases. CONCLUSIONS: Our study provides evidence regarding the gene expression and mutational profile in different clinical pathologic stages of TCGA HCC cases. Identifying four targeting genes, including GAS2L3, SNRPA, SNRPD2, and SEMA3F, offers insight into the molecular mechanisms associated with different prognoses of HCC.

Differential genetic mutations of ectoderm, mesoderm, and endoderm-derived tumors in TCGA database.

BACKGROUND: In terms of biological behavior, gene regulation, or signaling pathways, there is a certain similarity between tumorigenesis and embryonic development of humans. Three germ layer structure exhibits the distinct ability to form specific tissues and organs. METHODS: The present study set out to investigate the genetic mutation characteristics of germ layer differentiation-related genes using the tumor cases of the cancer genome atlas (TCGA) database. RESULTS: These tumor samples were divided into three groups, including the ectoderm, mesoderm, and endoderm. Children cases less than 9 years old accounted for a larger proportion for the cases in the ectoderm and mesoderm groups; whereas the middle-aged and elderly individuals (from 50 to 89 years old) were more susceptible to tumors of endoderm. There was a better prognosis for the cases of mesoderm, especially the male with the race of White, compared with the other groups. A missense mutation was frequently detected for the cases of ectoderm and endoderm, while deletion mutation was common for that of mesoderm. We could not identify the ectoderm, mesoderm, or endoderm-specific mutated genes or variants with high mutation frequency. However, there was a relatively higher mutation incidence of endoderm markers (GATA6, FOXA2, GATA4, AFP) in the endoderm group, compared with the groups of ectoderm and mesoderm. Additionally, four members (SMO, GLI1, GLI2, GLI3) within the Hedgehog signaling pathway genes showed a relatively higher mutation rate in the endoderm group than the other two groups. CONCLUSIONS: TCGA tumors of ectoderm, mesoderm, and endoderm groups exhibit the distinct subject distribution, survival status, and genomic alteration characteristics. The synergistic mutation effect of specific genes closely related to embryonic development may contribute to the tumorigenesis of tissues or organs derived from the specific germ layers. This study provides a novel reference for exploring the functional connection between embryogenesis and tumorigenesis.

Genomic evidence reveals a radiation of placental mammals uninterrupted by the KPg boundary.

The timing of the diversification of placental mammals relative to the Cretaceous-Paleogene (KPg) boundary mass extinction remains highly controversial. In particular, there have been seemingly irreconcilable differences in the dating of the early placental radiation not only between fossil-based and molecular datasets but also among molecular datasets. To help resolve this discrepancy, we performed genome-scale analyses using 4,388 loci from 90 taxa, including representatives of all extant placental orders and transcriptome data from flying lemurs (Dermoptera) and pangolins (Pholidota). Depending on the gene partitioning scheme, molecular clock model, and genic deviation from molecular clock assumptions, extensive sensitivity analyses recovered widely varying diversification scenarios for placental mammals from a given gene set, ranging from a deep Cretaceous origin and diversification to a scenario spanning the KPg boundary, suggesting that the use of suboptimal molecular clock markers and methodologies is a major cause of controversies regarding placental diversification timing. We demonstrate that reconciliation between molecular and paleontological estimates of placental divergence times can be achieved using the appropriate clock model and gene partitioning scheme while accounting for the degree to which individual genes violate molecular clock assumptions. A birth-death-shift analysis suggests that placental mammals underwent a continuous radiation across the KPg boundary without apparent interruption by the mass extinction, paralleling a genus-level radiation of multituberculates and ecomorphological diversification of both multituberculates and therians. These findings suggest that the KPg catastrophe evidently played a limited role in placental diversification, which, instead, was likely a delayed response to the slightly earlier radiation of angiosperms.

The effect of alignment uncertainty, substitution models and priors in building and dating the mammal tree of life.

BACKGROUND: The flood of genomic data to help build and date the tree of life requires automation at several critical junctures, most importantly during sequence assembly and alignment. It is widely appreciated that automated alignment protocols can yield inaccuracies, but the relative impact of various sources error on phylogenomic analysis is not yet known. This study employs an updated mammal data set of 5162 coding loci sampled from 90 species to evaluate the effects of alignment uncertainty, substitution models, and fossil priors on gene tree, species tree, and divergence time estimation. Additionally, a novel coalescent likelihood ratio test is introduced for comparing competing species trees against a given set of gene trees. RESULTS: The aligned DNA sequences of 5162 loci from 90 species were trimmed and filtered using trimAL and two filtering protocols. The final dataset contains 4 sets of alignments - before trimming, after trimming, filtered by a recently proposed pipeline, and further filtered by comparing ML gene trees for each locus with the concatenation tree. Our analyses suggest that the average discordance among the coalescent trees is significantly smaller than that among the concatenation trees estimated from the 4 sets of alignments or with different substitution models. There is no significant difference among the divergence times estimated with different substitution models. However, the divergence dates estimated from the alignments after trimming are more recent than those estimated from the alignments before trimming. CONCLUSIONS: Our results highlight that alignment uncertainty of the updated mammal data set and the choice of substitution models have little impact on tree topologies yielded by coalescent methods for species tree estimation, whereas they are more influential on the trees made by concatenation. Given the choice of calibration scheme and clock models, divergence time estimates are robust to the choice of substitution models, but removing alignments deemed problematic by trimming algorithms can lead to more recent dates. Although the fossil prior is important in divergence time estimation, Bayesian estimates of divergence times in this data set are driven primarily by the sequence data.

Identification and expression analysis of GRAS transcription factors in the wild relative of sweet potato Ipomoea trifida.

BACKGROUND: GRAS gene is an important transcription factor gene family that plays a crucial role in plant growth, development, adaptation to adverse environmental condition. Sweet potato is an important food, vegetable, industrial raw material, and biofuel crop in the world, which plays an essential role in food security in China. However, the function of sweet potato GRAS genes remains unknown. RESULTS: In this study, we identified and characterised 70 GRAS members from Ipomoea trifida, which is the progenitor of sweet potato. The chromosome distribution, phylogenetic tree, exon-intron structure and expression profiles were analysed. The distribution map showed that GRAS genes were randomly located in 15 chromosomes. In combination with phylogenetic analysis and previous reports in Arabidopsis and rice, the GRAS proteins from I. trifida were divided into 11 subfamilies. Gene structure showed that most of the GRAS genes in I. trifida lacked introns. The tissue-specific expression patterns and the patterns under abiotic stresses of ItfGRAS genes were investigated via RNA-seq and further tested by RT-qPCR. Results indicated the potential functions of ItfGRAS during plant development and stress responses. CONCLUSIONS: Our findings will further facilitate the functional study of GRAS gene and molecular breeding of sweet potato.

Genome-wide identification, structural and gene expression analysis of the bZIP transcription factor family in sweet potato wild relative Ipomoea trifida.

BACKGROUND: The basic leucine zipper (bZIP) transcription factor is one of the most abundant and conserved transcription factor families. In addition to being involved in growth and development, bZIP transcription factors also play an important role in plant adaption to abiotic stresses. RESULTS: A total of 41 bZIP genes that encode 66 proteins were identified in Ipomoea trifida. They were distributed on 14 chromosomes of Ipomoea trifida. Segmental and tandem duplication analysis showed that segmental duplication played an important role in the ItfbZIP gene amplification. ItfbZIPs were divided into ten groups (A, B, C, D, E, F, G, H, I and S groups) according to their phylogenetic relationships with Solanum lycopersicum and Arabidopsis thaliana. The regularity of the exon/intron numbers and distributions is consistent with the group classification in evolutionary tree. Prediction of the cis-acting elements found that promoter regions of ItfbZIPs harbored several stress responsive cis-acting elements. Protein three-dimensional structural analysis indicated that ItfbZIP proteins mainly consisted of α-helices and random coils. The gene expression pattern from transcriptome data and qRT-PCR analysis showed that ItfbZIP genes expressed with a tissue-specific manner and differently expressed under various abiotic stresses, suggesting that the ItfbZIPs were involved in stress response and adaption in Ipomoea trifida. CONCLUSIONS: Genome-wide identification, gene structure, phylogeny and expression analysis of bZIP gene in Ipomoea trifida supplied a solid theoretical foundation for the functional study of bZIP gene family and further facilitated the molecular breeding of sweet potato.

A Jurassic mammaliaform and the earliest mammalian evolutionary adaptations.

The earliest evolution of mammals and origins of mammalian features can be traced to the mammaliaforms of the Triassic and Jurassic periods that are extinct relatives to living mammals. Here we describe a new fossil from the Middle Jurassic that has a mandibular middle ear, a gradational transition of thoracolumbar vertebrae and primitive ankle features, but highly derived molars with a high crown and multiple roots that are partially fused. The upper molars have longitudinal cusp rows that occlude alternately with those of the lower molars. This specialization for masticating plants indicates that herbivory evolved among mammaliaforms, before the rise of crown mammals. The new species shares the distinctive dental features of the eleutherodontid clade, previously represented only by isolated teeth despite its extensive geographic distribution during the Jurassic. This eleutherodontid was terrestrial and had ambulatory gaits, analogous to extant terrestrial mammals such as armadillos or rock hyrax. Its fur corroborates that mammalian integument had originated well before the common ancestor of living mammals.

Genomic differentiation and patterns of gene flow between two long-tailed tit species (Aegithalos).

Patterns of heterogeneous genomic differentiation have been well documented between closely related species, with some highly differentiated genomic regions ("genomic differentiation islands") spread throughout the genome. Differential levels of gene flow are proposed to account for this pattern, as genomic differentiation islands are suggested to be resistant to gene flow. Recent studies have also suggested that genomic differentiation islands could be explained by linked selection acting on genomic regions with low recombination rates. Here, we investigate genomic differentiation and gene-flow patterns for autosomes using RAD-seq data between two closely related species of long-tailed tits (Aegithalos bonvaloti and A. fuliginosus) in both allopatric and contact zone populations. The results confirm recent or ongoing gene flow between these two species. However, there is little evidence that the genomic regions that were found to be highly differentiated between the contact zone populations are resistant to gene flow, suggesting that differential levels of gene flow is not the cause of the heterogeneous genomic differentiation. Linked selection may be the cause of genomic differentiation islands between the allopatric populations with no or very limited gene flow, but this could not account for the heterogeneous genomic differentiation between the contact zone populations, which show evidence of recent or ongoing gene flow.

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

BACKGROUND: Beak morphology exhibits considerable adaptive plasticity in birds, which results in highly varied or specialized forms in response to variations in ecology and life history. As the only parid species endemic to the Qinghai-Tibet Plateau, the Ground Tit (Parus humilis) has evolved a distinctly long and curved beak from other parids. An integration of morphometrics, phylogenetics, transcriptomics and embryology allows us to address the evolutionary and developmental mechanisms of the adaptive beak structure observed in the Ground Tit. RESULTS: A morphometric approach quantified that the Ground Tit has a comparatively longer and more decurved upper beaks than other parids. We estimated that the ancestor of the Ground Tit likely had a short straight upper beak similar to most current recognized parid species using an ancestral state reconstruction. This morphological specialization is considered an adaptation to its ground-oriented behavior on the high plateau. To identify genetic mechanisms behind this adaptive change, a comparative transcriptomic analysis was applied between the Ground Tit and its closely related species, the Great Tit (Parus major). We detected that 623 genes were significantly differentially expressed in embryonic upper beaks between the two species, 17 of which were functionally annotated to correlate with bone development and morphogenesis, although genes related to bone development were not found to undergo accelerated evolution in the Ground Tit. RT-qPCR validation confirmed differential expression of five out of eight genes that were selected from the 17 genes. Subsequent functional assays in chicken embryos demonstrated that two of these genes, FGF13 and ITGB3, may affect beak morphology by modulating levels of osteoblasts and osteoclasts. CONCLUSIONS: Our results provide preliminary evidence that development of the long decurved beak of the Ground Tit is likely regulated by transcriptional activities of multiple genes coordinating osteoblasts and osteoclasts. The integration of multiple approaches employed here sheds light on ecological and genetic mechanisms in the evolution of avian morphology.