Genomes of Subaerial Zygnematophyceae Provide Insights into Land Plant Evolution.

The transition to a terrestrial environment, termed terrestrialization, is generally regarded as a pivotal event in the evolution and diversification of the land plant flora that changed the surface of our planet. Through phylogenomic studies, a group of streptophyte algae, the Zygnematophyceae, have recently been recognized as the likely sister group to land plants (embryophytes). Here, we report genome sequences and analyses of two early diverging Zygnematophyceae (Spirogloea muscicola gen. nov. and Mesotaenium endlicherianum) that share the same subaerial/terrestrial habitat with the earliest-diverging embryophytes, the bryophytes. We provide evidence that genes (i.e., GRAS and PYR/PYL/RCAR) that increase resistance to biotic and abiotic stresses in land plants, in particular desiccation, originated or expanded in the common ancestor of Zygnematophyceae and embryophytes, and were gained by horizontal gene transfer (HGT) from soil bacteria. These two Zygnematophyceae genomes represent a cornerstone for future studies to understand the underlying molecular mechanism and process of plant terrestrialization.

Corallococcus senghenyddensis sp. nov., a myxobacterium with potent antimicrobial activity.

AIM: Corallococcus species are diverse in the natural environment with 10 new Corallococcus species having been characterized in just the last 5 years. As well as being an abundant myxobacterial genus, they produce several secondary metabolites, including Corallopyronin, Corramycin, Coralmycin, and Corallorazine. We isolated a novel strain Corallococcus spp RDP092CA from soil in South Wales, UK, using Candida albicans as prey bait and characterized its predatory activities against pathogenic bacteria and yeast. METHODS AND RESULTS: The size of the RDP092CA genome was 8.5 Mb with a G + C content of 71.4%. Phylogenetically, RDP092CA is closely related to Corallococcus interemptor, C. coralloides, and C. exiguus. However, genome average nucleotide identity and digital DNA-DNA hybridization values are lower than 95% and 70% when compared to those type strains, implying that it belongs to a novel species. The RDP092CA genome harbours seven types of biosynthetic gene clusters (BGCs) and 152 predicted antimicrobial peptides. In predation assays, RDP092CA showed good predatory activity against Escherichia coli, Pseudomonas aeruginosa, Citrobacter freundii, and Staphylococcus aureus but not against Enterococcus faecalis. It also showed good antibiofilm activity against all five bacteria in biofilm assays. Antifungal activity against eight Candida spp. was variable, with particularly good activity against Meyerozyma guillermondii DSM 6381. Antimicrobial peptide RDP092CA_120 exhibited potent antibiofilm activity with >50% inhibition and >60% dispersion of biofilms at concentrations down to 1 µg/ml. CONCLUSIONS: We propose that strain RDP092CA represents a novel species with promising antimicrobial activities, Corallococcus senghenyddensis sp. nov. (=NBRC 116490T =CCOS 2109T), based on morphological, biochemical, and genomic features.

Genomic Heterogeneity of Cryptosporidium parvum Isolates From Children in Bangladesh: Implications for Parasite Biology and Human Infection.

Cryptosporidium species are a major cause of diarrhea and associated with growth failure. There is currently only limited knowledge of the parasite's genomic variability. We report a genomic analysis of Cryptosporidium parvum isolated from Bangladeshi infants and reanalysis of sequences from the United Kingdom. Human isolates from both locations shared 154 variants not present in the cattle-derived reference genome, suggesting host-specific adaptation of the parasite. Remarkably 34.6% of single-nucleotide polymorphisms unique to human isolates were nonsynonymous and 8.2% of these were in secreted proteins. Linkage disequilibrium decay indicated frequent recombination. The genetic diversity of C. parvum has potential implications for vaccine and therapeutic design. Clinical Trials Registration. NCT02764918.

E. urophylla × E. grandis high-quality genome and comparative genomics provide insights on evolution and diversification of eucalyptus.

