Effect of dietary Agriophyllum squarrosum on average daily gain, meat quality and muscle fatty acids in growing Tan lambs.

The herb A. squarrosum is reputed to possess medicinal properties for humans, and has the potential to be a feed resource for livestock. We hypothesized that this herb would improve the meat quality of lambs. To test this hypothesis, 24 Tan ewe-lambs (27.7 ± 0.45 kg) were offered diets containing 0 (CON), 100 (AS100), 200 (AS200) and 300 (AS300) g A. squarrosum/kg DM, and average daily gain, carcass traits, blood metabolites, meat quality and fatty acid profiles were determined. Drip loss % and cooking loss % decreased with the AS100 and AS200 diets (P < 0.05). Dietary A. squarrosum reduced muscle fiber area and diameter and increased density of the meat (P < 0.05), which indicated that the meat was more tender. The concentrations of C10:0 and C18:1n-9 t were 1ower and of C17:0 and C18:3n-3 were greater in the AS200 and AS300 treatments than CON (P < 0.05). Our results suggest that feeding lambs up to 200 g/kg DM of A. squarrosum can increase the water-holding capacity and L* value of meat without compromising growth. Further research is needed to determine the optimal level.

A novel TF molecular switch-mechanism found in two contrasting ecotypes of a psammophyte, Agriophyllum squarrosum, in regulating transcriptional drought memory.

BACKGROUND: Prior drought stress may change plants response patterns and subsequently increase their tolerance to the same condition, which can be referred to as "drought memory" and proved essential for plants well-being. However, the mechanism of transcriptional drought memory in psammophytes remains unclear. Agriophyllum squarrosum, a pioneer species on mobile dunes, is widely spread in Northern China's vast desert areas with outstanding ability of water use efficiency. Here we conducted dehydration-rehydration treatment on A. squarrosum semi-arid land ecotype AEX and arid land ecotype WW to dissect the drought memory mechanism of A. squarrosum, and to determine the discrepancy in drought memory of two contrasting ecotypes that had long adapted to water heterogeneity. RESULT: Physiological traits monitoring unveiled the stronger ability and longer duration in drought memory of WW than that of AEX. A total of 1,642 and 1,339 drought memory genes (DMGs) were identified in ecotype AEX and WW, respectively. Furthermore, shared DMGs among A. squarrosum and the previously studied species depicted that drought memory commonalities in higher plants embraced pathways like primary and secondary metabolisms; while drought memory characteristics in A. squarrosum were mainly related to response to heat, high light intensity, hydrogen peroxide, and dehydration, which might be due to local adaptation to desert circumstances. Heat shock proteins (HSPs) occupied the center of the protein-protein interaction (PPI) network in drought memory transcription factors (TF), thus playing a key regulatory role in A. squarrosum drought memory. Co-expression analysis of drought memory TFs and DMGs uncovered a novel regulating module, whereby pairs of TFs might function as molecular switches in regulating DMG transforming between high and low expression levels, thus promoting drought memory reset. CONCLUSION: Based on the co-expression analysis, protein-protein interaction prediction, and drought memory metabolic network construction, a novel regulatory module of transcriptional drought memory in A. squarrosum was hypothesized here, whereby recurrent drought signal is activated by primary TF switches, then amplified by secondary amplifiers, and thus regulates downstream complicated metabolic networks. The present research provided valuable molecular resources on plants' stress-resistance basis and shed light on drought memory in A. squarrosum.

Corrigendum: Integrated metabolomics and transcriptomics insights on flavonoid biosynthesis of a medicinal functional forage, Agriophyllum squarrosum (L.), based on a common garden trial covering six ecotypes.

[This corrects the article DOI: 10.3389/fpls.2022.985572.].

Integrated metabolomics and transcriptomics insights on flavonoid biosynthesis of a medicinal functional forage, Agriophyllum squarrosum (L.), based on a common garden trial covering six ecotypes.

