Large-scale high-density brain-wide neural recording in nonhuman primates.

High-density, integrated silicon electrodes have begun to transform systems neuroscience, by enabling large-scale neural population recordings with single cell resolution. Existing technologies, however, have provided limited functionality in nonhuman primate species such as macaques, which offer close models of human cognition and behavior. Here, we report the design, fabrication, and performance of Neuropixels 1.0-NHP, a high channel count linear electrode array designed to enable large-scale simultaneous recording in superficial and deep structures within the macaque or other large animal brain. These devices were fabricated in two versions: 4416 electrodes along a 45 mm shank, and 2496 along a 25 mm shank. For both versions, users can programmatically select 384 channels, enabling simultaneous multi-area recording with a single probe. We demonstrate recording from over 3000 single neurons within a session, and simultaneous recordings from over 1000 neurons using multiple probes. This technology represents a significant increase in recording access and scalability relative to existing technologies, and enables new classes of experiments involving fine-grained electrophysiological characterization of brain areas, functional connectivity between cells, and simultaneous brain-wide recording at scale.

Comparative transcriptome analysis during tuberous stem formation in Kohlrabi (B. oleracea var. gongylodes) at early growth periods (seedling stages).

Tuberous stem of kohlrabi is an important agronomic trait, however, the molecular basis of tuberization is poorly understood. To elucidate the tuberization mechanism, we conducted a comparative transcriptomic analysis between kohlrabi and broccoli at 10 and 20 days after germination (DAG) as tuberous stem initiated between these time points. A total of 5580 and 2866 differentially expressed transcripts (DETs) were identified between genotypes (kohlrabi vs. broccoli) and growth stages (10 DAG vs. 20 DAG), respectively, and most of the DETs were down-regulated in kohlrabi. Gene ontology (GO) and KEGG pathway enrichment analyses showed that the DETs between genotypes are involved in cell wall loosening and expansion, cell cycle and division, carbohydrate metabolism, hormone transport, hormone signal transduction and in several transcription factors. The DETs identified in those categories may directly/indirectly relate to the initiation and development of tuberous stem in kohlrabi. In addition, the expression pattern of the hormone synthesis related DETs coincided with the endogenous IAA, IAAsp, GA, ABA, and tZ profiles in kohlrabi and broccoli seedlings, that were revealed in our phytohormone analysis. This is the first report on comparative transcriptome analysis for tuberous stem formation in kohlrabi at early growth periods. The resulting data could provide significant insights into the molecular mechanism underlying tuberous stem development in kohlrabi as well as in other tuberous organ forming crops.

High-performance pipeline for MutMap and QTL-seq.

