Potential Positive Selection Research Articles

Assembly and comparative analysis of the multichromosomal mitochondrial genome of globally endangered seagrass Halophila beccarii

BackgroundHalophila beccarii is one of the oldest two generations of seagrass plants and one of the 10 species of seagrass currently at risk of extinction worldwide. Therefore, how to effectively protect the H. beccarii resources from extinction is a huge challenge. Molecular biology research can provide a scientific basis for species conservation. So far, there has been no detailed analysis of the mitochondrial genome of the genus Halophila.ResultsThe mitochondrial genome of H. beccarii was assembled into 28 circular chromosomes, ranging in length from 41,738 bp to 104,744 bp, with a total length of 1,964,072 bp and a GC content of 46.71%. It contains 39 genes, including 26 protein coding genes, 10 tRNA genes, and 3 rRNA genes. Repeat sequence analysis and prediction of RNA editing sites revealed a total of 850 dispersed repeats, 1,205 simple repeats, 61 tandem repeats, and 120 RNA editing sites. Analysis of codon usage indicates that codons ending in A/U are preferred. Gene migration between the mitochondrial genome and the chloroplast genome was observed through homologous fragment detection. In addition, Ka/Ks analysis showed that most protein coding genes in the mitochondrial genome experienced negative selection, while only the nad3 gene experienced potential positive selection in most Alismatales. Nucleotide polymorphism analysis revealed variations in each gene, with rpl10 being the most significant. In addition, comparative analysis shows that the GC content is conserved, but there are significant differences in the size and structure of mitochondrial genomes among different species of Alismatales. The phylogenetic analysis based on the mitochondrial genome reflects the exact evolutionary and taxonomic status of H. beccarii.ConclusionIn this study, we sequenced and annotated the mitochondrial genome of H. beccarii, and compared it with the mitochondrial genomes of other plants in Alismatales. Our findings enrich the mitogenome database of seagrass plants and highlight the potential for mitochondrial genes to help decipher plant evolutionary history.

The Complete Chloroplast Genome Sequence of the Medicinal Moss Rhodobryum giganteum (Bryaceae, Bryophyta): Comparative Genomics and Phylogenetic Analyses.

Rhodobryum giganteum (Bryaceae, Bryophyta), a rare medicinal bryophyte, is valued for its cardiovascular therapeutic properties in traditional Chinese medicine. This study presents the first complete chloroplast genome sequence of R. giganteum, including its assembly and annotation. The circular chloroplast genome of R. giganteum is 124,315 bp in length, displaying a typical quadripartite structure with 128 genes: 83 protein-coding genes, 37 tRNAs, and 8 rRNAs. Analyses of codon usage bias, repetitive sequences, and simple sequence repeats (SSRs) revealed an A/U-ending codon preference, 96 repetitive sequences, and 385 SSRs in the R. giganteum chloroplast genome. Nucleotide diversity analysis identified 10 high mutational hotspots. Ka/Ks ratio analysis suggested potential positive selection in rpl20, rps18, petG, and psbM genes. Phylogenetic analysis of whole chloroplast genomes from 38 moss species positioned R. giganteum within Bryales, closely related to Rhodobryum laxelimbatum. This study augments the chloroplast genomic data for Bryales and provides a foundation for molecular marker development and genetic diversity analyses in medicinal bryophytes.

Evolutionary dynamics of canine kobuvirus in Vietnam and Thailand reveal the evidence of viral ability to evade host immunity

