Transgenic tools targeting the basal ganglia reveal both evolutionary conservation and specialization of neural circuits in zebrafish.

The cortico-basal ganglia circuit mediates decision making. Here, we generated transgenic tools for adult zebrafish targeting specific subpopulations of the components of this circuit and utilized them to identify evolutionary homologs of the mammalian direct- and indirect-pathway striatal neurons, which respectively project to the homologs of the internal and external segment of the globus pallidus (dorsal entopeduncular nucleus [dEN] and lateral nucleus of the ventral telencephalic area [Vl]) as in mammals. Unlike in mammals, the Vl mainly projects to the dEN directly, not by way of the subthalamic nucleus. Further single-cell RNA sequencing analysis reveals two pallidal output pathways: a major shortcut pathway directly connecting the dEN with the pallium and the evolutionarily conserved closed loop by way of the thalamus. Our resources and circuit map provide the common basis for the functional study of the basal ganglia in a small and optically tractable zebrafish brain for the comprehensive mechanistic understanding of the cortico-basal ganglia circuit.

The serotonergic neurons derived from rhombomere 2 are localized in the median raphe and project to the dorsal pallium in zebrafish.

The serotonergic neurons in the raphe nucleus are implicated in various cognitive functions such as learning and emotion. In vertebrates, the raphe nucleus is divided into the dorsal raphe and the median raphe. In contrast to the abundance of knowledge on the functions of the dorsal raphe, the roles of the serotonergic neurons in the median raphe are relatively unknown. The studies using zebrafish revealed that the median raphe serotonergic neurons receive input from the two distinct pathways from the habenula and the IPN. The use of zebrafish may reveal the function of the Hb-IPN-median raphe pathway. To clarify the functions of the median raphe serotonergic neurons, it is necessary to distinguish them from those in the dorsal raphe. Most median raphe serotonergic neurons originate from rhombomere 2 in mice, and we generated the transgenic zebrafish which can label the serotonergic neurons derived from rhombomere 2. In this study, we found the serotonergic neurons derived from rhombomere 2 are localized in the median raphe and project axons to the rostral dorsal pallium in zebrafish. This study suggests that this transgenic system has the potential to specifically reveal the function and information processing of the Hb-IPN-raphe-telencephalon circuit in learning.

Identification of a Functional Susceptibility Variant for Adolescent Idiopathic Scoliosis that Upregulates Early Growth Response 1 (EGR1)-Mediated UNCX Expression.

Adolescent idiopathic scoliosis (AIS) is a serious health problem affecting 3% of live births all over the world. Many loci associated with AIS have been identified by previous genome wide association studies, but their biological implication remains mostly unclear. In this study, we evaluated the AIS-associated variants in the 7p22.3 locus by combining in silico, in vitro, and in vivo analyses. rs78148157 was located in an enhancer of UNCX, a homeobox gene and its risk allele upregulated the UNCX expression. A transcription factor, early growth response 1 (EGR1), transactivated the rs78148157-located enhancer and showed a higher binding affinity for the risk allele of rs78148157. Furthermore, zebrafish larvae with UNCX messenger RNA (mRNA) injection developed body curvature and defective neurogenesis in a dose-dependent manner. rs78148157 confers the genetic susceptibility to AIS by enhancing the EGR1-regulated UNCX expression. © 2022 American Society for Bone and Mineral Research (ASBMR).

Genome sequence and analysis of the Japanese morning glory Ipomoea nil.

Ipomoea is the largest genus in the family Convolvulaceae. Ipomoea nil (Japanese morning glory) has been utilized as a model plant to study the genetic basis of floricultural traits, with over 1,500 mutant lines. In the present study, we have utilized second- and third-generation-sequencing platforms, and have reported a draft genome of I. nil with a scaffold N50 of 2.88 Mb (contig N50 of 1.87 Mb), covering 98% of the 750 Mb genome. Scaffolds covering 91.42% of the assembly are anchored to 15 pseudo-chromosomes. The draft genome has enabled the identification and cataloguing of the Tpn1 family transposons, known as the major mutagen of I. nil, and analysing the dwarf gene, CONTRACTED, located on the genetic map published in 1956. Comparative genomics has suggested that a whole genome duplication in Convolvulaceae, distinct from the recent Solanaceae event, has occurred after the divergence of the two sister families.

Social conflict resolution regulated by two dorsal habenular subregions in zebrafish.

