Cytonuclear Interactions and Subgenome Dominance Shape the Evolution of Organelle-Targeted Genes in the Brassica Triangle of U.

The interaction and coevolution between nuclear and cytoplasmic genomes are one of the fundamental hallmarks of eukaryotic genome evolution and, 2 billion yr later, are still major contributors to the formation of new species. Although many studies have investigated the role of cytonuclear interactions following allopolyploidization, the relative magnitude of the effect of subgenome dominance versus cytonuclear interaction on genome evolution remains unclear. The Brassica triangle of U features 3 diploid species that together have formed 3 separate allotetraploid species on similar evolutionary timescales, providing an ideal system for understanding the contribution of the cytoplasmic donor to hybrid polyploid. Here, we investigated the evolutionary pattern of organelle-targeted genes in Brassica carinata (BBCC) and 2 varieties of Brassica juncea (AABB) at the whole-genome level, with particular focus on cytonuclear enzyme complexes. We found partial evidence that plastid-targeted genes experience selection to match plastid genomes, but no obvious corresponding signal in mitochondria-targeted genes from these 2 separately formed allopolyploids. Interestingly, selection acting on plastid genomes always reduced the retention rate of plastid-targeted genes encoded by the B subgenome, regardless of whether the Brassica nigra (BB) subgenome was contributed by the paternal or maternal progenitor. More broadly, this study illustrates the distinct selective pressures experienced by plastid- and mitochondria-targeted genes, despite a shared pattern of inheritance and natural history. Our study also highlights an important role for subgenome dominance in allopolyploid genome evolution, even in genes whose function depends on separately inherited molecules.

Characterization of Firmiana danxiaensis plastomes and comparative analysis of Firmiana: insight into its phylogeny and evolution.

BACKGROUND: Firmiana danxiaensis is a critically endangered and ecologically important tree currently only found in four locations in Danxia or Karst habitats in northern Guangdong Province, China. The specialized habitat preference makes it an ideal model species for study of adaptive evolution. Meanwhile, the phylogenetic relationships of F. danxiaensis in four locations under two landforms are unclear. Therefore, we sequenced its complete chloroplast (cp.) genomes and conducted comprehensive interspecific and intrageneric plastome studies. RESULTS: The F. danxiaensis plastomes in four locations showed a typical quadripartite and circular structure that ranged from 160,832 to 161,206 bp in size, with 112 unique genes encoded. Comparative genomics showed that the plastomes of F. danxiaensis were relatively conserved with high similarity of genome organization, gene number, GC content and SSRs. While the genomes revealed higher biased codon preferences in Karst habitat than those in Danxia habitats. Eighteen and 11 divergent hotpots were identified at interspecific and intrageneric levels for species identification and further phylogenetic studies. Seven genes (clpP, accD, ccsA, ndhH, rpl20, rpoC2, and rps4) were under positive selection and may be related to adaptation. Phylogenetic analysis revealed that F. danxiaensis is sister to F. major and F. simplex. However, the interspecific relationships are not consistent with the habitat types. CONCLUSIONS: The characteristics and interspecific relationship of F. danxiaensis plastomes provide new insights into further integration of geographical factors, environmental factors, and genetic variations on the genomic study of F. danxiaensis. Together, our study will contribute to the study of species identification, population genetics, and conservation biology of F. danxiaensis.

Phylogeny, character evolution, and biogeography of the fern genus Bolbitis (Dryopteridaceae).

Bolbitis is a pantropical fern genus of Dryopteridaceae with ca. 80 species mainly in tropical Asia. Earlier studies confirmed the monophyly of Bolbitis when Mickelia is excluded and identified three major clades in Bolbitis. However, earlier studies are based on relatively small sampling and the majority of Asian species are not sampled. In this study, DNA sequences of three plastid markers of 169 accessions representing ca. 68 (85 % of total) species of Bolbitis in nine out of the 10 series recognized by Hennipman (1977), and 54 accessions representing the five remaining bolbitidoid genera are used to infer a global phylogeny with a focus on Asian species. The major results include: (1) Bolbitis is strongly supported as monophyletic; (2) species of Bolbitis are resolved into four major clades and their relationships are: the Malagasy/Mascarene clade is sister to the rest, followed by the African clade which is sister to the American clade + the Asian clade; (3) six well-supported subclades are identified in the most speciose Asian clade; (4) the free-veined Egenolfia is embedded in Bolbitis and is paraphyletic in relation to species with anastomosing venation; (5) three series sensu Hennipman (1977), B. ser. Alienae, B. ser. Egenolfianae, and B. ser. Heteroclitae, are paraphyletic or polyphyletic; (6) evolution of six morphological characters is analyzed and free venation is found to have evolved from anastomosing venation and reversed to free venation in Bolbitis; and (7) biogeographical implications are drawn and it is shown that a single recent dispersal from Asia resulted in continental disjunction of closely related ferns of Bolbitis between Africa and America.