Canine kobuvirus (CaKoV) is a pathogen associated with canine gastrointestinal disease (GID). This study examined 327 rectal swabs (RS), including 113 from Vietnam (46 healthy, 67 with GID) and 214 from Thailand (107 healthy and 107 with GID). CaKoV was detected in both countries, with prevalences of 28.3% (33/113) in Vietnam and 7.9% (17/214) in Thailand. Additionally, CaKoV was found in both dogs with diarrhea and healthy dogs. CaKoV was mainly found in puppies under six months of age (30.8%). Co-detection with other canine viruses were also observed. The complete coding sequence (CDS) of nine Vietnamese and four Thai CaKoV strains were characterized. Phylogenetic analysis revealed a close genetic relationship between Vietnamese and Thai CaKoV strains, which were related to the Chinese strains. CDS analysis indicated a distinct lineage for two Vietnamese CaKoV strains. Selective pressure analysis on the viral capsid (VP1) region showed negative selection, with potential positive selection sites on B-cell epitopes. This study, the first of its kind in Vietnam, provides insights into CaKoV prevalence in dogs of different ages and healthy statuses, updates CaKoV occurrence in Thailand, and sheds light on its molecular characteristics and immune evasion strategies.

Bona Fide Plant Steroid Receptors are Innovated in Seed Plants and Angiosperms through Successive Whole-Genome Duplication Events.

Whole-genome duplication (WGD) events are widespread in plants and animals, thus their long-term evolutionary contribution has long been speculated, yet a specific contribution is difficult to verify. Here, we show that ɛ-WGD and ζ-WGD contribute to the origin and evolution of bona fide brassinosteroid (BR) signaling through the innovation of active BR biosynthetic enzymes and active BR receptors from their respective ancestors. We found that BR receptors BRI1 (BR INSENSITIVE 1) and BRL1/3 (BRI1-LIKES 1/3) derived by ɛ-WGD and ζ-WGD, which occurred in the common ancestor of angiosperms and seed plants, respectively, while orphan BR receptor BRL2 first appeared in stomatophytes. Additionally, CYP85A enzymes synthesizing the bioactive BRs derived from a common ancestor of seed plants, while its sister enzymes CYP90 synthesizing BR precursors presented in all land plants, implying possible ligand-receptor coevolution. Consistently, the island domains (IDs) responsible for BR perception in BR receptors were most divergent among different receptor branches, supporting ligand-driven evolution. As a result, BRI1 was the most diversified BR receptor in angiosperms. Importantly, relative to the BR biosynthetic DET2 gene presented in all land plants, BRL2, BRL1/3 and BRI1 had high expression in vascular plants ferns, gymnosperms and angiosperms, respectively. Notably, BRI1 is the most diversified BR receptor with the most abundant expression in angiosperms, suggesting potential positive selection. Therefore, WGDs initiate a neofunctionalization process diverged by ligand-perception and transcriptional expression, which might optimize both BR biosynthetic enzymes and BR receptors, likely contributing to the evolution of land plants, especially seed plants and angiosperms.

Identification of genomic diversity and selection signatures in Luxi cattle using whole-genome sequencing data.

The objective of this study was to investigate the genetic diversity, population structure and whole-genome selection signatures of Luxi cattle to reveal its genomic characteristics in terms of meat and carcass traits, skeletal muscle development, body size, and other traits. To further analyze the genomic characteristics of Luxi cattle, this study sequenced the whole-genome of 16 individuals from the core conservation farm in Shandong region, and collected 174 published genomes of cattle for conjoint analysis. Furthermore, three different statistics (pi, Fst, and XP-EHH) were used to detect potential positive selection signatures related to selection in Luxi cattle. Moreover, gene ontology and Kyoto encyclopedia of genes and genomes pathway enrichment analyses were performed to reveal the potential biological function of candidate genes harbored in selected regions. The results showed that Luxi cattle had high genomic diversity and low inbreeding levels. Using three complementary methods (pi, Fst, and XP-EHH) to detect the signatures of selection in the Luxi cattle genome, there were 2,941, 2,221 and 1,304 potentially selected genes identified, respectively. Furthermore, there were 45 genes annotated in common overlapping genomic regions covered 0.723 Mb, including PLAG1 zinc finger (PLAG1), dedicator of cytokinesis 3 (DOCK3), ephrin A2 (EFNA2), DAZ associated protein 1 (DAZAP1), Ral GTPase activating protein catalytic subunit alpha 1 (RALGAPA1), mediator complex subunit 13 (MED13), and decaprenyl diphosphate synthase subunit 2 (PDSS2), most of which were enriched in pathways related to muscle growth and differentiation and immunity. In this study, we provided a series of genes associated with important economic traits were found in positive selection regions, and a scientific basis for the scientific conservation and genetic improvement of Luxi cattle.

