UNIFOLIATA regulates leaf and flower morphogenesis in pea

Background The vegetative phenotype of the pea mutant unifoliata (uni) is a simplification of the wild-type compound leaf to a single leaflet. Mutant uni plants are also self-sterile and the flowers resemble known floral meristem and organ identity mutants. In Antirrhinum and Arabidopsis, mutations in the floral meristem identity gene FLORICAULA/LEAFY (FLO/LFY) affect flower development alone, whereas the tobacco FLO/LFY homologue, NFL, is expressed in vegetative tissues, suggesting that NFL specifies determinacy in the progenitor cells for both flowers and leaves. In this paper, we characterised the pea homologue of FLO/LFY. Results The pea cDNA homologue of FLO/LFY, PEAFLO, mapped to the uni locus in recombinant-inbred mapping populations and markers based on PEAFLO cosegregated with uni in segregating sibling populations. The characterisation of two spontaneous uni mutant alleles, one containing a deletion and the other a point mutation in the PEAFLO coding sequences, predicted that PEAFLO corresponds to UNI and that the mutant vegetative phenotype was conferred by the defective PEAFLO gene. Conclusions The uni mutant demonstrates that there are shared regulatory processes in the morphogenesis of leaves and flowers and that floral meristem identity genes have an extended role in plant development. Pleiotropic regulatory genes such as UNI support the hypothesis that leaves and flowers derive from a common ancestral sporophyll-like structure. The regulation of indeterminacy during leaf and flower morphogenesis by UNI may reflect a primitive function for the gene in the pre-angiosperm era.

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.

Identification and characterization of flowering genes in kiwifruit : sequence conservation and role in kiwifruit flower development

Background Flower development in kiwifruit (Actinidia spp.) is initiated in the first growing season, when undifferentiated primordia are established in latent shoot buds. These primordia can differentiate into flowers in the second growing season, after the winter dormancy period and upon accumulation of adequate winter chilling. Kiwifruit is an important horticultural crop, yet little is known about the molecular regulation of flower development. Results To study kiwifruit flower development, nine MADS-box genes were identified and functionally characterized. Protein sequence alignment, phenotypes obtained upon overexpression in Arabidopsis and expression patterns suggest that the identified genes are required for floral meristem and floral organ specification. Their role during budbreak and flower development was studied. A spontaneous kiwifruit mutant was utilized to correlate the extended expression domains of these flowering genes with abnormal floral development. Conclusions This study provides a description of flower development in kiwifruit at the molecular level. It has identified markers for flower development, and candidates for manipulation of kiwifruit growth, phase change and time of flowering. The expression in normal and aberrant flowers provided a model for kiwifruit flower development.

The use of artificial MicroRNA technology to control gene expression in Arabidopsis thaliana

In plants, double-stranded RNA (dsRNA) is an effective trigger of RNA silencing, and several classes of endogenous small RNA (sRNA), processed from dsRNA substrates by DICER-like (DCL) endonucleases, are essential in controlling gene expression. One such sRNA class, the microRNAs (miRNAs) control the expression of closely related genes to regulate all aspects of plant development, including the determination of leaf shape, leaf polarity, flowering time, and floral identity. A single miRNA sRNA silencing signal is processed from a long precursor transcript of nonprotein-coding RNA, termed the primary miRNA (pri-miRNA). A region of the pri-miRNA is partially self-complementary allowing the transcript to fold back onto itself to form a stem-loop structure of imperfectly dsRNA. Artificial miRNA (amiRNA) technology uses endogenous pri-miRNAs, in which the miRNA and miRNA*(passenger strand of the miRNA duplex) sequences have been replaced with corresponding amiRNA/ amiRNA*sequences that direct highly efficient RNA silencing of the targeted gene. Here, we describe the rules for amiRNA design, as well as outline the PCR and bacterial cloning procedures involved in the construction of an amiRNA plant expression vector to control target gene expression in Arabidopsis thaliana. © 2014 Springer Science+Business Media New York.

Multiscale analysis of protein functions and stochastic modelling of gene transcriptional regulatory networks

