Comparative plastome analysis of Blumea,with implications for genome evolution and phylogeny of Asteroideae

The genus Blumea (Asteroideae, Asteraceae) comprises about 100 species, including herbs, shrubs, and small trees. Previous studies have been unable to resolve taxonomic issues and the phylogeny of the genus Blumea due to the low polymorphism of molecular markers. Therefore, suitable polymorphic regions need to be identified. Here, we de novo assembled plastomes of the three Blumea species B. oxyodonta, B. tenella, and B. balsamifera and compared them with 26 other species of Asteroideae after correction of annotations. These species have quadripartite plastomes with similar gene content, genome organization, and inverted repeat contraction and expansion comprising 113 genes, including 80 protein‐coding, 29 transfer RNA, and 4 ribosomal RNA genes. The comparative analysis of codon usage, amino acid frequency, microsatellite repeats, oligonucleotide repeats, and transition and transversion substitutions has revealed high resemblance among the newly assembled species of Blumea. We identified 10 highly polymorphic regions with nucleotide diversity above 0.02, including rps16‐trnQ, ycf1, ndhF‐rpl32, petN‐psbM, and rpl32‐trnL, and they may be suitable for the development of robust, authentic, and cost‐effective markers for barcoding and inference of the phylogeny of the genus Blumea. Among these highly polymorphic regions, five regions also co‐occurred with oligonucleotide repeats and support use of repeats as a proxy for the identification of polymorphic loci. The phylogenetic analysis revealed a close relationship between Blumea and Pluchea within the tribe Inuleae. At tribe level, our phylogeny supports a sister relationship between Astereae and Anthemideae rooted as Gnaphalieae, Calenduleae, and Senecioneae. These results are contradictory to recent studies which reported a sister relationship between “Senecioneae and Anthemideae” and “Astereae and Gnaphalieae” or a sister relationship between Astereae and Gnaphalieae rooted as Calenduleae, Anthemideae, and then Senecioneae using nuclear genome sequences. The conflicting phylogenetic signals observed at the tribal level between plastidt and nuclear genome data require further investigation.     

Georeferenced phylogenetic analysis of a global collection of wild and cultivated Citrullus species

The geographical origin of watermelon (Citrullus lanatus) remains debated. While a first hypothesis suggests the center of origin to be West Africa, where the endemic sister species C. mucosospermus thrives, a second hypothesis suggests northeastern Africa where the white‐fleshed Sudanese Kordophan melon is cultivated. In this study, we infer biogeographical and haplotype genealogy for C. lanatus, C. mucosospermus, C. amarus, and C. colocynthis using noncoding cpDNA sequences (trnT‐trnL and ndhF‐rpl32 regions) from a global collection of 135 accessions. In total, we identified 38 haplotypes in C. lanatus, C. mucosospermus, C. amarus, and C. colocynthis; of these, 21 were found in Africa and 17 appear endemic to the continent. The least diverse species was C. mucosospermus (5 haplotypes) and the most diverse was C. colocynthis (16 haplotypes). Some haplotypes of C. mucosospermus were nearly exclusive to West Africa, and C. lanatus and C. mucosospermus shared haplotypes that were distinct from those of both C. amarus and C. colocynthis. The results support previous findings that revealed C. mucosospermus to be the closest relative to C. lanatus (including subsp. cordophanus). West Africa, as a center of endemism of C. mucosospermus, is an area of interest in the search of the origin of C. lanatus. This calls for further historical and phylogeographical investigations and wider collection of samples in West and northeastern Africa.     

Plastid genome sequencing,identification of nuclear SNP markers,and quality assessment of medicinal rhizomatous herb Polygonatum odoratum (Asparagaceae) cultivars