BACKGROUND: Eucalyptus urophylla × Eucalyptus grandis, an economically important forest tree, provides important raw material for energy and reduces damage to native forests. However, the absence of a high-quality E. urophylla × E. grandis reference genome has significantly hindered its evolution and genetic analysis. RESULTS: We successfully presented a high-quality reference genome of E. urophylla × E. grandis (545.75 Mb; scaffold N50, 51.62 Mb) using a combination of the Illumina, PacBio HiFi, and Hi-C sequencing platforms. A total of 34,502 genes and 58.56% of the repetitive sequences in this genome were annotated. Using genome evolution analyses, we identified a recent whole-genome duplication (WGD) event in E. urophylla × E. grandis. We further found that gene families associated with starch and sucrose metabolism, flavonoid biosynthesis, and plant-pathogen interaction were significantly expanded in E. urophylla × E. grandis. Moreover, comparative genomic and evolutionary analyses showed large structural variations among the different chromosomes of the 34 Eucalyptus accessions, which were divided into six clades. CONCLUSIONS: Overall, our findings provide a valuable resource for expanding our understanding of the E. urophylla × E. grandis genome evolution, genetic improvement, and its comparative biology.

Exploring objective feature sets in constructing the evolution relationship of animal genome sequences.

BACKGROUND: Exploring evolution regularities of genome sequences and constructing more objective species evolution relationships at the genomic level are high-profile topics. Based on the evolution mechanism of genome sequences proposed in our previous research, we found that only the 8-mers containing CG or TA dinucleotides correlate directly with the evolution of genome sequences, and the relative frequency rather than the actual frequency of these 8-mers is more suitable to characterize the evolution of genome sequences. RESULT: Therefore, two types of feature sets were obtained, they are the relative frequency sets of CG1 + CG2 8-mers and TA1 + TA2 8-mers. The evolution relationships of mammals and reptiles were constructed by the relative frequency set of CG1 + CG2 8-mers, and two types of evolution relationships of insects were constructed by the relative frequency sets of CG1 + CG2 8-mers and TA1 + TA2 8-mers respectively. Through comparison and analysis, we found that evolution relationships are consistent with the known conclusions. According to the evolution mechanism, we considered that the evolution relationship constructed by CG1 + CG2 8-mers reflects the evolution state of genome sequences in current time, and the evolution relationship constructed by TA1 + TA2 8-mers reflects the evolution state in the early stage. CONCLUSION: Our study provides objective feature sets in constructing evolution relationships at the genomic level.

Whole-Genome Resequencing of Worldwide Wild and Domestic Sheep Elucidates Genetic Diversity, Introgression, and Agronomically Important Loci.

Domestic sheep and their wild relatives harbor substantial genetic variants that can form the backbone of molecular breeding, but their genome landscapes remain understudied. Here, we present a comprehensive genome resource for wild ovine species, landraces and improved breeds of domestic sheep, comprising high-coverage (∼16.10×) whole genomes of 810 samples from 7 wild species and 158 diverse domestic populations. We detected, in total, ∼121.2 million single nucleotide polymorphisms, ∼61 million of which are novel. Some display significant (P < 0.001) differences in frequency between wild and domestic species, or are private to continent-wide or individual sheep populations. Retained or introgressed wild gene variants in domestic populations have contributed to local adaptation, such as the variation in the HBB associated with plateau adaptation. We identified novel and previously reported targets of selection on morphological and agronomic traits such as stature, horn, tail configuration, and wool fineness. We explored the genetic basis of wool fineness and unveiled a novel mutation (chr25: T7,068,586C) in the 3'-UTR of IRF2BP2 as plausible causal variant for fleece fiber diameter. We reconstructed prehistorical migrations from the Near Eastern domestication center to South-and-Southeast Asia and found two main waves of migrations across the Eurasian Steppe and the Iranian Plateau in the Early and Late Bronze Ages. Our findings refine our understanding of genome variation as shaped by continental migrations, introgression, adaptation, and selection of sheep.

A Genetic History of the Near East from an aDNA Time Course Sampling Eight Points in the Past 4,000 Years.

