Intron Content Research Articles

Codon usage pattern of the ancestor of green plants revealed through Rhodophyta

Rhodophyta are among the closest known relatives of green plants. Studying the codons of their genomes can help us understand the codon usage pattern and characteristics of the ancestor of green plants. By studying the codon usage pattern of all available red algae, it was found that although there are some differences among species, high-bias genes in most red algae prefer codons ending with GC. Correlation analysis, Nc-GC3s plots, parity rule 2 plots, neutrality plot analysis, differential protein region analysis and comparison of the nucleotide content of introns and flanking sequences showed that the bias phenomenon is likely to be influenced by local mutation pressure and natural selection, the latter of which is the dominant factor in terms of translation accuracy and efficiency. It is worth noting that selection on translation accuracy could even be detected in the low-bias genes of individual species. In addition, we identified 15 common optimal codons in seven red algae except for G. sulphuraria for the first time, most of which were found to be complementary and bound to the tRNA genes with the highest copy number. Interestingly, tRNA modification was found for the highly degenerate amino acids of all multicellular red algae and individual unicellular red algae, which indicates that highly biased genes tend to use modified tRNA in translation. Our research not only lays a foundation for exploring the characteristics of codon usage of the red algae as green plant ancestors, but will also facilitate the design and performance of transgenic work in some economic red algae in the future.

(2955) Proposal to conserve the name Didymium against Mucilago and Spumaria (Physarales, Myxomycetes)

(2955) Proposal to conserve the name <i>Didymium</i> against <i>Mucilago</i> and <i>Spumaria</i> (<i>Physarales</i>, <i>Myxomycetes</i>)

Rate accelerations in plastid and mitochondrial genomes of Cyperaceae occur in the same clades

Cyperaceae, the second largest family in the monocot order Poales, comprises >5500 species and includes the genus Eleocharis with ∼ 250 species. A previous study of complete plastomes of two Eleocharis species documented extensive structural heteroplasmy, gene order changes, high frequency of dispersed repeats along with gene losses and duplications. To better understand the phylogenetic distribution of gene and intron content as well as rates and patterns of sequence evolution within and between mitochondrial and plastid genomes of Eleocharis and Cyperaceae, an additional 29 Eleocharis organelle genomes were sequenced and analyzed. Eleocharis experienced extensive gene loss in both genomes while loss of introns was mitochondria-specific. Eleocharis has higher rates of synonymous (dS) and nonsynonymous (dN) substitutions in the plastid and mitochondrion than most sampled angiosperms, and the pattern was distinct from other eudicot lineages with accelerated rates. Several clades showed higher dS and dN in mitochondrial genes than in plastid genes. Furthermore, nucleotide substitution rates of mitochondrial genes were significantly accelerated on the branch leading to Cyperaceae compared to most angiosperms. Mitochondrial genes of Cyperaceae exhibited dramatic loss of RNA editing sites and a negative correlation between RNA editing and dS values was detected among angiosperms. Mutagenic retroprocessing and dysfunction of DNA replication, repair and recombination genes are the most likely cause of striking rate accelerations and loss of edit sites and introns in Eleocharis and Cyperaceae organelle genomes.

Transposable elements drive intron gain in diverse eukaryotes

There is massive variation in intron numbers across eukaryotic genomes, yet the major drivers of intron content during evolution remain elusive. Rapid intron loss and gain in some lineages contrast with long-term evolutionary stasis in others. Episodic intron gain could be explained by recently discovered specialized transposons called Introners, but so far Introners are only known from a handful of species. Here, we performed a systematic search across 3,325 eukaryotic genomes and identified 27,563 Introner-derived introns in 175 genomes (5.2%). Species with Introners span remarkable phylogenetic diversity, from animals to basal protists, representing lineages whose last common ancestor dates to over 1.7 billion years ago. Aquatic organisms were 6.5 times more likely to contain Introners than terrestrial organisms. Introners exhibit mechanistic diversity but most are consistent with DNA transposition, indicating that Introners have evolved convergently hundreds of times from nonautonomous transposable elements. Transposable elements and aquatic taxa are associated with high rates of horizontal gene transfer, suggesting that this combination of factors may explain the punctuated and biased diversity of species containing Introners. More generally, our data suggest that Introners may explain the episodic nature of intron gain across the eukaryotic tree of life. These results illuminate the major source of ongoing intron creation in eukaryotic genomes.

