Compositional Bias Research Articles

The mitochondrial genome of Hercostomus brevipilosus (Diptera: Dolichopodidae)

The long-legged fly Hercostomus brevipilosus (Genbank accession number: MN840034) belongs to the subfamily Dolichopodinae of Dolichopodidae. The mitogenome of H. brevipilosus was sequenced, the first representative of the mitogenome of the genus Hercostomus. The mitogenome is 15,019 bp totally, consisting of 13 protein-coding genes, two rRNAs and 22 transfer RNAs. All genes have similar locations and strands with that of other published species of Dolichopodidae. The nucleotide composition biases toward A and T, which together made up 73.9％ of the entirety. Bayesian inference analysis strongly supported the monophyly of Dolichopodidae and Dolichopodinae. It suggested that Hercostomus is the sister group of Dolichopus.

Production of Intrinsically Disordered Proteins for Biophysical Studies: Tips and Tricks.

Intrinsically disordered proteins (IDPs) have no single, fixed tertiary structure, yet they take on many vital functions in biology. In recent years, considerable effort has been put into the structural characterization of their conformational ensembles, to understand the link between the transient, short- and long-range organizations of IDPs and their functions. Such biophysical studies require substantial amounts of pure protein, representing a major bottleneck in the studies of IDPs. However, the unique physicochemical properties resulting from their compositional bias may be exploited for simple yet effective purification strategies. In this chapter, we provide tips and tricks for IDP production and describe the most important analyses to carry out before bringing an IDP of interest to the laboratory. We outline four purification protocols utilizing the unique properties of IDPs as well as some commonly encountered challenges and pitfalls.

The mitochondrial genome of Medetera sp. (Diptera: Dolichopodidae)

The long-legged fly Medetera sp. (Genbank accession number: MN604695) belongs to the subfamily Medeterinae of Dolichopodidae. The mitogenome of Medetera sp. was sequenced, the first representative of the mitogenome of the subfamily. The mitogenome is 14,740 bp totally, consisting of 13 protein-coding genes, two rRNAs, and 22 transfer RNAs. All genes have similar locations and strands with that of other published species of Dolichopodidae. The nucleotide composition biases toward A and T, which together made up 67.6％ of the entirety. Bayesian inference analysis strongly supported the monophyly of Dolichopodidae. It suggested that subfamily Medeterinae is the sister group of subfamily Rhaphiinae.

Characterizing the interplay between gene nucleotide composition bias and splicing

BackgroundNucleotide composition bias plays an important role in the 1D and 3D organization of the human genome. Here, we investigate the potential interplay between nucleotide composition bias and the regulation of exon recognition during splicing.ResultsBy analyzing dozens of RNA-seq datasets, we identify two groups of splicing factors that activate either about 3200 GC-rich exons or about 4000 AT-rich exons. We show that splicing factor–dependent GC-rich exons have predicted RNA secondary structures at 5′ ss and are dependent on U1 snRNP–associated proteins. In contrast, splicing factor–dependent AT-rich exons have a large number of decoy branch points, SF1- or U2AF2-binding sites and are dependent on U2 snRNP–associated proteins. Nucleotide composition bias also influences local chromatin organization, with consequences for exon recognition during splicing. Interestingly, the GC content of exons correlates with that of their hosting genes, isochores, and topologically associated domains.ConclusionsWe propose that regional nucleotide composition bias over several dozens of kilobase pairs leaves a local footprint at the exon level and induces constraints during splicing that can be alleviated by local chromatin organization at the DNA level and recruitment of specific splicing factors at the RNA level. Therefore, nucleotide composition bias establishes a direct link between genome organization and local regulatory processes, like alternative splicing.

Human faecal collection methods demonstrate a bias in microbiome composition by cell wall structure

