MGA-seq: robust identification of extrachromosomal DNA and genetic variants using multiple genetic abnormality sequencing.

Genomic abnormalities are strongly associated with cancer and infertility. In this study, we develop a simple and efficient method - multiple genetic abnormality sequencing (MGA-Seq) - to simultaneously detect structural variation, copy number variation, single-nucleotide polymorphism, homogeneously staining regions, and extrachromosomal DNA (ecDNA) from a single tube. MGA-Seq directly sequences proximity-ligated genomic fragments, yielding a dataset with concurrent genome three-dimensional and whole-genome sequencing information, enabling approximate localization of genomic structural variations and facilitating breakpoint identification. Additionally, by utilizing MGA-Seq, we map focal amplification and oncogene coamplification, thus facilitating the exploration of ecDNA's transcriptional regulatory function.

Highly efficient and robust π-FISH rainbow for multiplexed in situ detection of diverse biomolecules.

In the unprecedented single-cell sequencing and spatial multiomics era of biology, fluorescence in situ hybridization (FISH) technologies with higher sensitivity and robustness, especially for detecting short RNAs and other biomolecules, are greatly desired. Here, we develop the robust multiplex π-FISH rainbow method to detect diverse biomolecules (DNA, RNA, proteins, and neurotransmitters) individually or simultaneously with high efficiency. This versatile method is successfully applied to detect gene expression in different species, from microorganisms to plants and animals. Furthermore, we delineate the landscape of diverse neuron subclusters by decoding the spatial distribution of 21 marker genes via only two rounds of hybridization. Significantly, we combine π-FISH rainbow with hybridization chain reaction to develop π-FISH+ technology for short nucleic acid fragments, such as microRNA and prostate cancer anti-androgen therapy-resistant marker ARV7 splicing variant in circulating tumour cells from patients. Our study provides a robust biomolecule in situ detection technology for spatial multiomics investigation and clinical diagnosis.

Transcriptomic and metabolomic insights into the variety of sperm storage in oviduct of egg layers.

In birds, the sperm storage tubules (SST) are dispersed in uterovaginal junction (UVJ) and highly correlated with differential capacity of sperm storage (SS) in and among species with unspecified mechanisms. Here, the SS duration of 252 egg layer breeders was evaluated in 5 rounds with 3 phenotypic traits to screen high- and low-SS individuals, respectively, followed with transcriptome of UVJ tissues and metabolome of serum (high-SS vs. low-SS) to decipher the candidate genes and biochemical markers correlated with differential SS capacity. Histological characterization suggested slightly higher density of SST in UVJ (high-SS vs. low-SS). Transcriptome analyses identified 596 differentially expressed genes (336 upregulated vs. 260 downregulated), which were mainly enriched in gene ontology terms of homeostasis, steroid and lipid metabolism and hormone activity, and 12 significant pathways (P < 0.05) represented by calcium, steroid, and lipid metabolism. Immunohistochemical staining of GNAQ, ST6GAL1, ADFP, and PCNA showed similar distribution in UVJ tissues between 2 groups. Several candidates (HSD11B2, DIO2, AQP3, GNAQ, NANS, ST6GAL1) combined with 4 (11ß-prostaglandin F2α, prostaglandin B1, 7α-hydroxytestosterone, and N-acetylneuraminic acid) of 40 differential metabolites enriched in serum metabolome were considered as regulators and biomarkers of SS duration in egg layer breeders. The integrated transcriptome and metabolome analyses of chicken breeder hens will provide novel insights for exploration and improvement of differential SS capacity in birds.

Deciphering the Forebrain Disorder in a Chicken Model of Cerebral Hernia.

Cerebral hernia in crested chicken has been characterized as the protrusion of cerebral hemispheres into the unsealed skull for hundreds of years, since Charles Darwin. The development of deformed forebrain (telencephalon) of cerebral hernia remains largely unknown. Here, the unsealed frontal skull combined with misplaced sphenoid bone was observed and potentially associated with brain protuberance. The shifted pallidum, elongated hippocampus, expanded mesopallium and nidopallium, and reduced hyperpallium were observed in seven regions of the malformed telencephalon. The neurons were detected with nuclear pyknosis and decreased density. Astrocytes showed uneven distribution and disordered protuberances in hyperpallium and hippocampus. Transcriptome analyses of chicken telencephalon (cerebral hernia vs. control) revealed 547 differentially expressed genes (DEGs), mainly related to nervous system development, and immune system processes, including astrocyte marker gene GFAP, and neuron and astrocyte developmental gene S100A6. The upregulation of GFAP and S100A6 genes in abnormal telencephalon was correlated with reduced DNA methylation levels in the promoter regions. The morphological, cellular, and molecular variations in the shape, regional specification, and cellular states of malformed telencephalon potentially participate in brain plasticity and previously reported behavior changes. Chickens with cerebral hernia might be an interesting and valuable disease model to further explore the recognition, diagnosis, and therapy of cerebral hernia development of crested chickens and other species.

Exploration of key regulators driving primary feather follicle induction in goose skin.