Agriophyllum squarrosum (L.) Moq., well known as sandrice, is an important wild forage in sandy areas and a promising edible and medicinal resource plant with great domestication potential. Previous studies showed flavonoids are one of the most abundant medicinal ingredients in sandrice, whereby isorhamnetin and isorhamnetin-3-glycoside were the top two flavonols with multiple health benefits. However, the molecular regulatory mechanisms of flavonoids in sandrice remain largely unclear. Based on a common garden trial, in this study, an integrated transcriptomic and flavonoids-targeted metabolomic analysis was performed on the vegetative and reproductive periods of six sandrice ecotypes, whose original habitats covered a variety of environmental factor gradients. Multiple linear stepwise regression analysis unveiled that flavonoid accumulation in sandrice was positively correlated with temperature and UVB and negatively affected by precipitation and sunshine duration, respectively. Weighted co-expression network analysis (WGCNA) indicated the bHLH and MYB transcription factor (TF) families might play key roles in sandrice flavonoid biosynthesis regulation. A total of 22,778 differentially expressed genes (DEGs) were identified between ecotype DL and ecotype AEX, the two extremes in most environmental factors, whereby 85 DEGs could be related to known flavonoid biosynthesis pathway. A sandrice flavonoid biosynthesis network embracing the detected 23 flavonoids in this research was constructed. Gene families Plant flavonoid O-methyltransferase (AsPFOMT) and UDP-glucuronosyltransferase (AsUGT78D2) were identified and characterized on the transcriptional level and believed to be synthases of isorhamnetin and isorhamnetin-3-glycoside in sandrice, respectively. A trade-off between biosynthesis of rutin and isorhamnetin was found in the DL ecotype, which might be due to the metabolic flux redirection when facing environmental changes. This research provides valuable information for understanding flavonoid biosynthesis in sandrice at the molecular level and laid the foundation for precise development and utilization of this functional resource forage.

Transcriptome profiling reveals that foliar water uptake occurs with C3 and crassulacean acid metabolism facultative photosynthesis in Tamarix ramosissima under extreme drought.

Tamarix ramosissima is a typical desert plant species that is widely distributed in the desert areas of Northwest China. It plays a significant role in sand fixation and soil water conservation. In particular, how it uses water to survive in the desert plays an important role in plant growth and ecosystem function. Previous studies have revealed that T. ramosissima can alleviate drought by absorbing water from its leaves under extreme drought conditions. To date, there is no clear molecular regulation mechanism to explain foliar water uptake (FWU). In the present study, we correlated diurnal meteorological data, sap flow and photosynthetic parameters to determine the physical and biological characteristics of FWU. Our results suggested that the lesser the groundwater, the easier it is for T. ramosissima to absorb water via the leaves. Gene ontology annotation and Kyoto Encyclopaedia of Genes and Genomes pathway analysis of the transcriptome profile of plants subjected to high humidity suggested that FWU was highly correlated to carbohydrate metabolism, energy transfer, pyruvate metabolism, hormone signal transduction and plant-pathogen interaction. Interestingly, as a C3 plant, genes such as PEPC, PPDK, MDH and RuBP, which are involved in crassulacean acid metabolism (CAM) photosynthesis, were highly upregulated and accompanied by FWU. Therefore, we proposed that in the case of sufficient water supply, C3 photosynthesis is used in T. ramosissima, whereas in cases of extreme drought, starch is degraded to provide CO2 for CAM photosynthesis to make full use of the water obtained via FWU and the water that was transported or stored to assimilating branches and stems. This study may provide not only an important theoretical foundation for FWU and conversion from C3 plants to CAM plants but also for engineering improved photosynthesis in high-yield drought-tolerant plants and mitigation of climate change-driven drought.

Variations in Flavonoid Metabolites Along Altitudinal Gradient in a Desert Medicinal Plant Agriophyllum squarrosum.