The Iron and Classical Ages in the Near East were marked by population expansions carrying cultural transformations that shaped human history, but the genetic impact of these events on the people who lived through them is little-known. Here, we sequenced the whole genomes of 19 individuals who each lived during one of four time periods between 800 BCE and 200 CE in Beirut on the Eastern Mediterranean coast at the center of the ancient world's great civilizations. We combined these data with published data to traverse eight archaeological periods and observed any genetic changes as they arose. During the Iron Age (∼1000 BCE), people with Anatolian and South-East European ancestry admixed with people in the Near East. The region was then conquered by the Persians (539 BCE), who facilitated movement exemplified in Beirut by an ancient family with Egyptian-Lebanese admixed members. But the genetic impact at a population level does not appear until the time of Alexander the Great (beginning 330 BCE), when a fusion of Asian and Near Easterner ancestry can be seen, paralleling the cultural fusion that appears in the archaeological records from this period. The Romans then conquered the region (31 BCE) but had little genetic impact over their 600 years of rule. Finally, during the Ottoman rule (beginning 1516 CE), Caucasus-related ancestry penetrated the Near East. Thus, in the past 4,000 years, three limited admixture events detectably impacted the population, complementing the historical records of this culturally complex region dominated by the elite with genetic insights from the general population.

Genetic determinants of Biofilm formation of Helicobacter pylori using whole-genome sequencing.

BACKGROUND: Infection with Helicobacter pylori as the cause of gastric cancer is a global public health concern. In addition to protecting germs from antibiotics, biofilms reduce the efficacy of H. pylori eradication therapy. The nucleotide polymorphisms (SNPs) related with the biofilm forming phenotype of Helicobacter pylori were studied. RESULTS: Fifty-six H. pylori isolate from Bangladeshi patients were included in this cross-sectional study. Crystal violet assay was used to quantify biofilm amount, and the strains were classified into high- and low-biofilm formers As a result, strains were classified as 19.6% high- and 81.4% low-biofilm formers. These phenotypes were not related to specific clades in the phylogenetic analysis. The accessories genes associated with biofilm from whole-genome sequences were extracted and analysed, and SNPs among the previously reported biofilm-related genes were analysed. Biofilm formation was significantly associated with SNPs of alpA, alpB, cagE, cgt, csd4, csd5, futB, gluP, homD, and murF (P < 0.05). Among the SNPs reported in alpB, strains encoding the N156K, G160S, and A223V mutations were high-biofilm formers. CONCLUSIONS: This study revealed the potential role of SNPs in biofilm formation and proposed a method to detect mutation in biofilm from whole-genome sequences.

Does selection occur at the intermediate zone of two insufficiently isolated populations? A whole-genome analysis along an altitudinal gradient.

Adaptive divergence occurs even between insufficiently isolated populations when there is a great difference in environments between their habitats. Individuals present in an intermediate zone of the two divergent populations are expected to have an admixed genetic structure due to gene flow. A selective pressure that acts on the genetically admixed individuals may limit the gene flow and maintain the adaptive divergence. Here, we addressed a question whether selection occurs in the genetically admixed individuals between two divergent populations. Arabidopsis halleri is a perennial montane plant, which has clear phenotypic dimorphisms between highland and lowland habitats in Mt. Ibuki, central Japan. We obtained the whole-genome sequences of Arabidopsis halleri plants along an altitudinal gradient of 359-1,317 m with a high spatial resolution (mean altitudinal interval of 20 m). We found a zone where the highland and lowland genes were mixing (intermediate subpopulation). In the intermediate subpopulation, we identified 5 and 13 genome regions, which included 3 and 8 genes, that had a high frequency of alleles that are accumulated in highland and lowland subpopulations, respectively. In addition, we also found that the frequency of highland alleles of these selected genome regions was smaller in the lowland subpopulation compared with that of the non-selected regions. These results suggest that the selection in the intermediate subpopulation might limit the gene flow and contribute to the adaptive divergence between altitudes. We also identified 7 genome regions that had low heterozygote frequencies in the intermediate subpopulation. We conclude that different types of selection in addition to gene flow occur at the intermediate altitude and shape the genetic structure across altitudes.

Genetic evaluations for endangered dual-purpose German Black Pied cattle using 50K SNPs, a breed-specific 200K chip, and whole-genome sequencing.

