<Emphasis Type="Italic">BcMF11</Emphasis> and its homologous sequences may form a lncRNA family in <Emphasis Type="Italic">Brassica</Emphasis> diploids

BcMF11 is a long non-coding RNA that has been identified in Brassica rapa and shown to be involved in pollen development. Here, when re-cloned the gene sequence, multiple paralogous copies of BcMF11 were identified in B. rapa (A genome). Multiple paralogous copies of BcMF11 were also found in B. nigra (B genome) and Brassica oleracea (C genome), the other two primary diploids of Brassica U triangle. While in the early diverging Brassicaceae lineage including Arabidopsis thaliana, no BcMF11 homolog was found. Phylogenetic analysis showed that the BcMF11 homologous sequences cloned from A genome or C genome could be clustered into a separate branch, respectively. However, there was no distinct cluster defined for BcMF11 homologous sequences cloned from B genome. The expression of BcMF11 in B. rapa was investigated and revealed a different result in the previous study. In addition, 12 expressed sequence tags from B. napus and B. rapa showing high similarities with BcMF11 were identified in the NCBI database, which further verified that rather than the useless repeat fragments in the genome, the BcMF11 homologous genes could transcribe. It is possible that BcMF11 and its homologous sequences may form a large gene family which might be originated in the recent ancestral lineage of Brassica.     

Natural genetic variation in <Emphasis Type="Italic">Brassica</Emphasis> homologs of <Emphasis Type="Italic">FLOWERING LOCUS T</Emphasis> and characterization of its expression domains

FLOWERING LOCUS T (FT), a major effect gene, regulates flowering time in Arabidopsis. We analyzed evolutionary changes distinguishing two FT homeologous loci in B. rapa, described genetic variation in homologs isolated and reported expression pattern of FT in B. juncea. Synteny analysis confirmed presence of two FT genomic copies in B. rapa ssp. pekinensis and resolved pre-existing anomalies regarding copy number in “AA” genome. Synteny analysis of B. rapa homeologous regions CR1 (129 kb) and CR2 (232 kb) revealed differential gene fractionation and wide-spread re-arrangements. Seven genomic DNA (gDNA) variants (2.1–2.2 kb) and 10 complementary DNA (cDNA) variants (528 bp) were isolated from 6 Brassica species. The gDNA variants shared 72–99 % similarity within Brassica and 58–60 % between Arabidopsis and Brassica. FT cDNA variants shared 92–100 % similarity within Brassica and 87 % between Arabidopsis and Brassica. Phylogenetic analysis of FT gDNA, cDNA and protein sequences revealed two major clades, differentiating homologs derived from species containing shared “BB” and “CC” genomes. Phylogram based on Brassica FT gDNA differentiated homeologs derived from AA-LF (Least fractioned) and AA-MF1 (Moderately fractioned) sub-genomes. Analysis of FT expression pattern in B. juncea revealed increasing levels correlating with attainment of physiological maturity; highest levels were detected in older leaves implying conservation in spatio-temporal expression pattern vis-à-vis Arabidopsis. In conclusion, our study reveals that polyploidy in Brassicas resulted in expansion of FT gene copies with homologs charting independent evolutionary course through accumulation of mutations. However, expression domains of FT remained conserved across Brassicaceae to preserve the critical function of FT in controlling flowering time.     

Expression analysis of lncRNA <Emphasis Type="Italic">AK370814</Emphasis> involved in the barley vitamin B6 salvage pathway under salinity

Long non-coding RNAs (lncRNAs), which are longer than >?200 nt, perform various functions in a variety of important biological processes. The aim of this study is the investigation of relative expression levels of AK372815 putative pyridoxal reductase (PLR) gene and sense lncRNA AK370814 on four barley genotypes (Hasat, Beysehir 99, Konevi 98 and Tarm 92) in response to 150 mM salinity application during 3 days post-germination. Seeds were placed randomly in petri dishes containing (a) only H₂O (control), (b) 150 mM NaCl, for 72 h. RNA isolation was carried out using TriPure® reagent from 150 mM salt-treated root and shoot samples. Relative expression levels of AK372815 PLR and sense lncRNA AK370814 were determined by qPCR. Results demonstrated that salinity affected the expression levels of both AK372815 PLR gene and sense lncRNA AK370814 during germination. Although expression levels of AK372815 PLR tended to be down-regulated under salinity, expression levels of sense lncRNA AK370814 were up-regulated. Another goal of this study is improvement of alternative approach to NGS technologies for determination of relative expression levels of sense lncRNAs under particular circumstances. This is the first report that demonstrates a relationship between lncRNA and vitamin B6 salvage pathway.     