Clinical trial faecal collections present challenges through geographical spread and inexperienced participants. Collection techniques have been developed and tested to overcome these challenges, but previous studies investigating these techniques have demonstrated a highly variable capacity for sample preservation. Furthermore, these studies typically only examine either preservation of genetic content or metabolites, not both. This study investigated the Stool Nucleic Acid Collection and Preservation Tube (Norgen BioTek Corp) for the preservation of both microbial DNA and microbial organic acid metabolites in human faecal samples when compared to frozen samples. Twenty six healthy adult participants were instructed to collect a bowel movement, subsample into collection tubes and immediately transfer the remaining bulk to −20 °C storage. Resulting organic acid concentrations remained comparable across methods when the preservation tubes were used correctly. The 16S rRNA gene sequencing data revealed twenty significantly different bacterial genera between the two collection methods. Ten Gram-negative genera were more abundant in the collection tubes, and ten Gram-positive genera were more abundant in the fresh frozen samples. This study has illustrated that faecal collection methods bias the microbial community profile according to Gram status and this should be considered when designing studies that collect and store human faecal samples.

A Puzzling Anomaly in the 4-Mer Composition of the Giant Pandoravirus Genomes Reveals a Stringent New Evolutionary Selection Process.

Pandoraviridae is a rapidly growing family of giant viruses, all of which have been isolated using laboratory strains of Acanthamoeba The genomes of 10 distinct strains have been fully characterized, reaching up to 2.5 Mb in size. These double-stranded DNA genomes encode the largest of all known viral proteomes and are propagated in oblate virions that are among the largest ever described (1.2 μm long and 0.5 μm wide). The evolutionary origin of these atypical viruses is the object of numerous speculations. Applying the chaos game representation to the pandoravirus genome sequences, we discovered that the tetranucleotide (4-mer) "AGCT" is totally absent from the genomes of 2 strains (Pandoravirus dulcis and Pandoravirus quercus) and strongly underrepresented in others. Given the amazingly low probability of such an observation in the corresponding randomized sequences, we investigated its biological significance through a comprehensive study of the 4-mer compositions of all viral genomes. Our results indicate that AGCT was specifically eliminated during the evolution of the Pandoraviridae and that none of the previously proposed host-virus antagonistic relationships could explain this phenomenon. Unlike the three other families of giant viruses (Mimiviridae, Pithoviridae, and Molliviridae) infecting the same Acanthamoeba host, the pandoraviruses exhibit a puzzling genomic anomaly suggesting a highly specific DNA editing in response to a new kind of strong evolutionary pressure.IMPORTANCE Recent years have seen the discovery of several families of giant DNA viruses infecting the ubiquitous amoebozoa of the genus Acanthamoeba With double-stranded DNA (dsDNA) genomes reaching 2.5 Mb in length packaged in oblate particles the size of a bacterium, the pandoraviruses are currently the most complex and largest viruses known. In addition to their spectacular dimensions, the pandoraviruses encode the largest proportion of proteins without homologs in other organisms, which is thought to result from a de novo gene creation process. While using comparative genomics to investigate the evolutionary forces responsible for the emergence of such an unusual giant virus family, we discovered a unique bias in the tetranucleotide composition of the pandoravirus genomes that can result only from an undescribed evolutionary process not encountered in any other microorganism.

Nucleotide composition and synonymous codon usage of open reading frames in Norovirus GII.4 variants