Gene loss, genome rearrangement, and accelerated substitution rates in plastid genome of Hypericum ascyron (Hypericaceae)

BackgroundComparative genomic analysis exhibits dynamic evolution of plastid genome (plastome) in the clusioid clade of Malpighiales, which comprise five families, including multiple inversions and gene losses. Little is known about the plastome evolution in Hypericaceae, a large family in the clade. Only the plastome of one species, Cratoxylum cochinchinense, has been published.ResultsWe generated a complete plastome sequence for Hypericum ascyron, providing the first complete plastome from the tribe Hypericeae (Hypericaceae). The H. ascyron plastome exhibits dynamic changes in gene and intron content, structure, and sequence divergence compared to the C. cochinchinense plastome from the tribe Cratoxyleae (Hypericaceae). Transcriptome data determined the evolutionary fate of the missing plastid genes infA, rps7, rps16, rpl23, and rpl32 in H. ascyron. Putative functional transfers of infA, rps7, and rpl32 were detected to the nucleus, whereas rps16 and rpl23 were substituted by nuclear-encoded homologs. The plastid rpl32 was integrated into the nuclear-encoded SODcp gene. Our findings suggested that the transferred rpl32 had undergone subfunctionalization by duplication rather than alternative splicing. The H. ascyron plastome rearrangements involved seven inversions, at least three inverted repeat (IR) boundary shifts, which generated gene relocations and duplications. Accelerated substitution rates of plastid genes were observed in the H. ascyron plastome compared with that of C. cochinchinense plastid genes. The higher substitution rates in the accD and clpP were correlated with structural change, including a large insertion of amino acids and losses of two introns, respectively. In addition, we found evidence of positive selection of the clpP, matK, and rps3 genes in the three branches related to H. ascyron. In particular, the matK gene was repeatedly under selection within the family Hypericaceae. Selective pressure in the H. ascyron matK gene was associated with the loss of trnK-UUU and relocation into the IR region.ConclusionsThe Hypericum ascyron plastome sequence provides valuable information for improving the understanding of plastome evolution among the clusioid of the Malpighiales. Evidence for intracellular gene transfer from the plastid to the nucleus was detected in the nuclear transcriptome, providing insight into the evolutionary fate of plastid genes in Hypericaceae.

GC content, but not nucleosome positioning, directly contributes to intron splicing efficiency in Paramecium.

Eukaryotic genes are interrupted by introns that must be accurately spliced from mRNA precursors. With an average length of 25 nt, the more than 90,000 introns of Paramecium tetraurelia stand among the shortest introns reported in eukaryotes. The mechanisms specifying the correct recognition of these tiny introns remain poorly understood. Splicing can occur cotranscriptionally, and it has been proposed that chromatin structure might influence splice site recognition. To investigate the roles of nucleosome positioning in intron recognition, we determined the nucleosome occupancy along the P. tetraurelia genome. We show that P. tetraurelia displays a regular nucleosome array with a nucleosome repeat length of ∼151 bp, among the smallest periodicities reported. Our analysis has revealed that introns are frequently associated with inter-nucleosomal DNA, pointing to an evolutionary constraint favoring introns at the AT-rich nucleosome edge sequences. Using accurate splicing efficiency data from cells depleted for nonsense-mediated decay effectors, we show that introns located at the edge of nucleosomes display higher splicing efficiency than those at the center. However, multiple regression analysis indicates that the low GC content of introns, rather than nucleosome positioning, is associated with high splicing efficiency. Our data reveal a complex link between GC content, nucleosome positioning, and intron evolution in Paramecium.

The complete mitochondrial genome of Cycas debaoensis revealed unexpected static evolution in gymnosperm species.