Downregulation of LncRNA GCLC-1 Promotes Microcystin-LR-Induced Malignant Transformation of Human Liver Cells by Regulating GCLC Expression.

Microcystin-LR (MCLR) is an aquatic toxin, which could lead to the development of hepatocellular carcinoma (HCC). Long non-coding RNAs (lncRNAs) are considered important regulatory elements in the occurrence and development of cancer. However, the roles and mechanisms of lncRNAs during the process of HCC, induced by MCLR, remain elusive. Here, we identified a novel lncRNA, namely lnc-GCLC-1 (lncGCLC), which is in close proximity to the chromosome location of glutamate-cysteine ligase catalytic subunit (GCLC). We then investigated the role of lncGCLC in MCLR-induced malignant transformation of WRL68, a human hepatic cell line. During MCLR-induced cell transformation, the expression of lncGCLC and GCLC decreased continuously, accompanied with a consistently high expression of miR-122-5p. Knockdown of lncGCLC promoted cell proliferation, migration and invasion, but reduced cell apoptosis. A xenograft nude mouse model demonstrated that knockdown of lncGCLC promoted tumor growth. Furthermore, knockdown of lncGCLC significantly upregulated miR-122-5p expression, suppressed GCLC expression and GSH levels, and enhanced oxidative DNA damages. More importantly, the expression of lncGCLC in human HCC tissues was significantly downregulated in the high-microcystin exposure group, and positively associated with GCLC level in HCC tissues. Together, these findings suggest that lncGCLC plays an anti-oncogenic role in MCLR-induced malignant transformation by regulating GCLC expression.

Comparative Plastomes of Curcuma alismatifolia (Zingiberaceae) Reveal Diversified Patterns among 56 Different Cut-Flower Cultivars.

Curcuma alismatifolia (Zingiberaceae) is an ornamental species with high economic value due to its recent rise in popularity among floriculturists. Cultivars within this species have mixed genetic backgrounds from multiple hybridization events and can be difficult to distinguish via morphological and histological methods alone. Given the need to improve identification resources, we carried out the first systematic study using plastomic data wherein genomic evolution and phylogenetic relationships from 56 accessions of C. alismatifolia were analyzed. The newly assembled plastomes were highly conserved and ranged from 162,139 bp to 164,111 bp, including 79 genes that code for proteins, 30 tRNA genes, and 4 rRNA genes. The A/T motif was the most common of SSRs in the assembled genomes. The Ka/Ks values of most genes were less than 1, and only two genes had Ka/Ks values above 1, which were rps15 (1.15), and ndhl (1.13) with petA equal to 1. The sequence divergence between different varieties of C. alismatifolia was large, and the percentage of variation in coding regions was lower than that in the non-coding regions. Such data will improve cultivar identification, marker assisted breeding, and preservation of germplasm resources.

The complete chloroplast genome of Eastern Asian fern Bolbitis laxireticulata (Dryopteridaceae).

Bolbitis laxireticulata is a potential ornamental plant, which is restricted to eastern Asia. Here, we sequenced the complete chloroplast (cp) genome of B. laxireticulata and constructed a phylogenetic cp tree of Dryopteridaceae to study their relationships. The cp genome of B. laxireticulata is 153,093 bp in length, being made up of large single-copy (LSC, 83,169 bp), small single-copy (SSC, 21,538 bp), and a pair of region inverted repeats (IRs, 24,193 bp). It has 124 genes including 83 protein-coding genes, 33 tRNA genes, and eight rRNA genes. With the maximum-likelihood tree indicating, B. laxireticulata is more closely related to B. subcordata.