When animals encounter conflict they initiate and escalate aggression to establish and maintain a social hierarchy. The neural mechanisms by which animals resolve fighting behaviors to determine such social hierarchies remain unknown. We identified two subregions of the dorsal habenula (dHb) in zebrafish that antagonistically regulate the outcome of conflict. The losing experience reduced neural transmission in the lateral subregion of dHb (dHbL)-dorsal/intermediate interpeduncular nucleus (d/iIPN) circuit. Silencing of the dHbL or medial subregion of dHb (dHbM) caused a stronger predisposition to lose or win a fight, respectively. These results demonstrate that the dHbL and dHbM comprise a dual control system for conflict resolution of social aggression.

The habenulo-raphe serotonergic circuit encodes an aversive expectation value essential for adaptive active avoidance of danger.

Anticipation of danger at first elicits panic in animals, but later it helps them to avoid the real threat adaptively. In zebrafish, as fish experience more and more danger, neurons in the ventral habenula (vHb) showed tonic increase in the activity to the presented cue and activated serotonergic neurons in the median raphe (MR). This neuronal activity could represent the expectation of a dangerous outcome and be used for comparison with a real outcome when the fish is learning how to escape from a dangerous to a safer environment. Indeed, inhibiting synaptic transmission from vHb to MR impaired adaptive avoidance learning, while panic behavior induced by classical fear conditioning remained intact. Furthermore, artificially triggering this negative outcome expectation signal by optogenetic stimulation of vHb neurons evoked place avoidance behavior. Thus, vHb-MR circuit is essential for representing the level of expected danger and behavioral programming to adaptively avoid potential hazard.

The habenula is crucial for experience-dependent modification of fear responses in zebrafish.

The zebrafish dorsal habenula (dHb) shows conspicuous asymmetry in its connection with the interpeduncular nucleus (IPN) and is equivalent to the mammalian medial habenula. Genetic inactivation of the lateral subnucleus of dHb (dHbL) biased fish towards freezing rather than the normal flight response to a conditioned fear stimulus, suggesting that the dHbL-IPN pathway is important for controlling experience-dependent modification of fear responses.

Novel mutations affecting axon guidance in zebrafish and a role for plexin signalling in the guidance of trigeminal and facial nerve axons.

In zebrafish embryos, the axons of the posterior trigeminal (Vp) and facial (VII) motoneurons project stereotypically to a small number of target muscles derived from the first and second branchial arches (BA1, BA2). Use of the Islet1 (Isl1)-GFP transgenic line enabled precise real-time observations of the growth cone behaviour of the Vp and VII motoneurons within BA1 and BA2. Screening for N-ethyl-N-nitrosourea-induced mutants identified seven distinct mutations affecting different steps in the axonal pathfinding of these motoneurons. The class 1 mutations caused severe defasciculation and abnormal pathfinding in both Vp and VII motor axons before they reached their target muscles in BA1. The class 2 mutations caused impaired axonal outgrowth of the Vp motoneurons at the BA1-BA2 boundary. The class 3 mutation caused impaired axonal outgrowth of the Vp motoneurons within the target muscles derived from BA1 and BA2. The class 4 mutation caused retraction of the Vp motor axons in BA1 and abnormal invasion of the VII motor axons in BA1 beyond the BA1-BA2 boundary. Time-lapse observations of the class 1 mutant, vermicelli (vmc), which has a defect in the plexin A3 (plxna3) gene, revealed that Plxna3 acts with its ligand Sema3a1 for fasciculation and correct target selection of the Vp and VII motor axons after separation from the common pathways shared with the sensory axons in BA1 and BA2, and for the proper exit and outgrowth of the axons of the primary motoneurons from the spinal cord.

A molecular neuroethological approach for identifying and characterizing a cascade of behaviorally regulated genes.

Songbirds have one of the most accessible neural systems for the study of brain mechanisms of behavior. However, neuroethological studies in songbirds have been limited by the lack of high-throughput molecular resources and gene-manipulation tools. To overcome these limitations, we constructed 21 regular, normalized, and subtracted full-length cDNA libraries from brains of zebra finches in 57 developmental and behavioral conditions in an attempt to clone as much of the brain transcriptome as possible. From these libraries, approximately 14,000 transcripts were isolated, representing an estimated 4,738 genes. With the cDNAs, we created a hierarchically organized transcriptome database and a large-scale songbird brain cDNA microarray. We used the arrays to reveal a set of 33 genes that are regulated in forebrain vocal nuclei by singing behavior. These genes clustered into four anatomical and six temporal expression patterns. Their functions spanned a large range of cellular and molecular categories, from signal transduction, trafficking, and structural, to synaptically released molecules. With the full-length cDNAs and a lentiviral vector system, we were able to overexpress, in vocal nuclei, proteins of representative singing-regulated genes in the absence of singing. This publicly accessible resource http://songbirdtranscriptome.net can now be used to study molecular neuroethological mechanisms of behavior.

Transcriptome analysis of the aphid bacteriocyte, the symbiotic host cell that harbors an endocellular mutualistic bacterium, Buchnera.