Mitochondrial genomes of vascular plants are well known for their liability in architecture evolution. However, the evolutionary features of mitogenomes at intra-generic level are seldom studied in vascular plants, especially among gymnosperms. Here we present the complete mitogenome of Cycas debaoensis, an endemic cycad species to the Guangxi region in southern China. In addition to assemblage of draft mitochondrial genome, we test the conservation of gene content and mitogenomic stability by comparing it to the previously published mitogenome of Cycas taitungensis. Furthermore, we explored the factors such as structural rearrangements and nuclear surveillance of double-strand break repair (DSBR) proteins in Cycas in comparison to other vascular plant groups. The C. debaoensis mitogenome is 413,715 bp in size and encodes 69 unique genes, including 40 protein coding genes, 26 tRNAs, and 3 rRNA genes, similar to that of C. taitungensis. Cycas mitogenomes maintained the ancestral intron content of seed plants (26 introns), which is reduced in other lineages of gymnosperms, such as Ginkgo biloba, Taxus cuspidata and Welwitschia mirabilis due to selective pressure or retroprocessing events. C. debaoensis mitogenome holds 1,569 repeated sequences (> 50 bp), which partially account for fairly large intron size (1200 bp in average) of Cycas mitogenome. The comparison of RNA-editing sites revealed 267 shared non-silent editing site among predicted vs. empirically observed editing events. Another 33 silent editing sites from empirical data increase the total number of editing sites in Cycas debaoensis mitochondrial protein coding genes to 300. Our study revealed unexpected conserved evolution between the two Cycas species. Furthermore, we found strict collinearity of the gene order along with the identical set of genomic content in Cycas mt genomes. The stability of Cycas mt genomes is surprising despite the existence of large number of repeats. This structural stability may be related to the relative expansion of three DSBR protein families (i.e., RecA, OSB, and RecG) in Cycas nuclear genome, which inhibit the homologous recombinations, by monitoring the accuracy of mitochondrial chromosome repair.

Evidence for a force favoring GC over AT at short intronic sites in Drosophila simulans and Drosophila melanogaster.

Population genetics studies often make use of a class of nucleotide site free from selective pressures, in order to make inferences about population size changes or natural selection at other sites. If such neutral sites can be identified, they offer the opportunity to avoid any confounding effects of selection. Here, we investigate evolution at putatively neutrally evolving short intronic sites in natural populations of Drosophila melanogaster and Drosophila simulans, in order to understand the properties of spontaneous mutations and the extent of GC-biased gene conversion in these species. Use of data on the genetics of natural populations is advantageous because it integrates information from large numbers of individuals over long timescales. In agreement with direct evidence from observations of spontaneous mutations in Drosophila, we find a bias in the spectrum of mutations toward AT basepairs. In addition, we find that this bias is stronger in the D. melanogaster lineage than in the D. simulans lineage. The evidence for GC-biased gene conversion in Drosophila has been equivocal. Here, we provide evidence for a weak force favoring GC in both species, which is correlated with the GC content of introns and is stronger in D. simulans than in D. melanogaster.

The mitochondrial genomes of bryophytes

In contrast to the highly variable mitogenomes of vascular plants, the composition and architecture of mitogenomes within the three bryophyte lineages appear stable and invariant. Currently, complete mitogenomes are available from 113 bryophyte accessions of 71 genera and 28 orders. Liverworts and mosses hold a rich mitochondrial (mt) gene repertoire among land plants with 40–42 protein-coding genes, whereas hornworts maintain the smallest functional gene set among land plants, of only around two dozen protein-coding genes, with the majority of ribosomal genes pseudogenized and all cytochrome c maturase genes lost. The rRNA and tRNA genes are also conserved and rich in mosses and liverworts, whereas subject to patchy losses in hornworts. In contrast to the conserved gene set, intron content varies significantly with only one intron shared among the three bryophyte lineages. Bryophytes hold relatively compact mitogenomes with narrow size fluctuations. Among the three bryophyte lineages, intergenic spacers and repeat content are smallest in mosses, largest in hornworts, and intermediate in liverworts, mirroring their size differences and levels of structural dynamics among the three lineages. Mosses, with the least repeated sequences, show the most static genome structure; whereas hornworts, with a relatively large set of repeated sequences, experience 1–4 rearrangements; liverworts, with intermediate repeat levels, see only one structural variant that requires two inversions to gain collinearity with the mitogenome of other liverworts. Repeat sequences were evoked to explain the mt gene order rearrangements in hornwort and liverwort mitogenomes; with the latter also supported by sequencing read evidence, which suggests that the conserved mitogenome structure observed in bryophyte lineages might be shaped by low repeat recombination level, and/or along with the intensified nucleus’ surveillance. Mitochondrial RNA editing is abundant in hornworts, with medium frequency and high variation in liverwort species, and generally limited in mosses, reflecting the diversity of nuclear encoded PPR proteins that are functionally related to RNA editing processes.