Norovirus GII.4 variants, a genotype in genogroup II belonging to the genus Norovirus, is a single-strand positive sense RNA containing three open reading frames (ORF1, ORF2 and ORF3) and is the most important pathogen causing nonbacterial gastroenteritis outbreaks. By using bioinformatic softwares such as Codon W, SPSS and so on, a total of 292 strains of the viruses isolated from 1974 to 2016 were analyzed for nucleotide composition and synonymous codon usage in each ORF. The result shows that it is enriched for A over the other bases in nucleotide composition, G behind the other bases in the 3rd site of all synonymous codons in the three ORFs. The patterns of nucleotide composition and codon bias of ORF2 are similar to those of ORF3 and different from those of ORF1. There are generally UpA motif and CpG motif in the codons with the lowest proportion. Correspondence analysis indicates that the codon usage may be changing over a certain time period for ORF1 in 2006 and 2012, ORF2 in 2012, and ORF3 in 2013. ENC (effective number of codons) plot and other analyses indicate that both natural selection and mutational pressure play partly roles in the ORFs, but natural selection is more important for ORF2 and ORF3. Besides, we also found all optimal codons in the ORFs. The study provides a basic understanding of the mechanism for norovirus GII.4 codon usage bias. AbbreviationsORFOpen Reading FrameENCEffective Number of CodonsCOAcorrespondence analysisRSCURelative Synonymous Codon UsageCAICodon Adaptation IndexCBICodon Bias IndexFopfrequency of optimal codonsL_symnumber of synonymous codonsL_aalength amino acidsGRAVYgrand average of hydropathicityAromaaromaticityCommunicated by Ramaswamy H. Sarma

IDDomainSpotter: Compositional bias reveals domains in long disordered protein regions-Insights from transcription factors.

Protein domains constitute regions of distinct structural properties and molecular functions that are retained when removed from the rest of the protein. However, due to the lack of tertiary structure, the identification of domains has been largely neglected for long (>50 residues) intrinsically disordered regions. Here we present a sequence-based approach to assess and visualize domain organization in long intrinsically disordered regions based on compositional sequence biases. An online tool to find putative intrinsically disordered domains (IDDomainSpotter) in any protein sequence or sequence alignment using any particular sequence trait is available at http://www.bio.ku.dk/sbinlab/IDDomainSpotter. Using this tool, we have identified a putative domain enriched in hydrophilic and disorder-promoting residues (Pro, Ser, and Thr) and depleted in positive charges (Arg and Lys) bordering the folded DNA-binding domains of several transcription factors (p53, GCR, NAC46, MYB28, and MYB29). This domain, from two different MYB transcription factors, was characterized biophysically to determine its properties. Our analyses show the domain to be extended, dynamic and highly disordered. It connects the DNA-binding domain to other disordered domains and is present and conserved in several transcription factors from different families and domains of life. This example illustrates the potential of IDDomainSpotter to predict, from sequence alone, putative domains of functional interest in otherwise uncharacterized disordered proteins.

Complete chloroplast genomes of all six Hosta species occurring in Korea: molecular structures, comparative, and phylogenetic analyses

BackgroundThe genus Hosta is a group of economically appreciated perennial herbs consisting of approximately 25 species that is endemic to eastern Asia. Due to considerable morphological variability, the genus has been well recognized as a group with taxonomic problems. Chloroplast is a cytoplasmic organelle with its own genome, which is the most commonly used for phylogenetic and genetic diversity analyses for land plants. To understand the genomic architecture of Hosta chloroplasts and examine the level of nucleotide and size variation, we newly sequenced four (H. clausa, H. jonesii, H. minor, and H. venusta) and analyzed six Hosta species (including the four, H. capitata and H. yingeri) distributed throughout South Korea.ResultsThe average size of complete chloroplast genomes for the Hosta taxa was 156,642 bp with a maximum size difference of ~ 300 bp. The overall gene content and organization across the six Hosta were nearly identical with a few exceptions. There was a single tRNA gene deletion in H. jonesii and four genes were pseudogenized in three taxa (H. capitata, H. minor, and H. jonesii). We did not find major structural variation, but there were a minor expansion and contractions in IR region for three species (H. capitata, H. minor, and H. venusta). Sequence variations were higher in non-coding regions than in coding regions. Four genic and intergenic regions including two coding genes (psbA and ndhD) exhibited the largest sequence divergence showing potential as phylogenetic markers. We found compositional codon usage bias toward A/T at the third position. The Hosta plastomes had a comparable number of dispersed and tandem repeats (simple sequence repeats) to the ones identified in other angiosperm taxa. The phylogeny of 20 Agavoideae (Asparagaceae) taxa including the six Hosta species inferred from complete plastome data showed well resolved monophyletic clades for closely related taxa with high node supports.ConclusionsOur study provides detailed information on the chloroplast genome of the Hosta taxa. We identified nucleotide diversity hotspots and characterized types of repeats, which can be used for developing molecular markers applicable in various research area.