Summary: Bulked segregant analysis implemented in MutMap and QTL-seq is a powerful and efficient method to identify loci contributing to important phenotypic traits. However, the previous pipelines were not user-friendly to install and run. Here, we describe new pipelines for MutMap and QTL-seq. These updated pipelines are approximately 5-8 times faster than the previous pipeline, are easier for novice users to use, and can be easily installed through bioconda with all dependencies. Availability: The new pipelines of MutMap and QTL-seq are written in Python and can be installed via bioconda. The source code and manuals are available online (MutMap: https://github.com/YuSugihara/MutMap, QTL-seq: https://github.com/YuSugihara/QTL-seq).

Cortical preparatory activity indexes learned motor memories.

The brain's remarkable ability to learn and execute various motor behaviours harnesses the capacity of neural populations to generate a variety of activity patterns. Here we explore systematic changes in preparatory activity in motor cortex that accompany motor learning. We trained rhesus monkeys to learn an arm-reaching task1 in a curl force field that elicited new muscle forces for some, but not all, movement directions2,3. We found that in a neural subspace predictive of hand forces, changes in preparatory activity tracked the learned behavioural modifications and reassociated4 existing activity patterns with updated movements. Along a neural population dimension orthogonal to the force-predictive subspace, we discovered that preparatory activity shifted uniformly for all movement directions, including those unaltered by learning. During a washout period when the curl field was removed, preparatory activity gradually reverted in the force-predictive subspace, but the uniform shift persisted. These persistent preparatory activity patterns may retain a motor memory of the learned field5,6 and support accelerated relearning of the same curl field. When a set of distinct curl fields was learned in sequence, we observed a corresponding set of field-specific uniform shifts which separated the associated motor memories in the neural state space7-9. The precise geometry of these uniform shifts in preparatory activity could serve to index motor memories, facilitating the acquisition, retention and retrieval of a broad motor repertoire.

Giant obscurin regulates migration and metastasis via RhoA-dependent cytoskeletal remodeling in pancreatic cancer.

Obscurins, encoded by the OBSCN gene, are giant cytoskeletal proteins with structural and regulatory roles. Large scale omics analyses reveal that OBSCN is highly mutated across different types of cancer, exhibiting a 5-8% mutation frequency in pancreatic cancer. Yet, the functional role of OBSCN in pancreatic cancer progression and metastasis has to be delineated. We herein show that giant obscurins are highly expressed in normal pancreatic tissues, but their levels are markedly reduced in pancreatic ductal adenocarcinomas. Silencing of giant obscurins in non-tumorigenic Human Pancreatic Ductal Epithelial (HPDE) cells and obscurin-expressing Panc5.04 pancreatic cancer cells induces an elongated, spindle-like morphology and faster cell migration via cytoskeletal remodeling. Specifically, depletion of giant obscurins downregulates RhoA activity, which in turn results in reduced focal adhesion density, increased microtubule growth rate and faster actin dynamics. Although OBSCN knockdown is not sufficient to induce de novo tumorigenesis, it potentiates tumor growth in a subcutaneous implantation model and exacerbates metastasis in a hemispleen murine model of pancreatic cancer metastasis, thereby shortening survival. Collectively, these findings reveal a critical role of giant obscurins as tumor suppressors in normal pancreatic epithelium whose loss of function induces RhoA-dependent cytoskeletal remodeling, and promotes cell migration, tumor growth and metastasis.

Dendritic calcium signals in rhesus macaque motor cortex drive an optical brain-computer interface.

Calcium imaging is a powerful tool for recording from large populations of neurons in vivo. Imaging in rhesus macaque motor cortex can enable the discovery of fundamental principles of motor cortical function and can inform the design of next generation brain-computer interfaces (BCIs). Surface two-photon imaging, however, cannot presently access somatic calcium signals of neurons from all layers of macaque motor cortex due to photon scattering. Here, we demonstrate an implant and imaging system capable of chronic, motion-stabilized two-photon imaging of neuronal calcium signals from macaques engaged in a motor task. By imaging apical dendrites, we achieved optical access to large populations of deep and superficial cortical neurons across dorsal premotor (PMd) and gyral primary motor (M1) cortices. Dendritic signals from individual neurons displayed tuning for different directions of arm movement. Combining several technical advances, we developed an optical BCI (oBCI) driven by these dendritic signalswhich successfully decoded movement direction online. By fusing two-photon functional imaging with CLARITY volumetric imaging, we verified that many imaged dendrites which contributed to oBCI decoding originated from layer 5 output neurons, including a putative Betz cell. This approach establishes new opportunities for studying motor control and designing BCIs via two photon imaging.

Genome-wide analysis of long noncoding RNAs, 24-nt siRNAs, DNA methylation and H3K27me3 marks in Brassica rapa.

Long noncoding RNAs (lncRNAs) are RNA fragments that generally do not code for a protein but are involved in epigenetic gene regulation. In this study, lncRNAs of Brassica rapa were classified into long intergenic noncoding RNAs, natural antisense RNAs, and intronic noncoding RNAs and their expression analyzed in relation to genome-wide 24-nt small interfering RNAs (siRNAs), DNA methylation, and histone H3 lysine 27 trimethylation marks (H3K27me3). More than 65% of the lncRNAs analyzed consisted of one exon, and more than 55% overlapped with inverted repeat regions (IRRs). Overlap of lncRNAs with IRRs or genomic regions encoding for 24-nt siRNAs resulted in increased DNA methylation levels when both were present. LncRNA did not overlap greatly with H3K27me3 marks, but the expression level of intronic noncoding RNAs that did coincide with H3K27me3 marks was higher than without H3K27me3 marks. The Brassica genus comprises important vegetables and oil seed crops grown across the world. B. rapa is a diploid (AA genome) thought to be one of the ancestral species of both B. juncea (AABB genome) and B. napus (AACC) through genome merging (allotetrapolyploidization). Complex genome restructuring and epigenetic alterations are thought to be involved in these allotetrapolyploidization events. Comparison of lncRNAs between B. rapa and B. nigra, B. oleracea, B. juncea, and B. napus showed the highest conservation with B. oleracea. This study presents a comprehensive analysis of the epigenome structure of B. rapa at multi-epigenetic levels (siRNAs, DNA methylation, H3K27me3, and lncRNAs) and identified a suite of candidate lncRNAs that may be epigenetically regulated in the Brassica genus.

Genetics of Clubroot and Fusarium Wilt Disease Resistance in Brassica Vegetables: The Application of Marker Assisted Breeding for Disease Resistance.

The genus Brassica contains important vegetable crops, which serve as a source of oil seed, condiments, and forages. However, their production is hampered by various diseases such as clubroot and Fusarium wilt, especially in Brassica vegetables. Soil-borne diseases are difficult to manage by traditional methods. Host resistance is an important tool for minimizing disease and many types of resistance (R) genes have been identified. More than 20 major clubroot (CR) disease-related loci have been identified in Brassica vegetables and several CR-resistant genes have been isolated by map-based cloning. Fusarium wilt resistant genes in Brassica vegetables have also been isolated. These isolated R genes encode the toll-interleukin-1 receptor/nucleotide-binding site/leucine-rice-repeat (TIR-NBS-LRR) protein. DNA markers that are linked with disease resistance allele have been successfully applied to improve disease resistance through marker-assisted selection (MAS). In this review, we focused on the recent status of identifying clubroot and Fusarium wilt R genes and the feasibility of using MAS for developing disease resistance cultivars in Brassica vegetables.

Causal Role of Motor Preparation during Error-Driven Learning.

Current theories suggest that an error-driven learning process updates trial-by-trial to facilitate motor adaptation. How this process interacts with motor cortical preparatory activity-which current models suggest plays a critical role in movement initiation-remains unknown. Here, we evaluated the role of motor preparation during visuomotor adaptation. We found that preparation time was inversely correlated to variance of errors on current trials and mean error on subsequent trials. We also found causal evidence that intracortical microstimulation during motor preparation was sufficient to disrupt learning. Surprisingly, stimulation did not affect current trials, but instead disrupted the update computation of a learning process, thereby affecting subsequent trials. This is consistent with a Bayesian estimation framework where the motor system reduces its learning rate by virtue of lowering error sensitivity when faced with uncertainty. This interaction between motor preparation and the error-driven learning system may facilitate new probes into mechanisms underlying trial-by-trial adaptation.

Publisher Correction: Development of an optogenetic toolkit for neural circuit dissection in squirrel monkeys.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The histone modification H3 lysine 27 tri-methylation has conserved gene regulatory roles in the triplicated genome of Brassica rapa L.

Brassica rapa L. is an important vegetable and oilseed crop. We investigated the distribution of the histone mark tri-methylation of H3K27 (H3K27me3) in B. rapa and its role in the control of gene expression at two stages of development (2-day cotyledons and 14-day leaves) and among paralogs in the triplicated genome. H3K27me3 has a similar distribution in two inbred lines, while there was variation of H3K27me3 sites between tissues. Sites that are specific to 2-day cotyledons have increased transcriptional activity, and low levels of H3K27me3 in the gene body region. In 14-day leaves, levels of H3K27me3 were associated with decreased gene expression. In the triplicated genome, H3K27me3 is associated with paralogs that have tissue-specific expression. Even though B. rapa and Arabidopsis thaliana are not closely related within the Brassicaceae, there is conservation of H3K27me3-marked sites in the two species. Both B. rapa and A. thaliana require vernalization for floral initiation with FLC being the major controlling locus. In all four BrFLC paralogs, low-temperature treatment increases H3K27me3 at the proximal nucleation site reducing BrFLC expression. Following return to normal temperature growth conditions, H3K27me3 spreads along all four BrFLC paralogs providing stable repression of the gene.

The role of FRIGIDA and FLOWERING LOCUS C genes in flowering time of Brassica rapa leafy vegetables.

There is a wide variation of flowering time among lines of Brassica rapa L. Most B. rapa leafy (Chinese cabbage etc.) or root (turnip) vegetables require prolonged cold exposure for flowering, known as vernalization. Premature bolting caused by low temperature leads to a reduction in the yield/quality of these B. rapa vegetables. Therefore, high bolting resistance is an important breeding trait, and understanding the molecular mechanism of vernalization is necessary to achieve this goal. In this study, we demonstrated that BrFRIb functions as an activator of BrFLC in B. rapa. We showed a positive correlation between the steady state expression levels of the sum of the BrFLC paralogs and the days to flowering after four weeks of cold treatment, suggesting that this is an indicator of the vernalization requirement. We indicate that BrFLCs are repressed by the accumulation of H3K27me3 and that the spreading of H3K27me3 promotes stable FLC repression. However, there was no clear relationship between the level of H3K27me3 in the BrFLC and the vernalization requirement. We also showed that if there was a high vernalization requirement, the rate of repression of BrFLC1 expression following prolonged cold treatments was lower.

Allele specific DNA marker for fusarium resistance gene FocBo1 in Brassica oleracea.

The fusarium yellows resistance (YR) gene FocBo1 was previously identified and the DNA markers were developed to assist the breeding of YR cultivars in Brassica oleracea. However, the further analysis revealed discrepancies between the phenotypes and the genotypes predicted by those DNA markers in cabbage commercial cultivars. Since this discrepancy seemed to be due to unknown susceptible alleles of focbo1, we sequenced the gene in 19 accessions to determine the sequence variations between alleles and found that there were two resistant FocBo1 alleles and six susceptible alleles in the investigated population. The newly designed PCR markers detected three mutations in the susceptible alleles that generate premature termination codons. These were shown to accurately distinguish resistant and susceptible alleles in more than 200 accessions of B. oleracea inbred lines and cultivars. The study revealed that the locus is represented by 37.2% resistant and 62.8% susceptible alleles within seventy-eight commercial cultivars. Structural analysis of the gene revealed that a part of the allelic variation comes from intragenic recombination between alleles. Our results enable a more precise prediction of the phenotype by marker assisted selection, promoting the production of YR cultivars in B. oleracea.

Accurate Estimation of Neural Population Dynamics without Spike Sorting.

A central goal of systems neuroscience is to relate an organism's neural activity to behavior. Neural population analyses often reduce the data dimensionality to focus on relevant activity patterns. A major hurdle to data analysis is spike sorting, and this problem is growing as the number of recorded neurons increases. Here, we investigate whether spike sorting is necessary to estimate neural population dynamics. The theory of random projections suggests that we can accurately estimate the geometry of low-dimensional manifolds from a small number of linear projections of the data. We recorded data using Neuropixels probes in motor cortex of nonhuman primates and reanalyzed data from three previous studies and found that neural dynamics and scientific conclusions are quite similar using multiunit threshold crossings rather than sorted neurons. This finding unlocks existing data for new analyses and informs the design and use of new electrode arrays for laboratory and clinical use.

Long noncoding RNAs in Brassica rapa L. following vernalization.

Brassica rapa L. is an important agricultural crop that requires a period of prolonged cold for flowering. This process is known as vernalization. Studies have shown that long noncoding RNAs (lncRNAs) play important roles in abiotic stress responses and several cold-responsive noncoding RNAs have been suggested to be involved in vernalization. We examined the transcriptome of the Chinese cabbage inbred line (B. rapa L. var. pekinensis) RJKB-T24, and identified 1,444 long intergenic noncoding RNAs (lincRNAs), 551 natural antisense transcripts (NATs), and 93 intronic noncoding RNAs (incRNAs); 549 of the 2,088 lncRNAs significantly altered their expression in response to four weeks of cold treatment. Most differentially expressed lncRNAs did not lead to a change of expression levels in mRNAs covering or near lncRNAs, suggesting that the transcriptional responses to four weeks of cold treatment in lncRNA and mRNA are independent. However, some differentially expressed mRNAs had NATs with expression altered in the same direction. These genes were categorized as having an abiotic stress response, suggesting that the paired-expression may play a role in the transcriptional response to vernalization or cold treatment. We also identified short-term cold treatment induced NATs in BrFLC and BrMAF genes, which are involved in vernalization. The lncRNAs we identified differed from those reported in Arabidopsis thaliana, suggesting the role of lncRNAs in vernalization differ between these two species.

A microfluidic assay for the quantification of the metastatic propensity of breast cancer specimens.

The challenge of predicting which patients with breast cancer will develop metastases leads to the overtreatment of patients with benign disease and to the inadequate treatment of aggressive cancers. Here, we report the development and testing of a microfluidic assay that quantifies the abundance and proliferative index of migratory cells in breast cancer specimens, for the assessment of their metastatic propensity and for the rapid screening of potential antimetastatic therapeutics. On the basis of the key roles of cell motility and proliferation in cancer metastasis, the device accurately predicts the metastatic potential of breast cancer cell lines and of patient-derived xenografts. Compared with unsorted cancer cells, highly motile cells isolated by the device exhibited similar tumourigenic potential but markedly increased metastatic propensity in vivo. RNA sequencing of the highly motile cells revealed an enrichment of motility-related and survival-related genes. The approach might be developed into a companion assay for the prediction of metastasis in patients and for the selection of effective therapeutic regimens.

A Direct Podocalyxin-Dynamin-2 Interaction Regulates Cytoskeletal Dynamics to Promote Migration and Metastasis in Pancreatic Cancer Cells.

The sialoglycoprotein podocalyxin is absent in normal pancreas but is overexpressed in pancreatic cancer and is associated with poor clinical outcome. Here, we investigate the role of podocalyxin in migration and metastasis of pancreatic adenocarcinomas using SW1990 and Pa03c as cell models. Although ezrin is regarded as a cytoplasmic binding partner of podocalyxin that regulates actin polymerization via Rac1 or RhoA, we did not detect podocalyxin-ezrin association in pancreatic cancer cells. Moreover, depletion of podocalyxin did not alter actin dynamics or modulate Rac1 and RhoA activities in pancreatic cancer cells. Using mass spectrometry, bioinformatics analysis, coimmunoprecipitation, and pull-down assays, we discovered a novel, direct binding interaction between the cytoplasmic tail of podocalyxin and the large GTPase dynamin-2 at its GTPase, middle, and pleckstrin homology domains. This podocalyxin-dynamin-2 interaction regulated microtubule growth rate, which in turn modulated focal adhesion dynamics and ultimately promoted efficient pancreatic cancer cell migration via microtubule- and Src-dependent pathways. Depletion of podocalyxin in a hemispleen mouse model of pancreatic cancer diminished liver metastasis without altering primary tumor size. Collectively, these findings reveal a novel mechanism by which podocalyxin facilitates pancreatic cancer cell migration and metastasis. SIGNIFICANCE: These findings reveal that a novel interaction between podocalyxin and dynamin-2 promotes migration and metastasis of pancreatic cancer cells by regulating microtubule and focal adhesion dynamics.

Inferring single-trial neural population dynamics using sequential auto-encoders.

Neuroscience is experiencing a revolution in which simultaneous recording of thousands of neurons is revealing population dynamics that are not apparent from single-neuron responses. This structure is typically extracted from data averaged across many trials, but deeper understanding requires studying phenomena detected in single trials, which is challenging due to incomplete sampling of the neural population, trial-to-trial variability, and fluctuations in action potential timing. We introduce latent factor analysis via dynamical systems, a deep learning method to infer latent dynamics from single-trial neural spiking data. When applied to a variety of macaque and human motor cortical datasets, latent factor analysis via dynamical systems accurately predicts observed behavioral variables, extracts precise firing rate estimates of neural dynamics on single trials, infers perturbations to those dynamics that correlate with behavioral choices, and combines data from non-overlapping recording sessions spanning months to improve inference of underlying dynamics.

The production and characterization of a BoFLC2 introgressed Brassica rapa by repeated backcrossing to an F1.

Flowering time is an important agronomic trait for Brassica rapa crops, and previous breeding work in Brassica has successfully transmitted other important agronomic traits from donor species. However, there has been no previous attempts to produce hybrids replacing the original Brassica FLC alleles with alien FLC alleles. In this paper, we introduce the creation of a chromosome substitution line (CSSL) containing a homozygous introgression of Flowering Locus C from Brassica oleracea (BoFLC2) into a B. rapa genomic background, and characterize the CSSL line with respect to the parental cultivars. The preferential transmission of alien chromosome inheritance and the pattern of transmission observed during the production of the CSSLs are also discussed.