Polygonatum odoratum (Mill.) Druce (Asparagaceae, Asparagales) is a widely cultivated medicinal herb in China. However, this useful herb is understudied despite being known as a medicinal resource with top grade medical and edible properties since long. In this study, P. odoratum and four cultivars were investigated. The variations in morphological characteristics and vegetative phases of each cultivar were observed. For genetic aspect, the plastid genome of P. odoratum varies in length from 154,569 bp to 155,491 bp, containing a large single‐copy region of 83,486–84,459 bp, a small single‐copy region of 18,292–18,471 bp, and two inverted repeats of 26,302–26,370 bp. A total of 131 genes were predicted, including 85 protein‐coding, 38 tRNA, and eight rRNA genes. Genome comparisons revealed a slight variation in the sequence across the five accessions, but two highly variable regions (trnC‐petN and rpl32‐trnL) were detected when comparing the four different cultivars. For the RAD‐seq markers, a total of 33.64 Gb of clean data, with an average value of 1.08 Gb per sample, were analyzed for the presence of single nucleotide polymorphisms (SNPs). Well‐resolved phylogenies of the P. odoratum cultivars are constructed; the nonmonophyletic relationship in the plastome‐based phylogenetic trees, yet monophyletic form in the RAD‐based linkage map suggested possibility of hybrid cultivar for P. odoratum “Dazhu” (GDDZ), which was further supported by morphological observations. Quality assessment based on the standards of the Chinese Pharmacopoeia on Polygonati Odorati Rhizoma (POR) on the four cultivars used in this study recorded that PORs from P. odoratum “Zhongzhu” (GDZZ) met the minimum criteria for the acceptance as raw material for medicinal drug production. This study has provided insights on the morphological variations, genetic background, and medicinal qualities of P. odoratum cultivars that could be explored for future genetic improvement as well as breeding programs of P. odoratum for POR production.     

Integrated analysis of three newly sequenced fern chloroplast genomes: Genome structure and comparative analysis

BackgroundSome ferns have medicinal properties and are used in therapeutic interventions. However, the classification and phylogenetic relationships of ferns remain incompletely reported. Considering that chloroplast genomes provide ideal information for species identification and evolution, in this study, three unpublished and one published ferns were sequenced and compared with other ferns to obtain comprehensive information on their classification and evolution.Materials and MethodsThe complete chloroplast genomes of Dryopteris goeringiana (Kunze) Koidz, D. crassirhizoma Nakai, Athyrium brevifrons Nakai ex Kitagawa, and Polystichum tripteron (Kunze) Presl were sequenced using the Illumina HiSeq 4,000 platform. Simple sequence repeats (SSRs), nucleotide diversity analysis, and RNA editing were investigated in all four species. Genome comparison and inverted repeats (IR) boundary expansion and contraction analyses were also performed. The relationships among the ferns were studied by phylogenetic analysis based on the whole chloroplast genomes.ResultsThe whole chloroplast genomes ranged from 148,539 to 151,341 bp in size and exhibited typical quadripartite structures. Ten highly variable loci with parsimony informative (Pi) values of > 0.02 were identified. A total of 75–108 SSRs were identified, and only six SSRs were present in all four ferns. The SSRs contained a higher number of A + T than G + C bases. C‐to‐U conversion was the most common type of RNA editing event. Genome comparison analysis revealed that single‐copy regions were more highly conserved than IR regions. IR boundary expansion and contraction varied among the four ferns. Phylogenetic analysis showed that species in the same genus tended to cluster together with and had relatively close relationships.ConclusionThe results provide valuable information on fern chloroplast genomes that will be useful to identify and classify ferns, and study their phylogenetic relationships and evolution.     

Phylogenetic relationships of the chloroplast genomes in the genus Glycine inferred from four intergenic spacer sequences

The nucleotide sequences of four intergenic spacer regions of chloroplast DNA, atpB-rbcL, trnS-trnG, rps11-rpl36, and rps3-rpl16, were analyzed in the genus Glycine. Phylogenetic analysis based on the sequence data using Neonotonia wightii as the outgroup generated trees supporting the classification of two subgenera, Soja and Glycine, and three plastome groups in the subgenus Glycine. The results were consistent with the presence of diversified chloroplast genomes within tetraploid plants of G. tabacina and G. tomentella, as well as with a close relationship between G. tomentella and G. dolichocarpa that had been suggested based on morphological analyses. Little sequence variation was found in the subgenus Soja, suggesting that G. soja rapidly expanded its distribution in East Asia. The analysis also showed that the differentiation into three plastome groups in the subgenus Glycine occurred in the early stages of its evolution, after the two subgenera diverged.     