The Glu69Asp Polymorphism of EME1 Gene is Associated with an Increased Risk of Hepatocellular Carcinoma in Guangxi Population, China.

Background: The dysfunction of Essential meiotic endonuclease 1 homolog 1 (EME1) can lead to genomic instability and tumorigenesis. Single nucleotide polymorphisms (SNPs) in the EME1 gene have been reported to be associated with the risk of several cancers, but its association with hepatocellular carcinoma (HCC) has not been investigated. This study aimed to determine the association between EME1 SNPs and the risk of HCC. Methods: This study included 645 HCC patients and 649 healthy controls from a Guangxi population of Southern China, and genotyped three functional SNPs (Glu69Asp: rs3760413A>C, Ile350Thr: rs12450550T>C, and rs11868055A>G) of the EME1 gene utilizing the Agena MassARRAY platform. Results: The rs3760413C variant genotypes (AC+CC: Glu/Asp+Asp/Asp) conferred a 1.419-fold risk of HCC compared to the AA (Glu/Glu) genotype (adjusted OR = 1.419, 95% CI = 1.017-1.980), and the allele C increased the risk of HCC in a dose-dependent manner (P trend = 0.017). Moreover, the effects of the rs3760413C variant genotypes were more pronounced in individuals who drank pond/ditch water (adjusted OR = 3.956, 95% CI = 1.413-11.076) than in those who never drank (P = 0.033). We further observed that a potential carcinogen microcystin-LR induced more DNA oxidative damages in peripheral blood mononuclear cells from the carriers of rs3760413C variant genotypes than those from the subjects with AA genotype (P = 0.006). A nomogram was also constructed combining the rs3760413A>C polymorphism and environmental risk factors for predicting HCC risk with a good discriminatory ability (concordance index = 0.892, 95% CI: 0.874-0.911) and good calibration (mean absolute error = 0.005). Conclusion: Our data suggest that the Glu69Asp missense polymorphism (rs3760413) of EME1 gene is associated with the risk of HCC, which may be a susceptible biomarker of HCC in the Guangxi population.

Genome size evolution of the extant lycophytes and ferns.

Ferns and lycophytes have remarkably large genomes. However, little is known about how their genome size evolved in fern lineages. To explore the origins and evolution of chromosome numbers and genome size in ferns, we used flow cytometry to measure the genomes of 240 species (255 samples) of extant ferns and lycophytes comprising 27 families and 72 genera, of which 228 species (242 samples) represent new reports. We analyzed correlations among genome size, spore size, chromosomal features, phylogeny, and habitat type preference within a phylogenetic framework. We also applied ANOVA and multinomial logistic regression analysis to preference of habitat type and genome size. Using the phylogeny, we conducted ancestral character reconstruction for habitat types and tested whether genome size changes simultaneously with shifts in habitat preference. We found that 2C values had weak phylogenetic signal, whereas the base number of chromosomes (x) had a strong phylogenetic signal. Furthermore, our analyses revealed a positive correlation between genome size and chromosome traits, indicating that the base number of chromosomes (x), chromosome size, and polyploidization may be primary contributors to genome expansion in ferns and lycophytes. Genome sizes in different habitat types varied significantly and were significantly correlated with habitat types; specifically, multinomial logistic regression indicated that species with larger 2C values were more likely to be epiphytes. Terrestrial habitat is inferred to be ancestral for both extant ferns and lycophytes, whereas transitions to other habitat types occurred as the major clades emerged. Shifts in habitat types appear be followed by periods of genomic stability. Based on these results, we inferred that habitat type changes and multiple whole-genome duplications have contributed to the formation of large genomes of ferns and their allies during their evolutionary history.

The complete chloroplast genome of Plateau herb Chesneya acaulis (Fabaceae).

Chesneya acaulis is a perennial herb, which restricts in Xizang (Tibet) of China, Afghanistan, and Pakistan. The complete chloroplast genome was sequenced using the Illumina Hiseq X-Ten platform. The genome lacks an inverted repeat (IR) region, containing 75 protein-coding genes, 29 tRNAs genes, and 4 rRNAs. The overall GC content is 34.6%. A phylogenetic tree based on the whole chloroplast genomes of 15 species indicated that C. acaulis had a close relationship with the genus Hedysarum, and it nested in the inverted repeat-lacking clade (IRLC), of the subfamily Papilionoideae (Leguminosae).