Mutational profiling of mitochondrial DNA reveals an epithelial ovarian cancer-specific evolutionary pattern contributing to high oxidative metabolism.

Epithelial ovarian cancer (EOC) heavily relies on oxidative phosphorylation (OXPHOS) and exhibits distinct mitochondrial metabolic reprogramming. Up to now, the evolutionary pattern of somatic mitochondrial DNA (mtDNA) mutations in EOC tissues and their potential roles in metabolic remodelling have not been systematically elucidated. Based on a large somatic mtDNA mutation dataset from private and public EOC cohorts (239 and 118 patients, respectively), we most comprehensively characterised the EOC-specific evolutionary pattern of mtDNA mutations and investigated its biological implication. Mutational profiling revealed that the mitochondrial genome of EOC tissues was highly unstable compared with non-cancerous ovary tissues. Furthermore, our data indicated the delayed heteroplasmy accumulation of mtDNA control region (mtCTR) mutations and near-complete absence of mtCTR non-hypervariable segment (non-HVS) mutations in EOC tissues, which is consistent with stringent negative selection against mtCTR mutation. Additionally, we observed a bidirectional and region-specific evolutionary pattern of mtDNA coding region mutations, manifested as significant negative selection against mutations in complex V (ATP6/ATP8) and tRNA loop regions, and potential positive selection on mutations in complex III (MT-CYB). Meanwhile, EOC tissues showed higher mitochondrial biogenesis compared with non-cancerous ovary tissues. Further analysis revealed the significant association between mtDNA mutations and both mitochondrial biogenesis and overall survival of EOC patients. Our study presents a comprehensive delineation of EOC-specific evolutionary patterns of mtDNA mutations that aligned well with the specific mitochondrial metabolic remodelling, conferring novel insights into the functional roles of mtDNA mutations in EOC tumourigenesis and progression.

Putative protective genomic variation in the Lithuanian population.

Genomic effect variants associated with survival and protection against complex diseases vary between populations due to microevolutionary processes. The aim of this study was to analyse diversity and distribution of effect variants in a context of potential positive selection. In total, 475 individuals of Lithuanian origin were genotyped using high-throughput scanning and/or sequencing technologies. Allele frequency analysis for the pre-selected effect variants was performed using the catalogue of single nucleotide polymorphisms. Comparison of the pre-selected effect variants with variants in primate species was carried out to ascertain which allele was derived and potentially of protective nature. Recent positive selection analysis was performed to verify this protective effect. Four variants having significantly different frequencies compared to European populations were identified while two other variants reached borderline significance. Effect variant in SLC30A8 gene may potentially protect against type 2 diabetes. The existing paradox of high rates of type 2 diabetes in the Lithuanian population and the relatively high frequencies of potentially protective genome variants against it indicate a lack of knowledge about the interactions between environmental factors, regulatory regions, and other genome variation. Identification of effect variants is a step towards better understanding of the microevolutionary processes, etiopathogenetic mechanisms, and personalised medicine.

A whole genome scan reveals distinct features of selection in Zhaotong cattle of Yunnan province.