Inheritance and gene mapping of the white flower trait in <Emphasis Type="Italic">Brassica juncea</Emphasis>

Despite being a unique marker trait, white flower inheritance in Brassica juncea remains poorly understood at the gene level. In this study, we investigated a B. juncea landrace with white petal in China. The white petal phenotype possessed defective chromoplasts with less plastoglobuli than the yellow petal phenotype. Genetic analysis confirmed that two independent recessive genes (Bjpc1 and Bjpc2) controlled the white flower trait. We then mapped the BjPC1 gene in a BC₄ population comprising 2295 individuals. We identified seven AFLP (amplified fragment length polymorphism) markers closely linked to the white flower gene. BLAST search revealed the sequence of AFLP fragments were highly homologous with the Scaffold000085 and Scaffold000031 sequences on the A02 chromosome in the Brassica rapa genome. Based on this sequence homology, we developed simple sequence repeat (SSR) primer pairs and identified 13 SSRs linked to the BjPC1 gene, including two that were co-segregated (SSR9 and SSR10). The two closest markers (SSR4 and SSR11) were respectively 0.9 and 0.4 cM on either side of BjPC1. BLAST analysis revealed that these marker sequences corresponded highly to A02 in B. juncea. They were mapped within a 33 kb genomic region on B. rapa A02 (corresponds to a 40 kb genomic region on B. juncea A02) that included three genes. Sequence BjuA008406, homologous to AtPES2 in Arabidopsis thaliana and Bra032956 in B. rapa, was the most likely candidate for BjPC1. These results should accelerate BjPC1 cloning and facilitate our understanding of the molecular mechanisms controlling B. juncea petal color.     

An Immune-Related Six-lncRNA Signature to Improve Prognosis Prediction of Glioblastoma Multiforme

Recent studies have demonstrated the utility and superiority of long non-coding RNAs (lncRNAs) as novel biomarkers for cancer diagnosis, prognosis, and therapy. In the present study, the prognostic value of lncRNAs in glioblastoma multiforme was systematically investigated by performing a genome-wide analysis of lncRNA expression profiles in 419 glioblastoma patients from The Cancer Genome Atlas (TCGA) project. Using survival analysis and Cox regression model, we identified a set of six lncRNAs (AC005013.5, UBE2R2-AS1, ENTPD1-AS1, RP11-89C21.2, AC073115.6, and XLOC_004803) demonstrating an ability to stratify patients into high- and low-risk groups with significantly different survival (median 0.899 vs. 1.611 years, p = 3.87e?09, log-rank test) in the training cohort. The six-lncRNA signature was successfully validated on independent test cohort of 219 patients with glioblastoma, and it revealed superior performance for risk stratification with respect to existing lncRNA-related signatures. Multivariate Cox and stratification analysis indicated that the six-lncRNA signature was an independent prognostic factor after adjusting for other clinical covariates. Further in silico functional analysis suggested that the six-lncRNA signature may be involved in the immune-related biological processes and pathways which are very well known in the context of glioblastoma tumorigenesis. The identified lncRNA signature had important clinical implication for improving outcome prediction and guiding the tailored therapy for glioblastoma patients with further prospective validation.     

Small RNA changes in synthetic <Emphasis Type="Italic">Brassica napus</Emphasis>

Main conclusion

Small RNAs and microRNAs were found to vary extensively in synthetic Brassica napus and subsequent generations, accompanied by the activation of transposable elements in response to hybridization and polyploidization.

Abstract

