Benchmarking spatial and single-cell transcriptomics integration methods for transcript distribution prediction and cell type deconvolution.

Spatial transcriptomics approaches have substantially advanced our capacity to detect the spatial distribution of RNA transcripts in tissues, yet it remains challenging to characterize whole-transcriptome-level data for single cells in space. Addressing this need, researchers have developed integration methods to combine spatial transcriptomic data with single-cell RNA-seq data to predict the spatial distribution of undetected transcripts and/or perform cell type deconvolution of spots in histological sections. However, to date, no independent studies have comparatively analyzed these integration methods to benchmark their performance. Here we present benchmarking of 16 integration methods using 45 paired datasets (comprising both spatial transcriptomics and scRNA-seq data) and 32 simulated datasets. We found that Tangram, gimVI, and SpaGE outperformed other integration methods for predicting the spatial distribution of RNA transcripts, whereas Cell2location, SpatialDWLS, and RCTD are the top-performing methods for the cell type deconvolution of spots. We provide a benchmark pipeline to help researchers select optimal integration methods to process their datasets.

Extended depth of field for Fresnel zone aperture camera via fast passive depth estimation.

The lensless camera with incoherent illumination has gained significant research interest for its thin and flexible structure. However, it faces challenges in resolving scenes with a wide depth of field (DoF) due to its depth-dependent point spread function (PSF). In this paper, we present a single-shot method for extending the DoF in Fresnel zone aperture (FZA) cameras at visible wavelengths through passive depth estimation. The improved ternary search method is utilized to determine the depth of targets rapidly by evaluating the sharpness of the back propagation reconstruction. Based on the depth estimation results, a set of reconstructed images focused on targets at varying depths are derived from the encoded image. After that, the DoF is extended through focus stacking. The experimental results demonstrate an 8-fold increase compared with the calibrated DoF at 130 mm depth. Moreover, our depth estimation method is five times faster than the traversal method, while maintaining the same level of accuracy. The proposed method facilitates the development of lensless imaging in practical applications such as photography, microscopy, and surveillance.

Address model mismatch and defocus in FZA lensless imaging via model-driven CycleGAN.

Mask-based lensless imaging systems suffer from model mismatch and defocus. In this Letter, we propose a model-driven CycleGAN, MDGAN, to reconstruct objects within a long distance. MDGAN includes two translation cycles for objects and measurements respectively, each consisting of a forward propagation and a backward reconstruction module. The backward module resembles the Wiener-U-Net, and the forward module consists of the estimated image formation model of a Fresnel zone aperture camera (FZACam), followed by CNN to compensate for the model mismatch. By imposing cycle consistency, the backward module can adaptively match the actual depth-varying imaging process. We demonstrate that MDGAN based on either a simulated or calibrated imaging model produces a higher-quality image compared to existing methods. Thus, it can be applied to other mask-based systems.

Role of Traditional Chinese Medicine for the Treatment of Lupus nephritis: Mechanisms and Applications.

Background: Lupus nephritis (LN), caused by Systemic lupus erythematosus, is a chronic autoimmune renal disease and a key risk factor for morbidity and fatality involving 50% damage to the kidney. LN is associated with aberrant functioning of the immune system, characterized by increased systemic inflammation, altered lymphocyte count, perivascular infiltration of inflammatory cells, and declined organ functioning. Current Therapies and Limitations: Conventional therapies to LN include high-dose glucocorticoids, immunosuppressants, calcineurin inhibitors, immune boosters, and targeted medicines that improve kidney functioning. However, these drugs triggered severe adverse side effects, and their prolonged usage resulted in drug resistance, accentuating LN complications. TCM in LN Treatment: Hence, safe and functional Traditional Chinese Medicines (TCM), with supporting clinical trials and observational studies, received significant recognition worldwide for the Treatment of LN. In the form of herbal extracts and preparations, TCM proved effective in treating immunodeficient disorders, including LN. Additionally, acupuncture as a TCM appeared promising in reducing LN-induced inflammation and joint pain. Mechanisms and Benefits: The therapeutic mechanisms included reduced antibody generation, pro-inflammatory cytokine release, immune complex formation, complement activation, extracellular matrix damage and proteinuria levels that played vital roles in chronic kidney diseases. They generated immunosuppressive effects by modulating apoptosis, oxidative stress, and inflammatory signaling pathways, such as JAK/STAT, NF-κB, AP-1 and MAPK and their cross-talk in LN and associated renal injury. These therapies improved blood circulation, alleviated renal pathological changes, restored glomerular capillary functioning, and regenerated renal tissues. However, an essential requisite for these therapies for LN included reduced side effects and improved hepatoprotection and detoxification. Clinical studies suggest that TCM formulations may demonstrate therapeutic benefits in alleviating the symptoms of LN, suggesting prospects of combined applications with Western medicine to enhance treatment efficiencies. Overall, TCM is beneficial for treating LN, and may serve as a potential alternative to conventional medicines.

