Evolutionary Dynamics of the Leucine-Rich Repeat Receptor-Like Kinase (LRR-RLK) Subfamily in Angiosperms

Gene duplications are an important factor in plant evolution, and lineage-specific expanded (LSE) genes are of particular interest. Receptor-like kinases expanded massively in land plants, and leucine-rich repeat receptor-like kinases (LRR-RLK) constitute the largest receptor-like kinases family. Based on the phylogeny of 7,554 LRR-RLK genes from 31 fully sequenced flowering plant genomes, the complex evolutionary dynamics of this family was characterized in depth. We studied the involvement of selection during the expansion of this family among angiosperms. LRR-RLK subgroups harbor extremely contrasting rates of duplication, retention, or loss, and LSE copies are predominantly found in subgroups involved in environmental interactions. Expansion rates also differ significantly depending on the time when rounds of expansion or loss occurred on the angiosperm phylogenetic tree. Finally, using a d_N/d_S-based test in a phylogenetic framework, we searched for selection footprints on LSE and single-copy LRR-RLK genes. Selective constraint appeared to be globally relaxed at LSE genes, and codons under positive selection were detected in 50% of them. Moreover, the leucine-rich repeat domains, and specifically four amino acids in them, were found to be the main targets of positive selection. Here, we provide an extensive overview of the expansion and evolution of this very large gene family.Receptor-like kinases (RLKs) constitute one of the largest gene families in plants and expanded massively in land plants (Embryophyta; Lehti-Shiu et al., 2009, 2012). For plant RLK gene families, the functions of most members are often not known (especially in recently expanded families), but some described functions include innate immunity (Albert et al., 2010), pathogen response (Dodds and Rathjen, 2010), abiotic stress (Yang et al., 2010), development (De Smet et al., 2009), and sometimes multiple functions (Lehti-Shiu et al., 2012). The RLKs usually consist of three domains: an N-terminal extracellular domain, a transmembrane domain, and a C-terminal kinase domain (KD). In plants, the KD usually has a Ser/Thr specificity (Shiu and Bleecker, 2001), but Tyr-specific RLKs were also described (e.g. BRASSINOSTEROID INSENSITIVE1; Oh et al., 2009). Interestingly, it was estimated that approximately 20% of RLKs contain a catalytically inactive KD (e.g. STRUBBELIG and CORYNE; Chevalier et al., 2005; Castells and Casacuberta, 2007; Gish and Clark, 2011). In Arabidopsis (Arabidopsis thaliana), 44 RLK subgroups (SGs) were defined by inferring the phylogenetic relationships between the KDs (Shiu and Bleecker, 2001). Interestingly, different SGs show different duplication/retention rates (Lehti-Shiu et al., 2009). Specifically, RLKs involved in stress responses show a high number of tandemly duplicated genes whereas those involved in development do not (Shiu et al., 2004), which suggests that some RLK genes are important for the responses of land plants to a changing environment (Lehti-Shiu et al., 2012). There seem to be relatively few RLK pseudogenes compared with other large gene families, and copy retention was argued to be driven by both drift and selection (Zou et al., 2009; Lehti-Shiu et al., 2012). As most SGs are relatively old and RLK subfamilies expanded independently in several plant lineages, duplicate retention cannot be explained by drift alone, and natural selection is expected to be an important driving factor in RLK gene family retention (Lehti-Shiu et al., 2009).Leucine-rich repeat-receptor-like kinases (LRR-RLKs), which contain up to 30 leucine-rich repeat (LRRs) in their extracellular domain, constitute the largest RLK family (Shiu and Bleecker, 2001). Based on the KD, 15 LRR-RLK SGs have been established in Arabidopsis (Shiu et al., 2004; Lehti-Shiu et al., 2009). So far, two major functions have been attributed to them: defense against pathogens and development (Tang et al., 2010b). LRR-RLKs involved in defense are predominantly found in lineage-specific expanded (LSE) gene clusters, whereas LRR-RLKs involved in development are mostly found in nonexpanded groups (Tang et al., 2010b). It was also discovered that the LRR domains are significantly less conserved than the remaining domains of the LRR-RLK genes (Tang et al., 2010b). In addition, a study of four plant genomes (Arabidopsis, grape [Vitis vinifera], poplar [Populus trichocarpa], and rice [Oryza sativa]) showed that LRR-RLK genes from LSE gene clusters show significantly more indications of positive selection or relaxed constraint than LRR-RLKs from nonexpanded groups (Tang et al., 2010b).The genomes of flowering plants (angiosperms) have been shown to be highly dynamic compared with most other groups of land plants (Leitch and Leitch, 2012). This dynamic is mostly caused by the frequent multiplication of genetic material, followed by a complex pattern of differential losses (i.e. the fragmentation process) and chromosomal rearrangements (Langham et al., 2004; Leitch and Leitch, 2012). Most angiosperm genomes sequenced so far show evidence for at least one whole-genome multiplication event during their evolution (Jaillon et al., 2007; D’Hont et al., 2012; Tomato Genome Consortium, 2012). At a smaller scale, tandem and segmental duplications are also very common in angiosperms (Arabidopsis Genome Initiative, 2000; International Rice Genome Sequencing Project, 2005; Rizzon et al., 2006). Although the most common fate of duplicated genes is to be progressively lost, in some cases they can be retained in the genome, and adaptive as well as nonadaptive scenarios have been discussed to play a role in this preservation process (for review, see Moore and Purugganan, 2005; Hahn, 2009; Innan, 2009; Innan and Kondrashov, 2010). Whole-genome sequences also revealed that the same gene may undergo several rounds of duplication and retention. These LSE genes were shown to evolve under positive selection more frequently than single-copy genes in angiosperms (Fischer et al., 2014). That study analyzed general trends over whole genomes. Here, we ask if, and to what extent, this trend is observable at LRR-RLK genes. As this gene family is very dynamic and large, and in accordance with the results of Tang et al. (2010b), we expect the effect of positive selection to be even more pronounced than in the whole-genome average.We analyzed 33 Embryophyta genomes to investigate the evolutionary history of the LRR-RLK gene family in a phylogenetic framework. Twenty LRR-RLK SGs were identified, and from this data set, we deciphered the evolutionary dynamics of this family within angiosperms. The expansion/reduction rates were contrasted between SGs and species as well as in ancestral branches of the angiosperm phylogeny. We then focused on genes whose number increased dramatically in an SG- and/or species-specific manner (i.e. LSE genes). Those genes are likely to be involved in species-specific cellular processes or adaptive interactions and were used as a template to infer the potential occurrence of positive selection. This led to the identification of sites at which positive selection likely acted. We discuss our results in the light of angiosperm genome evolution and current knowledge of LRR-RLK functions. Positive selection footprints identified in LSE genes highlight the importance of combining evolutionary analysis and functional knowledge to guide further investigations.     