The primary feather follicles are universal skin appendages widely distributed in the skin of feathered birds. The morphogenesis and development of the primary feather follicles in goose skin remain largely unknown. Here, the induction of primary feather follicles in goose embryonic skin (pre-induction vs induction) was investigated by de novo transcriptome analyses to reveal 409 differentially expressed genes (DEGs). The DEGs were characterized to potentially regulate the de novo formation of feather follicle primordia consisting of placode (4 genes) and dermal condensate (12 genes), and the thickening of epidermis (5 genes) and dermal fibroblasts (17 genes), respectively. Further analyses enriched DEGs into GO terms represented as cell adhesion and KEGG pathways including Wnt and Hedgehog signaling pathways that are highly correlated with cell communication and molecular regulation. Six selected Wnt pathway genes were detected by qPCR with up-regulation in goose skin during the induction of primary feather follicles. The localization of WNT16, SFRP1 and FRZB by in situ hybridization showed weak expression in the primary feather primordia, whereas FZD1, LEF1 and DKK1 were expressed initially in the inter-follicular skin and feather follicle primordia, then mainly restricted in the feather primordia. The spatial-temporal expression patterns indicate that Wnt pathway genes DKK1, FZD1 and LEF1 are the important regulators functioned in the induction of primary feather follicle in goose skin. The dynamic molecular changes and specific gene expression patterns revealed in this report provide the general knowledge of primary feather follicle and skin development in waterfowl, and contribute to further understand the diversity of hair and feather development beyond the mouse and chicken models.

Expression patterns of three JAK-STAT pathway genes in feather follicle development during chicken embryogenesis.

The Janus kinase (JAK)-signal transducer and activator of transcription (STAT) (JAK-STAT) pathway is shown to restrain the hair follicles in catagen and telogen and prevent anagen reentry in murine hair follicle cycling. The early roles of JAK-STAT pathway genes in skin development remain uncharacterized in mouse and chicken models. Here, we revealed the expression patterns of three JAK-STAT pathway genes (JAK1, JAK2, and TYK2) in chicken embryonic skin at E6-E10 stages which are key to feather follicle morphogenesis. Multiple sequence alignment of the three genes from chicken and other species all showed a closely related homology with birds like quail and goose. Whole mount in situ hybridization (WISH) revealed weak expression of JAK1, JAK2, and TYK2 in chicken skin at E6 and E7, and followed with the focally restricted signals in the feather follicles of neck and body skin located dorsally at E8 for JAK1, E9 for TYK2 and E10 for JAK2 gene. All three genes displayed stronger expression in feather follicles of neck skin than that of body skin. The expression levels of JAK1 and TYK2 were much stronger than those of JAK2. Quantitative real-time PCR (qRT-PCR) analysis revealed the increased expression tendency for JAK2 both in the neck and body skin from E6 to E10, and the much stronger expression in neck and body skin at later stages (E8-E10) than earlier stages (E6 and E7) for JAK1 and TYK2. Overall, these findings suggest that JAK1 and TYK2, not JAK2 are important to specify the feather follicle primordia, and to arrange the proximal-distal axis of feather follicles, respectively, during the morphogenesis of feather follicles in embryonic chicken skin.

Defining Key Genes Regulating Morphogenesis of Apocrine Sweat Gland in Sheepskin.

The apocrine sweat gland is a unique skin appendage in humans compared to mouse and chicken models. The absence of apocrine sweat glands in chicken and murine skin largely restrains further understanding of the complexity of human skin biology and skin diseases, like hircismus. Sheep may serve as an additional system for skin appendage investigation owing to the distributions and histological similarities between the apocrine sweat glands of sheep trunk skin and human armpit skin. To understand the molecular mechanisms underlying morphogenesis of apocrine sweat glands in sheepskin, transcriptome analyses were conducted to reveal 1631 differentially expressed genes that were mainly enriched in three functional groups (cellular component, molecular function and biological process), particularly in gland, epithelial, hair follicle and skin development. There were 7 Gene Ontology (GO) terms enriched in epithelial cell migration and morphogenesis of branching epithelium that were potentially correlated with the wool follicle peg elongation. An additional 5 GO terms were enriched in gland morphogenesis (20 genes), gland development (42 genes), salivary gland morphogenesis and development (8 genes), branching involved in salivary gland morphogenesis (6 genes) and mammary gland epithelial cell differentiation (4 genes). The enriched gland-related genes and two Kyoto Encyclopedia of Genes and Genomes pathway genes (WNT and TGF-ß) were potentially involved in the induction of apocrine sweat glands. Genes named BMPR1A, BMP7, SMAD4, TGFB3, WIF1, and WNT10B were selected to validate transcript expression by qRT-PCR. Immunohistochemistry was performed to localize markers for hair follicle (SOX2), skin fibroblast (PDGFRB), stem cells (SOX9) and BMP signaling (SMAD5) in sheepskin. SOX2 and PDGFRB were absent in apocrine sweat glands. SOX9 and SMAD5 were both observed in precursor cells of apocrine sweat glands and later in gland ducts. These results combined with the upregulation of BMP signaling genes indicate that apocrine sweat glands were originated from outer root sheath of primary wool follicle and positively regulated by BMP signaling. This report established the primary network regulating early development of apocrine sweat glands in sheepskin and will facilitate the further understanding of histology and pathology of apocrine sweat glands in human and companion animal skin.