Genomic and proteomic analyses have attracted a great deal of interests in biological research in recent years. Many methods have been applied to discover useful information contained in the enormous databases of genomic sequences and amino acid sequences. The results of these investigations inspire further research in biological fields in return. These biological sequences, which may be considered as multiscale sequences, have some specific features which need further efforts to characterise using more refined methods. This project aims to study some of these biological challenges with multiscale analysis methods and stochastic modelling approach. The first part of the thesis aims to cluster some unknown proteins, and classify their families as well as their structural classes. A development in proteomic analysis is concerned with the determination of protein functions. The first step in this development is to classify proteins and predict their families. This motives us to study some unknown proteins from specific families, and to cluster them into families and structural classes. We select a large number of proteins from the same families or superfamilies, and link them to simulate some unknown large proteins from these families. We use multifractal analysis and the wavelet method to capture the characteristics of these linked proteins. The simulation results show that the method is valid for the classification of large proteins. The second part of the thesis aims to explore the relationship of proteins based on a layered comparison with their components. Many methods are based on homology of proteins because the resemblance at the protein sequence level normally indicates the similarity of functions and structures. However, some proteins may have similar functions with low sequential identity. We consider protein sequences at detail level to investigate the problem of comparison of proteins. The comparison is based on the empirical mode decomposition (EMD), and protein sequences are detected with the intrinsic mode functions. A measure of similarity is introduced with a new cross-correlation formula. The similarity results show that the EMD is useful for detection of functional relationships of proteins. The third part of the thesis aims to investigate the transcriptional regulatory network of yeast cell cycle via stochastic differential equations. As the investigation of genome-wide gene expressions has become a focus in genomic analysis, researchers have tried to understand the mechanisms of the yeast genome for many years. How cells control gene expressions still needs further investigation. We use a stochastic differential equation to model the expression profile of a target gene. We modify the model with a Gaussian membership function. For each target gene, a transcriptional rate is obtained, and the estimated transcriptional rate is also calculated with the information from five possible transcriptional regulators. Some regulators of these target genes are verified with the related references. With these results, we construct a transcriptional regulatory network for the genes from the yeast Saccharomyces cerevisiae. The construction of transcriptional regulatory network is useful for detecting more mechanisms of the yeast cell cycle.

Cloning and characterization of a novel human gene related to vascular endothelial growth factor

This paper describes the cloning and characterization of a new member of the vascular endothelial growth factor (VEGF) gene family, which we have designated VRF for VEGF-related-factor. Sequencing of cDNAs from a human fetal brain library and RT-PCR products from normal and tumor tissue cDNA pools indicate two alternatively spliced messages with open reading frames of 621 and 564 bp, respectively. The predicted proteins differ at their carboxyl ends resulting from a shift in the open reading frame. Both isoforms show strong homology to VEGF at their amino termini, but only the shorter isoform maintains homology to VEGF at its carboxyl terminus and conserves all 16 cysteine residues of VEGF165. Similarity comparisons of this isoform revealed overall protein identity of 48% and conservative substitution of 69% with VEGF189. VRF is predicted to contain a signal peptide, suggesting that it may be a secreted factor. The VRF gene maps to the D11S750 locus at chromosome band 11q13, and the protein coding region, spanning approximately 5 kb, is comprised of 8 exons that range in size from 36 to 431 bp. Exons 6 and 7 are contiguous and the two isoforms of VRF arise through alternate splicing of exon 6. VRF appears to be ubiquitously expressed as two transcripts of 2.0 and 5.5 kb; the level of expression is similar among normal and malignant tissues.

Prediction of multilocus identity-by-descent

Previous studies have enabled exact prediction of probabilities of identity-by-descent (IBD) in randommating populations for a few loci (up to four or so), with extension to more using approximate regression methods. Here we present a precise predictor of multiple-locus IBD using simple formulas based on exact results for two loci. In particular, the probability of non-IBD X ABC at each of ordered loci A, B, and C can be well approximated by XABC = XABXBC/XB and generalizes to X123. . .k = X12X23. . .Xk-1,k/ Xk-2, where X is the probability of non-IBD at each locus. Predictions from this chain rule are very precise with population bottlenecks and migration, but are rather poorer in the presence of mutation. From these coefficients, the probabilities of multilocus IBD and non-IBD can also be computed for genomic regions as functions of population size, time, and map distances. An approximate but simple recurrence formula is also developed, which generally is less accurate than the chain rule but is more robust with mutation. Used together with the chain rule it leads to explicit equations for non-IBD in a region. The results can be applied to detection of quantitative trait loci (QTL) by computing the probability of IBD at candidate loci in terms of identity-by-state at neighboring markers.

Prediction of IBD based on population history for fine gene mapping

A novel multiple regression method (RM) is developed to predict identity-by-descent probabilities at a locus L (IBDL), among individuals without pedigree, given information on surrounding markers and population history. These IBDL probabilities are a function of the increase in linkage disequilibrium (LD) generated by drift in a homogeneous population over generations. Three parameters are sufficient to describe population history: effective population size (Ne), number of generations since foundation (T), and marker allele frequencies among founders (p). IBD L are used in a simulation study to map a quantitative trait locus (QTL) via variance component estimation. RM is compared to a coalescent method (CM) in terms of power and robustness of QTL detection. Differences between RM and CM are small but significant. For example, RM is more powerful than CM in dioecious populations, but not in monoecious populations. Moreover, RM is more robust than CM when marker phases are unknown or when there is complete LD among founders or Ne is wrong, and less robust when p is wrong. CM utilises all marker haplotype information, whereas RM utilises information contained in each individual marker and all possible marker pairs but not in higher order interactions. RM consists of a family of models encompassing four different population structures, and two ways of using marker information, which contrasts with the single model that must cater for all possible evolutionary scenarios in CM.