Over the years, indigenous cattle have not only played an essential role in securing primary food sources but have also been utilized for labor by humans, making them invaluable genetic resources. The Zhaotong cattle, a native Chinese breed from the Yunnan province, possess excellent meat quality and resistance to heat and humidity. Here we used whole genome sequencing data of 104 animals to delve into the population structure, genomic diversity and potential positive selection signals in Zhaotong cattle. The findings of this study demonstrate that the genetic composition of Zhaotong cattle was primarily derived from Chinese indicine cattle and East Asian cattle. The nucleotide diversity of Zhaotong cattle was only lower than that of Chinese indicine cattle, which was much higher than that of other taurine cattle. Genome-wide selection scans detected a series of positive candidate regions containing multiple key genes related to bone development and metabolism (CA10, GABRG3, GLDN and NOTUM), meat quality traits (ALG8, LINGO2, MYO5B, PRKG1 and GABRB1), immune response (ADA2, BMF, LEF1 and PAK6) and heat resistance (EIF2AK4 and LEF1). In summary, this study supplies essential genetic insights into the genome diversity within Zhaotong cattle and provides a foundational framework for comprehending the genetic basis of indigenous cattle breeds.

Assembly and comparative analysis of the complete mitochondrial genome of Ilex metabaptista (Aquifoliaceae), a Chinese endemic species with a narrow distribution

BackgroundIlex metabaptista is a woody tree species with strong waterlogging tolerance and is also admired as a landscape plant with high development prospects and scientific research value. Unfortunately, populations of this species have declined due to habitat loss. Thus, it is a great challenge for us to efficiently protect I. metabaptista resources from extinction. Molecular biology research can provide the scientific basis for the conservation of species. However, the study of I. metabaptista genetics is still in its infancy. To date, no mitochondrial genome (mitogenome) in the genus Ilex has been analysed in detail.ResultsThe mitogenome of I. metabaptista was assembled based on the reads from Illumina and Nanopore sequencing platforms; it was a typical circular DNA molecule of 529,560 bp with a GC content of 45.61% and contained 67 genes, including 42 protein-coding genes, 22 tRNA genes, and 3 rRNA genes. Repeat sequence analysis and prediction of RNA editing sites revealed a total of 286 dispersed repeats, 140 simple repeats, 18 tandem repeats, and 543 RNA editing sites. Analysis of codon usage showed that codons ending in A/T were preferred. Gene migration was observed to occur between the mitogenome and chloroplast genome via the detection of homologous fragments. In addition, Ka/Ks analysis revealed that most of the protein-coding genes in the mitogenome had undergone negative selection, and only the ccmB gene had undergone potential positive selection in most asterids. Nucleotide polymorphism analysis revealed the variation in each gene, with atp9 being the most notable. Furthermore, comparative analysis showed that the GC contents were conserved, but the sizes and structure of mitogenomes varied greatly among asterids. Phylogenetic analysis based on the mitogenomes reflected the exact evolutionary and taxonomic status of I. metabaptista.ConclusionIn this study, we sequenced and annotated the mitogenome of I. metabaptista and compared it with the mitogenomes of other asterids, which provided essential background information for further understanding of the genetics of this plant and helped lay the foundation for future studies on molecular breeding of I. metabaptista.

Hybridization in late stages of speciation: Strong but incomplete genome-wide reproductive isolation and 'large Z-effect' in a moving hybrid zone.

In organisms reproducing sexually, speciation occurs when increasing divergence results in pre- or post-zygotic reproductive isolation between lineages. Studies focusing on reproductive isolation origin in early stages of speciation are common and many rely on genomic scans to infer introgression providing limited information on the genomic architecture of reproductive isolation long-term maintenance. This study analyses a natural hybrid zone between two species in a late stage of speciation. We used ddRADseq genotyping in the contact between Podarcis bocagei and P. carbonelli to examine admixture extent, analyse hybrid zone stability and assess genome-wide variation in selection against introgression. We confirmed strong but incomplete reproductive isolation in a bimodal hybrid zone. New findings revealed population genetic structure within P. carbonelli in the contact zone; geographical and genomic clines analysis suggested strong selection against gene flow, but a relatively small proportion of the loci can introgress, mostly within the narrow contact zone. However, geographical clines revealed that a few introgressed loci show signs of potential positive selection, particularly into P. bocagei. Geographical clines also detected a signal of hybrid zone movement towards P. bocagei distribution. Genomic cline analysis revealed heterogeneous patterns of introgression among loci within the syntopy zone, but the majority maintain a strong association with the genomic background of origin. However, incongruences between both cline approaches were found, potentially driven by confounding effects on genomic clines. Last, an important role of the Z chromosome in reproductive isolation is suggested. Importantly, overall patterns of restricted introgression seem to result from numerous strong intrinsic barriers across the genome.

