Hypervariable genes--experimental error or hidden dynamics

In a homogeneous group of samples, not all genes of high variability stem from experimental errors in microarray experiments. These expression variations can be attributed to many factors including natural biological oscillations or metabolic processes. The behavior of these genes can tease out important clues about naturally occurring dynamic processes in the organism or experimental system under study. We developed a statistical procedure for the selection of genes with high variability denoted hypervariable (HV) genes. After the exclusion of low expressed genes and a stabilizing log-transformation, the majority of genes have comparable residual variability. Based on an F-test, HV genes are selected as having a statistically significant difference from the majority of variability stabilized genes measured by the 'reference group'. A novel F-test clustering technique, further noted as 'F-means clustering', groups HV genes with similar variability patterns, presumably from their participation in a common dynamic biological process. F-means clustering establishes, for the first time, groups of co-expressed HV genes and is illustrated with microarray data from patients with juvenile rheumatoid arthritis and healthy controls.     

Application of independent component analysis to microarrays

We apply linear and nonlinear independent component analysis (ICA) to project microarray data into statistically independent components that correspond to putative biological processes, and to cluster genes according to over- or under-expression in each component. We test the statistical significance of enrichment of gene annotations within clusters. ICA outperforms other leading methods, such as principal component analysis, k-means clustering and the Plaid model, in constructing functionally coherent clusters on microarray datasets from Saccharomyces cerevisiae, Caenorhabditis elegans and human.     

Heritable,De Novo Resistance to Leaf Rust and Other Novel Traits in Selfed Descendants of Wheat Responding to Inoculation with Wheat Streak Mosaic Virus

Stable resistance to infection with Wheat streak mosaic virus (WSMV) can be evolved de novo in selfing bread wheat lines subjected to cycles of WSMV inoculation and selection of best-performing plants or tillers. To learn whether this phenomenon might be applied to evolve resistance de novo to pathogens unrelated to WSMV, we examined the responses to leaf rust of succeeding generations of the rust- and WSMV-susceptible cultivar ‘Lakin’ following WSMV inoculation and derived rust-resistant sublines. After three cycles of the iterative protocol five plants, in contrast to all others, expressed resistance to leaf and stripe rust. A subset of descendant sublines of one of these, ‘R1’, heritably and uniformly expressed the new trait of resistance to leaf rust. Such sublines, into which no genes from a known source of resistance had been introgressed, conferred resistance to progeny of crosses with susceptible parents. The F₁ populations produced from crosses between, respectively, susceptible and resistant ‘Lakin’ sublines 4-3-3 and 4-12-3 were not all uniform in their response to seedling inoculation with race TDBG. In seedling tests against TDBG and MKPS races the F₂s from F₁ populations that were uniformly resistant had 3∶1 ratios of resistant to susceptible individuals but the F₂s from susceptible F₁ progenitors were uniformly susceptible. True-breeding lines derived from resistant individuals in F₂ populations were resistant to natural stripe and leaf rust inoculum in the field, while the ‘Lakin’ progenitor was susceptible. The next generation of six of the ‘Lakin’-derived lines exhibited moderate to strong de novo resistance to stem rust races TPMK, QFCS and RKQQ in seedling tests while the ‘Lakin’ progenitor was susceptible. These apparently epigenetic effects in response to virus infection may help researchers fashion a new tool that expands the range of genetic resources already available in adapted germplasm.     

Exploring Candidate Genes for Pericarp Russet Pigmentation of Sand Pear (Pyrus pyrifolia) via RNA-Seq Data in Two Genotypes Contrasting for Pericarp Color

