The zebrafish intestine and liver, as in other vertebrates, are derived from the endoderm. Great effort has been devoted to deciphering the molecular mechanisms controlling the specification and development of the zebrafish intestine and liver; however, genome-wide comparison of the transcriptomes between these two organs at the larval stage remains unexplored. There is a lack of extensive identification of feature genes marking specific cell types in the zebrafish intestine and liver at 5 days post-fertilization, when the larval fish starts food intake. In this report, through RNA sequencing and single-cell RNA sequencing of intestines and livers separately dissected from wild-type zebrafish larvae at 5 days post-fertilization, together with the experimental validation of 47 genes through RNA whole-mount in situ hybridization, we identified not only distinctive transcriptomes for the larval intestine and liver, but also a considerable number of feature genes for marking the intestinal bulb, mid-intestine and hindgut, and for marking hepatocytes and cholangiocytes. Meanwhile, we identified 135 intestine- and 97 liver-enriched transcription factor genes in zebrafish larvae at 5 days post-fertilization. Our findings provide rich molecular and cellular resources for studying cell patterning and specification during the early development of the zebrafish intestine and liver.

The hepatopancreatic duct (HPD) system links the liver and pancreas to the intestinal tube and is composed of the extrahepatic biliary duct, gallbladder and pancreatic duct. Haematopoietically-expressed-homeobox (Hhex) protein plays an essential role in the establishment of HPD, however, the molecular mechanism remains elusive. Here we show that zebrafish hhex -null mutants fail to develop the HPD system characterized by lacking the biliary marker Annexin A4 and the HPD marker sox9b . The mutant HPD system is replaced by an intrahepatic intestinal tube characterized by expressing the intestinal marker fatty-acid-binding–protein 2a ( fabp2a ). Cell lineage analysis showed that this intrahepatic intestinal tube is not originated from hepatocytes or cholangiocytes. Further analysis revealed that cdx1b and pdx1 were expressed ectopically in the intrahepatic intestinal tube and knockdown of cdx1b and pdx1 restored the expression of sox9b in the mutant. Chromatin-immunoprecipitation analysis shows that Hhex binds to the promoters of pdx1 and cdx1b genes to repress their expression. We therefore propose that Hhex, Cdx1b and Pdx1 form a genetic network governing the patterning and morphogenesis of the HPD and digestive tract systems in zebrafish.