Complete mitochondrial genome of Thuja sutchuenensis and its implications on evolutionary analysis of complex mitogenome architecture in Cupressaceae

BackgroundThe complex physical structure and abundant repeat sequences make it difficult to assemble the mitogenomes of seed plants, especially gymnosperms. Only approximately 33 mitogenomes of gymnosperms have been reported. However, as the most widely distributed and the second largest family among gymnosperms, Cupressaceae has only six assembled mitogenomes, including five draft mitogenomes and one complete mitogenome, which has greatly hindered the understanding of mitogenome evolution within this large family, even gymnosperms.ResultsIn this study, we assembled and validated the complete mitogenome of Thuja sutchuenensis, with a size of 2.4 Mb. Multiple sequence units constituted its complex structure, which can be reduced to three linear contigs and one small circular contig. The analysis of repeat sequences indicated that the numbers of simple sequence repeats increased during the evolutionary history of gymnosperms, and the mitogenome of Thuja sutchuenensis harboured abundant extra-long repeats (more than 5 kb). Additionally, the longest repeat sequence identified in these seven gymnosperms also came from the mitogenome of Thuja sutchuenensis, with a length of up to 47 kb. The analysis of colinear blocks and gene clusters both revealed that the orders of mitochondrial genes within gymnosperms was not conserved. The comparative analysis showed that only four tRNAs were shared by seven gymnosperms, namely, trnD-GUC, trnE-UUC, trnI-CAU and trnY-GUA. Furthermore, four genes have undergone potential positive selection in most gymnosperm species, namely, atp8, ccmB, mttB and sdh4.ConclusionWe successfully assembled the second complete mitogenome within Cupressaceae and verified that it consisted of multiple sequence units. Our study also indicated that abundant long repeats may contribute to the generation of the complex conformation of the mitogenome of Thuja sutchuenensis. The investigation of Thuja sutchuenensis’s mitogenome in our study provides new insight into further understanding the complex mitogenome architecture within gymnosperms.

SARS-CoV-2 Pandemic Tracing in Italy Highlights Lineages with Mutational Burden in Growing Subsets.

Tracing the appearance and evolution of virus variants is essential in the management of the COVID-19 pandemic. Here, we focus on SARS-CoV-2 spread in Italian patients by using viral sequences deposited in public databases and a tracing procedure which is used to monitor the evolution of the pandemic and detect the spreading, within the infected population of emergent sub-clades with a potential positive selection. Analyses of a collection of monthly samples focused on Italy highlighted the appearance and evolution of all the main viral sub-trees emerging at the end of the first year of the pandemic. It also identified additional expanding subpopulations which spread during the second year (i.e., 2021). Three-dimensional (3D) modelling of the main amino acid changes in mutated viral proteins, including ORF1ab (nsp3, nsp4, 2’-o-ribose methyltransferase, nsp6, helicase, nsp12 [RdRp]), N, ORF3a, ORF8, and spike proteins, shows the potential of the analysed structural variations to result in epistatic modulation and positive/negative selection pressure. These analyzes will be of importance to the early identification of emerging clades, which can develop into new “variants of concern” (i.e., VOC). These analyses and settings will also help SARS-CoV-2 coronet genomic centers in other countries to trace emerging worldwide variants.

Genomic Insight Into the Population Structure and Admixture History of Tai-Kadai-Speaking Sui People in Southwest China.