Evidence of horizontal gene transfer between land plant plastids has surprising conservation implications

Background and AimsHorizontal gene transfer (HGT) is an important evolutionary mechanism because it transfers genetic material that may code for traits or functions between species or genomes. It is frequent in mitochondrial and nuclear genomes but has not been demonstrated between plastid genomes of different green land plant species.MethodsWe Sanger-sequenced the nuclear internal transcribed spacers (ITS1 and 2) and the plastid rpl16 G2 intron (rpl16). In five individuals with foreign rpl16 we also sequenced atpB-rbcL and trnL_UAA-trnF_GAA.Key ResultsWe discovered 14 individuals of a moss species with typical nuclear ITSs but foreign plastid rpl16 from a species of a distant lineage. None of the individuals with three plastid markers sequenced contained all foreign markers, demonstrating the transfer of plastid fragments rather than the entire plastid genome, i.e. entire plastids were not transferred. The two lineages diverged 165–185 Myr BP. The extended time interval since lineage divergence suggests that the foreign rpl16 is more likely explained by HGT than by hybridization or incomplete lineage sorting.ConclusionsWe provide the first conclusive evidence of interspecific plastid-to-plastid HGT among land plants. Two aspects are critical: it occurred at several localities during the massive colonization of recently disturbed open habitats that were created by large-scale liming as a freshwater biodiversity conservation measure; and it involved mosses whose unique life cycle includes spores that first develop a filamentous protonema phase. We hypothesize that gene transfer is facilitated when protonema filaments of different species intermix intimately when colonizing disturbed early succession habitats.     

Phylogeography of sugar kelp: Northern ice‐age refugia in the Gulf of Alaska

Many Northeast (NE) Pacific fishes and invertebrates survived Pleistocene glaciations in northern refugia, but the extent that kelps survived in northern areas is uncertain. Here, we test the hypothesis that populations of sugar kelp (Saccharina latissima) persisted in the Gulf of Alaska during ice‐age maxima when the western margin of the Cordilleran ice sheet covered coastal areas around the NE Pacific Ocean. We estimated genetic diversities within and phylogeographical relationships among 14 populations along 2,800 km in the NE Pacific and Bering Sea with partial sequences of mitochondrial DNA 5′‐cytochrome oxidase subunit I (COI, bp = 624, n = 543), chloroplast DNA ribulose‐1,5‐bisphosphate carboxylase large subunit‐3′ (rbcL, bp = 735, n = 514), and 11 microsatellite loci. Concatenated sequences of rbcL and COI showed moderate levels of within‐population genetic diversity (mean h = 0.200) but substantial differences among populations (Φ_ST = 0.834, p < .0001). Microsatellites showed moderate levels of heterozygosity within populations (mean H _E = 0.391). Kelps in the same organellar lineage tended to cluster together, regardless of geographic origins, as indicated in a principal coordinate analysis (PCoA) of microsatellite genotypes. The PCoA also showed evidence of nuclear hybridizations between co‐occurring organellar lineages. Individual admixture plots with population clusters of K = 2, 6, and 9 showed increasing complexity with considerable historical admixture between some clusters. A time‐calibrated phylogeny placed divergences between rbcL‐COI lineages at 1.4 million years at most. The time frames of mutation in the rbcL‐COI lineages and microsatellite population clusters differed among locations. The existence of ancient lineages in the Gulf of Alaska, moderate levels of genetic diversity, and the absence of departures from neutrality are consistent with northern refugia during multiple Croll‐Milankovitch climate cycles in the Pleistocene Epoch.     

Dynamics of NK,CD8 and Tfh cell mediated the production of cytokines and antiviral antibodies in Chinese patients with moderate COVID‐19