Resynthesizing B. napus by hybridization and chromosome doubling provides an approach to create novel polyploids and increases the usable genetic variability in oilseed rape. Although many studies have shown that small RNAs (sRNAs) act as important factor during hybridization and polyploidization in plants, much less is known on how sRNAs change in synthetic B. napus, particularly in subsequent generations after formation. We performed high-throughput sequencing of sRNAs in S₁–S₄ generations of synthetic B. napus and in the homozygous B. oleracea and B. rapa parent lines. We found that the number of small RNAs (sRNAs) and microRNAs (miRNAs) doubled in synthetic B. napus relative to the parents. The proportions of common sRNAs detected varied from the S₁ to S₄ generations, suggesting sRNAs are unstable in synthetic B. napus. The majority of miRNAs (67.2 %) were non-additively expressed in the synthesized Brassica allotetraploid, and 33.3 % of miRNAs were novel in the resynthesized B. napus. The percentage of miRNAs derived from transposable elements (TEs) also increased, indicating transposon activation and increased transposon-associated miRNA production in response to hybridization and polyploidization. The number of target genes for each miRNA in the synthesized Brassica allotetraploid was doubled relative to the parents, enhancing the complexity of gene expression regulation. The potential roles of miRNAs and their targets are discussed. Our data demonstrate generational changes in sRNAs and miRNAs in synthesized B. napus.

     

Patterns of subspecies genetic diversity among oilseed <Emphasis Type="Italic">Brassica rapa</Emphasis> as revealed by agro-morphological traits and SSR markers

Brassica rapa (2n = 20, AA genome) is an important oil yielding species of the family Brassicaceae and characterized by wide range of genetic and morphological subtypes suitable for cultivation under diverse agro-climatic regions of India. In this study, genetic diversity among three subspecies of B. rapa including yellow sarson, toria and outlier brown sarson was estimated using various agro-morphological traits and simple sequence repeat (SSR) markers. Maximum variability was recorded for siliqua angle (Coefficient of variation = 30.9%), followed by seeds/siliqua (CV = 18.8%), leaf length (CV = 10%) and plant height (CV = 16.8%). Principal component analysis explained more than 50% of the total observed morphological variability for first two components. Of the 107 SSR markers tested, 80 generated reproducible, clear and distinct amplicons of which, 65 (81.25%) were found polymorphic. The number of alleles at each locus ranged from 2 to 7, with an average of 3.03 alleles per marker. A total of 197 alleles were detected at 65 SSR loci with average PIC value of 0.457 and a mean resolving power of 3.04. Neighbor-Joining cluster analysis based on morphological traits and SSR markers separately classified all the 28 genotypes into five major groups. The population structure analysis resulted into three sub-populations with certain extent of admixture among the earlier established taxonomic sub-groups. Twenty-three unique alleles were detected in thirteen B. rapa varieties. The clustering analysis and principal coordinate analysis outlined the genetic relationships among different varieties belonging to the three subspecies of B. rapa. Genetically diverse genotypes as illustrated by score plots and from the clustering patterns brought out the wide range of diversity present among B. rapa genotypes and the underlying options available for selecting parental genotypes for hybridization and developing high yielding cultivars suitable for Indian conditions.     

Microarray analysis reveals a potential role of LncRNAs expression in cardiac cell proliferation

Background

Long non-coding RNAs (LncRNAs) have been identified to play important roles in epigenetic processes that underpin organogenesis. However, the role of LncRNAs in the regulation of transition from fetal to adult life of human heart has not been evaluated.

Methods

Immunofiuorescent staining was used to determine the extent of cardiac cell proliferation. Human LncRNA microarrays were applied to define gene expression signatures of the fetal (13–17 weeks of gestation, n?=?4) and adult hearts (30–40 years old, n?=?4). Pathway analysis was performed to predict the function of differentially expressed mRNAs (DEM). DEM related to cell proliferation were selected to construct a lncRNA-mRNA co-expression network. Eight lncRNAs were confirmed by quantificational real-time polymerase chain reaction (n?=?6).

Results

