A petavoxel fragment of human cerebral cortex reconstructed at nanoscale resolution.

To fully understand how the human brain works, knowledge of its structure at high resolution is needed. Presented here is a computationally intensive reconstruction of the ultrastructure of a cubic millimeter of human temporal cortex that was surgically removed to gain access to an underlying epileptic focus. It contains about 57,000 cells, about 230 millimeters of blood vessels, and about 150 million synapses and comprises 1.4 petabytes. Our analysis showed that glia outnumber neurons 2:1, oligodendrocytes were the most common cell, deep layer excitatory neurons could be classified on the basis of dendritic orientation, and among thousands of weak connections to each neuron, there exist rare powerful axonal inputs of up to 50 synapses. Further studies using this resource may bring valuable insights into the mysteries of the human brain.

Randomized Controlled Trial of a Dichoptic Digital Therapeutic for Amblyopia.

PURPOSE: Digital therapeutics are a new class of interventions that are software driven and are intended to treat various conditions. We developed and evaluated a dichoptic digital therapeutic for amblyopia, a neurodevelopmental disorder for which current treatments may be limited by poor adherence and residual vision deficits. DESIGN: Randomized controlled trial. PARTICIPANTS: One hundred five children 4 to 7 years of age with amblyopia were enrolled at 21 academic and community sites in the United States. Participants were randomized 1:1 to the treatment or comparison group, stratified by site. METHODS: We conducted a phase 3 randomized controlled trial to evaluate the safety and efficacy of a dichoptic digital therapeutic for amblyopia. Participants in the treatment group used the therapeutic at home for 1 hour per day, 6 days per week and wore glasses full-time. Participants in the comparison group continued wearing glasses full-time alone. MAIN OUTCOME MEASURES: The primary efficacy outcome was change in amblyopic eye visual acuity (VA) from baseline at 12 weeks, and VA was measured by masked examiners. Safety was evaluated using the frequency and severity of study-related adverse events. Primary analyses were conducted using the intention-to-treat population. RESULTS: Between January 16, 2019, and January 15, 2020, 105 participants were enrolled; 51 were randomized to the treatment group and 54 were randomized to the comparison group. At 12 weeks, amblyopic eye VA improved by 1.8 lines (95% confidence interval [CI], 1.4-2.3 lines; n = 45) in the treatment group and by 0.8 lines (95% CI, 0.4-1.3 lines; n = 45) in the comparison group. At the planned interim analysis (adjusted α = 0.0193), the difference between groups was significant (1.0 lines; P = 0.0011; 96.14% CI, 0.33-1.63 lines) and the study was stopped early for success, according to the protocol. No serious adverse events were reported. CONCLUSIONS: Our findings support the value of the therapeutic in clinical practice as an effective treatment. Future studies should evaluate the therapeutic compared with other methods and in additional patient populations.

Digital therapeutic improves visual acuity and encourages high adherence in amblyopic children in open-label pilot study.

BACKGROUND: The effectiveness of amblyopia therapy can be limited by poor adherence. Dichoptic therapies are a new approach, but recent trials have demonstrated difficulty maintaining high adherence over extended periods of at-home treatment. We evaluated the efficacy and adherence of Luminopia One-a dichoptic treatment that applies therapeutic modifications to streaming content chosen by the patient. METHODS: This single-arm, multicenter prospective pilot study enrolled children aged 4-12 with anisometropic, strabismic, or mixed amblyopia at 10 pediatric ophthalmic and optometric practices across the United States. The therapeutic was prescribed for 1 hour/day, 6 days/week for 12 weeks of at-home use. The primary endpoint was best-corrected visual acuity (BCVA) at the 12-week follow-up visit. RESULTS: In total, 90 participants (mean age, 6.7 ± 2.0 years) were enrolled, and 73/90 participants (81%) had prior treatment beyond refractive correction. For those who completed the 12-week visit, mean amblyopic eye BCVA improved from 0.50 logMAR to 0.35 logMAR (1.5 logMAR lines; 95% CI, 1.2-1.8 lines; P < 0.0001). Mean stereoacuity improved by 0.28 log arcsec (95% CI, 0.14-0.42 log arcsec; P < 0.0001). Median adherence was 86% (interquartile range, 70%-97%). CONCLUSIONS: In our study cohort, adherence over the 12-week study period was high, and participants demonstrated clinically and statistically significant improvements in visual acuity and stereoacuity.