Amino acid substitution scoring matrices specific to intrinsically disordered regions in proteins

An amino acid substitution scoring matrix encapsulates the rates at which various amino acid residues in proteins are substituted by other amino acid residues, over time. Database search methods make use of substitution scoring matrices to identify sequences with homologous relationships. However, widely used substitution scoring matrices, such as BLOSUM series, have been developed using aligned blocks that are mostly devoid of disordered regions in proteins. Hence, these substitution-scoring matrices are mostly inappropriate for homology searches involving proteins enriched with disordered regions as the disordered regions have distinct amino acid compositional bias, and therefore expected to have undergone amino acid substitutions that are distinct from those in the ordered regions. We, therefore, developed a novel series of substitution scoring matrices referred to as EDSSMat by exclusively considering the substitution frequencies of amino acids in the disordered regions of the eukaryotic proteins. The newly developed matrices were tested for their ability to detect homologs of proteins enriched with disordered regions by means of SSEARCH tool. The results unequivocally demonstrate that EDSSMat matrices detect more number of homologs than the widely used BLOSUM, PAM and other standard matrices, indicating their utility value for homology searches of intrinsically disordered proteins.

Structural and phylogenetic implications of the complete mitochondrial genome of Ledra auditura

We sequenced and annotated the first complete mitochondrial genome (mitogenome) of Ledra auditura (Hemiptera: Cicadellidae: Ledrinae) and reconstructed phylogenetic relationships among 47 species (including 2 outgroup species) on the basis of 3 datasets using maximum likelihood (ML) and Bayesian inference (BI) analyses. The complete L. auditura mitogenome (length, 16,094 bp) comprises 37 genes [13 protein-coding genes (PCGs), 22 tRNAs, and 2 rRNAs], 1 control region, and 2 long non-coding regions. The first long non-coding region (length, 211 bp) is located between tRNA-I and tRNA-Q and the second region (length, 994 bp) between tRNA-S2 and ND1. All PCGs show ATN (Met/Ile) as their start codon and TAR as their stop codon. Except tRNA-S1 (AGN), which lacks the dihydrouridine arm, all tRNAs can fold into the typical cloverleaf secondary structure. The complete L. auditura mitogenome shows a base composition bias of 76.3% A + T (A = 29.9%, T = 46.4%, G = 13.3%, and C = 10.5%), negative AT skew of −0.22, and positive GC skew of 0.12. In ML and BI analyses, L. auditura was clustered with Evacanthus heimianus (Hemiptera: Cicadellidae: Evacanthinae) with strong branch support.

The Complete Mitochondrial Genome of the Caecal Fluke of Poultry, Postharmostomum commutatum, as the First Representative from the Superfamily Brachylaimoidea.

Postharmostomum commutatum (Platyhelminthes: Brachylaimoidea), a parasite of the caeca of poultry, has been frequently reported from many countries and regions, including China. However, the molecular epidemiology, population genetics and phylogenetics of this parasite are poorly understood. In the present study, we determined and characterized the complete mitochondrial (mt) genome of P. commutatum, as the first representative from the superfamily Brachylaimoidea. The mt genome of P. commutatum is a circular DNA molecule of 13,799 bp in size and encodes the complete set of 36 genes (12 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes) as well as a typical control region. The mt genome of P. commutatum presents a clear bias in nucleotide composition with a negative AT-skew on average (−0.306) and a positive GC-skew on average (0.466). Phylogenetic analyses showed that P. commutatum (superfamily Brachylaimoidea) and other ten members of the order Diplostomida were recovered as sister groups of the order Plagiorchiida, indicating that the order Diplostomida is paraphyletic. This is the first mt genome of any member of the superfamily Brachylaimoidea and should represent a rich source of genetic markers for molecular epidemiological, population genetic and phylogenetic studies of parasitic flukes of socio-economic importance in poultry.