Aphids possess bacteriocytes, cells specifically differentiated to harbor obligatory mutualistic bacteria of the genus Buchnera, which have lost many genes that are essential for common bacterial functions. To understand the host's role in maintaining the symbiotic relationship, bacteriocytes were isolated from the pea aphid, Acyrthosiphon pisum, and the host transcriptome was investigated by using EST analysis and real-time quantitative RT-PCR. A number of genes were highly expressed specifically in the bacteriocyte, including (i) genes for amino acid metabolism, including those for biosynthesis of amino acids that Buchnera cannot produce, and those for utilization of amino acids that Buchnera can synthesize; (ii) genes related to transport, including genes for mitochondrial transporters and a gene encoding Rab, a G protein that regulates vesicular transport; and (iii) genes for putative lysozymes that degrade bacterial cell walls. Significant up-regulation of i clearly indicated that the bacteriocyte is involved in the exchange of amino acids between the host aphid and Buchnera, the key metabolic process in the symbiotic system. Conspicuously high expression of ii and iii shed light on previously unknown aspects of the host-Buchnera interactions in the symbiotic system.

Identification of unique transcripts from a mouse full-length, subtracted inner ear cDNA library.

A small-scale full-length library construction approach was developed to facilitate production of a mouse full-length cDNA encyclopedia representing approximately 250 enriched, normalized, and/or subtracted cDNA libraries. One library produced using this approach was a subtracted adult mouse inner ear cDNA library (sIEa). The average size of the inserts was approximately 2.5 kb, with the majority ranging from 0.5 to 7.0 kb. From this library 22,574 sequence reads were obtained from 15,958 independent clones. Sequencing and chromosomal localization established 5240 clusters, with 1302 clusters being unique and 359 representing new ESTs. Our sIEa library contributed 56.1% of the 7773 nonredundant Unigene clusters associated with the four mouse inner ear libraries in the NCBI dbEST. Based on homologous chromosomal regions between human and mouse, we identified 1018 UniGene clusters associated with the deafness locus critical regions. Of these, 59 clusters were found only in our sIEa library and represented approximately 50% of the identified critical regions.

Absolute expression values for mouse transcripts: re-annotation of the READ expression database by the use of CAGE and EST sequence tags.

The RIKEN expression array database (READ) provides comprehensive gene expression data for the mouse, which were obtained as relative values from microarray double-staining experiments with E17.5 mRNA as common reference. To assign absolute expression values for mouse transcripts within READ, we applied the E17.5 reference sample to CAGE (cap analysis of gene expression) and expressed sequence tag (EST) high-throughput tag sequencing. Newly assigned values within the READ database were validated by comparison to expression data from serial analysis of gene expression, CAGE and EST experiments. These experiments confirmed the great significance of the absolute expression values within the improved READ database. The new Absolute READ database on absolute expression data is available under.

Libraries enriched for alternatively spliced exons reveal splicing patterns in melanocytes and melanomas.

It is becoming increasingly clear that alternative splicing enables the complex development and homeostasis of higher organisms. To gain a better understanding of how splicing contributes to regulatory pathways, we have developed an alternative splicing library approach for the identification of alternatively spliced exons and their flanking regions by alternative splicing sequence enriched tags sequencing. Here, we have applied our approach to mouse melan-c melanocyte and B16-F10Y melanoma cell lines, in which 5,401 genes were found to be alternatively spliced. These genes include those encoding important regulatory factors such as cyclin D2, Ilk, MAPK12, MAPK14, RAB4, melastatin 1 and previously unidentified splicing events for 436 genes. Real-time PCR further identified cell line-specific exons for Tmc6, Abi1, Sorbs1, Ndel1 and Snx16. Thus, the ASL approach proved effective in identifying splicing events, which suggest that alternative splicing is important in melanoma development.

Cap analysis gene expression for high-throughput analysis of transcriptional starting point and identification of promoter usage.

We introduce cap analysis gene expression (CAGE), which is based on preparation and sequencing of concatamers of DNA tags deriving from the initial 20 nucleotides from 5' end mRNAs. CAGE allows high-throughout gene expression analysis and the profiling of transcriptional start points (TSP), including promoter usage analysis. By analyzing four libraries (brain, cortex, hippocampus, and cerebellum), we redefined more accurately the TSPs of 11-27% of the analyzed transcriptional units that were hit. The frequency of CAGE tags correlates well with results from other analyses, such as serial analysis of gene expression, and furthermore maps the TSPs more accurately, including in tissue-specific cases. The high-throughput nature of this technology paves the way for understanding gene networks via correlation of promoter usage and gene transcriptional factor expression.