Plastomes from tribe Plantagineae (Plantaginaceae) reveal infrageneric structural synapormorphies and localized hypermutation for Plantago and functional loss of ndh genes from Littorella

Tribe Plantagineae (Plantaginaceae) comprises ~ 270 species in three currently recognized genera (Aragoa, Littorella, Plantago), of which Plantago is most speciose. Plantago plastomes exhibit several atypical features including large inversions, expansions of the inverted repeat, increased repetitiveness, intron losses, and gene-specific increases in substitution rate, but the prevalence of these plastid features among species and subgenera is unknown. To assess phylogenetic relationships and plastomic evolutionary dynamics among Plantagineae genera and Plantago subgenera, we generated 25 complete plastome sequences and compared them with existing plastome sequences from Plantaginaceae. Using whole plastome and partitioned alignments, our phylogenomic analyses provided strong support for relationships among major Plantagineae lineages. General plastid features—including size, GC content, intron content, and indels—provided additional support that reinforced major Plantagineae subdivisions. Plastomes from Plantago subgenera Plantago and Coronopus have synapomorphic expansions and inversions affecting the size and gene order of the inverted repeats, and particular genes near the inversion breakpoints exhibit accelerated nucleotide substitution rates, suggesting localized hypermutation associated with rearrangements. The Littorella plastome lacks functional copies of ndh genes, which may be related to an amphibious lifestyle and partial reliance on CAM photosynthesis.

Multiple Intraspecific Variations of Mitochondrial Genomes in the Green-Tide Forming Alga, Ulva compressa Linnaeus (Ulvophyceae, Chlorophyta)

To gain further insights into intraspecific evolution of Ulva mitochondrial genomes, the mitogenomes of three morphotypes of the green-tide forming alga, Ulva compressa Linnaeus from China and USA, were sequenced and compared with the available data from Ulvophyceae. The U. compressa mitogenomes displayed significantly variable size at intraspecific level ranging from 61,700 to &gt;66,587 bp, due to difference in acquisition of foreign DNA fragments, gain or loss of both group I and II introns, and non-coding intergenic spacer regions. The U. compressa mitogenomes harbored variable gene content ranging from 69 genes (including orfs) in Uco1 to 76 in Uco5, and contained different intron content from 4 introns in Uco3 and Uco4 to 7 in Uco1. A total of 63 genes and only two group IB introns (intron cox1-1107 and cox1-1125) were shared by these five mitogenomes. The U. compressa mitogenomes accumulated much more inverted repeat (IR) elements, ranging from 45 in Uco1 to 88 in Uco2, than that of the other Ulva species (3 - 34). A locally collinear block of eight genes (rps11-rps19-rps4-rpl16-trnR-trnQ-trnE-trnS) with the size of 3,631 bp has been inverted only in Uco1 indicating that the rearrangement event happened after its divergence from Uco2-5 and might be relate to a specific IR element. A majority of the common genes (≈76%) displayed high identity (&gt;98%) among these five mitogenomes, while some low values were observed in six genes mainly due to duplication and insertion/deletion mutations of small DNA fragments. Our study presented the first case of multiple intraspecific variations in ulvophycean mitogenomes, and indicated that the mitogenome will be a valuable tool for understanding the native or nonindigenous nature of the cosmopolitan Ulva species.