Recent studies have demonstrated a marked decrease in peripheral lymphocyte levels in patients with coronavirus disease 2019 (COVID‐19) caused by severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2). Few studies have focused on the changes of NK, T‐ and B‐cell subsets, inflammatory cytokines and virus‐specific antibodies in patients with moderate COVID‐19. A total of 11 RT‐PCR‐confirmed convalescent patients with COVID‐19 and 11 patients with non‐SARS‐CoV‐2 pneumonia (control patients) were enrolled in this study. NK, CD8⁺ T, CD4⁺ T, Tfh‐like and B‐cell subsets were analysed using flow cytometry. Cytokines and SARS‐CoV‐2‐specific antibodies were analysed using an electrochemiluminescence immunoassay. NK cell counts were significantly higher in patients with COVID‐19 than in control patients (P = 0.017). Effector memory CD8⁺ T‐cell counts significantly increased in patients with COVID‐19 during a convalescent period of 1 week (P = 0.041). TIM‐3⁺ Tfh‐like cell and CD226⁺ Tfh‐like cell counts significantly increased (P = 0.027) and decreased (P = 0.022), respectively, during the same period. Moreover, ICOS⁺ Tfh‐like cell counts tended to decrease (P = 0.074). No abnormal increase in cytokine levels was observed. The high expression of NK cells is important in innate immune response against SARS‐CoV‐2. The increase in effector memory CD8⁺ T‐cell counts, the up‐regulation of inhibitory molecules and the down‐regulation of active molecules on CD4⁺ T cells and Tfh‐like cells in patients with COVID‐19 would benefit the maintenance of balanced cellular and humoural immune responses, may prevent the development of severe cases and contribute to the recovery of patients with COVID‐19.     

Quorum sensing governs a transmissive Legionella subpopulation at the pathogen vacuole periphery

The Gram‐negative bacterium Legionella pneumophila is the causative agent of Legionnaires'' disease and replicates in amoebae and macrophages within a distinct compartment, the Legionella‐containing vacuole (LCV). The facultative intracellular pathogen switches between a replicative, non‐virulent and a non‐replicating, virulent/transmissive phase. Here, we show on a single‐cell level that at late stages of infection, individual motile (P_flaA‐GFP‐positive) and virulent (P_ralF‐ and P_sidC‐GFP‐positive) L. pneumophila emerge in the cluster of non‐growing bacteria within an LCV. Comparative proteomics of P_flaA‐GFP‐positive and P_flaA‐GFP‐negative L. pneumophila subpopulations reveals distinct proteomes with flagellar proteins or cell division proteins being preferentially produced by the former or the latter, respectively. Toward the end of an infection cycle (˜ 48 h), the P_flaA‐GFP‐positive L. pneumophila subpopulation emerges at the cluster periphery, predominantly escapes the LCV, and spreads from the bursting host cell. These processes are mediated by the Legionella quorum sensing (Lqs) system. Thus, quorum sensing regulates the emergence of a subpopulation of transmissive L. pneumophila at the LCV periphery, and phenotypic heterogeneity underlies the intravacuolar bi‐phasic life cycle of L. pneumophila.     

IFN‐γ+IL‐17+Th17 cells regulate fibrosis through secreting IL‐21 in systemic scleroderma

This study aimed to explore the function of IFN‐γ⁺IL‐17⁺Th17 cells on fibrosis in systemic scleroderma (SSc). Blood and skin samples were collected from 20 SSc cases and 10 healthy individuals. The percentage of IFN‐γ⁺IL‐17⁺Th17 cells was detected using flow cytometry. The in vitro induction of IFN‐γ⁺IL‐17⁺Th17 cells was performed adopting PHA and rIL‐12. Gene expression was detected via quantitative real‐time polymerase chain reaction (qRT‐PCR), whereas western blot analysis was adopted for protein analysis. The distribution of IFN‐γ⁺IL‐17⁺Th17 cells was significantly increased in SSc cases and positively correlated with SSc stages (P = .031), disease duration (P = .016), activity (P = .025) and skin scores (P < .001). In vitro, IFN‐γ⁺IL‐17⁺Th17 cells could promote the expressions of α‐SMA and COL1A1, revealing increased fibroblasts’ proliferation and enhanced collagen‐secreting capacity. In addition, IL‐21 expression was significantly increased in co‐culture medium of IFN‐γ⁺IL‐17⁺Th17 cells and fibroblasts (P < .001). IL‐21 neutralizer treatment resulted in the down‐regulation of α‐SMA and COL1A1. IL‐21 was confirmed as an effector of IFN‐γ⁺IL‐17⁺Th17 cells in fibrosis process. The distribution of IFN‐γ⁺IL‐17⁺Th17 cells was significantly increased in SSc cases and positively correlated with disease activity. IFN‐γ⁺IL‐17⁺Th17 cells could promote fibroblast proliferation and enhance collagen‐secreting ability via producing IL‐21, thus contributing to fibrosis in SSc.     