Improved adherence and treatment outcomes with an engaging, personalized digital therapeutic in amblyopia.

Given the prevalence of poor adherence to therapy and the biases of self-reporting across healthcare, we hypothesized that an engaging, personalized therapy may improve adherence and treatment outcomes in the home. We tested this hypothesis in the initial indication of amblyopia, a neurodevelopmental disorder for which available treatments are limited by low adherence. We designed a novel digital therapeutic that modifies patient-selected cinematic content in real-time into therapeutic visual input, while objectively monitoring adherence. The therapeutic design integrated a custom-designed headset that delivers precise visual input to each eye, computational algorithms that apply real-time therapeutic modifications to source content, a cloud-based content management system that enables treatment in the home, and a broad library of licensed content. In a proof-of-concept human study on the therapeutic, we found that amblyopic eye vision improved significantly after 12 weeks of treatment, with higher adherence than that of available treatments. These initial results support the utility of personalized therapy in amblyopia and may have broader relevance for improving treatment outcomes in additional indications.

High throughput generation of promoter reporter (GFP) transgenic lines of low expressing genes in Arabidopsis and analysis of their expression patterns.

BACKGROUND: Although the complete genome sequence and annotation of Arabidopsis were released at the end of year 2000, it is still a great challenge to understand the function of each gene in the Arabidopsis genome. One way to understand the function of genes on a genome-wide scale is expression profiling by microarrays. However, the expression level of many genes in Arabidopsis genome cannot be detected by microarray experiments. In addition, there are many more novel genes that have been discovered by experiments or predicted by new gene prediction programs. Another way to understand the function of individual genes is to investigate their in vivo expression patterns by reporter constructs in transgenic plants which can provide basic information on the patterns of gene expression. RESULTS: A high throughput pipeline was developed to generate promoter-reporter (GFP) transgenic lines for Arabidopsis genes expressed at very low levels and to examine their expression patterns in vivo. The promoter region from a total of 627 non- or low-expressed genes in Arabidopsis based on Arabidopsis annotation release 5 were amplified and cloned into a Gateway vector. A total of 353 promoter-reporter (GFP) constructs were successfully transferred into Agrobacterium (GV3101) by triparental mating and subsequently used for Arabidopsis transformation. Kanamycin-resistant transgenic lines were obtained from 266 constructs and among them positive GFP expression was detected from 150 constructs. Of these 150 constructs, multiple transgenic lines exhibiting consistent expression patterns were obtained for 112 constructs. A total 81 different regions of expression were discovered during our screening of positive transgenic plants and assigned Plant Ontology (PO) codes. CONCLUSIONS: Many of the genes tested for which expression data were lacking previously are indeed expressed in Arabidopsis during the developmental stages screened. More importantly, our study provides plant researchers with another resource of gene expression information in Arabidopsis. The results of this study are captured in a MySQL database and can be searched at http://www.jcvi.org/arabidopsis/qpcr/index.shtml. Transgenic seeds and constructs are also available for the research community.

A community resource for high-throughput quantitative RT-PCR analysis of transcription factor gene expression in Medicago truncatula.

BACKGROUND: Medicago truncatula is a model legume species that is currently the focus of an international genome sequencing effort. Although several different oligonucleotide and cDNA arrays have been produced for genome-wide transcript analysis of this species, intrinsic limitations in the sensitivity of hybridization-based technologies mean that transcripts of genes expressed at low-levels cannot be measured accurately with these tools. Amongst such genes are many encoding transcription factors (TFs), which are arguably the most important class of regulatory proteins. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) is the most sensitive method currently available for transcript quantification, and one that can be scaled up to analyze transcripts of thousands of genes in parallel. Thus, qRT-PCR is an ideal method to tackle the problem of TF transcript quantification in Medicago and other plants. RESULTS: We established a bioinformatics pipeline to identify putative TF genes in Medicago truncatula and to design gene-specific oligonucleotide primers for qRT-PCR analysis of TF transcripts. We validated the efficacy and gene-specificity of over 1000 TF primer pairs and utilized these to identify sets of organ-enhanced TF genes that may play important roles in organ development or differentiation in this species. This community resource will be developed further as more genome sequence becomes available, with the ultimate goal of producing validated, gene-specific primers for all Medicago TF genes. CONCLUSION: High-throughput qRT-PCR using a 384-well plate format enables rapid, flexible, and sensitive quantification of all predicted Medicago transcription factor mRNAs. This resource has been utilized recently by several groups in Europe, Australia, and the USA, and we expect that it will become the 'gold-standard' for TF transcript profiling in Medicago truncatula.