Sui people, which belong to the Tai-Kadai-speaking family, remain poorly characterized due to a lack of genome-wide data. To infer the fine-scale population genetic structure and putative genetic sources of the Sui people, we genotyped 498,655 genome-wide single-nucleotide polymorphisms (SNPs) using SNP arrays in 68 Sui individuals from seven indigenous populations in Guizhou province and Guangxi Zhuang Autonomous Region in Southwest China and co-analyzed with available East Asians via a series of population genetic methods including principal component analysis (PCA), ADMIXTURE, pairwise Fst genetic distance, f-statistics, qpWave, and qpAdm. Our results revealed that Guangxi and Guizhou Sui people showed a strong genetic affinity with populations from southern China and Southeast Asia, especially Tai-Kadai- and Hmong-Mien-speaking populations as well as ancient Iron Age Taiwan Hanben, Gongguan individuals supporting the hypothesis that Sui people came from southern China originally. The indigenous Tai-Kadai-related ancestry (represented by Li), Northern East Asian-related ancestry, and Hmong-Mien-related lineage contributed to the formation processes of the Sui people. We identified the genetic substructure within Sui groups: Guizhou Sui people were relatively homogeneous and possessed similar genetic profiles with neighboring Tai-Kadai-related populations, such as Maonan. While Sui people in Yizhou and Huanjiang of Guangxi might receive unique, additional gene flow from Hmong-Mien-speaking populations and Northern East Asians, respectively, after the divergence within other Sui populations. Sui people could be modeled as the admixture of ancient Yellow River Basin farmer-related ancestry (36.2–54.7%) and ancient coastal Southeast Asian-related ancestry (45.3–63.8%). We also identified the potential positive selection signals related to the disease susceptibility in Sui people via integrated haplotype score (iHS) and number of segregating sites by length (nSL) scores. These genomic findings provided new insights into the demographic history of Tai-Kadai-speaking Sui people and their interaction with neighboring populations in Southern China.

Detection of differential selection pressure and functional-specific sites in subunits of vertebrate neuronal nicotinic acetylcholine receptors

The nicotinic acetylcholine receptors (nAChR) are made of subunits evolved from a common ancestor. Despite the similarity in their sequences and structures, the properties of these subunits vary significantly. Thus, identifying the evolution features and function-related sites specific to each subunit is essential for understanding the characteristics of the subunits and the receptors assembled by them. In this study, we examined the sequence features of the nine neuronal nAChRs subunits from representative vertebrate species. Analysis revealed that all the subunits were subject to strong purifying selection in evolution, and each was under a unique pattern of selection pressures. At the same time, the functional constraints were not uniform within each subunit, with different domains in the molecule being subject to different selection pressures. We also detected potential positive selection events in the subunits or subunit clusters, and identified the sites might be associated with the function specificity of each subunit. Furthermore, positive selection at some domains might contribute to the diversity of subunit function; for example, the β9 strand might be related to the agonist specificity of α subunit in heteromeric receptor and β4-β5 linker could be involved in Ca2+ permeability. Subunits α7, α4 and β2 subunits possess a strong adaptability in vertebrates. Our results highlighted the importance of tracking functional differentiation in protein sequence underlying functional properties of nAChRs. In summary, our work may provide clues on understanding the diversity and the function specificity of the nAChR subunits, as well as the receptors co-assembled by them. Communicated by Ramaswamy H. Sarma

Genotyping and Molecular Characterization of Classical Swine Fever Virus Isolated in China during 2016–2018

Classical swine fever (CSF) is a highly contagious disease of swine caused by classical swine fever virus (CSFV). For decades the disease has been controlled in China by a modified live vaccine (C-strain) of genotype 1. The emergent genotype 2 strains have become predominant in China in the past years that are genetically distant from the vaccine strain. Here, we aimed to evaluate the current infectious status of CSF, and for this purpose 24 isolates of CSFV were identified from different areas of China during 2016–2018. Phylogenetic analysis of NS5B, E2 and full genome revealed that the new isolates were clustered into subgenotype 2.1d and 2.1b, while subgenotype 2.1d was predominant. Moreover, E2 and Erns displayed multiple variations in neutralizing epitope regions. Furthermore, the new isolates exhibited capacity to escape C-strain-derived antibody neutralization compared with the Shimen strain (genotype 1). Potential positive selection sites were identified in antigenic regions of E2 and Erns, which are related with antibody binding affinity. Recombination events were predicted in the new isolates with vaccine strains in the E2 gene region. In conclusion, the new isolates showed molecular variations and antigenic alterations, which provide evidence for the emergence of vaccine-escaping mutants and emphasize the need of updated strategies for CSF control.