Sand pear (Pyrus pyrifolia) russet pericarp is an important trait affecting both the quality and stress tolerance of fruits. This trait is controlled by a relative complex genetic process, with some fundamental biological questions such as how many and which genes are involved in the process remaining elusive. In this study, we explored differentially expressed genes between the russet- and green-pericarp offspring from the sand pear (Pyrus pyrifolia) cv. ‘Qingxiang’ × ‘Cuiguan’ F1 group by RNA-seq-based bulked segregant analysis (BSA). A total of 29,100 unigenes were identified and 206 of which showed significant differences in expression level (log₂fold values>1) between the two types of pericarp pools. Gene Ontology (GO) analyses detected 123 unigenes in GO terms related to ‘cellular_component’ and ‘biological_process’, suggesting developmental and growth differentiations between the two types. GO categories associated with various aspects of ‘lipid metabolic processes’, ‘transport’, ‘response to stress’, ‘oxidation-reduction process’ and more were enriched with genes with divergent expressions between the two libraries. Detailed examination of a selected set of these categories revealed repressed expressions of candidate genes for suberin, cutin and wax biosynthesis in the russet pericarps.Genes encoding putative cinnamoyl-CoA reductase (CCR), cinnamyl alcohol dehydrogenase (CAD) and peroxidase (POD) that are involved in the lignin biosynthesis were suggested to be candidates for pigmentation of sand pear russet pericarps. Nine differentially expressed genes were analyzed for their expressions using qRT-PCR and the results were consistent with those obtained from Illumina RNA-sequencing. This study provides a comprehensive molecular biology insight into the sand pear pericarp pigmentation and appearance quality formation.     

Cinnamic Aldehyde,the main monomer component of Cinnamon,exhibits anti‐inflammatory property in OA synovial fibroblasts via TLR4/MyD88 pathway

Cinnamon is a wildly used traditional Chinese herbal medicine for osteoarthritis (OA) treatment, but the underlying mechanism remains ambiguous. The purpose of this study is to explore the mechanism of cinnamic aldehyde (CA), a bioactive substance extracted from Cinnamon, on synovial inflammation in OA. A total of 144 CA‐OA co‐targeted genes were identified by detect databases (PubChem, HIT, TCMSP, TTD, DrugBank and GeneCards). The results of GO enrichment analysis indicated that these co‐targeted genes have participated in many biological processes including ‘inflammatory response’, ‘cellular response to lipopolysaccharide’, ‘response to drug’, ‘immune response’, ‘lipopolysaccharide‐mediated signalling pathway’, etc. KEGG pathway analysis showed these co‐targeted genes were mainly enriched in ‘Toll‐like receptor signalling pathway’, ‘TNF signalling pathway’, ‘NF‐kappa B signalling pathway’, etc. Molecular docking demonstrated that CA could successfully bind to TLR2 and TLR4. The results of in vitro experiments showed no potential toxicity of 10, 20 and 50 μM/L CA on human OA FLS, and CA can significantly inhibit the inflammation in LPS‐induced human FLS. Further experimental mechanism evidence confirmed CA can inhibited the inflammation in LPS‐induced human OA FLS via blocking the TLR4/MyD88 signalling pathway. Our results demonstrated that CA exhibited strong anti‐inflammation effect in OA FLS through blocking the activation of TLR4/MyD88 signalling pathway, suggesting its potential as a hopeful candidate for the development of novel agents for the treatment of OA.     

Genome sequencing and analysis of two early-flowering cherry (Cerasus × kanzakura) varieties, ‘Kawazu-zakura’ and ‘Atami-zakura’

To gain genetic insights into the early-flowering phenotype of ornamental cherry, also known as sakura, we determined the genome sequences of two early-flowering cherry (Cerasus × kanzakura) varieties, ‘Kawazu-zakura’ and ‘Atami-zakura’. Because the two varieties are interspecific hybrids, likely derived from crosses between Cerasus campanulata (early-flowering species) and Cerasus speciosa, we employed the haplotype-resolved sequence assembly strategy. Genome sequence reads obtained from each variety by single-molecule real-time sequencing (SMRT) were split into two subsets, based on the genome sequence information of the two probable ancestors, and assembled to obtain haplotype-phased genome sequences. The resultant genome assembly of ‘Kawazu-zakura’ spanned 519.8 Mb with 1,544 contigs and an N50 value of 1,220.5 kb, while that of ‘Atami-zakura’ totalled 509.6 Mb with 2,180 contigs and an N50 value of 709.1 kb. A total of 72,702 and 69,528 potential protein-coding genes were predicted in the genome assemblies of ‘Kawazu-zakura’ and ‘Atami-zakura’, respectively. Gene clustering analysis identified 2,634 clusters uniquely presented in the C. campanulata haplotype sequences, which might contribute to its early-flowering phenotype. Genome sequences determined in this study provide fundamental information for elucidating the molecular and genetic mechanisms underlying the early-flowering phenotype of ornamental cherry tree varieties and their relatives.     