A ∼ 4.1 kb deletion in IRX1 gene upstream is completely associated with rumplessness in Piao chicken.

Piao chicken, a Chinese indigenous rumpless chicken breed, lacks pygostyle, caudal vertebra, uropygial gland and tail feathers. The rumplessness in Piao chicken presents an autosomal dominant inheritance pattern. However, the molecular genetic mechanisms underlying the rumplessness in Piao chicken remains unclear. In this study, whole-genome resequencing was performed for 146 individuals from 10 chicken breeds, including 9 tailed chicken breeds and Piao rumpless breed. Tailbone CT scan for Piao chickens and WL chickens, revealed that some Piao chicken tails were normal in number, and for a few Piao chickens tail length and tail bone numbers were between the rumpless and the normal tailed chickens. The results showed that the rumpless phenotype has not been completely fixed in Piao chicken breed. Using selection signature analysis and structural variation detection, we found a 4174 bp deletion located in the upstream region of IRX1 gene on chromosome 2 related to rumpless phenotype. Structural variation genotyping showed that the deletion was present in all 32 rumpless Piao chickens (del/del, wild/del) and absent from all 112 tailed chickens included in the dataset for the other 9 breeds and 2 tailed Piao chickens (wild/wild). In summary, all rumpless Piao chickens tested here carry this deletion mutation, to show a complete linkage association with rumplessness trait. We suggested that the 4174 bp deletion could be causative for rumpless phenotype in Piao chicken since this is the only mutation to show the complete linkage disequilibrium with rumplessness on whole genome level across all of 146 chickens from the 10 breeds. This study could facilitate a better understanding of the genetic characteristics of Piao chicken.

Genomic and transcriptomic analyses reveal genetic adaptation to cold conditions in the chickens.

Under the pressure of natural and artificial selection, domestic animals, including chickens, have evolved unique mechanisms of genetic adaptations such as high-altitude adaptation, hot and arid climate adaptation, and desert adaptation. Here, we investigated the genetic basis of cold tolerance in chicken by integrating whole-genome and transcriptome sequencing technologies. Genome-wide comparative analyses of 118 chickens living in different latitudes showed 46 genes and several pathways that may be involved in cold adaptation. The results of the functional enrichment analysis of differentially expressed genes proved the important role of metabolic pathways and immune-related pathways in cold tolerance in chickens. The subsequent integration of whole genome and transcriptome sequencing technology further identified six genes - dnah5 (dynein axonemal heavy chain 5), ptgs2 (prostaglandin-endoperoxide synthase 2), inhba (inhibin beta A subunit), irx2 (iroquois homeobox 2), ensgalg00000054917, and ensgalg00000046652 - requiring more detailed studies. In addition, we also discovered different allele frequency distributions of five SNPs (single nucleotide polymorphisms) within ptgs2 and nine SNPs within dnah5 in chickens in different latitudes, suggesting strong selective pressure of these two genes in chickens. We provide a novel insight into the genetic adaptation in chickens to cold environments, and provide a reference for evaluating and developing adaptive chicken breeds in cold environments.

Systematic Selection Signature Analysis of Chinese Gamecocks Based on Genomic and Transcriptomic Data.

Selection pressures driven by natural causes or human interference are key factors causing genome variants and signatures of selection in specific regions of the genome. Gamecocks were bred for cockfighting, presenting pea-combs, larger body sizes, stronger limbs, and higher levels of aggression than other chickens. In this study, we aimed to explore the genomic differences between Chinese gamecocks and commercial, indigenous, foreign, and cultivated breeds by detecting the regions or sites under natural or artificial selection using genome-wide association studies (GWAS), genome-wide selective sweeps based on the genetic differentiation index (FST), and transcriptome analyses. Ten genes were identified using GWAS and FST: gga-mir-6608-1, SOX5, DGKB, ISPD, IGF2BP1, AGMO, MEOX2, GIP, DLG5, and KCNMA1. The ten candidate genes were mainly associated with muscle and skeletal development, glucose metabolism, and the pea-comb phenotype. Enrichment analysis results showed that the differentially expressed genes between the Luxi (LX) gamecock and Rhode Island Red (RIR) chicken were mainly related to muscle development and neuroactive-related pathways. This study will help to understand the genetic basis and evolution of Chinese gamecocks and support the further use of gamecocks as an excellent breeding material from a genetic perspective.