Telomere length and metabolic syndrome traits: A Mendelian randomisation study

Observational studies have revealed associations between short leucocyte telomere length (LTL), a TL marker in somatic tissues and multiple Metabolic Syndrome (MetS) traits. Animal studies have supported these findings by showing that increased telomere attrition leads to adipose tissue dysfunction and insulin resistance. We investigated the associations between genetically instrumented LTL and MetS traits using Mendelian Randomisation (MR). Fifty‐two independent variants identified at FDR<0.05 from a genome‐wide association study (GWAS) including 78,592 Europeans and collectively accounting for 2.93% of LTL variance were selected as genetic instruments for LTL. Summary‐level data for MetS traits and for the MetS as a binary phenotype were obtained from the largest publicly available GWAS and two‐sample MR analyses were used to estimate the associations of LTL with these traits. The combined effect of the genetic instruments was modelled using inverse variance weighted regression and sensitivity analyses with MR‐Egger, weighted‐median and MR‐PRESSO were performed to test for and correct horizonal pleiotropy. Genetically instrumented longer LTL was associated with higher waist‐to‐hip ratio adjusted for body mass index (β = 0.045 SD, SE = 0.018, p = 0.01), raised systolic (β = 1.529 mmHg, SE = 0.332, p = 4x10⁻⁶) and diastolic (β = 0.633 mmHg, SE = 0.222, p = 0.004) blood pressure, and increased MetS risk (OR = 1.133, 95% CI 1.057–1.215). Consistent results were obtained in sensitivity analyses, which provided no evidence of unbalanced horizontal pleiotropy. Telomere shortening might not be a major driver of cellular senescence and dysfunction in human adipose tissue. Future experimental studies should examine the mechanistic bases for the links between longer LTL and increased upper‐body fat distribution and raised blood pressure.     

Low base‐substitution mutation rate and predominance of insertion‐deletion events in the acidophilic bacterium Acidobacterium capsulatum

Analyses of spontaneous mutation have shown that total genome‐wide mutation rates are quantitatively similar for most prokaryotic organisms. However, this view is mainly based on organisms that grow best around neutral pH values (6.0–8.0). In particular, the whole‐genome mutation rate has not been determined for an acidophilic organism. Here, we have determined the genome‐wide rate of spontaneous mutation in the acidophilic Acidobacterium capsulatum using a direct and unbiased method: a mutation‐accumulation experiment followed by whole‐genome sequencing. Evaluation of 69 mutation accumulation lines of A. capsulatum after an average of ~2900 cell divisions yielded a base‐substitution mutation rate of 1.22 × 10⁻¹⁰ per site per generation or 4 × 10⁻⁴ per genome per generation, which is significantly lower than the consensus value (2.5−4.6 × 10⁻³) of mesothermophilic (~15–40°C) and neutrophilic (pH 6–8) prokaryotic organisms. However, the insertion‐deletion rate (0.43 × 10⁻¹⁰ per site per generation) is high relative to the base‐substitution mutation rate. Organisms with a similar effective population size and a similar expected effect of genetic drift should have similar mutation rates. Because selection operates on the total mutation rate, it is suggested that the relatively high insertion‐deletion rate may be balanced by a low base‐substitution rate in A. capsulatum, with selection operating on the total mutation rate.     

Contributions of PARP‐1 rs1136410 C>T polymorphism to the development of cancer