Agriophyllum squarrosum (L.) Moq., a pioneer plant endemic to the temperate deserts of Asia, could be domesticated into an ideal crop with outstanding ecological and medicinal characteristics. A previous study showed differential flavonoid accumulation between two in situ altitudinal ecotypes. To verify whether this accumulation was determined by environmental or genetic factors, we conducted flavonoid-targeted metabolic profiling among 14 populations of A. squarrosum collected from regions with different altitudes based on a common garden experiment. Results showed that the most abundant flavonoid in A. squarrosum was isorhamnetin (48.40%, 557.45 µg/g), followed by quercetin (13.04%, 150.15 µg/g), tricin (11.17%, 128.70 µg/g), isoquercitrin (7.59%, 87.42 µg/g), isovitexin (7.20%, 82.94 µg/g), and rutin (7.00%, 80.62 µg/g). However, based on a common garden at middle-altitude environment, almost none of the flavonoids was enriched in the high-altitude populations, and even some flavonoids, such as quercetin, tricin, and rutin, were significantly enriched in low-altitude populations. This phenomenon indicated that the accumulation of flavonoids was not a result of local adaptation to high altitude. Furthermore, association analysis with in situ environmental variables showed that the contents of quercetin, tricin, and rutin were strongly positively correlated with latitude, longitude, and precipitation gradients and negatively correlated with temperature gradients. Thus, we could conclude that the accumulations of flavonoids in A. squarrosum were more likely as a result of local adaption to environmental heterogeneity combined with precipitation and temperature other than high altitude. This study not only provides an example to understand the molecular ecological basis of pharmacognosy, but also supplies methodologies for developing a new industrial crop with ecological and agricultural importance.

Organic acid metabolites involved in local adaptation to altitudinal gradient in Agriophyllum squarrosum, a desert medicinal plant.

Agriophyllum squarrosum (L.) Moq., a pioneer plant endemic to the temperate deserts of Asia, could be domesticated into an ideal crop with outstanding ecological and medicinal characteristics. A previous study showed differential organic acid accumulation between two in situ altitudinal ecotypes. To verify whether this accumulation was determined by environmental or genetic factors, we conducted organic acid targeted metabolic profiling among 14 populations of A. squarrosum collected from regions with different altitudes based on a common garden experiment. Results showed that the most abundant organic acid in A. squarrosum was citric acid (96.03%, 2322.90 µg g-1). Association analysis with in situ environmental variables showed that salicylic acid content was positively correlated with altitudinal gradient. Considering the enrichment of salicylic acid and protocatechualdehyde in high-altitude populations based on the common garden experiment, and the high expression of their biosynthesis relative genes (i.e., PAL and C4H) in the in situ high-altitude ecotype, we propose that organic acid accumulation could be involved in local adaptation to high altitudes. This study not only addresses the molecular basis of local adaptation involving the accumulation of organic acids in the desert plant A. squarrosum but also provides a method to screen wild germplasms to mitigate the impact of global climate change.

Genomic Adaptive Evolution of Sand Rice (Agriophyllum squarrosum) and Its Implications for Desert Ecosystem Restoration.

Natural selection is a significant driver of population divergence and speciation of plants. Due to local adaptation to geographic regions with ecological gradients, plant populations harbored a wide range of adaptive genetic variation to enable them to survive the heterogeneous habitats. This is all the more necessary for desert plants, as they must tolerant more striking gradients of abiotic stresses. However, the genomic mechanism by which desert plants adapt to ecological heterogeneity remains unclear, which could help to guide the sustainability of desert ecosystems. Here, using restriction-site-associated DNA sequencing in 38 natural populations, we investigated the genomic divergence and environmental adaptation of sand rice, Agriophyllum squarrosum, an annual pioneer species that covers sand dunes in northern China. Population genetic structure analyses showed that sand rice could be divided into three geographically distinct lineages, namely, Northwest, Central, and East. Phylogeographic analyses revealed that the plant might originate locally in Bergen County and further differentiated into the East lineage and then the Central lineage. Ecological niche modeling found that different lineages occupied distinct ecological niches, suggesting that the ecological gradient would have triggered genomic differentiation among sand rice lineages. Ecological association study supported that the three SNPs under divergent selection were closely correlated with precipitation gradients, indicating that precipitation might be the most important stress trigger for lineage diversity in sand rice. These adaptive SNPs could be used to genotype suitable germplasms for the ecological restoration of specific desertified lands. Further analyses found that genetic structure could significantly overestimate the signals for balancing selection. Within the Central lineage, we still found that 175 SNPs could be subject to balancing selection, which could be the means by which sand rice maintains genetic diversity and adapts to multiple stresses across heterogeneous deserts and sandy lands. From a genomic point of view, this study highlighted the local and global adaptation patterns of a desert plant to extreme and heterogeneous habitats. Our data provide molecular guidance for the restoration of desertified lands in the arid and semi-arid regions of China and could facilitate the marker assistant breeding of this potential crop to mitigate climate change.