A biofilm model to understand the onset of sulfate reduction in denitrifying membrane biofilm reactors

This work presents a multispecies biofilm model that describes the co‐existence of nitrate‐ and sulfate‐reducing bacteria in the H₂‐based membrane biofilm reactor (MBfR). The new model adapts the framework of a biofilm model for simultaneous nitrate and perchlorate removal by considering the unique metabolic and physiological characteristics of autotrophic sulfate‐reducing bacteria that use H₂ as their electron donor. To evaluate the model, the simulated effluent H₂, UAP (substrate‐utilization‐associated products), and BAP (biomass‐associated products) concentrations are compared to experimental results, and the simulated biomass distributions are compared to real‐time quantitative polymerase chain reaction (qPCR) data in the experiments for parameter optimization. Model outputs and experimental results match for all major trends and explain when sulfate reduction does or does not occur in parallel with denitrification. The onset of sulfate reduction occurs only when the nitrate concentration at the fiber's outer surface is low enough so that the growth rate of the denitrifying bacteria is equal to that of the sulfate‐reducing bacteria. An example shows how to use the model to design an MBfR that achieves satisfactory nitrate reduction, but suppresses sulfate reduction. Biotechnol. Bioeng. 2013; 110: 763–772. © 2012 Wiley Periodicals, Inc.     

Balancing redox cofactor generation and ATP synthesis: Key microaerobic responses in thermophilic fermentations

Geobacillus thermoglucosidasius is a Gram‐positive, thermophilic bacterium capable of ethanologenic fermentation of both C5 and C6 sugars and may have possible use for commercial bioethanol production [Tang et al., 2009; Taylor et al. (2009) Trends Biotechnol 27(7): 398–405]. Little is known about the physiological changes that accompany a switch from aerobic (high redox) to microaerobic/fermentative (low redox) conditions in thermophilic organisms. The changes in the central metabolic pathways in response to a switch in redox potential were analyzed using quantitative real‐time PCR and proteomics. During low redox (fermentative) states, results indicated that glycolysis was uniformly up‐regulated, the Krebs (tricarboxylic acid or TCA) cycle non‐uniformly down‐regulated and that there was little to no change in the pentose phosphate pathway. Acetate accumulation was accounted for by strong down‐regulation of the acetate CoA ligase gene (acs) in addition to up‐regulation of the pta and ackA genes (involved in acetate production), thus conserving ATP while reducing flux through the TCA cycle. Substitution of an NADH dehydrogenase (down‐regulated) by an up‐regulated NADH:FAD oxidoreductase and up‐regulation of an ATP synthase subunit, alongside the observed shifts in the TCA cycle, suggested that an oxygen‐scavenging electron transport chain likely remained active during low redox conditions. Together with the observed up‐regulation of a glyoxalase and down‐regulation of superoxide dismutase, thought to provide protection against the accumulation of toxic phosphorylated glycolytic intermediates and reactive oxygen species, respectively, the changes observed in G. thermoglucosidasius NCIMB 11955 under conditions of aerobic‐to‐microaerobic switching were consistent with responses to low pO₂ stress. Biotechnol. Bioeng. 2013; 110: 1057–1065. © 2012 Wiley Periodicals, Inc.     

A thin‐walled polydimethylsiloxane bioreactor for high‐density hepatocyte sandwich culture