Gene Silencing in Arabidopsis Spreads from the Root to the Shoot, through a Gating Barrier, by Template-Dependent, Nonvascular, Cell-to-Cell Movement

Upward long-distance mobile silencing has been shown to be phloem mediated in several different solanaceous species. We show that the Arabidopsis (Arabidopsis thaliana) seedling grafting system and a counterpart inducible system generate upwardly spreading long-distance silencing that travels not in the phloem but by template-dependent reiterated short-distance cell-to-cell spread through the cells of the central stele. Examining the movement of the silencing front revealed a largely unrecognized zone of tissue, below the apical meristem, that is resistant to the silencing signal and that may provide a gating or protective barrier against small RNA signals. Using a range of auxin and actin transport inhibitors revealed that, in this zone, alteration of vesicular transport together with cytoskeleton dynamics prevented or retarded the spread of the silencing signal. This suggests that small RNAs are transported from cell to cell via plasmodesmata rather than diffusing from their source in the phloem.

Comparison of the coat protein, movement protein and RNA polymerase gene sequences of Australian, Chinese, and American isolates of barley yellow dwarf virus transmitted by Rhopalosiphum padi

Barley yellow dwarf luteovirus-GPV (BYDV-GPV) is a common problem in Chinese wheat crops but is unrecorded elsewhere. A defining characteristic of GPV is its capacity to be transmitted efficiently by both Schizaphis graminum and Rhopaloshiphum padi. This dual aphid species transmission contrasts with those of BYDV-RPV and BYDV-SGV, globally distributed viruses, which are efficiently transmitted only by Rhopaloshiphum padi and Schizaphis graminum respectively. The viral RNA sequences encoding the coat protein (22K) gene, the movement protein (17K) gene, the region surrounding the conserved GDD motif of the polymerase gene and the intergenic sequences between these genes were determined for GPV and an Australian isolate of BYDV-RPV (RPVa). In all three genes, the sequences of GPV and RPVa were more similar to those of an American isolate of BYDV-RPV (RPVu) than to any other luteovirus for which there is data available. RPVa and RPVu were very similar, especially their coat proteins which had 97% identity at the amino acid level. The coat protein of GPV had 76% and 78% amino acid identity with RPVa and RPVu respectively. The data suggest that RPVu and RPVa are correctly named as strains of the same serotype and that GPV is sufficiently different from either RPV strain to be considered a distinct BYDV type. The coat protein and movement protein genes of GPV are very dissimilar to SGV. The polymerase sequences of RPVu, RPVa and GPV show close affinities with those of the sobemo-like luteoviruses and little similarity with those of the carmo-like luteoviruses. The sequences of the coat proteins, movement proteins and the polymerase segments of BYDV serotypes, other than RPV and GPV, form a cluster that is separate from their counterpart sequences from dicot-infecting luteoviruses. The RPV and GPV isolates consistently fall within a dicot-infecting cluster. This suggests that RPV and GPV evolved from within this group of viruses. Since these other viruses all infect dicots it seems likely that their common ancestor infected a dicot and that RPV and GPV evolved from a virus that switched hosts from a dicot to a monocot.

Nucleotide sequence of coat protein gene for GPV isolate of barley yellow dwarf virus and construction of expression plasmid for plant

GPV is a Chinese serotype isolate of barley yellow dwarf virus (BYDV) that has no reaction with antiserum of MAV, PAV, SGV, RPV and RMV The sequence of the coat protein (CP) of GPV isolate of BYDV was identified and its amino acid sequence was deduced. The coding region for the putative GPV CP is 603 bases nucleotides and encodes a Mr 22 218 (22 ku) protein. The same as MAV, PAV and RPV, GPV contained a second ORF within the coat protein coding region. This protein of 17 024 Mr (17 ku) is thought to correspond to the Virion protein genome linked (Vpg). Sequence comparisons of the CP coding region between the GPV isolate of BYDV and other isolates of BYDV have been done. The nucleotide and amino acid sequence homology of GPV has a greater identity to the sequence of RPV than those of PAV and MAV. The GPV CP sequence stored 83.7% of nucleotide similarity and 77.5% of deduced amino acid similarity, whereas that of the PAV and MAV shared 56.9%, 53.2% and 44.1%, 43.8% respectively. According to BYDV-GPV CP sequence, two primers were designed. The cDNA of CP was produced by RT-PCR. Full-length cDNA of CP was inserted into plasmid to construct expression plasmids named pPPI1, pPPI2 and pPPI5 based on different promoters. The recombinant plasmids were identified by using α-32P-dATP labelled CP probe, α-32P-ATP labelled GPV RNA probe and sequencing to confirm real GPV CP gene cDNA in plasmids.