Leaf size variations in a dominant desert shrub, Reaumuria soongarica, adapted to heterogeneous environments.

The climate in arid Central Asia (ACA) has changed rapidly in recent decades, but the ecological consequences of this are far from clear. To predict the impacts of climate change on ecosystem functioning, greater attention should be given to the relationships between leaf functional traits and environmental heterogeneity. As a dominant constructive shrub widely distributed in ACA, Reaumuria soongarica provided us with an ideal model to understand how leaf functional traits of desert ecosystems responded to the heterogeneous environments of ACA. Here, to determine the influences of genetic and ecological factors, we characterized species-wide variations in leaf traits among 30 wild populations of R. soongarica and 16 populations grown in a common garden. We found that the leaf length, width, and leaf length to width ratio (L/W) of the northern lineage were significantly larger than those of other genetic lineages, and principal component analysis based on the in situ environmental factors distinguished the northern lineage from the other lineages studied. With increasing latitude, leaf length, width, and L/W in the wild populations increased significantly. Leaf length and L/W were negatively correlated with altitude, and first increased and then decreased with increasing mean annual temperature (MAT) and mean annual precipitation (MAP). Stepwise regression analyses further indicated that leaf length variation was mainly affected by latitude. However, leaf width was uncorrelated with altitude, MAT, or MAP. The common garden trial showed that leaf width variation among the eastern populations was caused by both local adaptation and phenotypic plasticity. Our findings suggest that R. soongarica preferentially changes leaf length to adjust leaf size to cope with environmental change. We also reveal phenotypic evidence for ecological speciation of R. soongarica. These results will help us better understand and predict the consequences of climate change for desert ecosystem functioning.

The East Asian Winter Monsoon Acts as a Major Selective Factor in the Intraspecific Differentiation of Drought-Tolerant Nitraria tangutorum in Northwest China.

The influence of Quaternary climate fluctuation on the geographical structure and genetic diversity of species distributed in the regions of the Qinghai-Tibet Plateau (QTP) has been well established. However, the underlying role of the East Asian monsoon system (EAMS) in shaping the genetic structure of the population and the demography of plants located in the arid northwest of China has not been explored. In the present study, Nitraria tangutorum, a drought-tolerant desert shrub that is distributed in the EAMS zone and has substantial ecological and economic value, was profiled to better understand the influence of EAMS evolution on its biogeographical patterns and demographic history. Thus, the phylogeographical structure and historical dynamics of this plant species were elucidated using its five chloroplast DNA (cpDNA) fragments. Hierarchical structure analysis revealed three distinct, divergent lineages: West, East-A, and East-B. The molecular dating was carried out using a Bayesian approach to estimate the time of intraspecies divergence. Notably, the eastern region, which included East-A and East-B lineages, was revealed to be the original center of distribution and was characterized by a high level of genetic diversity, with the intraspecific divergence time dated to be around 2.53 million years ago (Ma). These findings, combined with the data obtained by ecological niche modeling analysis, indicated that the East lineages have undergone population expansion and differentiation, which were closely correlated with the development of the EAMS, especially the East Asian winter monsoon (EAWM). The West lineage appears to have originated from the migration of N. tangutorum across the Hexi corridor at around 1.85 Ma, and subsequent colonization of the western region. These results suggest that the EAWM accelerated the population expansion of N. tangutorum and subsequent intraspecific differentiation. These findings collectively provide new information on the impact of the evolution of the EAMS on intraspecific diversification and population demography of drought-tolerant plant species in northwest China.

Adaptive signals of flowering time pathways in wild barley from Israel over 28 generations.

