Supplementary material : large scale read classification for next generation sequencing

Next Generation Sequencing (NGS) has revolutionised molecular biology, resulting in an explosion of data sets and an increasing role in clinical practice. Such applications necessarily require rapid identification of the organism as a prelude to annotation and further analysis. NGS data consist of a substantial number of short sequence reads, given context through downstream assembly and annotation, a process requiring reads consistent with the assumed species or species group. Highly accurate results have been obtained for restricted sets using SVM classifiers, but such methods are difficult to parallelise and success depends on careful attention to feature selection. This work examines the problem at very large scale, using a mix of synthetic and real data with a view to determining the overall structure of the problem and the effectiveness of parallel ensembles of simpler classifiers (principally random forests) in addressing the challenges of large scale genomics.

Analysis of the NF-κB p50 dimer interface by diversity screening

An in vivo screen has been devised for NF-κB p50 activity in Escherichia coli exploiting the ability of the mammalian transcription factor to emulate a prokaryotic repressor. Active intracellular p50 was shown to repress the expression of a green fluorescent protein reporter gene allowing for visual screening of colonies expressing active p50 on agar plates. A library of mutants was constructed in which the residues Y267, L269, A308 and V310 of the dimer interface were simultaneously randomised and twenty-five novel functional interfaces were selected which repressed the reporter gene to similar levels as the wild-type protein. The leucine-269 alanine-308 core was repeatedly, but not exclusively, selected from the library whilst a diversity of predominantly non-polar residues were selected at positions 267 and 310. These results indicate that L269 and A308 may form a hot spot of interaction and allow an insight into the processes of dimer selectivity and evolution within this family of transcription factors.

An ancient duplication of apple MYB transcription factors is responsible for novel red fruit-flesh phenotyp

Anthocyanin accumulation is coordinated in plants by a number of conserved transcription factors. In apple (Malus × domestica), an R2R3 MYB transcription factor has been shown to control fruit flesh and foliage anthocyanin pigmentation (MYB10) and fruit skin color (MYB1). However, the pattern of expression and allelic variation at these loci does not explain all anthocyanin-related apple phenotypes. One such example is an open-pollinated seedling of cv Sangrado that has green foliage and develops red flesh in the fruit cortex late in maturity. We used methods that combine plant breeding, molecular biology, and genomics to identify duplicated MYB transcription factors that could control this phenotype. We then demonstrated that the red-flesh cortex phenotype is associated with enhanced expression of MYB110a, a paralog of MYB10. Functional characterization of MYB110a showed that it was able to up-regulate anthocyanin biosynthesis in tobacco (Nicotiana tabacum). The chromosomal location of MYB110a is consistent with a whole-genome duplication event that occurred during the evolution of apple within the Maloideae family. Both MYB10 and MYB110a have conserved function in some cultivars, but they differ in their expression pattern and response to fruit maturity.

pGreen : a versatile and flexible binary Ti vector for Agrobacterium-mediated plant transformation