Paleovirology of ‘syncytins’, retroviral env genes exapted for a role in placentation

The development of the emerging field of ‘paleovirology’ allows biologists to reconstruct the evolutionary history of fossil endogenous retroviral sequences integrated within the genome of living organisms and has led to the retrieval of conserved, ancient retroviral genes ‘exapted’ by ancestral hosts to fulfil essential physiological roles, syncytin genes being undoubtedly among the most remarkable examples of such a phenomenon. Indeed, syncytins are ‘new’ genes encoding proteins derived from the envelope protein of endogenous retroviral elements that have been captured and domesticated on multiple occasions and independently in diverse mammalian species, through a process of convergent evolution. Knockout of syncytin genes in mice provided evidence for their absolute requirement for placenta development and embryo survival, via formation by cell–cell fusion of syncytial cell layers at the fetal–maternal interface. These genes of exogenous origin, acquired ‘by chance’ and yet still ‘necessary’ to carry out a basic function in placental mammals, may have been pivotal in the emergence of mammalian ancestors with a placenta from egg-laying animals via the capture of a founding retroviral env gene, subsequently replaced in the diverse mammalian lineages by new env-derived syncytin genes, each providing its host with a positive selective advantage.     

Molecular Characterization of Barley 3H Semi-Dwarf Genes

The barley chromosome 3H accommodates many semi-dwarfing genes. To characterize these genes, the two-rowed semi-dwarf Chinese barley landrace ‘TX9425’ was crossed with the Australian barley variety ‘Franklin’ to generate a doubled haploid (DH) population, and major QTLs controlling plant height have been identified in our previous study. The major QTL derived from ‘TX9425’ was targeted to investigate the allelism of the semi-dwarf gene uzu in barley. Twelve sets of near-isogenic lines and a large NILF₂ fine mapping population segregating only for the dwarfing gene from ‘TX9425’ were developed. The semi-dwarfing gene in ‘TX9425’ was located within a 2.8 cM region close to the centromere on chromosome 3H by fine mapping. Molecular cloning and sequence analyses showed that the ‘TX9425’-derived allele contained a single nucleotide substitution from A to G at position 2612 of the HvBRI1 gene. This was apparently the same mutation as that reported in six-rowed uzu barley. Markers co-segregating with the QTL were developed from the sequence of the HvBRI1 gene and were validated in the ‘TX9425’/‘Franklin’ DH population. The other major dwarfing QTL derived from the Franklin variety was distally located on chromosome 3HL and co-segregated with the sdw1 diagnostic marker hv20ox2. A third dwarfing gene, expressed only in winter-sown trials, was identified and located on chromosome 3HS. The effects and interactions of these dwarfing genes under different growing conditions are discussed. These results improve our understanding of the genetic mechanisms controlling semi-dwarf stature in barley and provide diagnostic markers for the selection of semi-dwarfness in barley breeding programs.     

Differential Responses of CO2 Assimilation,Carbohydrate Allocation and Gene Expression to NaCl Stress in Perennial Ryegrass with Different Salt Tolerance