Flowering time is one of the most critical traits for plants' life cycles, which is influenced by various environment changes, such as global warming. Previous studies have suggested that to guarantee reproductive success, plants have shifted flowering times to adapt to global warming. Although many studies focused on the molecular mechanisms of early flowering, little was supported by the repeated sampling at different time points through the changing climate. To fully dissect the temporal and spatial evolutionary genetics of flowering time, we investigated nucleotide variation in ten flowering time candidate genes and nine reference genes for the same ten wild-barley populations sampled 28 years apart (1980-2008). The overall genetic differentiation was significantly greater in the descendant populations (2008) compared with the ancestral populations (1980); however, local adaptation tests failed to detect any single-nucleotide polymorphism (SNP)/indel under spatial-diversifying selection at either time point. By contrast, the WFABC (Wright-Fisher ABC-based approach) that detected 54 SNPs/indels was under strong selection during the past 28 generations. Moreover, all these 54 alleles were segregated in the ancestral populations, but fixed in the descendent populations. Among the top ten SNPs/indels, seven were located in genes of FT1 (FLOWERING TIME LOCUS T 1), CO1 (CONSTANS-LIKE PROTEIN 1), and VRN-H2 (VERNALIZATION-H2), which have been documented to be associated with flowering time regulation in barley cultivars. This study might suggest that all ten populations have undergone parallel evolution over the past few decades in response to global warming, and even an overwhelming local adaptation and ecological differentiation.

Transcriptomes Divergence of Ricotia lunaria Between the Two Micro-Climatic Divergent Slopes at "Evolution Canyon" I, Israel.

As one of the hotspot regions for sympatric speciation studies, Evolution Canyon (EC) became an ideal place for its high level of microclimatic divergence interslopes. In this study, to highlight the genetic mechanisms of sympatric speciation, phenotypic variation on flowering time and transcriptomic divergence were investigated between two ecotypes of Ricotia lunaria, which inhabit the opposite temperate and tropical slopes of EC I (Lower Nahal Oren, Mount Carmel, Israel) separated by 100 m at the bottom of the slopes. Growth chamber results showed that flowering time of the ecotype from south-facing slope population # 3 (SFS 3) was significantly 3 months ahead of the north-facing slope population # 5 (NFS 5). At the same floral development stage, transcriptome analysis showed that 1,064 unigenes were differentially expressed between the two ecotypes, which enriched in the four main pathways involved in abiotic and/or biotic stresses responses, including flavonoid biosynthesis, α-linolenic acid metabolism, plant-pathogen interaction and linoleic acid metabolism. Furthermore, based on Ka/Ks analysis, nine genes were suggested to be involved in the ecological divergence between the two ecotypes, whose homologs functioned in RNA editing, ABA signaling, photoprotective response, chloroplasts protein-conducting channel, and carbohydrate metabolism in Arabidopsis thaliana. Among them, four genes, namely, SPDS1, FCLY, Tic21 and BGLU25, also showed adaptive divergence between R. lunaria and A. thaliana, suggesting that these genes could play an important role in plant speciation, at least in Brassicaceae. Based on results of both the phenotype of flowering time and comparative transcriptome, we hypothesize that, after long-time local adaptations to their interslope microclimatic environments, the molecular functions of these nine genes could have been diverged between the two ecotypes. They might differentially regulate the expression of the downstream genes and pathways that are involved in the interslope abiotic stresses, which could further diverge the flowering time between the two ecotypes, and finally induce the reproductive isolation establishment by natural selection overruling interslope gene flow, promoting sympatric speciation.

A whole-transcriptome approach to evaluate reference genes for quantitative diurnal gene expression studies under natural field conditions in Tamarix ramosissima leaves

Background: Tamarix ramosissima is a desert forest tree species that is widely distributed in the drought-stricken areas to sustain the fragile ecosystem. Owing to its wide usage in the desert restoration of Asia, it can be used as an ecophysiological model plant. To obtain reliable and accurate results, a set of reference genes should be screened before gene expression. However, up to date, systematical evaluation of reference genes has not been conducted in T. ramosissima. Results: In this study, we used eigenvalues derived from principal component analysis to identify stable expressed genes from 72,035 unigenes from diurnal transcriptomes under natural field conditions. With combined criteria of read counts above 900 and CV of FPKM below 0.3, a total of 7385 unigenes could be qualified as candidate reference genes in T. ramosissima. By using three statistical algorithm packages, geNorm, NormFinder, and BestKeeper, the stabilities of these novel reference genes were further compared with a panel of traditional reference genes. The expression patterns of three aquaporins (AQPs) suggested that at least UBQ (high expression), EIF4A2 (low expression), and GAPDH (moderate expression) could be qualified as ideal reference genes in both RT-PCR and RNA-seq analysis of T. ramosissima. Conclusions: This work will not only facilitate future studies on gene expression and functional analysis of genetic resources of desert plants but also improve our understanding of the molecular regulation of water transport in this plant, which could provide a new clue to further investigate the drought adaptation mechanism of desert plant species under harsh environments.