MPIDB: the microbial protein interaction database.

SUMMARY: The microbial protein interaction database (MPIDB) aims to collect and provide all known physical microbial interactions. Currently, 22,530 experimentally determined interactions among proteins of 191 bacterial species/strains can be browsed and downloaded. These microbial interactions have been manually curated from the literature or imported from other databases (IntAct, DIP, BIND, MINT) and are linked to 24,060 experimental evidences (PubMed ID, PSI-MI methods). In contrast to these databases, interactions in MPIDB are further supported by 8150 additional evidences based on interaction conservation, co-purification and 3D domain contacts (iPfam, 3did). AVAILABILITY: http://www.jcvi.org/mpidb/

Small cysteine-rich peptides resembling antimicrobial peptides have been under-predicted in plants.

Multicellular organisms produce small cysteine-rich antimicrobial peptides as an innate defense against pathogens. While defensins, a well-known class of such peptides, are common among eukaryotes, there are other classes restricted to the plant kingdom. These include thionins, lipid transfer proteins and snakins. In earlier work, we identified several divergent classes of small putatively secreted cysteine-rich peptides (CRPs) in legumes [Graham et al. (2004)Plant Physiol. 135, 1179-97]. Here, we built sequence motif models for each of these classes of peptides, and iteratively searched for related sequences within the comprehensive UniProt protein dataset, the Institute for Genomic Research's 33 plant gene indices, and the entire genomes of the model dicot, Arabidopsis thaliana, and the model monocot and crop species, Oryza sativa (rice). Using this search strategy, we identified approximately 13,000 plant genes encoding peptides with common features: (i) an N-terminal signal peptide, (ii) a small divergent charged or polar mature peptide with conserved cysteines, (iii) a similar intron/exon structure, (iv) spatial clustering in the genomes studied, and (v) overrepresentation in expressed sequences from reproductive structures of specific taxa. The identified genes include classes of defensins, thionins, lipid transfer proteins, and snakins, plus other protease inhibitors, pollen allergens, and uncharacterized gene families. We estimate that these classes of genes account for approximately 2-3% of the gene repertoire of each model species. Although 24% of the genes identified were not annotated in the latest Arabidopsis genome releases (TIGR5, TAIR6), we confirmed expression via RT-PCR for 59% of the sequences attempted. These findings highlight limitations in current annotation procedures for small divergent peptide classes.

Phosphorus stress in common bean: root transcript and metabolic responses.

Phosphorus (P) is an essential element for plant growth. Crop production of common bean (Phaseolus vulgaris), the most important legume for human consumption, is often limited by low P in the soil. Functional genomics were used to investigate global gene expression and metabolic responses of bean plants grown under P-deficient and P-sufficient conditions. P-deficient plants showed enhanced root to shoot ratio accompanied by reduced leaf area and net photosynthesis rates. Transcript profiling was performed through hybridization of nylon filter arrays spotted with cDNAs of 2,212 unigenes from a P deficiency root cDNA library. A total of 126 genes, representing different functional categories, showed significant differential expression in response to P: 62% of these were induced in P-deficient roots. A set of 372 bean transcription factor (TF) genes, coding for proteins with Inter-Pro domains characteristic or diagnostic for TF, were identified from The Institute of Genomic Research/Dana Farber Cancer Institute Common Bean Gene Index. Using real-time reverse transcription-polymerase chain reaction analysis, 17 TF genes were differentially expressed in P-deficient roots; four TF genes, including MYB TFs, were induced. Nonbiased metabolite profiling was used to assess the degree to which changes in gene expression in P-deficient roots affect overall metabolism. Stress-related metabolites such as polyols accumulated in P-deficient roots as well as sugars, which are known to be essential for P stress gene induction. Candidate genes have been identified that may contribute to root adaptation to P deficiency and be useful for improvement of common bean.