Binary Ti vectors are the plasmid vectors of choice in Agrobacterium-mediated plant transformation protocols. The pGreen series of binary Ti vectors are configured for ease-of-use and to meet the demands of a wide range of transformation procedures for many plant species. This plasmid system allows any arrangement of selectable marker and reporter gene at the right and left T-DNA borders without compromising the choice of restriction sites for cloning, since the pGreen cloning sites are based on the well-known pBluescript general vector plasmids. Its size and copy number in Escherichia coli offers increased efficiencies in routine in vitro recombination procedures. pGreen can replicate in Agrobacterium only if another plasmid, pSoup, is co-resident in the same strain. pSoup provides replication functions in trans for pGreen. The removal of RepA and Mob functions has enabled the size of pGreen to be kept to a minimum. Versions of pGreen have been used to transform several plant species with the same efficiencies as other binary Ti vectors. Information on the pGreen plasmid system is supplemented by an Internet site (http://www.pgreen.ac.uk) through which comprehensive information, protocols, order forms and lists of different pGreen marker gene permutations can be found.

Repeated sequences as genetic markers in pooled tissue samples

We show, using the PDR1 element of pea, that dispersed repeated sequences of moderate copy number can be used simply and efficiently to generate markers linked to a trait of interest. Inspection of hybridization patterns of repeated sequences to DNA mixtures of pooled genotypes is a sensitive way of detecting such markers. The large number of bands in tracks of digests of these mixtures allows the simultaneous sampling of loci at many places in the genome, and the many unlinked loci serve as internal controls. It is also shown that intensity ratios calculated from these band differences can be used to give a rough estimate of linkage distance.

Coincident sequence-specific RNA degradation of linked transgenes in the plant genome

The expression of transgenes in plant genomes can be inhibited by either transcriptional gene silencing or posttranscriptional gene silencing (PTGS). Overexpression of the chalcone synthase-A (CHS-A) transgene triggers PTGS of CHS-A and thus results in loss of flower pigmentation in petunia. We previously demonstrated that epigenetic inactivation of CHS-A transgene transcription leads to a reversion of the PTGS phenotype. Although neomycin phosphotransferase II (nptII), a marker gene co-introduced into the genome with the CHS-A transgene, is not normally silenced in petunia, even when CHS-A is silenced, here we found that nptII was silenced in a petunia line in which CHS-A PTGS was induced, but not in the revertant plants that had no PTGS of CHS-A. Transcriptional activity, accumulation of short interfering RNAs, and restoration of mRNA level after infection with viruses that had suppressor proteins of gene silencing indicated that the mechanism for nptII silencing was posttranscriptional. Read-through transcripts of the CHS-A gene toward the nptII gene were detected. Deep-sequencing analysis revealed a striking difference between the predominant size class of small RNAs produced from the read-through transcripts (22 nt) and that from the CHS-A RNAs (21 nt). These results implicate the involvement of read-through transcription and distinct phases of RNA degradation in the coincident PTGS of linked transgenes and provide new insights into the destabilization of transgene expression.

A copia-like element in Pisum demonstrates the uses of dispersed repeated sequences in genetic analysis

A DNA sequence between two legumin genes in Pisum is a member of the copia-like class of retrotransposons and represents one member of a polymorphic and heterogeneous dispersed repeated sequence family in Pisum. This sequence can be exploited in genetic studies either by RFLP analysis where several markers can be scored together, or the segregation of individual elements can be followed after PCR amplification of specific members.

Genetic aspects of the organization of legumin genes in pea

We have compared physical and genetic maps of the region around the legJ gene in pea. In this vicinity there are four B-type legumin genes, arranged as two close pairs. The detection of a recombination event within this gene cluster allows the orientation of this group of genes within the surrounding linkage group to be determined. The relationship between physical and genetic distances in this region is discussed, as are the implications of this for relating physical and genetic maps elsewhere in the pea genome.

Kiwifruit floral gene APETALA2 is alternatively spliced and accumulates in aberrant indeterminate flowers in the absence of miR172

In Arabidopsis, the identity of perianth and reproductive organs are specified by antagonistic action of two floral homeotic genes, APETALA2 (AP2) and AGAMOUS (AG). AP2 is also negatively regulated by an evolutionary conserved interaction with a microRNA, miR172, and has additional roles in general plant development. A kiwifruit gene with high levels of homology to AP2 and AP2-like genes from other plant species was identified. The transcript was abundant in the kiwifruit flower, particularly petal, suggesting a role in floral organ identity. Splice variants were identified, all containing both AP2 domains, including a variant that potentially produces a shorter transcript without the miRNA172 targeting site. Increased AP2 transcript accumulation was detected in the aberrant flowers of the mutant 'Pukekohe dwarf' with multiple perianth whorls and extended petaloid features. In contrast to normal kiwifruit flowers, the aberrant flowers failed to accumulate miR172 in the developing whorls, although accumulation was detected at the base of the flower. An additional role during dormancy in kiwifruit was proposed based on AP2 transcript accumulation in axillary buds before and after budbreak.

Mini-scale method for nuclear run-on transcription assay in plants

We present a mini-scale method for nuclear run-on transcription assay. In our method, all the centrifuge steps can be carried out by using micro-tubes for short time (5 min each) throughout the process, including isolation of transcriptionally active nuclei and purification of labeled RNA after synthesis of RNA in isolated nuclei. The assay can be performed using a small amount of plant tissue, which enables analysis of developmental changes in transcriptional status of given genes in a single individual plant. Successful results were obtained using the tissues of flower and leaf of petunia and embryo of pea, suggesting that the method is potentially applicable to a variety of plant tissues.

Quantitative stem-loop RT-PCR for detection of microRNAs

Plant microRNAs (miRNAs) are a class of endogenous small RNAs that are essential for plant development and survival. They arise from larger precursor RNAs with a characteristic hairpin structure and regulate gene activity by targeting mRNA transcripts for cleavage or translational repression. Efficient and reliable detection and quantification of miRNA expression has become an essential step in understanding their specific roles. The expression levels of miRNAs can vary dramatically between samples and they often escape detection by conventional technologies such as cloning, northern hybridization and microarray analysis. The stem-loop RT-PCR method described here is designed to detect and quantify mature miRNAs in a fast, specific, accurate and reliable manner. First, a miRNA-specific stem-loop RT primer is hybridized to the miRNA and then reverse transcribed. Next, the RT product is amplified and monitored in real time using a miRNA-specific forward primer and the universal reverse primer. This method enables miRNA expression profiling from as little as 10 pg of total RNA and is suitable for high-throughput miRNA expression analysis.

qRT-PCR of Small RNAs

Plant small RNAs are a class of 19- to 25-nucleotide (nt) RNA molecules that are essential for genome stability, development and differentiation, disease, cellular communication, signaling, and adaptive responses to biotic and abiotic stress. Small RNAs comprise two major RNA classes, short interfering RNAs (siRNAs) and microRNAs (miRNAs). Efficient and reliable detection and quantification of small RNA expression has become an essential step in understanding their roles in specific cells and tissues. Here we provide protocols for the detection of miRNAs by stem-loop RT-PCR. This method enables fast and reliable miRNA expression profiling from as little as 20 pg of total RNA extracted from plant tissue and is suitable for high-throughput miRNA expression analysis. In addition, this method can be used to detect other classes of small RNAs, provided the sequence is known and their GC contents are similar to those specific for miRNAs.

Sphingolipid distribution changes with age in the human lens

The formation of an internal barrier to the diffusion of small molecules in the lens during middle age is hypothesized to be a key event in the development of age-related nuclear (ARN) cataract. Changes in membrane lipids with age may be responsible. In this study, we investigated the effect of age on the distribution of sphingomyelins, the most abundant lens phospholipids. Human lens sections were initially analyzed by MALDI mass spectrometry imaging. A distinct annular distribution of the dihydrosphingomyelin, DHSM (d18:0/16:0), in the barrier region was observed in 64- and 70-year-old lenses but not in a 23-year-old lens. An increase in the dihydroceramide, DHCer (d18:0/16:0), in the lens nucleus was also observed in the older lenses. These findings were supported by ESI mass spectrometry analysis of lipid extracts from lenses dissected into outer, barrier, and nuclear regions. A subsequent analysis of 18 lenses ages 20-72 years revealed that sphingomyelin levels increased with age in the barrier region until reaching a plateau at approximately 40 years of age. Such changes in lipid composition will have a significant impact on the physical properties of the fiber cell membranes and may be associated with the formation of a barrier.-Deeley, J. M., J. A. Hankin, M. G. Friedrich, R. C. Murphy, R. J. W. Truscott, T. W. Mitchell, and S. J. Blanksby. Sphingolipid distribution changes with age in the human lens. J. Lipid Res. 2010. 51: 2753-2760.

Factors regulating basement membrane invasion by tumor cells

Basement membranes serve as significant barriers to the passage of tumor cells but ones which metastatic cells can pass. This involves the production of a cascade of proteases leading to the activation of a specific collagenase that degrades the unique collagen network in basement membrane. Breast cancer cells, when estrogen dependent, show a requirement for estrogen for invasive activity. However, when these cells progress to an estrogen independent state and increased malignancy, they express an invasive phenotype constitutively. Studies with various anti-estrogens suggest that these responses are mediated via the estrogen receptor. Anti-estrogens lacking agonist activity suppress invasiveness as well as growth of the breast cancer cells.

Cell adhesion molecule uvomorulin expression in human breast cancer cell lines : relationship to morphology and invasive capacities

Loss of cell-cell adhesion in carcinoma cells may be an important step in the acquisition of an invasive, metastatic phenotype. We have examined the expression of the epithelial-specific cell adhesion molecule uvomorulin (E-cadherin, cell-CAM 120/80, L-CAM) in human breast cancer cell lines. We find that fibroblastoid, highly invasive, vimentin-expressing breast cancer cell lines do not express uvomorulin. Of the more epithelial-appearing, less invasive, keratin-expressing breast cancer cell lines, some express uvomorulin, and some do not. We examined the morphologies of the cell lines in the reconstituted basement membrane matrix Matrigel and measured the ability of the cells to traverse a Matrigel-coated filter as in vitro models for detachment of carcinoma cells from neighboring cells and invasion through basement membrane into surrounding tissue. Colonies of uvomorulin-positive cells have a characteristic fused appearance in Matrigel, whereas uvomorulin-negative cells appear detached. Cells which are uvomorulin negative and vimentin positive have a stellate morphology in Matrigel. We show that uvomorulin is responsible for the fused colony morphology in Matrigel since treatment of uvomorulin-positive MCF-7 cells with an antibody to uvomorulin caused the cells to detach from one another but did not induce invasiveness in these cells, as measured by their ability to cross a Matrigel-coated polycarbonate filter in a modified Boyden chamber assay. Two uvomorulin-negative, vimentin-negative cell lines are also not highly invasive as measured by this assay. We suggest that loss of uvomorulin-mediated cell-cell adhesion may be one of many changes involved in the progression of a carcinoma cell to an invasive phenotype.