Genetic evaluations of local cattle breeds are hampered due to small reference groups or biased due to the utilization of SNP effects estimated in other large populations. Against this background, there is a lack of studies addressing the possible advantage of whole-genome sequences (WGS) or consideration of specific variants from WGS data in genomic predictions for local breeds with small population size. Consequently, the aim of this study was to compare genetic parameters and accuracies of genomic estimated breeding values (GEBV) for 305-d production traits, fat-to protein ratio (FPR), and somatic cell score (SCS) at the first test date after calving and confirmation traits of the endangered German Black Pied cattle (DSN) breed using 4 different marker panels: (1) the commercial 50K Illumina BovineSNP50 BeadChip, (2) a customized 200K chip designed for DSN (DSN200K) which considers the most important variants for DSN from WGS, (3) randomly generated 200K chips based on WGS data, and (4) a WGS panel. The same number of animals was considered for all marker panel analyses (i.e., 1,811 genotyped or sequenced cows for conformation traits, 2,383 cows for lactation production traits, and 2,420 cows for FPR and SCS). Mixed models for the estimation of genetic parameters directly included the respective genomic relationship matrix from the different marker panels plus the trait-specific fixed effects. For the calculation of GEBV accuracies, we applied repeated random subsampling validation. In the process of separate cross-validations per trait, we created a validation set including 20% of cows with masked phenotypes, and a training set comprising 80% of the cows. The cows were selected randomly in a procedure with 10 replicates considering replacements in the different scenarios. The accuracy was defined as the correlation between the direct GEBV and the phenotypes with subtracted corresponding fixed effects for the cows in the validation set. For FPR and SCS, as well as for lactation production traits, heritabilities were largest based on WGS data, but the increase compared with the 50K or DSN200K applications was quite small in the range from 0.01 to 0.03. Also, for most of the conformation traits, heritabilities were largest based on WGS and DSN200K data, but the increase was in the range of the corresponding standard error. Accordingly, GEBV accuracies for most of the studied traits were highest based on WGS data or when utilizing the DSN200K chip, but the accuracy differences across the marker panels were quite small and nonsignificant. In conclusion, WGS data and the DSN200K chip only contributed to minor improvements in genomic predictions, still justifying the use of the commercial 50K chip. Nevertheless, WGS and the 200KDSN chip harbor breed-specific variants, which are valuable for studying causal genetic mechanisms in the endangered DSN population.

Genome Resources of Verticillium dahliae VdGL16: The Causal Agent of Vascular Wilt on the Invasive Species Ailanthus altissima.

Verticillium species are known as plant pathogens responsible for wilt diseases in a large variety of dicotyledon plants and crops in many parts of the world. Here we present the draft genome sequence of Verticillium dahliae Kleb. (strain VdGL16) isolated in Italy from the invasive alien species Ailanthus altissima (Mill.; commonly known as tree-of-heaven) showing Verticillium wilt symptoms. The comparison between the newly sequenced genome with those publicly available revealed candidate genes putatively involved in pathogenicity. The genome represents a new useful source for future research on Verticillium genetics and biology as well as research on novel approaches in the control of A. altissima.

getSequenceInfo: a suite of tools allowing to get genome sequence information from public repositories.

BACKGROUND: Biological sequences are increasing rapidly and exponentially worldwide. Nucleotide sequence databases play an important role in providing meaningful genomic information on a variety of biological organisms. RESULTS: The getSequenceInfo software tool allows to access sequence information from various public repositories (GenBank, RefSeq, and the European Nucleotide Archive), and is compatible with different operating systems (Linux, MacOS, and Microsoft Windows) in a programmatic way (command line) or as a graphical user interface. getSequenceInfo or gSeqI v1.0 should help users to get some information on queried sequences that could be useful for specific studies (e.g. the country of origin/isolation or the release date of queried sequences). Queries can be made to retrieve sequence data based on a given kingdom and species, or from a given date. This program allows the separation between chromosomes and plasmids (or other genetic elements/components) by arranging each component in a given folder. Some basic statistics are also performed by the program (such as the calculation of GC content for queried assemblies). An empirically designed nucleotide ratio is calculated using nucleotide information in order to tentatively provide a "NucleScore" for studied genome assemblies. Besides the main gSeqI tool, other additional tools have been developed to perform various tasks related to sequence analysis. CONCLUSION: The aim of this study is to democratize the use of public repositories in programmatic ways, and to facilitate sequence data analysis in a pedagogical perspective. Output results are available in FASTA, FASTQ, Excel/TSV or HTML formats. The program is freely available at: https://github.com/karubiotools/getSequenceInfo . getSequenceInfo and supplementary tools are partly available through the recently released Galaxy KaruBioNet platform ( http://calamar.univ-ag.fr/c3i/galaxy_karubionet.html ).