Patterns and Rates of Plastid rps12 Gene Evolution Inferred in a Phylogenetic Context using Plastomic Data of Ferns

The trans-splicing rps12 gene of fern plastomes (plastid genomes) exhibits a unique structure owing to its variations in intragenic exon location and intron content, and thus, it provides an excellent model system for examining the effect of plastid gene structure on rates and patterns of molecular evolution. In this study, 16 complete fern plastome sequences were newly generated via the Illumina HiSeq sequencing platform. We reconstructed the phylogeny of ferns and inferred the patterns and rates of plastid rps12 gene evolution in a phylogenetic context by combining these plastome data with those of previously published fern species. We uncovered the diversity of fern plastome evolution by characterizing the structures of these genomes and obtained a highly supported phylogenetic framework for ferns. Furthermore, our results revealed molecular evolutionary patterns that were completely different from the patterns revealed in previous studies. There were significant differences in the patterns and rates of nucleotide substitutions in both intron-containing and intron-less rps12 alleles. Rate heterogeneity between single-copy (SC) and inverted repeat (IR) exons was evident. Unexpectedly, however, IR exons exhibited significantly higher synonymous substitution rates (dS) than SC exons, a pattern that contrasts the regional effect responsible for decreased rates of nucleotide substitutions in IRs. Our results reveal that structural changes in plastid genes have important effects on evolutionary rates, and we propose possible mechanisms to explain the variations in the nucleotide substitution rates of this unusual gene.

Distinct effects on mRNA export factor GANP underlie neurological disease phenotypes and alter gene expression depending on intron content.

Defects in the mRNA export scaffold protein GANP, encoded by the MCM3AP gene, cause autosomal recessive early-onset peripheral neuropathy with or without intellectual disability. We extend here the phenotypic range associated with MCM3AP variants, by describing a severely hypotonic child and a sibling pair with a progressive encephalopathic syndrome. In addition, our analysis of skin fibroblasts from affected individuals from seven unrelated families indicates that disease variants result in depletion of GANP except when they alter critical residues in the Sac3 mRNA binding domain. GANP depletion was associated with more severe phenotypes compared with the Sac3 variants. Patient fibroblasts showed transcriptome alterations that suggested intron content-dependent regulation of gene expression. For example, all differentially expressed intronless genes were downregulated, including ATXN7L3B, which couples mRNA export to transcription activation by association with the TREX-2 and SAGA complexes. Our results provide insight into the molecular basis behind genotype-phenotype correlations in MCM3AP-associated disease and suggest mechanisms by which GANP defects might alter RNA metabolism.

The complete mitogenomes of the green algae Jenufa minuta and Jenufa perforata (Chlorophyceae, incertae sedis) reveal a variant genetic code previously unrecognized in the Chlorophyceae

We sampled the first mitogenomes from a novel lineage of the Chlorophyceae that is sister to Sphaeropleales. Although the Jenufa perforata (27,180 bp) and Jenufa minuta (41,488 bp) mitogenomes differ markedly in size, gene order and intron content, their repertoires of conserved genes are essentially identical. Their rRNA genes display the same fragmentation patterns as their counterparts in Sphaeropleales. Their tRNA gene complements include two ancestral genes shared with Oedogoniales and Chaetophorales (Chlorophyceae) as well as a derived gene for a tRNA species decoding UGA as tryptophan. This noncanonical genetic code differs from those previously reported for Sphaeropleales.

Pseudoexons of the DMD Gene.

The DMD gene is the largest in the human genome, with a total intron content exceeding 2.2Mb. In the decades since DMD was discovered there have been numerous reported cases of pseudoexons (PEs) arising in the mature DMD transcripts of some individuals, either as the result of mutations or as low-frequency errors of the spliceosome. In this review, I collate from the literature 58 examples of DMD PEs and examine the diversity and commonalities of their features. In particular, I note the high frequency of PEs that arise from deep intronic SNVs and discuss a possible link between PEs induced by distal mutations and the regulation of recursive splicing.