Molecular and phylogenetic analysis of Vietnamese canine parvovirus 2C originated from dogs reveals a new Asia-IV clade.

Canine parvovirus type 2 (CPV-2) is a small, single-stranded DNA virus causing fatal haemorrhagic enteritis in dogs. Currently, CPV-2 is classified into CPV-2a, CPV-2b and CPV-2c based on genetic variation in the VP2 gene. The CPV-2c variant has become ubiquitous worldwide and gained attention for monitoring parvoviral evolution. In this study, we characterized the full-length genome sequences of CPV-2c strains obtained from 59 dogs in Vietnam. Molecular analysis revealed that Vietnamese CPV-2c shared a common evolutionary pattern with the Asian CPV-2 clade, which is marked by genetic signature patterns in the structural and nonstructural proteins. In addition, these Vietnamese CPV-2c strains exhibited unique Thr112Ile and Ile447Met mutations in the VP1 and VP2 sequence, respectively. Interestingly, phylogenetic analysis indicated that the mutations of amino acid residues in both the structural and nonstructural genes have contributed to the emergence of a new clade, designated here as the Asia-IV clade. The substitution rates, estimated from a dataset containing 199 sequences over the last 42years, confirmed that CPV-2 showed a high rate of nucleotide substitution, at about 2.49×10-4 nucleotide substitutions per site per year (nt/s/y), with VP1/2 and NS1/2 estimates of 3.06×10-4 and 3.16×10-4 nt/s/y, respectively. Even though no evidence of genetic recombination in these Vietnamese CPV-2c strains was established, potential positive selection sites were observed in both the structural and nonstructural genes, suggesting the viral evolutionary process has occurred in both the structural and nonstructural proteins. Genetic and evolutionary analysis of the full-length genome sequence is necessary to gain evolutionary insight of CPV-2.

Mitochondrial genome diversity and population mitogenomics of polar cod (Boreogadus saida) and Arctic dwelling gadoids

High-latitude fish typically exhibit a narrow thermal tolerance window, which may pose challenges when coping with temperatures that shift outside of a species’ range of tolerance. Due to its role in aerobic metabolism and energy balance, the mitochondrial genome is likely critical for the acclimation and adaptation to differing temperature regimes in marine ectotherms. As oceans continue to warm, there is growing need to understand the ability of organisms to respond to changing environmental conditions given evidence that some species, in particular cold-water species, may already be experiencing difficulties. To assess how Arctic gadids in Alaska have responded to differential thermal preferences in the past and how regions are interconnected, we sequenced complete mitochondrial genomes for four Arctic gadids to determine the distribution of mitochondrial diversity and population-level structure as well as to detect signatures of selection acting on the mitochondrial genome. We found little population-level structure within all four species with the clear exception of Gulf of Alaska saffron cod (Eleginus gracilis). Northern localities exhibited higher levels of genetic diversity and primarily northern lineages were observed within polar cod (Boreogadus saida) and saffron cod, likely reflecting asymmetrical dispersal and potentially admixture of distinct lineages via ocean currents. The main evolutionary force shaping the evolution of the mitogenome appears to be purifying selection, but we also identified potential positive selection of candidate amino acid replacements primarily in complex I (ND genes) in polar cod. The high levels of mitochondrial diversity observed in our study and large population size may provide this species with the ability to respond evolutionarily (i.e. long-term) to a changing environment.