Fine mapping of a quantitative resistance gene for gray leaf spot of maize (Zea mays L.) derived from teosinte (Z. mays ssp. parviglumis).

KEY MESSAGE: In this study we mapped the QTL Qgls8 for gray leaf spot (GLS) resistance in maize to a ~130 kb region on chromosome 8 including five predicted genes. In previous work, using near isogenic line (NIL) populations in which segments of the teosinte (Zea mays ssp. parviglumis) genome had been introgressed into the background of the maize line B73, we had identified a QTL on chromosome 8, here called Qgls8, for gray leaf spot (GLS) resistance. We identified alternate teosinte alleles at this QTL, one conferring increased GLS resistance and one increased susceptibility relative to the B73 allele. Using segregating populations derived from NIL parents carrying these contrasting alleles, we were able to delimit the QTL region to a ~130 kb (based on the B73 genome) which encompassed five predicted genes.

Genetic dissection of the maize (Zea mays L.) MAMP response.

KEY MESSAGE: Loci associated with variation in maize responses to two microbe-associated molecular patterns (MAMPs) were identified. MAMP responses were correlated. No relationship between MAMP responses and quantitative disease resistance was identified. Microbe-associated molecular patterns (MAMPs) are highly conserved molecules commonly found in microbes which can be recognized by plant pattern recognition receptors. Recognition triggers a suite of responses including production of reactive oxygen species (ROS) and nitric oxide (NO) and expression changes of defense-related genes. In this study, we used two well-studied MAMPs (flg22 and chitooctaose) to challenge different maize lines to determine whether there was variation in the level of responses to these MAMPs, to dissect the genetic basis underlying that variation and to understand the relationship between MAMP response and quantitative disease resistance (QDR). Naturally occurring quantitative variation in ROS, NO production, and defense genes expression levels triggered by MAMPs was observed. A major quantitative traits locus (QTL) associated with variation in the ROS production response to both flg22 and chitooctaose was identified on chromosome 2 in a recombinant inbred line (RIL) population derived from the maize inbred lines B73 and CML228. Minor QTL associated with variation in the flg22 ROS response was identified on chromosomes 1 and 4. Comparison of these results with data previously obtained for variation in QDR and the defense response in the same RIL population did not provide any evidence for a common genetic basis controlling variation in these traits.

Unconditional and conditional QTL analysis of kernel weight related traits in wheat (Triticum aestivum L.) in multiple genetic backgrounds.

Wheat thousand kernel weight (TKW) is a complex trait, and is largely controlled by several kernel traits, including kernel length (KL) and kernel width (KW). In order to reveal the genetic relationship between TKW and these kernel traits (KW and KL) as accurate as possible, we applied both unconditional and conditional mapping analyses to three distinct genetic populations, one DH population and two RIL populations. This report describes the identifications of 36 unconditional and conditional additive QTLs and 30 pairs of unconditional and conditional epistatic QTLs, all of which are closely associated with TKW. While the conditional additive locus Qtkw1B, detected in the RIL2 population, exhibited the largest contribution, explaining 14.12 % of TKW variance, the unconditional epistatic QTLs Qtkw3A-2/Qtkw5B.1, detected in the DH population, accounted for 11.95 % of phenotypic variance. This study also showed that, compared with unconditional mapping, conditional mapping resulted in very different numbers and different extent of effects of additive and epistatic QTLs that were associated with TKW when TKW was conditioned on kernel traits (KW and KL). These data strongly suggest that KW and KL indeed play a significant role in determining TKW. Furthermore, we demonstrated that the effects of the 25 additive QTLs for TKW were either entirely or largely determined by KW, while the effects of the other 25 additive QTLs for TKW were either entirely or largely affected by KL. We conclude that the conditional mapping can be useful for a better understanding of the interrelationship between the yield contributing traits at the QTL level.

Injectable Tannin-Containing Hydroxypropyl Chitin Hydrogel as Novel Bioactive Pulp Capping Material Accelerates Repair of Inflamed Dental Pulp.