In vitro drug testing requires long‐term maintenance of hepatocyte liver specific functions. Hepatocytes cultured at a higher seeding density in a sandwich configuration exhibit an increased level of liver specific functions when compared to low density cultures due to the better cell to cell contacts that promote long term maintenance of polarity and liver specific functions. However, culturing hepatocytes at high seeding densities in a standard 24‐well plate poses problems in terms of the mass transport of nutrients and oxygen to the cells. In view of this drawback, we have developed a polydimethylsiloxane (PDMS) bioreactor that was able to maintain the long‐term liver specific functions of a hepatocyte sandwich culture at a high seeding density. The bioreactor was fabricated with PDMS, an oxygen permeable material, which allowed direct oxygenation and perfusion to take place simultaneously. The mass transport of oxygen and the level of shear stress acting on the cells were analyzed by computational fluid dynamics (CFD). The combination of both direct oxygenation and perfusion has a synergistic effect on the liver specific function of a high density hepatocyte sandwich culture over a period of 9 days. Biotechnol. Bioeng. 2013; 110: 1663–1673. © 2012 Wiley Periodicals, Inc.     

Chemical Diversity of the Natural Populations of Dalmatian Pyrethrum (Tanacetum cinerariifolium (Trevir.) Sch.Bip.) in Croatia

Dalmatian pyrethrum (Tanacetum cinerariifolium (Trevir. ) Sch.Bip. ) is a plant species endemic to the east Adriatic coast. The bioactive substance of Dalmatian pyrethrum is a natural insecticide, pyrethrin, a mixture of six active components (pyrethrins I and II, cinerins I and II, and jasmolins I and II). The insecticidal potential of pyrethrin was recognized decades ago, and dried and ground flowers have traditionally been used in Croatian agriculture and households. A total of 25 Dalmatian pyrethrum populations from Croatia were studied to determine the pyrethrin content and composition, and to identify distinct chemotypes. The total pyrethrin content ranged from 0.36 to 1.30% (dry flower weight; DW) and the pyrethrin I/pyrethrin II ratio ranged from 0.64 to 3.33%. The statistical analyses revealed that the correlations between the percentage of pyrethrin I and of all the other components were significant and negative. The total pyrethrin content was positively correlated with the percentage of pyrethrin I and negatively correlated with cinerin II. The multivariate analysis of the chemical variability enabled the identification of five chemotypes among 25 Dalmatian pyrethrum populations. The chemical characterization of indigenous Dalmatian pyrethrum populations may serve as a good background for future breeding and agricultural exploitation.     

Embryonic lethality leads to hybrid male inviability in hybrids between Drosophila melanogaster and D. santomea

The study of the morphological defects unique to interspecific hybrids can reveal which developmental pathways have diverged between species. Drosophila melanogaster and D. santomea diverged more than 10 million years ago, and when crossed produce sterile adult females. Adult hybrid males are absent from all interspecific crosses. We aimed to determine the fate of these hybrid males. To do so, we tracked the development of hybrid females and males using classic genetic markers and techniques. We found that hybrid males die predominantly as embryos with severe segment‐specification defects while a large proportion of hybrid females embryos hatch and survive to adulthood. In particular, we show that most male embryos show a characteristic abdominal ablation phenotype, not observed in either parental species. This suggests that sex‐specific embryonic developmental defects eliminate hybrid males in this interspecific cross. The study of the developmental abnormalities that occur in hybrids can lead to the understanding of cryptic molecular divergence between species sharing a conserved body plan.     

Role of evolutionary and ecological factors in the reproductive success and the spatial genetic structure of the temperate gorgonian Paramuricea clavata

Dispersal and mating features strongly influence the evolutionary dynamics and the spatial genetic structure (SGS) of marine populations. For the first time in a marine invertebrate, we examined individual reproductive success, by conducting larval paternity assignments after a natural spawning event, combined with a small‐scale SGS analysis within a population of the gorgonian Paramuricea clavata. Thirty four percent of the larvae were sired by male colonies surrounding the brooding female colonies, revealing that the bulk of the mating was accomplished by males from outside the studied area. Male success increased with male height and decreased with increasing male to female distance. The parentage analyses, with a strong level of self‐recruitment (25%), unveiled the occurrence of a complex family structure at a small spatial scale, consistent with the limited larval dispersal of this species. However, no evidence of small scale SGS was revealed despite this family structure. Furthermore, temporal genetic structure was not observed, which appears to be related to the rather large effective population size. The low level of inbreeding found suggests a pattern of random mating in this species, which disagrees with expectations that limited larval dispersal should lead to biparental inbreeding. Surface brooding and investment in sexual reproduction in P. clavata contribute to multiple paternity (on average 6.4 fathers were assigned per brood), which enhance genetic diversity of the brood. Several factors may have contributed to the lack of biparental inbreeding in our study such as (i) the lack of sperm limitation at a small scale, (ii) multiple paternity, and (iii) the large effective population size. Thus, our results indicate that limited larval dispersal and complex family structure do not necessarily lead to biparental inbreeding and SGS. In the framework of conservation purposes, our results suggested that colony size, proximity among colonies and the population size should be taken into consideration for restoration projects.