Experimental validation of novel genes predicted in the un-annotated regions of the Arabidopsis genome.

BACKGROUND: Several lines of evidence support the existence of novel genes and other transcribed units which have not yet been annotated in the Arabidopsis genome. Two gene prediction programs which make use of comparative genomic analysis, Twinscan and EuGene, have recently been deployed on the Arabidopsis genome. The ability of these programs to make use of sequence data from other species has allowed both Twinscan and EuGene to predict over 1000 genes that are intergenic with respect to the most recent annotation release. A high throughput RACE pipeline was utilized in an attempt to verify the structure and expression of these novel genes. RESULTS: 1,071 un-annotated loci were targeted by RACE, and full length sequence coverage was obtained for 35% of the targeted genes. We have verified the structure and expression of 378 genes that were not present within the most recent release of the Arabidopsis genome annotation. These 378 genes represent a structurally diverse set of transcripts and encode a functionally diverse set of proteins. CONCLUSION: We have investigated the accuracy of the Twinscan and EuGene gene prediction programs and found them to be reliable predictors of gene structure in Arabidopsis. Several hundred previously un-annotated genes were validated by this work. Based upon this information derived from these efforts it is likely that the Arabidopsis genome annotation continues to overlook several hundred protein coding genes.

The TIGR Plant Transcript Assemblies database.