Conventional pulp capping materials have limited anti-inflammatory capacity. It is necessary to develop more effective pulp capping material for the treatment of inflamed pulps. Tannic acid (TA) is a natural, water-soluble polyphenol with antimicrobial and anti-inflammatory properties. This study aimed to investigate the effects of a tannin-containing hydroxypropyl chitin hydrogel (HPCH/TA hydrogel) as an innovative pulp capping material. The physicochemical properties of the composite hydrogels were characterized. The effects of HPCH/TA hydrogel as a pulp capping material were evaluated in vitro and in vivo. The underlying mechanism of the anti-inflammatory effects of HPCH/TA hydrogel was explored. The HPCH/TA hydrogel demonstrated favorable temperature sensitivity, injectability, and antibacterial properties. In vitro, the HPCH/TA hydrogel effectively promoted the proliferation of human dental pulp cells and inhibited interleukin-1ß, interleukin-6, and tumor necrosis factor-α expression, possibly by suppressing the nuclear factor kappa-B pathway. In vivo, on the fourth day after capping, the HPCH/TA hydrogel group showed lower inflammatory scores compared to the control and iRoot BP Plus (commercial pulp capping material) group. By the sixth week, complete reparative dentin formation was observed in the HPCH/TA hydrogel group, with no difference in thickness compared to the iRoot BP Plus group. Collectively, the HPCH/TA hydrogel holds promise as a bioactive pulp capping material for promoting the repair of inflamed pulp in vital pulp therapy.

An Optimized Human Erythroblast Differentiation System Reveals Cholesterol-Dependency of Robust Production of Cultured Red Blood Cells Ex Vivo.

The generation of cultured red blood cells (cRBCs) ex vivo represents a potentially unlimited source for RBC transfusion and other cell therapies. Human cRBCs can be generated from the terminal differentiation of proliferating erythroblasts derived from hematopoietic stem/progenitor cells or erythroid precursors in peripheral blood mononuclear cells. Efficient differentiation and maturation into cRBCs highly depend on replenishing human plasma, which exhibits variable potency across donors or batches and complicates the consistent cRBC production required for clinical translation. Hence, the role of human plasma in erythroblast terminal maturation is investigated and uncovered that 1) a newly developed cell culture basal medium mimicking the metabolic profile of human plasma enhances cell growth and increases cRBC yield upon erythroblast terminal differentiation and 2) LDL-carried cholesterol, as a substitute for human plasma, is sufficient to support erythroid survival and terminal differentiation ex vivo. Consequently, a chemically-defined optimized medium (COM) is developed, enabling robust generation of cRBCs from erythroblasts of multiple origins, with improved enucleation efficiency and higher reticulocyte yield, without the need for supplementing human plasma or serum. In addition, the results reveal the crucial role of lipid metabolism during human terminal erythropoiesis.

Engineered extracellular vesicles enable high-efficient delivery of intracellular therapeutic proteins.

Developing an intracellular delivery system is of key importance in the expansion of protein-based therapeutics acting on cytosolic or nuclear targets. Recently, extracellular vesicles (EVs) have been exploited as next-generation delivery modalities due to their natural role in intercellular communication and biocompatibility. However, fusion of protein of interest to a scaffold represents a widely used strategy for cargo enrichment in EVs, which could compromise the stability and functionality of cargo. Herein, we report intracellular delivery via EV-based approach (IDEA) that efficiently packages and delivers native proteins both in vitro and in vivo without the use of a scaffold. As a proof-of-concept, we applied the IDEA to deliver cyclic GMP-AMP synthase (cGAS), an innate immune sensor. The results showed that cGAS-carrying EVs activated interferon signaling and elicited enhanced antitumor immunity in multiple syngeneic tumor models. Combining cGAS EVs with immune checkpoint inhibition further synergistically boosted antitumor efficacy in vivo. Mechanistically, scRNA-seq demonstrated that cGAS EVs mediated significant remodeling of intratumoral microenvironment, revealing a pivotal role of infiltrating neutrophils in the antitumor immune milieu. Collectively, IDEA, as a universal and facile strategy, can be applied to expand and advance the development of protein-based therapeutics.

Molecular genetic foundation of a sex-linked tailless trait in Hongshan chicken by whole genome data analysis.