SNP Discovery and Genetic Variation of Candidate Genes Relevant to Heat Tolerance and Agronomic Traits in Natural Populations of Sand Rice (Agriophyllum squarrosum).

The extreme stress tolerance and high nutritional value of sand rice (Agriophyllum squarrosum) make it attractive for use as an alternative crop in response to concerns about ongoing climate change and future food security. However, a lack of genetic information hinders understanding of the mechanisms underpinning the morphological and physiological adaptations of sand rice. In the present study, we sequenced and analyzed the transcriptomes of two individuals representing semi-arid [Naiman (NM)] and arid [Shapotou (SPT)] sand rice genotypes. A total of 105,868 pairwise single nucleotide polymorphisms (SNPs) distributed in 24,712 Unigenes were identified among SPT and NM samples; the average SNP frequency was 0.3% (one SNP per 333 base pair). Characterization of gene annotation demonstrated that variations in genes involved in DNA recombination were associated with the survival of the NM population in the semi-arid environment. A set of genes predicted to be relevant to heat stress response and agronomic traits was functionally annotated using the accumulated knowledge from Arabidopsis and several crop plants, including rice, barley, maize, and sorghum. Four candidate genes related to heat tolerance (heat-shock transcription factor, HsfA1d), seed size (DA1-Related, DAR1), and flowering (early flowering 3, ELF3 and late elongated hypocotyl, LHY) were subjected to analysis of the genetic diversity in 10 natural populations, representing the core germplasm resource across the area of sand rice distribution in China. Only one SNP was detected in each of HsfA1d and DAR1, among 60 genotypes, with two in ELF3 and four in LHY. Nucleotide diversity ranged from 0.00032 to 0.00118. Haplotype analysis indicated that the NM population carried a specific allele for all four genes, suggesting that divergence has occurred between NM and other populations. These four genes could be further analyzed to determine whether they are associated with phenotype variation and identify alleles favorable for sand rice breeding.

Population dynamics of Agriophyllum squarrosum, a pioneer annual plant endemic to mobile sand dunes, in response to global climate change.

Climate change plays an important role in the transition of ecosystems. Stratigraphic investigations have suggested that the Asian interior experienced frequent transitions between grassland and desert ecosystems as a consequence of global climate change. Using maternally and bi-parentally inherited markers, we investigated the population dynamics of Agriophyllum squarrosum (Chenopodiaceae), an annual pioneer plant endemic to mobile sand dunes. Phylogeographic analysis revealed that A. squarrosum could originate from Gurbantunggut desert since ~1.6 Ma, and subsequently underwent three waves of colonisation into other deserts and sandy lands corresponding to several glaciations. The rapid population expansion and distribution range shifts of A. squarrosum from monsoonal climate zones suggested that the development of the monsoonal climate significantly enhanced the population growth and gene flow of A. squarrosum. These data also suggested that desertification of the fragile grassland ecosystems in the Qinghai-Tibetan Plateau was more ancient than previously suggested and will be aggravated under global warming in the future. This study provides new molecular phylogeographic insights into how pioneer annual plant species in desert ecosystems respond to global climate change, and facilitates evaluation of the ecological potential and genetic resources of future crops for non-arable dry lands to mitigate climate change.

Characterization of the complete mitochondrial genome of Collared Finchbill, Spizixos semitorques (Passeriformes: Pycnonotidae).

The complete mitochondrial genome of Spizixos semitorques is 17,045 bp long with an A+T content of 55.52%, which composed of 13 typical protein-coding genes, 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and 2 extensive heteroplasmic control regions (CR1 and CR2). Genome organization of S. semitorques is similar to that of other birds. All of the PCGs started with ATG, except for cox1, which began with GTG. Six types of termination codons were used, canonical TAA and TAG termination codons are found in seven protein-coding genes, while AGA for nad1 and nad5, AGG for cox1, and the remaining genes (nad2, cox3, nad4) have incomplete termination codons (TA- or T- -). All tRNA genes have a typical clover-leaf structure. Furthermore, a minisatellite of 24 nucleotides (5'-CACAAAAATCAAACAAAAATACAA-3') with 2.9 tandem repeats at the end of CR2 were detected.