Comprehensive analysis of synonymous codon usage patterns in orf3 gene of porcine epidemic diarrhea virus in China

The ORF3 protein of porcine epidemic diarrhea virus (PEDV) is found to function as an ion channel which influences virus virulence and production. Taking consideration of the importance of PEDV orf3 gene, we have performed comprehensive analysis to investigate its synonymous codon usage patterns. In this study, the results of base composition analysis showed A/T rich and G/C poor in PEDV orf3 genes, and the most abundant base was nucleotide T. The relative synonymous codon usage value in each codon revealed that codon usage bias existed. The mean ENC value of each gene was 48.75, indicating a low codon usage bias, as well as a relatively instable change in PEDV orf3 genes. The general correlation analysis between base composition and codon usage bias indicated that mutational bias has an impact on the PEDV codon usage bias. Neutral analysis suggested that natural selection pressure takes a more important influence than mutational bias in shaping codon usage bias. Moreover, other factors including hydrophobicity and aromaticity have been also found to influence the codon usage variation among the PEDV orf3 genes. This study not only represents the most systematic analysis of codon usage patterns in PEDV orf3 genes, but also provides a basic shaping mechanism of the codon usage bias.

SOCIAL CHANGES IN IMPRESSIONABLE YEARS AND ADULT POLITICAL ATTITUDES: EVIDENCE FROM JEWISH EXPULSIONS IN NAZI GERMANY

This study provides evidence that individuals who grew up during the 1930's Jewish expulsions are less likely to show interest and participate in politics. The estimates imply that, at the mean, individuals in their impressionable ages at the time of the expulsions are about 13% less likely to be interested in politics and 26% less likely to participate in politics. These results are not found for individuals who were older at the time of the expulsions nor for those growing up during world War (WWII). Results are robust to fixed region and birth‐year characteristics, various definitions of impressionable ages, and composition bias induced by differential migration and mortality rates across regions and cohorts. The estimates are also not driven by other regional differences in 1930's political participation, party support, Catholic share, exposure and destruction during WWII, urbanization, and other regional characteristics. We provide evidence that the adverse effects on political attitudes we find are explained by a model of political participation emphasizing the role of civic skills and socioeconomic status acquired at younger ages. Exposure to the expulsions when young is associated with lower adult volunteerism, trust, church attendance, and socioeconomic status. (JELO12, D72, D74, J15)

SURVEYING TO GAUGE STAKEHOLDERS' PARTICIPATION IN URBAN HERITAGE CONSERVATION

Abstract. Attention to the conservation of urban areas has been already raised in international legal instruments since 1962 and recently gained momentum with the adoption of the Recommendation on Historic Urban Landscapes (HUL) in 2011 (henceforth referred as the 2011 Recommendation). Considering heritage areas as drivers of economic growth but also targets of development pressures, many cities have adopted the tools provided in the 2011 Recommendation depending on their specific geographic and socio-economic contexts testing various systems to innovate in matters of urban heritage conservation and sustainable development. The present research targets the reporting of the process of implementation of the 2011 Recommendation through the creation of a survey to gauge the participation of stakeholders in worldwide case studies to establish a comparative framework. The methodology used for the survey originates from a specific case, Cuenca in Ecuador, and the secondary data collection from 101 additional case studies. The hypothesis is that the questionnaire can compile data as a tool to compare the processes for the implementation of the 2011 Recommendation and to find trends in the cases where stakeholders have become skilled at conserving urban heritage. The design of the survey demonstrates the feasibility to carry out an exploratory validation method (COSMIN) to use significant evidence coming from active stakeholders. The survey as developed becomes a powerful tool to be adapted in different contexts and as a scientific method to reach a proper statistic confidential level, margin error and to control bias in the sample composition.