Subtraction of cap-trapped full-length cDNA libraries to select rare transcripts.

The normalization and subtraction of highly expressed cDNAs from relatively large tissues before cloning dramatically enhanced the gene discovery by sequencing for the mouse full-length cDNA encyclopedia, but these methods have not been suitable for limited RNA materials. To normalize and subtract full-length cDNA libraries derived from limited quantities of total RNA, here we report a method to subtract plasmid libraries excised from size-unbiased amplified lambda phage cDNA libraries that avoids heavily biasing steps such as PCR and plasmid library amplification. The proportion of full-length cDNAs and the gene discovery rate are high, and library diversity can be validated by in silico randomization.

Targeting a complex transcriptome: the construction of the mouse full-length cDNA encyclopedia.

We report the construction of the mouse full-length cDNA encyclopedia,the most extensive view of a complex transcriptome,on the basis of preparing and sequencing 246 libraries. Before cloning,cDNAs were enriched in full-length by Cap-Trapper,and in most cases,aggressively subtracted/normalized. We have produced 1,442,236 successful 3'-end sequences clustered into 171,144 groups, from which 60,770 clones were fully sequenced cDNAs annotated in the FANTOM-2 annotation. We have also produced 547,149 5' end reads,which clustered into 124,258 groups. Altogether, these cDNAs were further grouped in 70,000 transcriptional units (TU),which represent the best coverage of a transcriptome so far. By monitoring the extent of normalization/subtraction, we define the tentative equivalent coverage (TEC),which was estimated to be equivalent to >12,000,000 ESTs derived from standard libraries. High coverage explains discrepancies between the very large numbers of clusters (and TUs) of this project,which also include non-protein-coding RNAs,and the lower gene number estimation of genome annotations. Altogether,5'-end clusters identify regions that are potential promoters for 8637 known genes and 5'-end clusters suggest the presence of almost 63,000 transcriptional starting points. An estimate of the frequency of polyadenylation signals suggests that at least half of the singletons in the EST set represent real mRNAs. Clones accounting for about half of the predicted TUs await further sequencing. The continued high-discovery rate suggests that the task of transcriptome discovery is not yet complete.

The Drosophila gene collection: identification of putative full-length cDNAs for 70% of D. melanogaster genes.

Collections of full-length nonredundant cDNA clones are critical reagents for functional genomics. The first step toward these resources is the generation and single-pass sequencing of cDNA libraries that contain a high proportion of full-length clones. The first release of the Drosophila Gene Collection Release 1 (DGCr1) was produced from six libraries representing various tissues, developmental stages, and the cultured S2 cell line. Nearly 80,000 random 5' expressed sequence tags (5' expressed sequence tags [ESTs]from these libraries were collapsed into a nonredundant set of 5849 cDNAs, corresponding to ~40% of the 13,474 predicted genes in Drosophila. To obtain cDNA clones representing the remaining genes, we have generated an additional 157,835 5' ESTs from two previously existing and three new libraries. One new library is derived from adult testis, a tissue we previously did not exploit for gene discovery; two new cap-trapped normalized libraries are derived from 0-22-h embryos and adult heads. Taking advantage of the annotated D. melanogaster genome sequence, we clustered the ESTs by aligning them to the genome. Clusters that overlap genes not already represented by cDNA clones in the DGCr1 were analyzed further, and putative full-length clones were selected for inclusion in the new DGC. This second release of the DGC (DGCr2) contains 5061 additional clones, extending the collection to 10,910 cDNAs representing >70% of the predicted genes in Drosophila.

A computer-based method of selecting clones for a full-length cDNA project: simultaneous collection of negligibly redundant and variant cDNAs.

We describe a computer-based method that selects representative clones for full-length sequencing in a full-length cDNA project. Our method classifies end sequences using two kinds of criteria, grouping, and clustering. Grouping places together variant cDNAs, family genes, and cDNAs with sequencing errors. Clustering separates those cDNA clones into distinct clusters. The full-length sequences of the clones selected by grouping are determined preferentially, and then the sequences selected by clustering are determined. Grouping reduced the number of rice cDNA clones for full-length sequencing to 21% and mouse cDNA clones to 25%. Rice full-length sequences selected by grouping showed a 1.07-fold redundancy. Mouse full-length sequences showed a 1.04-fold redundancy, which can be reduced by approximately 30% from the selection using our previous method. To estimate the coverage of unique genes, we used FANTOM (Functional Annotation of RIKEN Mouse cDNA Clones) clusters (). Grouping covered almost all unique genes (93% of FANTOM clusters), and clustering covered all genes. Therefore, our method is useful for the selection of appropriate representative clones for full-length sequencing, thereby greatly reducing the cost, labor, and time necessary for this process.