Poly(ADP‐ribose) polymerase‐1 (PARP‐1) is a nuclear chromatin‐associated enzyme involved in the DNA damage response. SNP rs1136410 C>T, the most studied polymorphism in PARP‐1 gene, is highly implicated in the susceptibility of cancer. However, the roles of PARP‐1 rs1136410 C>T on cancer risk vary from different studies. We comprehensively screened all qualified publications from several databases, including PubMed, EMBASE, MEDLINE, CNKI and Wanfang. The searching was updated to April 2020. Our meta‐analysis included 60 articles with 65 studies, comprised of a total of 23 996 cases with cancer and 33 015 controls. Overall, pooled data showed that the PARP‐1 rs1136410 C>T polymorphism was significantly but a border‐line associated with an increased risk of overall cancer (CC vs. TT/TC: OR = 1.11, 95% CI = 1.00‐1.24; C vs T: OR = 1.07, 95% CI = 1.01‐1.14). Subgroup analysis indicated that rs1136410 C allele contributed to high risk among gastric, thyroid, and cervical cancer, but lower risk among brain cancer. Furthermore, increased cancer risk was detected in the subgroups of Asian, controls from population‐based design studies, and HWE ≤ 0.05 studies. Sensitivity analysis and Egger''s test showed that results of the meta‐analysis were fairly stable. The current study indicated that PARP1 rs1136410 C>T polymorphism may have an impact on certain types of cancer susceptibility.     

Identification of novel susceptibility loci for non‐syndromic cleft lip with or without cleft palate

Although several genome‐wide association studies (GWAS) of non‐syndromic cleft lip with or without cleft palate (NSCL/P) have been reported, more novel association signals are remained to be exploited. Here, we performed an in‐depth analysis of our previously published Chinese GWAS cohort study with replication in an extra dbGaP case‐parent trios and another in‐house Nanjing cohort, and finally identified five novel significant association signals (rs11119445: 3’ of SERTAD4, P = 6.44 × 10⁻¹⁴; rs227227 and rs12561877: intron of SYT14, P = 5.02 × 10⁻¹³ and 2.80 × 10⁻¹¹, respectively; rs643118: intron of TRAF3IP3, P = 4.45 × 10⁻⁶; rs2095293: intron of NR6A1, P = 2.98 × 10⁻⁵). The mean (standard deviation) of the weighted genetic risk score (wGRS) from these SNPs was 1.83 (0.65) for NSCL/P cases and 1.58 (0.68) for controls, respectively (P = 2.67 × 10⁻¹⁶). Rs643118 was identified as a shared susceptible factor of NSCL/P among Asians and Europeans, while rs227227 may contribute to the risk of NSCL/P as well as NSCPO. In addition, sertad4 knockdown zebrafish models resulted in down‐regulation of sox2 and caused oedema around the heart and mandibular deficiency, compared with control embryos. Taken together, this study has improved our understanding of the genetic susceptibility to NSCL/P and provided further clues to its aetiology in the Chinese population.     

Temporal and spatial variation in population structure among brooding sea stars in the genus Leptasterias

Temporal genetic studies of low‐dispersing organisms are rare. Marine invertebrates lacking a planktonic larval stage are expected to have lower dispersal, low gene flow, and a higher potential for local adaptation than organisms with planktonic dispersal. Leptasterias is a genus of brooding sea stars containing several cryptic species complexes. Population genetic methods were used to resolve patterns of fine‐scale population structure in central California Leptasterias species using three loci from nuclear and mitochondrial genomes. Historic samples (collected between 1897 and 1998) were compared to contemporary samples (collected between 2008 and 2014) to delineate changes in species distributions in space and time. Phylogenetic analysis of contemporary samples confirmed the presence of a bay‐localized clade and revealed the presence of an additional bay‐localized and previously undescribed clade of Leptasterias. Analysis of contemporary and historic samples indicates two clades are experiencing a constriction in their southern range limit and suggests a decrease in clade‐specific abundance at sites at which they were once prevalent. Historic sampling revealed a dramatically different distribution of diversity along the California coastline compared to contemporary sampling and illustrates the importance of temporal genetic sampling in phylogeographic studies. These samples were collected prior to significant impacts of Sea Star Wasting Disease (SSWD) and represent an in‐depth analysis of genetic structure over 117 years prior to the SSWD‐associated mass die‐off of Leptasterias.     