Overexpression of the kiwifruit SVP3 gene affects reproductive development and suppresses anthocyanin biosynthesis in petals, but has no effect on vegetative growth, dormancy, or flowering time

SVP-like MADS domain transcription factors have been shown to regulate flowering time and both inflorescence and flower development in annual plants, while having effects on growth cessation and terminal bud formation in perennial species. Previously, four SVP genes were described in woody perennial vine kiwifruit (Actinidia spp.), with possible distinct roles in bud dormancy and flowering. Kiwifruit SVP3 transcript was confined to vegetative tissues and acted as a repressor of flowering as it was able to rescue the Arabidopsis svp41 mutant. To characterize kiwifruit SVP3 further, ectopic expression in kiwifruit species was performed. Ectopic expression of SVP3 in A. deliciosa did not affect general plant growth or the duration of endodormancy. Ectopic expression of SVP3 in A. eriantha also resulted in plants with normal vegetative growth, bud break, and flowering time. However, significantly prolonged and abnormal flower, fruit, and seed development were observed, arising from SVP3 interactions with kiwifruit floral homeotic MADS-domain proteins. Petal pigmentation was reduced as a result of SVP3-mediated interference with transcription of the kiwifruit flower tissue-specific R2R3 MYB regulator, MYB110a, and the gene encoding the key anthocyanin biosynthetic step, F3GT1. Constitutive expression of SVP3 had a similar impact on reproductive development in transgenic tobacco. The flowering time was not affected in day-neutral and photoperiod-responsive Nicotiana tabacum cultivars, but anthesis and seed germination were significantly delayed. The accumulation of anthocyanin in petals was reduced and the same underlying mechanism of R2R3 MYB NtAN2 transcript reduction was demonstrated.

Identity-by-Descent Mapping to Detect Rare Variants Conferring Susceptibility to Multiple Sclerosis

Genome-wide association studies (GWAS) have identified around 60 common variants associated with multiple sclerosis (MS), but these loci only explain a fraction of the heritability of MS. Some missing heritability may be caused by rare variants that have been suggested to play an important role in the aetiology of complex diseases such as MS. However current genetic and statistical methods for detecting rare variants are expensive and time consuming. 'Population-based linkage analysis' (PBLA) or so called identity-by-descent (IBD) mapping is a novel way to detect rare variants in extant GWAS datasets. We employed BEAGLE fastIBD to search for rare MS variants utilising IBD mapping in a large GWAS dataset of 3,543 cases and 5,898 controls. We identified a genome-wide significant linkage signal on chromosome 19 (LOD = 4.65; p = 1.9×10-6). Network analysis of cases and controls sharing haplotypes on chromosome 19 further strengthened the association as there are more large networks of cases sharing haplotypes than controls. This linkage region includes a cluster of zinc finger genes of unknown function. Analysis of genome wide transcriptome data suggests that genes in this zinc finger cluster may be involved in very early developmental regulation of the CNS. Our study also indicates that BEAGLE fastIBD allowed identification of rare variants in large unrelated population with moderate computational intensity. Even with the development of whole-genome sequencing, IBD mapping still may be a promising way to narrow down the region of interest for sequencing priority. © 2013 Lin et al.

Mobile gene silencing in Arabidopsis is regulated by hydrogen peroxide

In plants and nematodes, RNAi can spread from cells from which it is initiated to other cells in the organism. The underlying mechanism controlling the mobility of RNAi signals is not known, especially in the case of plants. A genetic screen designed to recover plants impaired in the movement but not the production or effectiveness of the RNAi signal identified RCI3, which encodes a hydrogen peroxide (H2O2)-producing type III peroxidase, as a key regulator of silencing mobility in Arabidopsis thaliana. Silencing initiated in the roots of rci3 plants failed to spread into leaf tissue or floral tissue. Application of exogenous H2O2 reinstated the spread in rci3 plants and accelerated it in wild-type plants. The addition of catalase or MnO2, which breaks down H2O2, slowed the spread of silencing in wild-type plants. We propose that endogenous H2O2, under the control of peroxidases, regulates the spread of gene silencing by altering plasmodesmata permeability through remodelling of local cell wall structure, and may play a role in regulating systemic viral defence.