As a Chinese local chicken breed, Hongshan chickens have 2 kinds of tail feather phenotypes, normal and taillessness. Our previous studies showed that taillessness was a sex-linked dominant trait. Abnormal development of the tail vertebrae could be explained this phenomenon in some chicken breeds. However, the number of caudal vertebrae in rumpless Hongshan chickens was normal, so rumplessness in Hongshan chicken was not related to the development of the caudal vertebrae. Afterwards, we found that rumplessness in Hongshan was due to abnormal development of tail feather rather than abnormal development of caudal vertebrae. In order to understand the genetic foundation of the rumplessness of Hongshan chickens, we compared and reanalyzed 2 sets of data in normal and rumpless Hongshan chickens from our previous studies. By joint analysis of genome-wide selection signature analysis and genome-wide association approach, we found that 1 overlapping gene (EDIL3) and 16 peak genes (ENSGALG00000051843, ENSGALG00000053498, ENSGALG00000054800, KIF27, PTPRD, ENSGALG00000047579, ENSGALG00000041052, ARHGEF28, CAMK4, SERINC5, ENSGALG00000050776, ERCC8, MCC, ADAMTS19, ENSGALG00000053322, CHRNA8) located on the Z chromosome was associated with the rumpless trait. The results of this study furtherly revealed the molecular mechanism of the rumpless trait in Hongshan chickens, and identified the candidate genes associated with this trait. Our results will help to improve the shape of chicken tail feathers and to rise individual economic value in some specific market in China.

Whole genome resequencing reveals genomic regions related to red plumage in ducks.

Plumage color is a characteristic trait of ducks that originates as a result of natural and artificial selection. As a conspicuous phenotypic feature, it is a breed characteristic. Previous studies have identified some genes associated with the formation of black and white plumage in ducks. However, studies on the genetic basis underlying the red plumage phenotype in ducks are limited. Here, genome-wide association analysis (GWAS) and selection signal detection (Fst, θπ ratio, and cross-population composite likelihood ratio [XP-CLR]) were conducted to identify candidate regions and genes underlying duck plumage color phenotype. Selection signal detection revealed 29 overlapping genes (including ENPP1 and ULK1) significantly associated with red plumage color in Ji'an Red ducks. ENSAPLG00000012679, ESRRG, and SPATA5 were identified as candidate genes associated with red plumage using GWAS. Selection signal detection revealed that 19 overlapping genes (including GMDS, PDIA6, and ODC1) significantly correlated with light brown plumage in Brown Tsaiya ducks. GWAS to narrow down the significant regions further revealed nine candidate genes (AKT1, ATP6V1C2, GMDS, LRP4, MAML3, PDIA6, PLD5, TMEM63B, and TSPAN8). Notably, in Brown Tsaiya ducks, GMDS, ODC1, and PDIA6 exhibit significantly differentiated allele frequencies among other feather-colored ducks, while in Ji'an Red ducks, ENSAPLG00000012679 has different allele frequency distributions compared with that in other feather-colored ducks. This study offers new insights into the variation and selection of the red plumage phenotype using GWAS and selective signals.

Genetic Foundation of Male Spur Length and Its Correlation with Female Egg Production in Chickens.

Spurs, which mainly appear in roosters, are protrusions near the tarsometatarsus on both sides of the calves of chickens, and are connected to the tarsometatarsus by a bony core. As a male-biased morphological characteristic, the diameter and length of spurs vary significantly between different individuals, mainly related to genetics and age. As a specific behavior of hens, egg-laying also varies greatly between individuals in terms of traits such as age at first egg (AFE), egg weight (EW), and so on. At present, there are few studies on chicken spurs. In this study, we investigated the inheritance pattern of the spur trait in roosters with different phenotypes and the correlations between spur length, body weight at 18 weeks of age (BW18), shank length at 18 weeks of age (SL18), and the egg-laying trait in hens (both hens and roosters were from the same population and were grouped according to their family). These traits related to egg production included AFE, body weight at first egg (BWA), and first egg weight (FEW). We estimated genetic parameters based on pedigree and phenotype data, and used variance analysis to calculate broad-sense heritability for correcting the parameter estimation results. The results showed that the heritability of male left and right spurs ranged from 0.6 to 0.7. There were significant positive correlations between left and right spur length, BW18, SL18, and BWA, as well as between left and right spur length and AFE. We selected 35 males with the longest spurs and 35 males with the shortest spurs in the population, and pooled them into two sets to obtain the pooled genome sequencing data. After genome-wide association and genome divergency analysis by FST, allele frequency differences (AFDs), and XPEHH methods, we identified 7 overlapping genes (CENPE, FAT1, FAM149A, MANBA, NFKB1, SORBS2, UBE2D3) and 14 peak genes (SAMD12, TSPAN5, ENSGALG00000050071, ENSGALG00000053133, ENSGALG00000050348, CNTN5, TRPC6, ENSGALG00000047655,TMSB4X, LIX1, CKB, NEBL, PRTFDC1, MLLT10) related to left and right spur length through genome-wide selection signature analysis and a genome-wide association approach. Our results identified candidate genes associated with chicken spurs, which helps to understand the genetic mechanism of this trait and carry out subsequent research around it.