The influence of phylogeny and life history on telomere lengths and telomere rate of change among bird species: A meta‐analysis

Longevity is highly variable among animal species and has coevolved with other life‐history traits, such as body size and rates of reproduction. Telomeres, through their erosion over time, are one of the cell mechanisms that produce senescence at the cell level and might even have an influence on the rate of aging in whole organisms. However, uneroded telomeres are also risk factors of cell immortalization. The associations of telomere lengths, their rate of change, and life‐history traits independent of body size are largely underexplored for birds. To test associations of life‐history traits and telomere dynamics, we conducted a phylogenetic meta‐analysis using studies of 53 species of birds. We restricted analyses to studies that applied the telomere restriction fragment length (TRF) method, and examined relationships between mean telomere length at the chick (Chick TL) and adult (Adult TL) stages, the mean rate of change in telomere length during life (TROC), and life‐history traits. We examined 3 principal components of 12 life‐history variables that represented: body size (PC1), the slow–fast continuum of pace of life (PC2), and postfledging parental care (PC3). Phylogeny had at best a small‐to‐medium influence on Adult and Chick TL (r ² = .190 and .138, respectively), but a substantial influence on TROC (r ² = .688). Phylogeny strongly influenced life histories: PC1 (r ² = .828), PC2 (.838), and PC3 (.613). Adult TL and Chick TL were poorly associated with the life‐history variables. TROC, however, was negatively and moderate‐to‐strongly associated with PC2 (unadjusted r = −.340; with phylogenetic correction, r = −.490). Independent of body size, long‐lived species with smaller clutches, and slower embryonic rate of growth may exhibit less change in telomere length over their lifetimes. We suggest that telomere lengths may have diverged, even among closely avian‐related species, yet telomere dynamics are strongly linked to the pace of life.     

High‐altitude adaptation in vertebrates as revealed by mitochondrial genome analyses

The high‐altitude environment may drive vertebrate evolution in a certain way, and vertebrates living in different altitude environments might have different energy requirements. We hypothesized that the high‐altitude environment might impose different influences on vertebrate mitochondrial genomes (mtDNA). We used selection pressure analyses and PIC (phylogenetic independent contrasts) analysis to detect the evolutionary rate of vertebrate mtDNA protein‐coding genes (PCGs) from different altitudes. The results showed that the ratio of nonsynonymous/synonymous substitutions (dN/dS) in the mtDNA PCGs was significantly higher in high‐altitude vertebrates than in low‐altitude vertebrates. The seven rapidly evolving genes were shared by the high‐altitude vertebrates, and only one positive selection gene (ND5 gene) was detected in the high‐altitude vertebrates. Our results suggest the mtDNA evolutionary rate in high‐altitude vertebrates was higher than in low‐altitude vertebrates as their evolution requires more energy in a high‐altitude environment. Our study demonstrates the high‐altitude environment (low atmospheric O₂ levels) drives vertebrate evolution in mtDNA PCGs.     

Evolutionary patterns of nucleotide substitution rates in plastid genomes of Quercus

Molecular evolution, including nucleotide substitutions, plays an important role in understanding the dynamics and mechanisms of species evolution. Here, we sequenced whole plastid genomes (plastomes) of Quercus fabri, Quercus semecarpifolia, Quercus engleriana, and Quercus phellos and compared them with 14 other Quercus plastomes to explore their evolutionary relationships using 67 shared protein‐coding sequences. While many previously identified evolutionary relationships were found, our findings do not support previous research which retrieve Quercus subg. Cerris sect. Ilex as a monophyletic group, with sect. Ilex found to be polyphyletic and composed of three strongly supported lineages inserted between sections Cerris and Cyclobalanposis. Compared with gymnosperms, Quercus plastomes showed higher evolutionary rates (Dn/Ds = 0.3793). Most protein‐coding genes experienced relaxed purifying selection, and the high Dn value (0.1927) indicated that gene functions adjusted to environmental changes effectively. Our findings suggest that gene interval regions play an important role in Quercus evolution. We detected greater variation in the intergenic regions (trnH‐psbA, trnK_UUU‐rps16, trnfM_CAU‐rps14, trnS_GCU‐trnG_GCC, and atpF‐atpH), intron losses (petB and petD), and pseudogene loss and degradation (ycf15). Additionally, the loss of some genes suggested the existence of gene exchanges between plastid and nuclear genomes, which affects the evolutionary rate of the former. However, the connective mechanism between these two genomes is still unclear.     