Historical Introgression from Wild Relatives Enhanced Climatic Adaptation and Resistance to Pneumonia in Sheep.

How animals, particularly livestock, adapt to various climates and environments over short evolutionary time is of fundamental biological interest. Further, understanding the genetic mechanisms of adaptation in indigenous livestock populations is important for designing appropriate breeding programs to cope with the impacts of changing climate. Here, we conducted a comprehensive genomic analysis of diversity, interspecies introgression, and climate-mediated selective signatures in a global sample of sheep and their wild relatives. By examining 600K and 50K genome-wide single nucleotide polymorphism data from 3,447 samples representing 111 domestic sheep populations and 403 samples from all their seven wild relatives (argali, Asiatic mouflon, European mouflon, urial, snow sheep, bighorn, and thinhorn sheep), coupled with 88 whole-genome sequences, we detected clear signals of common introgression from wild relatives into sympatric domestic populations, thereby increasing their genomic diversities. The introgressions provided beneficial genetic variants in native populations, which were significantly associated with local climatic adaptation. We observed common introgression signals of alleles in olfactory-related genes (e.g., ADCY3 and TRPV1) and the PADI gene family including in particular PADI2, which is associated with antibacterial innate immunity. Further analyses of whole-genome sequences showed that the introgressed alleles in a specific region of PADI2 (chr2: 248,302,667-248,306,614) correlate with resistance to pneumonia. We conclude that wild introgression enhanced climatic adaptation and resistance to pneumonia in sheep. This has enabled them to adapt to varying climatic and environmental conditions after domestication.

A Transient Pulse of Genetic Admixture from the Crusaders in the Near East Identified from Ancient Genome Sequences.

During the medieval period, hundreds of thousands of Europeans migrated to the Near East to take part in the Crusades, and many of them settled in the newly established Christian states along the Eastern Mediterranean coast. Here, we present a genetic snapshot of these events and their aftermath by sequencing the whole genomes of 13 individuals who lived in what is today known as Lebanon between the 3rd and 13th centuries CE. These include nine individuals from the "Crusaders' pit" in Sidon, a mass burial in South Lebanon identified from the archaeology as the grave of Crusaders killed during a battle in the 13th century CE. We show that all of the Crusaders' pit individuals were males; some were Western Europeans from diverse origins, some were locals (genetically indistinguishable from present-day Lebanese), and two individuals were a mixture of European and Near Eastern ancestries, providing direct evidence that the Crusaders admixed with the local population. However, these mixtures appear to have had limited genetic consequences since signals of admixture with Europeans are not significant in any Lebanese group today-in particular, Lebanese Christians are today genetically similar to local people who lived during the Roman period which preceded the Crusades by more than four centuries.

OrchidBase 5.0: updates of the orchid genome knowledgebase.

Containing the largest number of species, the orchid family provides not only materials for studying plant evolution and environmental adaptation, but economically and culturally important ornamental plants for human society. Previously, we collected genome and transcriptome information of Dendrobium catenatum, Phalaenopsis equestris, and Apostasia shenzhenica which belong to two different subfamilies of Orchidaceae, and developed user-friendly tools to explore the orchid genetic sequences in the OrchidBase 4.0. The OrchidBase 4.0 offers the opportunity for plant science community to compare orchid genomes and transcriptomes and retrieve orchid sequences for further study.In the year 2022, two whole-genome sequences of Orchidoideae species, Platanthera zijinensis and Platanthera guangdongensis, were de novo sequenced, assembled and analyzed. In addition, systemic transcriptomes from these two species were also established. Therefore, we included these datasets to develop the new version of OrchidBase 5.0. In addition, three new functions including synteny, gene order, and miRNA information were also developed for orchid genome comparisons and miRNA characterization.OrchidBase 5.0 extended the genetic information to three orchid subfamilies (including five orchid species) and provided new tools for orchid researchers to analyze orchid genomes and transcriptomes. The online resources can be accessed at https://cosbi.ee.ncku.edu.tw/orchidbase5/.