Significant genomic introgression from grey junglefowl (Gallus sonneratii) to domestic chickens (Gallus gallus domesticus).

BACKGROUND: Chicken is one of the most numerous and widely distributed species around the world, and many studies support the multiple ancestral origins of domestic chickens. The research regarding the yellow skin phenotype in domestic chickens (regulated by BCO2) likely originating from the grey junglefowl serves as crucial evidence for demonstrating the multiple origins of chickens. However, beyond the BCO2 gene region, much remains unknown about the introgression from the grey junglefowl into domestic chickens. Therefore, in this study, based on whole-genome data of 149 samples including 4 species of wild junglefowls and 13 local domestic chicken breeds, we explored the introgression events from the grey junglefowl to domestic chickens. RESULTS: We successfully detected introgression regions besides BCO2, including two associated with growth trait (IGFBP2 and TKT), one associated with angiogenesis (TIMP3) and two members of the heat shock protein family (HSPB2 and CRYAB). Our findings suggest that the introgression from the grey junglefowl may impact the growth performance of chickens. Furthermore, we revealed introgression events from grey junglefowl at the BCO2 region in multiple domestic chicken breeds, indicating a phenomenon where the yellow skin phenotype likely underwent strong selection and was retained. Additionally, our haplotype analysis shed light on BCO2 introgression event from different sources of grey junglefowl into domestic chickens, possibly suggesting multiple genetic flows between the grey junglefowl and domestic chickens. CONCLUSIONS: In summary, our findings provide evidences of the grey junglefowl contributing to the genetic diversity of domestic chickens, laying the foundation for a deeper understanding of the genetic composition within domestic chickens, and offering new perspectives on the impact of introgression on domestic chickens.

Genetic parameter estimation and molecular foundation of chicken egg-laying trait.

The age of first egg (AFE) in chicken can affect early and even life-time egg production performance to some extent, and therefore is an important economic trait that affects production efficiency. To better understand the genetic patterns of AFE and other production traits including body weight at first egg (BWA), first egg weight (FEW), and total egg number from AFE to 58 wk of age (total-EN), we recorded the production performance of 2 widely used layer breeds, white leghorn (WL) and Rhode Island Red (RIR) and estimated genetic parameters based on pedigree and production data. The results showed that the heritability of AFE in both breeds ranged from 0.4 to 0.6, and AFE showed strong positive genetic and phenotypic correlations to BWA as well as FEW, while showing strong negative genetic and phenotypic correlations with total-EN. Furtherly, by genome-wide association analysis study (GWAS), we identified 12 and 26 significant SNPs to be related to AFE in the 2-layer breeds, respectively. A total of 18 genes were identified that could affect AFE based on the significant SNP annotations obtained, but there were no gene overlapped in the 2 breeds indicating the genetic foundation of AFE could differ from breed to breed. Our results provided a deeper understanding of genetic patterns and molecular basement of AFE in different breeds and could help in the selection of egg production traits.

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.

Evolution, Inheritance, and Strata Formation of the W Chromosome in Duck (Anas platyrhynchos).

The nonrecombining female-limited W chromosome is predicted to experience unique evolutionary processes. Difficulties in assembling W chromosome sequences have hindered the identification of duck W-linked sequences and their evolutionary footprint. To address this, we conducted three initial contig-level genome assemblies and developed a rigorous pipeline by which to successfully expand the W-linked data set, including 11 known genes and 24 newly identified genes. Our results indicate that the W chromosome expression may not be subject to female-specific selection; a significant convergent pattern of upregulation associated with increased female-specific selection was not detected. The genetic stability of the W chromosome is also reflected in the strong evolutionary correlation between it and the mitochondria; the complete consistency of the cladogram topology constructed from their gene sequences proves the shared maternal coevolution. By detecting the evolutionary trajectories of W-linked sequences, we have found that recombination suppression started in four distinct strata, of which three were conserved across Neognathae. Taken together, our results have revealed a unique evolutionary pattern and an independent stratum evolutionary pattern for sex chromosomes.