The TIGR Plant Transcript Assemblies (TA) database (http://plantta.tigr.org) uses expressed sequences collected from the NCBI GenBank Nucleotide database for the construction of transcript assemblies. The sequences collected include expressed sequence tags (ESTs) and full-length and partial cDNAs, but exclude computationally predicted gene sequences. The TA database includes all plant species for which more than 1000 EST or cDNA sequences are publicly available. The EST and cDNA sequences are first clustered based on an all-versus-all pairwise sequence comparison, followed by the generation of consensus sequences (TAs) from individual clusters. The clustering and assembly procedures use the TGICL tool, Megablast and the CAP3 assembler. The UniProt Reference Clusters (UniRef100) protein database is used as the reference database for the functional annotation of the assemblies. The transcription orientation of each TA is determined based on the orientation of the alignment with the best protein hit. The TA sequences and annotation are available via web interfaces and FTP downloads. Assemblies can be retrieved by a text-based keyword search or a sequence-based BLAST search. The current version of the TA database is Release 2 (July 17, 2006) and includes a total of 215 plant species.

Development of Arabidopsis whole-genome microarrays and their application to the discovery of binding sites for the TGA2 transcription factor in salicylic acid-treated plants.

We have developed two long-oligonucleotide microarrays for the analysis of genome features in Arabidopsis thaliana, in particular for the high-throughput identification of transcription factor-binding sites. The first platform contains 190,000 probes representing the 2-kb regions upstream of all annotated genes at a density of seven probes per promoter. The second platform is divided into three chips, each of over 390,000 features, and represents the entire Arabidopsis genome at a density of one probe per 90 bases. Protein-DNA complexes resulting from the formaldehyde fixation of leaves of plants 2 h after exposure to 1 mm salicylic acid (SA) were immunoprecipitated using antibodies against the TGA2 transcription factor. After reversal of the cross-links and amplification, the resulting ChIP sample was hybridized to both platforms. High signal ratios of the ChIP sample versus raw chromatin for clusters of neighboring probes provided evidence for 51 putative binding sites for TGA2, including the only previously confirmed site in the promoter of PR-1 (At2g14610). Enrichment of several regions was confirmed by quantitative real-time PCR. Motif search revealed that the palindromic octamer TGACGTCA was found in 55% of the enriched regions. Interestingly, 15 of the putative binding sites for TGA2 lie outside the presumptive promoter regions. The effect of the 2-h SA treatment on gene expression was measured using Affymetrix ATH1 arrays, and SA-induced genes were found to be significantly over-represented among genes neighboring putative TGA2-binding sites.

Analysis of the cDNAs of hypothetical genes on Arabidopsis chromosome 2 reveals numerous transcript variants.

In the fully sequenced Arabidopsis (Arabidopsis thaliana) genome, many gene models are annotated as "hypothetical protein," whose gene structures are predicted solely by computer algorithms with no support from either expressed sequence matches from Arabidopsis, or nucleic acid or protein homologs from other species. In order to confirm their existence and predicted gene structures, a high-throughput method of rapid amplification of cDNA ends (RACE) was used to obtain their cDNA sequences from 11 cDNA populations. Primers from all of the 797 hypothetical genes on chromosome 2 were designed, and, through 5' and 3' RACE, clones from 506 genes were sequenced and cDNA sequences from 399 target genes were recovered. The cDNA sequences were obtained by assembling their 5' and 3' RACE polymerase chain reaction products. These sequences revealed that (1) the structures of 151 hypothetical genes were different from their predictions; (2) 116 hypothetical genes had alternatively spliced transcripts and 187 genes displayed polyadenylation sites; and (3) there were transcripts arising from both strands, from the strand opposite to that of the prediction and possible dicistronic transcripts. Promoters from five randomly chosen hypothetical genes (At2g02540, At2g31270, At2g33640, At2g35550, and At2g36340) were cloned into report constructs, and their expressions are tissue or development stage specific. Our results indicate at least 50% of hypothetical genes on chromosome 2 are expressed in the cDNA populations with about 38% of the gene structures differing from their predictions. Thus, by using this targeted approach, high-throughput RACE, we revealed numerous transcripts including many uncharacterized variants from these hypothetical genes.

Whole genome shotgun sequencing of Brassica oleracea and its application to gene discovery and annotation in Arabidopsis.

Through comparative studies of the model organism Arabidopsis thaliana and its close relative Brassica oleracea, we have identified conserved regions that represent potentially functional sequences overlooked by previous Arabidopsis genome annotation methods. A total of 454,274 whole genome shotgun sequences covering 283 Mb (0.44 x) of the estimated 650 Mb Brassica genome were searched against the Arabidopsis genome, and conserved Arabidopsis genome sequences (CAGSs) were identified. Of these 229,735 conserved regions, 167,357 fell within or intersected existing gene models, while 60,378 were located in previously unannotated regions. After removal of sequences matching known proteins, CAGSs that were close to one another were chained together as potentially comprising portions of the same functional unit. This resulted in 27,347 chains of which 15,686 were sufficiently distant from existing gene annotations to be considered a novel conserved unit. Of 192 conserved regions examined, 58 were found to be expressed in our cDNA populations. Rapid amplification of cDNA ends (RACE) was used to obtain potentially full-length transcripts from these 58 regions. The resulting sequences led to the creation of 21 gene models at 17 new Arabidopsis loci and the addition of splice variants or updates to another 19 gene structures. In addition, CAGSs overlapping already annotated genes in Arabidopsis can provide guidance for manual improvement of existing gene models. Published genome-wide expression data based on whole genome tiling arrays and massively parallel signature sequencing were overlaid on the Brassica-Arabidopsis conserved sequences, and 1399 regions of intersection were identified. Collectively our results and these data sets suggest that several thousand new Arabidopsis genes remain to be identified and annotated.

Empirical analysis of transcriptional activity in the Arabidopsis genome.

Functional analysis of a genome requires accurate gene structure information and a complete gene inventory. A dual experimental strategy was used to verify and correct the initial genome sequence annotation of the reference plant Arabidopsis. Sequencing full-length cDNAs and hybridizations using RNA populations from various tissues to a set of high-density oligonucleotide arrays spanning the entire genome allowed the accurate annotation of thousands of gene structures. We identified 5817 novel transcription units, including a substantial amount of antisense gene transcription, and 40 genes within the genetically defined centromeres. This approach resulted in completion of approximately 30% of the Arabidopsis ORFeome as a resource for global functional experimentation of the plant proteome.