Complete mitochondrial genome of the Black-capped Bulbul, Pycnonotus melanicterus (Passeriformes: Pycnonotidae).

The first complete maternally-inherited mitochondrial genome of Pycnonotus melanicterus has been sequenced and annotated in this study. The mitogenome is 17,011 bp in size with an A+T content of 55.15%, consisting of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (srRNA and lrRNA) genes, and 2 extensive heteroplasmic control regions. Genome organization and gene arrangements of P. melanicterus are similar to those of other bulbuls. All 13 PCGs initiated with ATG as the start codon, while terminated by seven types of stop codons. All tRNA genes sequences have the potential to fold into typical cloverleaf secondary structures. Gene overlapping was observed in the P. melanicterus mitogenome. Moreover, our results of phylogenetic analysis well support the monophyly of Pycnonotidae as previously studies.

The role of East Asian monsoon system in shaping population divergence and dynamics of a constructive desert shrub Reaumuria soongarica.

Both of the uplift of Qinghai-Tibet Plateau (QTP) and the development of East Asian monsoon system (EAMS) could have comprehensively impacted the formation and evolution of Arid Central Asia (ACA). To understand how desert plants endemic to ACA responded to these two factors, we profiled the historical population dynamics and distribution range shift of a constructive desert shrub Reaumuria soongarica (Tamaricaceae) based on species wide investigation of sequence variation of chloroplast DNA and nuclear ribosomal ITS. Phylogenetic analysis uncovered a deep divergence occurring at ca. 2.96 Mya between the western and eastern lineages of R. soongarica, and ecological niche modeling analysis strongly supported that the monsoonal climate could have fragmented its habitats in both glacial and interglacial periods and impelled its intraspecific divergence. Additionally, the population from the east monsoonal zone expanded rapidly, suggesting that the local monsoonal climate significantly impacted its population dynamics. The isolation by distance tests supported strong maternal gene flow along the direction of the East Asian winter monsoon, whose intensification induced the genetic admixture along the latitudinal populations of R. soongarica. Our results presented a new case that the development of EAMS had prominently impacted the intraspecific divergence and population dynamics of this desert plant.

The Mitochondrial Genomes of Aquila fasciata and Buteo lagopus (Aves, Accipitriformes): Sequence, Structure and Phylogenetic Analyses.

The family Accipitridae is one of the largest groups of non-passerine birds, including 68 genera and 243 species globally distributed. In the present study, we determined the complete mitochondrial sequences of two species of accipitrid, namely Aquila fasciata and Buteo lagopus, and conducted a comparative mitogenome analysis across the family. The mitogenome length of A. fasciata and B. lagopus are 18,513 and 18,559 bp with an A + T content of 54.2% and 55.0%, respectively. For both the two accipitrid birds mtDNAs, obvious positive AT-skew and negative GC-skew biases were detected for all 12 PCGs encoded by the H strand, whereas the reverse was found in MT-ND6 encoded by the L strand. One extra nucleotide'C'is present at the position 174 of MT-ND3 gene of A. fasciata, which is not observed at that of B. lagopus. Six conserved sequence boxes in the Domain II, named boxes F, E, D, C, CSBa, and CSBb, respectively, were recognized in the CRs of A. fasciata and B. lagopus. Rates and patterns of mitochondrial gene evolution within Accipitridae were also estimated. The highest dN/dS was detected for the MT-ATP8 gene (0.32493) among Accipitridae, while the lowest for the MT-CO1 gene (0.01415). Mitophylogenetic analysis supported the robust monophyly of Accipitriformes, and Cathartidae was basal to the balance of the order. Moreover, we performed phylogenetic analyses using two other data sets (two mitochondrial loci, and combined nuclear and mitochondrial loci). Our results indicate that the subfamily Aquilinae and all currently polytypic genera of this subfamily are monophyletic. These two novel mtDNA data will be useful in refining the phylogenetic relationships and evolutionary processes of Accipitriformes.