Little is known about the effects of NaCl stress on perennial ryegrass (Lolium perenne L.) photosynthesis and carbohydrate flux. The objective of this study was to understand the carbohydrate metabolism and identify the gene expression affected by salinity stress. Seventy-four days old seedlings of two perennial ryegrass accessions (salt-sensitive ‘PI 538976’ and salt-tolerant ‘Overdrive’) were subjected to three levels of salinity stress for 5 days. Turf quality in all tissues (leaves, stems and roots) of both grass accessions negatively and significantly correlated with GFS (Glu+Fru+Suc) content, except for ‘Overdrive’ stems. Relative growth rate (RGR) in leaves negatively and significantly correlated with GFS content in ‘Overdrive’ (P<0.01) and ‘PI 538976’ (P<0.05) under salt stress. ‘Overdrive’ had higher CO₂ assimilation and F_v/F_m than ‘PI 538976’. Intercellular CO₂ concentration, however, was higher in ‘PI 538976’ treated with 400 mM NaCl relative to that with 200 mM NaCl. GFS content negatively and significantly correlated with RGR in ‘Overdrive’ and ‘PI 538976’ leaves and in ‘PI 538976’ stems and roots under salt stress. In leaves, carbohydrate allocation negatively and significantly correlated with RGR (r² = 0.83, P<0.01) and turf quality (r² = 0.88, P<0.01) in salt-tolerant ‘Overdrive’, however, the opposite trend for salt-sensitive ‘PI 538976’ (r² = 0.71, P<0.05 for RGR; r² = 0.62, P>0.05 for turf quality). A greater up-regulation in the expression of SPS, SS, SI, 6-SFT gene was observed in ‘Overdrive’ than ‘PI 538976’. A higher level of SPS and SS expression in leaves was found in ‘PI 538976’ relative to ‘Overdrive’. Accumulation of hexoses in roots, stems and leaves can induce a feedback repression to photosynthesis in salt-stressed perennial ryegrass and the salt tolerance may be changed with the carbohydrate allocation in leaves and stems.     

Unique Features of a Japanese ‘Candidatus Liberibacter asiaticus’ Strain Revealed by Whole Genome Sequencing

Citrus greening (huanglongbing) is the most destructive disease of citrus worldwide. It is spread by citrus psyllids and is associated with phloem-limited bacteria of three species of α-Proteobacteria, namely, ‘Candidatus Liberibacter asiaticus’, ‘Ca. L. americanus’, and ‘Ca. L. africanus’. Recent findings suggested that some Japanese strains lack the bacteriophage-type DNA polymerase region (DNA pol), in contrast to the Floridian psy62 strain. The whole genome sequence of the pol-negative ‘Ca. L. asiaticus’ Japanese isolate Ishi-1 was determined by metagenomic analysis of DNA extracted from ‘Ca. L. asiaticus’-infected psyllids and leaf midribs. The 1.19-Mb genome has an average 36.32% GC content. Annotation revealed 13 operons encoding rRNA and 44 tRNA genes, but no typical bacterial pathogenesis-related genes were located within the genome, similar to the Floridian psy62 and Chinese gxpsy. In contrast to other ‘Ca. L. asiaticus’ strains, the genome of the Japanese Ishi-1 strain lacks a prophage-related region.     

Bayesian Network Expansion Identifies New ROS and Biofilm Regulators

Signaling and regulatory pathways that guide gene expression have only been partially defined for most organisms. However, given the increasing number of microarray measurements, it may be possible to reconstruct such pathways and uncover missing connections directly from experimental data. Using a compendium of microarray gene expression data obtained from Escherichia coli, we constructed a series of Bayesian network models for the reactive oxygen species (ROS) pathway as defined by EcoCyc. A consensus Bayesian network model was generated using those networks sharing the top recovered score. This microarray-based network only partially agreed with the known ROS pathway curated from the literature and databases. A top network was then expanded to predict genes that could enhance the Bayesian network model using an algorithm we termed ‘BN+1’. This expansion procedure predicted many stress-related genes (e.g., dusB and uspE), and their possible interactions with other ROS pathway genes. A term enrichment method discovered that biofilm-associated microarray data usually contained high expression levels of both uspE and gadX. The predicted involvement of gene uspE in the ROS pathway and interactions between uspE and gadX were confirmed experimentally using E. coli reporter strains. Genes gadX and uspE showed a feedback relationship in regulating each other''s expression. Both genes were verified to regulate biofilm formation through gene knockout experiments. These data suggest that the BN+1 expansion method can faithfully uncover hidden or unknown genes for a selected pathway with significant biological roles. The presently reported BN+1 expansion method is a generalized approach applicable to the characterization and expansion of other biological pathways and living systems.