Protein Abundance Biases the Amino Acid Composition of Disordered Regions to Minimize Non-functional Interactions

In eukaryotes, disordered regions cover up to 50% of proteomes and mediate fundamental cellular processes. In contrast to globular domains, where about half of the amino acids are buried in the protein interior, disordered regions show higher solvent accessibility, which makes them prone to engage in non-functional interactions. Such interactions are exacerbated by the law of mass action, prompting the question of how they are minimized in abundant proteins. We find that interaction propensity or “stickiness” of disordered regions negatively correlates with their cellular abundance, both in yeast and human. Strikingly, considering yeast proteins where a large fraction of the sequence is disordered, the correlation between stickiness and abundance reaches R=−0.55. Beyond this global amino-acid composition bias, we identify three rules by which amino-acid composition of disordered regions adjusts with high abundance. First, lysines are preferred over arginines, consistent with the latter amino acid being stickier than the former. Second, compensatory effects exist, whereby a sticky region can be tolerated if it is compensated by a distal non-sticky region. Third, such compensation requires a lower average stickiness at the same abundance when compared to a scenario where stickiness is homogeneous throughout the sequence. We validate these rules experimentally, employing them as different strategies to rescue an otherwise sticky protein fragment from aggregation. Our results highlight that non-functional interactions represent a significant constraint in cellular systems and reveal simple rules by which protein sequences adapt to that constraint. Data from this work are deposited in Figshare, at https://doi.org/10.6084/m9.figshare.8068937.v3.

Repeatability in protein sequences