Characterization and analysis of the transcriptome in Gymnocypris selincuoensis on the Qinghai-Tibetan Plateau using single-molecule long-read sequencing and RNA-seq.

The lakes on the Qinghai-Tibet Plateau (QTP) are the largest and highest lake group in the world. Gymnocypris selincuoensis is the only cyprinid fish living in lake Selincuo, the largest lake on QTP. However, its genetic resource is still blank, limiting studies on molecular and genetic analysis. In this study, the transcriptome of G. selincuoensis was first generated by using PacBio Iso-Seq and Illumina RNA-seq. A full-length (FL) transcriptome with 75,435 transcripts was obtained by Iso-Seq with N50 length of 3,870 bp. Among all transcripts, 75,016 were annotated to public databases, 64,710 contain complete open reading frames and 2,811 were long non-coding RNAs. Based on all- vs.-all BLAST, 2,069 alternative splicing events were detected, and 80% of them were validated by reverse transcription polymerase chain reaction (RT-PCR). Tissue gene expression atlas showed that the number of detected expressed transcripts ranged from 37,397 in brain to 19,914 in muscle, with 10,488 transcripts detected in all seven tissues. Comparative genomic analysis with other cyprinid fishes identified 77 orthologous genes with potential positive selection (Ka/Ks > 0.3). A total of 56,696 perfect simple sequence repeats were identified from FL transcripts. Our results provide valuable genetic resources for further studies on adaptive evolution, gene expression and population genetics in G. selincuoensis and other congeneric fishes.

Time matters! Developmental shift in gene expression between the head and the trunk region of the cichlid fish Astatotilapia burtoni

BackgroundDifferential gene expression can be translated into differing phenotypic traits. Especially during embryogenesis, specific gene expression networks regulate the development of different body structures. Cichlid fishes, with their impressive phenotypic diversity and propensity to radiate, are an emerging model system in the genomics era. Here we set out to investigate gene expression throughout development in the well-studied cichlid fish Astatotilapia burtoni, native to Lake Tanganyika and its affluent rivers.ResultsCombining RNA-sequencing from different developmental time points as well as integrating adult gene expression data, we constructed a new genome annotation for A. burtoni comprising 103,253 transcripts (stemming from 52,584 genomic loci) as well as a new reference transcriptome set. We compared our transcriptome to the available reference genome, redefining transcripts and adding new annotations. We show that about half of these transcripts have coding potential. We also characterize transcripts that are not present in the genome assembly. Next, using our newly constructed comprehensive reference transcriptome, we characterized differential gene expression through time and showed that gene expression is shifted between different body parts. We constructed a gene expression network that identified connected genes responsible for particular phenotypes and made use of it to focus on genes under potential positive selection in A. burtoni, which were implicated in fin development and vision.ConclusionsWe provide new genomic resources for the cichlid fish Astatotilapia burtoni, which will contribute to its further establishment as a model system. Tracing gene expression through time, we identified gene networks underlying particular functions, which will help to understand the genetic basis of phenotypic diversity in cichlids.

Identification of Caspase-6 and Caspase-7 from miiuy croaker and evolution analysis in fish

Apoptosis is a basic biological phenomenon of cells, which is an important component in the evolution of organisms, the stabilization of the internal environment and the development of multiple systems. In addition, the caspase protein family plays an important role in these pathways of apoptosis. Among them, apoptotic executors can directly act on specific substrates to complete the apoptotic response. In this study, we identified the Caspase-6 and Caspase-7 genes of miiuy croaker, and then analyzed the evolution of the whole Caspase family, furthermore described the evolutionary selection sites of the caspase-6 and caspase-7 genes in fish. The results showed that Caspase-6 gene appeared earlier than Caspase-7 in species evolution and gene duplication in teleost fish. Moreover, we also found that caspase-6 gene had no potential positive selection sites in the evolution of fish. Unlike the caspase-6 gene, the caspase-7 gene did not appear to be missed or replicated during the evolution of the species, while, it to be found two potential positive selection sites.