Differential Gene Expression between Fungal Mating Types Is Associated with Sequence Degeneration.

Degenerative mutations in non-recombining regions, such as in sex chromosomes, may lead to differential expression between alleles if mutations occur stochastically in one or the other allele. Reduced allelic expression due to degeneration has indeed been suggested to occur in various sex-chromosome systems. However, whether an association occurs between specific signatures of degeneration and differential expression between alleles has not been extensively tested, and sexual antagonism can also cause differential expression on sex chromosomes. The anther-smut fungus Microbotryum lychnidis-dioicae is ideal for testing associations between specific degenerative signatures and differential expression because 1) there are multiple evolutionary strata on the mating-type chromosomes, reflecting successive recombination suppression linked to mating-type loci; 2) separate haploid cultures of opposite mating types help identify differential expression between alleles; and 3) there is no sexual antagonism as a confounding factor accounting for differential expression. We found that differentially expressed genes were enriched in the four oldest evolutionary strata compared with other genomic compartments, and that, within compartments, several signatures of sequence degeneration were greater for differentially expressed than non-differentially expressed genes. Two particular degenerative signatures were significantly associated with lower expression levels within differentially expressed allele pairs: upstream insertion of transposable elements and mutations truncating the protein length. Other degenerative mutations associated with differential expression included nonsynonymous substitutions and altered intron or GC content. The association between differential expression and allele degeneration is relevant for a broad range of taxa where mating compatibility or sex is determined by genes located in large regions where recombination is suppressed.

The complete mitochondrial genome of Taxus cuspidata (Taxaceae): eight protein-coding genes have transferred to the nuclear genome

BackgroundGymnosperms represent five of the six lineages of seed plants. However, most sequenced plant mitochondrial genomes (mitogenomes) have been generated for angiosperms, whereas mitogenomic sequences have been generated for only six gymnosperms. In particular, complete mitogenomes are available for all major seed plant lineages except Conifer II (non-Pinaceae conifers or Cupressophyta), an important lineage including six families, which impedes a comprehensive understanding of the mitogenomic diversity and evolution in gymnosperms.ResultsHere, we report the complete mitogenome of Taxus cuspidata in Conifer II. In comparison with previously released gymnosperm mitogenomes, we found that the mitogenomes of Taxus and Welwitschia have lost many genes individually, whereas all genes were identified in the mitogenomes of Cycas, Ginkgo and Pinaceae. Multiple tRNA genes and introns also have been lost in some lineages of gymnosperms, similar to the pattern observed in angiosperms. In general, gene clusters could be less conserved in gymnosperms than in angiosperms. Moreover, fewer RNA editing sites were identified in the Taxus and Welwitschia mitogenomes than in other mitogenomes, which could be correlated with fewer introns and frequent gene losses in these two species.ConclusionsWe have sequenced the Taxus cuspidata mitogenome, and compared it with mitogenomes from the other four gymnosperm lineages. The results revealed the diversity in size, structure, gene and intron contents, foreign sequences, and mutation rates of gymnosperm mitogenomes, which are different from angiosperm mitogenomes.

Mitochondrial genomes of the early land plant lineage liverworts (Marchantiophyta): conserved genome structure, and ongoing low frequency recombination

BackgroundIn contrast to the highly labile mitochondrial (mt) genomes of vascular plants, the architecture and composition of mt genomes within the main lineages of bryophytes appear stable and invariant. The available mt genomes of 18 liverwort accessions representing nine genera and five orders are syntenous except for Gymnomitrion concinnatum whose genome is characterized by two rearrangements. Here, we expanded the number of assembled liverwort mt genomes to 47, broadening the sampling to 31 genera and 10 orders spanning much of the phylogenetic breadth of liverworts to further test whether the evolution of the liverwort mitogenome is overall static.ResultsLiverwort mt genomes range in size from 147 Kb in Jungermanniales (clade B) to 185 Kb in Marchantiopsida, mainly due to the size variation of intergenic spacers and number of introns. All newly assembled liverwort mt genomes hold a conserved set of genes, but vary considerably in their intron content. The loss of introns in liverwort mt genomes might be explained by localized retroprocessing events. Liverwort mt genomes are strictly syntenous in genome structure with no structural variant detected in our newly assembled mt genomes. However, by screening the paired-end reads, we do find rare cases of recombination, which means multiple concurrent genome structures may exist in the vegetative tissues of liverworts. Our phylogenetic analyses of the nuclear encoded double stand break repair protein families revealed liverwort-specific subfamilies expansions.ConclusionsThe low repeat recombination level, selection, along with the intensified nuclear surveillance, might together shape the structural evolution of liverwort mt genomes.