Genotyping validates the efficacy of photographic identification in a capture‐mark‐recapture study based on the head scale patterns of the prairie lizard (Sceloporus consobrinus)

Population studies often incorporate capture‐mark‐recapture (CMR) techniques to gather information on long‐term biological and demographic characteristics. A fundamental requirement for CMR studies is that an individual must be uniquely and permanently marked to ensure reliable reidentification throughout its lifespan. Photographic identification involving automated photographic identification software has become a popular and efficient noninvasive method for identifying individuals based on natural markings. However, few studies have (a) robustly assessed the performance of automated programs by using a double‐marking system or (b) determined their efficacy for long‐term studies by incorporating multi‐year data. Here, we evaluated the performance of the program Interactive Individual Identification System (I³S) by cross‐validating photographic identifications based on the head scale pattern of the prairie lizard (Sceloporus consobrinus) with individual microsatellite genotyping (N = 863). Further, we assessed the efficacy of the program to identify individuals over time by comparing error rates between within‐year and between‐year recaptures. Recaptured lizards were correctly identified by I³S in 94.1% of cases. We estimated a false rejection rate (FRR) of 5.9% and a false acceptance rate (FAR) of 0%. By using I³S, we correctly identified 97.8% of within‐year recaptures (FRR = 2.2%; FAR = 0%) and 91.1% of between‐year recaptures (FRR = 8.9%; FAR = 0%). Misidentifications were primarily due to poor photograph quality (N = 4). However, two misidentifications were caused by indistinct scale configuration due to scale damage (N = 1) and ontogenetic changes in head scalation between capture events (N = 1). We conclude that automated photographic identification based on head scale patterns is a reliable and accurate method for identifying individuals over time. Because many lizard or reptilian species possess variable head squamation, this method has potential for successful application in many species.     

Inhibition of endoplasmic reticulum stress prevents high‐fat diet mediated atrial fibrosis and fibrillation

Obesity is a significant risk factor for atrial fibrillation (AF), which is the most common sustained arrhythmia with increased mortality and morbidity. High‐fat diet (HFD)‐induced obesity is associated with the activation of endoplasmic reticulum stress (ERS). However, the role of ERS in HFD‐induced AF remains elusive. Human atrium samples were examined for the ERS activation test. C57BL/6J mice were divided into four groups, including the control group, the HFD group, the 4‐phenylbutyric acid (4‐PBA) group, and the HFD + 4‐PBA group. At the age of 4 weeks, the HFD group and the HFD + 4‐PBA group were given HFD to construct the obesity model, while the other two groups were given a normal diet (ND). Transesophageal programmed electrical stimulation was conducted to evaluate the AF inducibility and duration. Atrial fibrosis and ERS activation were also investigated.We found that CHOP and GRP‐78 protein were significantly higher in overweight patients than the controls (both P < 0.05). AF inducibility and duration of the HFD group were significantly higher than the other groups (both P < 0.05), while there was no difference between those groups (P > 0.05). The mice of the HFD group had significantly higher collagen volume fraction (CVF%) than the other groups (P < 0.05). ERS marker protein of GRP78, p‐PERK, ATF6 and CHOP protein expression level was increased in the HFD group, which were significantly mitigated in the HFD + 4‐PBA group. In summary, HFD‐induced ERS activation facilitates atrial fibrosis and AF. The inhibition of ERS might alleviate atrial fibrosis and reduce the incidence of AF‐associated obesity.