First report of Ranunculus mild mosaic virus in Ranunculus asiaticus in Korea.

Ranunculus (Ranunculus asiaticus L.) is a popular ornamental plant mainly cultivated for cut flowers and flowering potted plants. In January 2021, a leaf sample of R. asiaticus that showed virus-like symptoms including mosaic, yellowing and malformation on leaves was collected from a greenhouse in Jangheung, South Korea for disease diagnosis (Fig. S1). Disease incidence was greater than 30% in the greenhouse (~1,000 m2). Transmission electron microscopy (TEM) of symptomatic leaves identified potyvirus-like filamentous virus particles of about 800 nm. To confirm the TEM results, a symptomatic leaf sample was further analyzed by reverse-transcription polymerase chain reaction (RT-PCR) using species-specific detection primers for six potyviruses that infect R. asiaticus (Sacco et al., 2018). The sample was positive only for ranunculus mild mosaic virus (RanMMV). Additional analysis of nine symptomatic R. asiaticus plants from the infected greenhouse found that all samples were positive for RanMMV. To exclude the presence of the other viruses, next generation sequencing (NGS) was carried out. Total RNA was extracted from symptomatic leaves using the RNeasy Plant Mini Kit (Qiagen, Germany) and a transcriptome library was generated using the TruSeq Stranded Total RNA LT Sample Prep kit (Illumina, San Diego, CA) acccording to the recommended protocol. NGS was performed using an Illumina NovaSeq 6000 system (Macrogen Inc., Korea). A total of 75.58 million reads were obtained, and the reads were de novo assembled to contigs using Trinity software (Grabherr et al., 2011). BLASTn and BLASTx analysis of the contigs against the NCBI viral reference database identified the assembled large contig of 9,539 nt (5,321 mapped reads, mean read coverage of 84.2 times) as RanMMV. This sequence shared 98% nt identity (99% coverage) with the RanMMV NL isolate (acc. no. LC604020) isolated from an anemone plant (A. blanda cv. Charmer) from Netherlands. To obtain the complete genome sequence, the termini sequences were determined by 5' and 3' rapid amplification of cDNA ends (RACE) methods as reported recently (Imamura et al., 2021). The assembled full-length genome sequence of RanMMV-JH is 9,574 nt in length, excluding the poly(A) tail, and encoding a polyprotein of 3,074aa. The sequence was deposited in GenBank under the accession no. OL742438. RanMMV is transmitted by aphids in a nonpersistent manner and has very narrow host range. RanMMV, one of causative agents of ranunculus mosaic disease, has been problematic in ranunculus production area of Japan (Hayahi et al., 2018; Kamikawa et al., 2022). Recently, some perennial weeds from the Ranunculaceae family (e.g. R. japonicus, R. silerifolius and R. tachiroei) are known to may act as a virus reservoir (Kamikawa et al., 2022). As R. asiaticus is cultivated by vegetative propagation, there is need to develop certification system for producing virus-free R. asiaticus. To our knowledge, this is the first report of RanMMV infection in R. asiaticus in Korea.

Introgression contributes to distribution of structural variations in cattle.

Structural variations (SVs) are an important source of phenotypic diversity in cattle. Here, 72 whole genome sequences representing taurine and zebu cattle were used to identify SVs. Applying multiple approaches, 16,738 SVs were identified. A comparison against the Database of Genomic Variants archives revealed that 1575 SVs were novel in our data. A novel duplication covering the entire GALNT15 gene, was observed only in N'Dama. A duplication, which was previously reported only in zebu and associated with navel length, was also observed in N'Dama. Investigation of a novel deletion located upstream of CAST13 gene and identified only in Italian cattle and zebu, revealed its introgressed origin in the former. Overall, our data highlights how the SVs distribution in cattle is also shaped by forces such as demographical differences and gene flow. The cattle SVs of this study and its meta-data can be visualized on an interactive genome browser at https://tinyurl.com/svCowArs.

Assembly and Annotation of Red Spruce (Picea rubens) Chloroplast Genome, Identification of Simple Sequence Repeats, and Phylogenetic Analysis in Picea.