Low complexity regions (LCRs) in protein sequences have special properties that are very different from those of globular proteins. The rules that define secondary structure elements do not apply when the distribution of amino acids becomes biased. While there is a tendency towards structural disorder in LCRs, various examples, and particularly homorepeats of single amino acids, suggest that very short repeats could adopt structures very difficult to predict. These structures are possibly variable and dependant on the context of intra- or inter-molecular interactions. In general, short repeats in LCRs can induce structure. This could explain the observation that very short (non-perfect) repeats are widespread and many define regions with a function in protein interactions. For these reasons, we have developed an algorithm to quickly analyze local repeatability along protein sequences, that is, how close a protein fragment is from a perfect repeat. Using this algorithm we identified that the proteins of the yeast Saccharomyces cerevisiae are depleted in short repeats (approximate or not) of odd-length, while the human proteins are not, that the fish Danio rerio has many proteins with repeats of length two and that the plant Arabidopsis thaliana has an unusually large amount of repeats of length seven. Our method (REpeatability Scanner, RES, accessible at http://cbdm-01.zdv.uni-mainz.de/~munoz/res/) allows to find regions with approximate short repeats in protein sequences, and helps to characterize the variable use of LCRs and compositional bias in different organisms.

Information theoretic measures for quantifying sequence-ensemble relationships of intrinsically disordered proteins.

Intrinsically disordered proteins (IDPs) contribute to a multitude of functions. De novo design of IDPs should open the door to modulating functions and phenotypes controlled by these systems. Recent design efforts have focused on compositional biases and specific sequence patterns as the design features. Analysis of the impact of these designs on sequence-function relationships indicates that individual sequence/compositional parameters are insufficient for describing sequence-function relationships in IDPs. To remedy this problem, we have developed information theoretic measures for sequence–ensemble relationships (SERs) of IDPs. These measures rely on prior availability of statistically robust conformational ensembles derived from all atom simulations. We show that the measures we have developed are useful for comparing sequence-ensemble relationships even when sequence is poorly conserved. Based on our results, we propose that de novo designs of IDPs, guided by knowledge of their SERs, should provide improved insights into their sequence–ensemble–function relationships.

Genomic Evidence for Simultaneous Optimization of Transcription and Translation through Codon Variants in the pmoCAB Operon of Type Ia Methanotrophs.

Understanding the interplay between genotype and phenotype is a fundamental goal of functional genomics. Methane oxidation is a microbial phenotype with global-scale significance as part of the carbon biogeochemical cycle and a sink for greenhouse gas. Microorganisms that oxidize methane (methanotrophs) are taxonomically diverse and widespread around the globe. In methanotrophic bacteria, enzymes in the methane oxidation metabolic module (KEGG module M00174, conversion of methane to formaldehyde) are encoded in four operons (pmoCAB, mmoXYZBCD, mxaFI, and xoxF). Recent reports have suggested that methanotrophs in Proteobacteria acquired methane monooxygenases through horizontal gene transfer. Here, we used a genomic meta-analysis to infer the transcriptional and translational advantages of coding sequences from the methane oxidation metabolic modules of different types of methanotrophs. By analyzing isolate and metagenome-assembled genomes from phylogenetically and geographically diverse sources, we detected an anomalous nucleotide composition bias in the coding sequences of particulate methane monooxygenase genes (pmoCAB) from type Ia methanotrophs. We found that this nucleotide bias increases the level of codon bias by decreasing the GC content in the third base of codons, a strategy that contrasts with that of other coding sequences in the module. Further codon usage analyses uncovered that codon variants of the type Ia pmoCAB coding sequences deviate from the genomic signature to match ribosomal protein-coding sequences. Subsequently, computation of transcription and translation metrics revealed that the pmoCAB coding sequences of type Ia methanotrophs optimize the usage of codon variants to maximize translation efficiency and accuracy, while minimizing the synthesis cost of transcripts and proteins.IMPORTANCE Microbial methane oxidation plays a fundamental role in the biogeochemical cycle of Earth's system. Recent reports have provided evidence for the acquisition of methane monooxygenases by horizontal gene transfer in methane-oxidizing bacteria from different environments, but how evolution has shaped the coding sequences to execute methanotrophy efficiently remains unexplored. In this work, we provide genomic evidence that among the different types of methanotrophs, type Ia methanotrophs possess a unique coding sequence of the pmoCAB operon that is under positive selection for optimal resource allocation and efficient synthesis of transcripts and proteins. This adaptive trait possibly enables type Ia methanotrophs to respond robustly to fluctuating methane availability and explains their global prevalence.

Chloroplast phylogenomics of liverworts: a reappraisal of the backbone phylogeny of liverworts with emphasis on Ptilidiales.

As one of the four main lineages diverging from the early diversification of land plants, the phylogeny of liverworts holds the information about nearly 500Myr of independent adaptation to changing environments. Thus, resolving the phylogenetic history of liverworts will provide unique insights into the successful diversification of early land plants in terrestrial ecosystems. However, the deep diverging events of this group remain incompletely resolved, such as the definite position of Ptilidiales. Here, we aimed to reconstruct the backbone relationships of liverworts using 84 protein-coding chloroplast genes, a dataset comprising 35 representatives from all major lineages of liverworts, and three phylogenetic analyses, namely maximum parsimony, maximum likelihood and Bayesian inference. To test the impact of composition biases, the phylogenetic analyses were carried out using three alignments representing the same dataset either as: (i) nucleotides, (ii) amino acids, or (iii) recoded nucleotides applying ambiguity base code. Chloroplast genome data consistently supported the monophyletic origin of three major lineages in liverworts, as well as the majority of backbone relationships. Ptilidiales were found to be sister to Jungermanniales. The rapid accumulation of G/C tracks as a consequence of increased GC content is an important cause for the long branches inferred in this group. Our study not only provides empirical evidence to support the significance of plastid genome sequencing to reconstruct the phylogeny of this important plant lineage, but also suggests that the GC content has played a critical role in the evolutionary dynamics of plastid genomes in land plants.