Cardiac cell proliferation was significant in the fetal heart. Two thousand six hundred six lncRNAs and 3079 mRNAs were found to be differentially expressed. Cell cycle was the most enriched pathway in down-regulated genes in the adult heart. Eight lncRNAs (RP11-119 F7.5, AX747860, HBBP1, LINC00304, TPTE2P6, AC034193.5, XLOC_006934 and AL833346) were predicted to play a central role in cardiac cell proliferation.

Conclusions

We discovered a profile of lncRNAs differentially expressed between the human fetal and adult heart. Several meaningful lncRNAs involved in cardiac cell proliferation were disclosed.

     

Genome-wide characterization and stress-responsive expression profiling of <Emphasis Type="Italic">MCM</Emphasis> genes in <Emphasis Type="Italic">Brassica oleracea</Emphasis> and <Emphasis Type="Italic">Brassica rapa</Emphasis>

The Minichromosome maintenance protein [MCM (2-7)] complex is associated with helicase activity for replication fork formation during DNA replication. We identified and characterized each 12 putative MCM genes from Brassica oleracea and Brassica rapa. MCM genes were classified into nine groups according to their evolutionary relationships. A high number of syntenic regions were present on chromosomes C03 and A03 in B. oleracea and B. rapa, respectively, compared to the other chromosomes. Expression analysis showed that most of the MCM(2-7) helicase-subunit genes and their coregulating MCM genes were upregulated during hydroxyurea (HU) induced stress in B. oleracea. In B. rapa, MCM(2-7) helicase genes BrMCM2_2, BrMCM7_1, BrMCM7_2 and their co-regulating genes were upregulated during replication stress. During cold stress, BoMCM6 in B. oleracea and BrMCM5 in B. rapa were remarkably upregulated. During salt stress, BoMCM6_2, BoMCM7_1, BoMCM8, BoMCM9, and BoMCM10 were markedly upregulated in B. oleracea. Hence, our study identified the candidate MCM family genes those possess abiotic stress-responsive behavior and DNA replication stress tolerance. As the first genome-wide analysis of MCM genes in B. oleracea and B. rapa, this work provides a foundation to develop stress responsive plants. Further functional and molecular studies on MCM genes will be helpful to enhance stress tolerance in plants.     

Identification of candidate genes of QTLs for seed weight in <Emphasis Type="Italic">Brassica napus</Emphasis> through comparative mapping among <Emphasis Type="Italic">Arabidopsis</Emphasis> and <Emphasis Type="Italic">Brassica</Emphasis>species

Background

Map-based cloning of quantitative trait loci (QTLs) in polyploidy crop species remains a challenge due to the complexity of their genome structures. QTLs for seed weight in B. napus have been identified, but information on candidate genes for identified QTLs of this important trait is still rare.

Results

In this study, a whole genome genetic linkage map for B. napus was constructed using simple sequence repeat (SSR) markers that covered a genetic distance of 2,126.4 cM with an average distance of 5.36 cM between markers. A procedure was developed to establish colinearity of SSR loci on B. napus with its two progenitor diploid species B. rapa and B. oleracea through extensive bioinformatics analysis. With the aid of B. rapa and B. oleracea genome sequences, the 421 homologous colinear loci deduced from the SSR loci of B. napus were shown to correspond to 398 homologous loci in Arabidopsis thaliana. Through comparative mapping of Arabidopsis and the three Brassica species, 227 homologous genes for seed size/weight were mapped on the B. napus genetic map, establishing the genetic bases for the important agronomic trait in this amphidiploid species. Furthermore, 12 candidate genes underlying 8 QTLs for seed weight were identified, and a gene-specific marker for BnAP2 was developed through molecular cloning using the seed weight/size gene distribution map in B. napus.

Conclusions

Our study showed that it is feasible to identify candidate genes of QTLs using a SSR-based B. napus genetic map through comparative mapping among Arabidopsis and B. napus and its two progenitor species B. rapa and B. oleracea. Identification of candidate genes for seed weight in amphidiploid B. napus will accelerate the process of isolating the mapped QTLs for this important trait, and this approach may be useful for QTL identification of other traits of agronomic significance.