We have sequenced the chloroplast genome of red spruce (Picea rubens) for the first time using the single-end, short-reads (44 bp) Illumina sequences, assembled and functionally annotated it, and identified simple sequence repeats (SSRs). The contigs were assembled using SOAPdenovo2 following the retrieval of chloroplast genome sequences using the black spruce (Picea mariana) chloroplast genome as the reference. The assembled genome length was 122,115 bp (gaps included). Comparatively, the P. rubens chloroplast genome reported here may be considered a near-complete draft. Global genome alignment and phylogenetic analysis based on the whole chloroplast genome sequences of Picea rubens and 10 other Picea species revealed high sequence synteny and conservation among 11 Picea species and phylogenetic relationships consistent with their known classical interrelationships and published molecular phylogeny. The P. rubens chloroplast genome sequence showed the highest similarity with that of P. mariana and the lowest with that of P. sitchensis. We have annotated 107 genes including 69 protein-coding genes, 28 tRNAs, 4 rRNAs, few pseudogenes, identified 42 SSRs, and successfully designed primers for 26 SSRs. Mononucleotide A/T repeats were the most common followed by dinucleotide AT repeats. A similar pattern of microsatellite repeats occurrence was found in the chloroplast genomes of 11 Picea species.

Investigation of ancestral alleles in the Bovinae subfamily.

BACKGROUND: In evolutionary theory, divergence and speciation can arise from long periods of reproductive isolation, genetic mutation, selection and environmental adaptation. After divergence, alleles can either persist in their initial state (ancestral allele - AA), co-exist or be replaced by a mutated state (derived alleles -DA). In this study, we aligned whole genome sequences of individuals from the Bovinae subfamily to the cattle reference genome (ARS.UCD-1.2) for defining ancestral alleles necessary for selection signatures study. RESULTS: Accommodating independent divergent of each lineage from the initial ancestral state, AA were defined based on fixed alleles on at least two groups of yak, bison and gayal-gaur-banteng resulting in ~ 32.4 million variants. Using non-overlapping scanning windows of 10 Kb, we counted the AA observed within taurine and zebu cattle. We focused on the extreme points, regions with top 0. 1% (high count) and regions without any occurrence of AA (null count). High count regions preserved gene functions from ancestral states that are still beneficial in the current condition, while null counts regions were linked to mutated ones. For both cattle, high count regions were associated with basal lipid metabolism, essential for survival of various environmental pressures. Mutated regions were associated to productive traits in taurine, i.e. higher metabolism, cell development and behaviors and in immune response domain for zebu. CONCLUSIONS: Our findings suggest that retaining and losing AA in some regions are varied and made it species-specific with possibility of overlapping as it depends on the selective pressure they had to experience.

OrchidBase 4.0: a database for orchid genomics and molecular biology.

BACKGROUND: The Orchid family is the largest families of the monocotyledons and an economically important ornamental plant worldwide. Given the pivotal role of this plant to humans, botanical researchers and breeding communities should have access to valuable genomic and transcriptomic information of this plant. Previously, we established OrchidBase, which contains expressed sequence tags (ESTs) from different tissues and developmental stages of Phalaenopsis as well as biotic and abiotic stress-treated Phalaenopsis. The database includes floral transcriptomic sequences from 10 orchid species across all the five subfamilies of Orchidaceae. DESCRIPTION: Recently, the whole-genome sequences of Apostasia shenzhenica, Dendrobium catenatum, and Phalaenopsis equestris were de novo assembled and analyzed. These datasets were used to develop OrchidBase 4.0, including genomic and transcriptomic data for these three orchid species. OrchidBase 4.0 offers information for gene annotation, gene expression with fragments per kilobase of transcript per millions mapped reads (FPKM), KEGG pathways and BLAST search. In addition, assembled genome sequences and location of genes and miRNAs could be visualized by the genome browser. The online resources in OrchidBase 4.0 can be accessed by browsing or using BLAST. Users can also download the assembled scaffold sequences and the predicted gene and protein sequences of these three orchid species. CONCLUSIONS: OrchidBase 4.0 is the first database that contain the whole-genome sequences and annotations of multiple orchid species. OrchidBase 4.0 is available at http://orchidbase.itps.ncku.edu.tw/.