Plastid genome evolution in tribe Desmodieae (Fabaceae: Papilionoideae).

Recent plastid genome (plastome) studies of legumes (family Fabaceae) have shown that this family has undergone multiple atypical plastome evolutions from each of the major clades. The tribe Desmodieae belongs to the Phaseoloids, an important but systematically puzzling clade within Fabaceae. In this study, we investigated the plastome evolution of Desmodieae and analyzed its phylogenetic signaling. We sequenced six complete plastomes from representative members of Desmodieae and from its putative sister Phaseoloid genus Mucuna. Those genomes contain 128 genes and range in size from 148,450 to 153,826 bp. Analyses of gene and intron content revealed similar characters among the members of Desmodieae and Mucuna. However, there were also several distinct characters identified. The loss of the rpl2 intron was a feature shared between Desmodieae and Mucuna, whereas the loss of the rps12 intron was specific to Desmodieae. Likewise, gene loss of rps16 was observed in Mucuna but not in Desmodieae. Substantial sequence variation of ycf4 was detected from all the sequenced plastomes, but pseudogenization was restricted to the genus Desmodium. Comparative analysis of gene order revealed a distinct plastome conformation of Desmodieae compared with other Phaseoloid legumes, i.e., an inversion of an approximately 1.5-kb gene cluster (trnD-GUC, trnY-GUA, and trnE-UUC). The inversion breakpoint suggests that this event was mediated by the recombination of an 11-bp repeat motif. A phylogenetic analysis based on the plastome-scale data set found the tribe Desmodieae is a highly supported monophyletic group nested within the paraphyletic Phaseoleae, as has been found in previous phylogenetic studies. Two subtribes (Desmodiinae and Lespedezinae) of Desmodieae were also supported as monophyletic groups. Within the subtribe Lespedezinae, Lespedeza is closer to Kummerowia than Campylotropis.

Single-cell correlations of mRNA and protein content in a human monocytic cell line after LPS stimulation.

The heterogeneity of mRNA and protein expression at the single-cell level can reveal fundamental information about cellular response to external stimuli, including the sensitivity, timing, and regulatory interactions of genes. Here we describe a fully automated system to digitally count the intron, mRNA, and protein content of up to five genes of interest simultaneously in single-cells. Full system automation of 3D microscope scans and custom image analysis routines allows hundreds of individual cells to be automatically segmented and the mRNA-protein content to be digitally counted. Single-molecule intron and mRNA content is measured by single-molecule fluorescence in-situ hybridization (smFISH), while protein content is quantified though the use of antibody probes. To mimic immune response to bacterial infection, human monocytic leukemia cells (THP-1) were stimulated with lipopolysaccharide (LPS), and the expression of two inflammatory genes, IL1β (interleukin 1β) and TNF-α (tumor necrosis factor α), were simultaneously quantified by monitoring the intron, mRNA, and protein levels over time. The simultaneous labeling of cellular content allowed for a series of correlations at the single-cell level to be explored, both in the progressive maturation of a single gene (intron-mRNA-protein) and comparative analysis between the two immune response genes. In the absence of LPS stimulation, mRNA expression of IL1β and TNF-α were uncorrelated. Following LPS stimulation, mRNA expression of the two genes became more correlated, consistent with a model in which IL1β and TNF-α upregulation occurs in parallel through independent mechanistic pathways. This smFISH methodology can be applied to different complex biological systems to provide valuable insight into highly dynamic gene mechanisms that determine cell plasticity and heterogeneity of cellular response.