Complete Mitogenome of "Pumpo" (Bos taurus), a Top Bull from a Peruvian Genetic Nucleus, and Its Phylogenetic Analysis.

The mitochondrial genome of Pumpo (Bos taurus), a prominent breed contributing to livestock farming, was sequenced using the Illumina HiSeq 2500 platform. Assembly and annotation of the mitochondrial genome were achieved through a multifaceted approach employing bioinformatics tools such as Trim Galore, SPAdes, and Geseq, followed by meticulous manual inspection. Additionally, analyses covering tRNA secondary structure and codon usage bias were conducted for comprehensive characterization. The 16,341 base pair mitochondrial genome comprises 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes. Phylogenetic analysis places Pumpo within a clade predominantly composed of European cattle, reflecting its prevalence in Europe. This comprehensive study underscores the importance of mitochondrial genome analysis in understanding cattle evolution and highlights the potential of genetic improvement programs in livestock farming, thus contributing to enhanced livestock practices.

Greatly reduced phylogenetic structure in the cultivated potato clade (Solanum section Petota pro parte).

PREMISE OF THE STUDY: The species boundaries of wild and cultivated potatoes are controversial, with most of the taxonomic problems in the cultivated potato clade. We here provide the first in-depth phylogenetic study of the cultivated potato clade to explore possible causes of these problems. METHODS: We examined 131 diploid accessions, using 12 nuclear orthologs, producing an aligned data set of 14,072 DNA characters, 2171 of which are parsimony-informative. We analyzed the data to produce phylogenies and perform concordance analysis and goodness-of-fit tests. KEY RESULTS: There is good phylogenetic structure in clades traditionally referred to as clade 1+2 (North and Central American diploid potatoes exclusive of Solanum verrucosum), clade 3, and a newly discovered basal clade, but drastically reduced phylogenetic structure in clade 4, the cultivated potato clade. The results highlight a clade of species in South America not shown before, 'neocardenasii', sister to clade 1+2, that possesses key morphological traits typical of diploids in Mexico and Central America. Goodness-of-fit tests suggest potential hybridization between some species of the cultivated potato clade. However, we do not have enough phylogenetic signal with the data at hand to explicitly estimate such hybridization events with species networks methods. CONCLUSIONS: We document the close relationships of many of the species in the cultivated potato clade, provide insight into the cause of their taxonomic problems, and support the recent reduction of species in this clade. The discovery of the neocardenasii clade forces a reevaluation of a hypothesis that section Petota originated in Mexico and Central America.

Genotyping-by-sequencing provides the discriminating power to investigate the subspecies of Daucus carota (Apiaceae).

BACKGROUND: The majority of the subspecies of Daucus carota have not yet been discriminated clearly by various molecular or morphological methods and hence their phylogeny and classification remains unresolved. Recent studies using 94 nuclear orthologs and morphological characters, and studies employing other molecular approaches were unable to distinguish clearly many of the subspecies. Fertile intercrosses among traditionally recognized subspecies are well documented. We here explore the utility of single nucleotide polymorphisms (SNPs) generated by genotyping-by-sequencing (GBS) to serve as an effective molecular method to discriminate the subspecies of the D. carota complex. RESULTS: We used GBS to obtain SNPs covering all nine Daucus carota chromosomes from 162 accessions of Daucus and two related genera. To study Daucus phylogeny, we scored a total of 10,814 or 38,920 SNPs with a maximum of 10 or 30 % missing data, respectively. To investigate the subspecies of D. carota, we employed two data sets including 150 accessions: (i) rate of missing data 10 % with a total of 18,565 SNPs, and (ii) rate of missing data 30 %, totaling 43,713 SNPs. Consistent with prior results, the topology of both data sets separated species with 2n = 18 chromosome from all other species. Our results place all cultivated carrots (D. carota subsp. sativus) in a single clade. The wild members of D. carota from central Asia were on a clade with eastern members of subsp. sativus. The other subspecies of D. carota were in four clades associated with geographic groups: (1) the Balkan Peninsula and the Middle East, (2) North America and Europe, (3) North Africa exclusive of Morocco, and (4) the Iberian Peninsula and Morocco. Daucus carota subsp. maximus was discriminated, but neither it, nor subsp. gummifer (defined in a broad sense) are monophyletic. CONCLUSIONS: Our study suggests that (1) the morphotypes identified as D. carota subspecies gummifer (as currently broadly circumscribed), all confined to areas near the Atlantic Ocean and the western Mediterranean Sea, have separate origins from sympatric members of other subspecies of D. carota, (2) D. carota subsp. maximus, on two clades with some accessions of subsp. carota, can be distinguished from each other but only with poor morphological support, (3) D. carota subsp. capillifolius, well distinguished morphologically, is an apospecies relative to North African populations of D. carota subsp. carota, (4) the eastern cultivated carrots have origins closer to wild carrots from central Asia than to western cultivated carrots, and (5) large SNP data sets are suitable for species-level phylogenetic studies in Daucus.

Morphometrics of Daucus (Apiaceae): a counterpart to a phylogenomic study.

PREMISE OF STUDY: Molecular phylogenetics of genome-scale data sets (phylogenomics) often produces phylogenetic trees with unprecedented resolution. A companion phylogenomics analysis of Daucus using 94 conserved nuclear orthologs supported many of the traditional species but showed unexpected results that require morphological analyses to help interpret them in a practical taxonomic context. METHODS: We evaluated character state distributions, stepwise discriminant analyses, canonical variate analyses, and hierarchical cluster analyses from 40 morphological characters from 81 accessions of 14 taxa of Daucus and eight species in related genera in an experimental plot. KEY RESULTS: Most characters showed tremendous variation with character state overlap across many taxa. Multivariate analyses separated the outgroup taxa easily from the Daucus ingroup. Concordant with molecular analyses, most species form phenetic groups, except the same taxa that are problematical in the molecular results: (1) the subspecies of D. carota, (2) D. sahariensis and D. syrticus, and (3) D. broteri and D. guttatus. CONCLUSIONS: Phenetic analyses, in combination with molecular data, support many Daucus species, but mostly by overlapping ranges of size and meristic variation. The subspecies of D. carota are poorly separated morphologically, are paraphyletic, and all could be recognized at the subspecies rank under D. carota. Daucus sahariensis and D. syrticus are so similar morphologically that they could be placed in synonymy. Combined molecular and morphological data support three species in accessions previously identified as D. broteri and D. guttatus. Molecular and morphological results support the new combination Daucus carota subsp. capillifolius.

Phylogenomics of the carrot genus (Daucus, Apiaceae).

UNLABELLED: • PREMISE OF THE STUDY: We explored the utility of multiple nuclear orthologs for the taxonomic resolution of wild and cultivated carrot, Daucus species.• METHODS: We studied the phylogeny of 92 accessions of 13 species and two subspecies of Daucus and 15 accessions of related genera (107 accessions total) with DNA sequences of 94 nuclear orthologs. Reiterative analyses examined data of both alleles using ambiguity codes or a single allele with the highest coverage, trimmed vs. untrimmed homopolymers; pure exonic vs. pure intronic data; the use of all 94 markers vs. a reduced subset of markers; and analysis of a concatenated data set vs. a coalescent (species tree) approach.• KEY RESULTS: Our maximum parsimony and maximum likelihood trees were highly resolved, with 100% bootstrap support for most of the external and many of the internal clades. They resolved multiple accessions of many different species as monophyletic with strong support, but failed to support other species. The single allele analysis gave slightly better topological resolution; trimming homopolymers failed to increase taxonomic resolution; the exonic data had a smaller proportion of parsimony-informative characters. Similar results demonstrating the same dominant topology can be obtained with many fewer markers. A Bayesian concordance analysis provided an overall similar phylogeny, but the coalescent analysis provided drastic changes in topology to all the above.• CONCLUSIONS: Our research highlights some difficult species groups in Daucus and misidentifications in germplasm collections. It highlights a useful subset of markers and approaches for future studies of dominant topologies in Daucus.

Changes in Bulk and Rhizosphere Soil Microbial Diversity Communities of Native Quinoa Due to the Monocropping in the Peruvian Central Andes.

Quinoa (Chenopodium quinoa) is a highly nutritious crop that is resistant to adverse conditions. Due to the considerable increase in its commercial production in Andean soils, the plant is suffering the negative effects of monocropping, which reduces its yield. We used for the first time a high-throughput Illumina MiSeq sequencing approach to explore the composition, diversity, and functions of fungal and bacterial communities of the bulk and rhizosphere in soils of native C. quinoa affected by monocropping in the central Andes of Peru. The results showed that the bacterial and fungal community structure among the treatments was significantly changed by the monocropping and the types of soil (rhizosphere and bulk). Also, in soils subjected to monocropping, there was an increase in Actinobacteria and a decrease in Proteobacteria, and the reduction in the presence of Ascomycota and the increase in Basidiomycota. By alpha-diversity indices, lower values of bacteria and fungi were observed in the monoculture option compared to the soil not affected by monocropping, and sometimes significant differences were found between both. We detected differentially abundant phytopathogenic fungi and bacteria with growth-stimulating effects on plants. Also, we denoted a decrease in the abundance of the functional predictions in bacteria in the monocropped soils. This research will serve as a starting point to explore the importance and effects of microorganisms in degraded soils and their impact on the growth and quality of quinoa crops.

Genetic diversity and population structure of a Peruvian cattle herd using SNP data.

New-generation sequencing technologies, among them SNP chips for massive genotyping, are useful for the effective management of genetic resources. To date, molecular studies in Peruvian cattle are still scarce. For the first time, the genetic diversity and population structure of a reproductive nucleus cattle herd of four commercial breeds from a Peruvian institution were determined. This nucleus comprises Brahman (N = 9), Braunvieh (N = 9), Gyr (N = 5), and Simmental (N = 15) breeds. Additionally, samples from a locally adapted creole cattle, the Arequipa Fighting Bull (AFB, N = 9), were incorporated. Female individuals were genotyped with the GGPBovine100K and males with the BovineHD. Quality control, and the proportion of polymorphic SNPs, minor allele frequency, expected heterozygosity, observed heterozygosity, and inbreeding coefficient were estimated for the five breeds. Admixture, principal component analysis (PCA), and discriminant analysis of principal components (DAPC) were performed. Also, a dendrogram was constructed using the Neighbor-Joining clustering algorithm. The genetic diversity indices in all breeds showed a high proportion of polymorphic SNPs, varying from 51.42% in Gyr to 97.58% in AFB. Also, AFB showed the highest expected heterozygosity estimate (0.41 ± 0.01), while Brahman the lowest (0.33 ± 0.01). Besides, Braunvieh possessed the highest observed heterozygosity (0.43 ± 0.01), while Brahman the lowest (0.37 ± 0.02), indicating that Brahman was less diverse. According to the molecular variance analysis, 75.71% of the variance occurs within individuals, whereas 24.29% occurs among populations. The pairwise genetic differentiation estimates (FST) between breeds showed values that ranged from 0.08 (Braunvieh vs. AFB) to 0.37 (Brahman vs. Braunvieh). Similarly, pairwise Reynold's distance ranged from 0.09 (Braunvieh vs. AFB) to 0.46 (Brahman vs. Braunvieh). The dendrogram, similar to the PCA, identified two groups, showing a clear separation between Bos indicus (Brahman and Gyr) and B. taurus breeds (Braunvieh, Simmental, and AFB). Simmental and Braunvieh grouped closely with the AFB cattle. Similar results were obtained for the population structure analysis with K = 2. The results from this study would contribute to the appropriate management, avoiding loss of genetic variability in these breeds and for future improvements in this nucleus. Additional work is needed to speed up the breeding process in the Peruvian cattle system.

Microsatellite-Based Genetic Diversity and Population Structure of Huacaya Alpacas (Vicugna pacos) in Southern Peru.

The alpaca population mostly consists of the Huacaya phenotype and is widely distributed in Southern Peru. This study aimed to estimate the genetic diversity and population structure of two Huacaya alpaca populations (Ajoyani and Quimsachata) using fourteen and twelve microsatellite markers for each population, respectively. A total of 168 alpaca biological samples were outsourced to Peruvian laboratories for DNA extraction and genotyping. For genetic diversity, observed heterozygosity (Ho), expected heterozygosity (He), polymorphism information content (PIC), and fixation indices values were estimated. An admixture analysis was performed for the population structure analysis. Different programs were used for these estimations. In total, 133 (Ajoyani) and 129 (Quimsachata) alleles were found, with a range of 4 to 17 by locus. The mean HO, HE, and PIC per marker for Ajoyani were 0.764 ± 0.112, 0.771 ± 0.1, and 0.736; for Quimsachata, they were 0.783 ± 0.087, 0.773 ± 0.095, and 0.738, respectively. The population structure showed no structure with K = 2. This study provides useful indicators for the creation of appropriate alpaca conservation programs.

Draft Genome Sequence Resource of Erwinia sp. Strain INIA01, a Phytopathogen Isolated from a Diseased Stalk of Peruvian Maize.

Here, we report the complete genome sequence of Erwinia sp. strain INIA01, a bacterium isolated from lesions of Zea mays from northern Peru. This genome possesses two circular replicons, a 4.2-Mb chromosome, and a 438-kb plasmid.

Unlocking the Complete Chloroplast Genome of a Native Tree Species from the Amazon Basin, Capirona (Calycophyllum Spruceanum, Rubiaceae), and Its Comparative Analysis with Other Ixoroideae Species.

Capirona (Calycophyllum spruceanum Benth.) belongs to subfamily Ixoroideae, one of the major lineages in the Rubiaceae family, and is an important timber tree. It originated in the Amazon Basin and has widespread distribution in Bolivia, Peru, Colombia, and Brazil. In this study, we obtained the first complete chloroplast (cp) genome of capirona from the department of Madre de Dios located in the Peruvian Amazon. High-quality genomic DNA was used to construct libraries. Pair-end clean reads were obtained by PE 150 library and the Illumina HiSeq 2500 platform. The complete cp genome of C. spruceanum has a 154,480 bp in length with typical quadripartite structure, containing a large single copy (LSC) region (84,813 bp) and a small single-copy (SSC) region (18,101 bp), separated by two inverted repeat (IR) regions (25,783 bp). The annotation of C. spruceanum cp genome predicted 87 protein-coding genes (CDS), 8 ribosomal RNA (rRNA) genes, 37 transfer RNA (tRNA) genes, and one pseudogene. A total of 41 simple sequence repeats (SSR) of this cp genome were divided into mononucleotides (29), dinucleotides (5), trinucleotides (3), and tetranucleotides (4). Most of these repeats were distributed in the noncoding regions. Whole chloroplast genome comparison with the other six Ixoroideae species revealed that the small single copy and large single copy regions showed more divergence than inverted regions. Finally, phylogenetic analyses resolved that C. spruceanum is a sister species to Emmenopterys henryi and confirms its position within the subfamily Ixoroideae. This study reports for the first time the genome organization, gene content, and structural features of the chloroplast genome of C. spruceanum, providing valuable information for genetic and evolutionary studies in the genus Calycophyllum and beyond.

Characterization of the complete chloroplast genome of a Peruvian landrace of Capsicum chinense Jacq. (Solanaceae), arnaucho chili pepper.

In this study, we sequenced the first complete chloroplast (cp) genome of a Peruvian chili pepper landrace, "arnacucho" (Capsicum chinense). This cp genome has a 156,931 bp in length with typical quadripartite structure, containing a large single copy (LSC) region (87,325 bp) and a 17,912 bp small single-copy (SSC) region, separated by two inverted repeat (IR) regions (25,847 bp); and the percentage of GC content was 37.71%. Arnaucho chili pepper chloroplast genome possesses 133 genes that consists of 86 protein-coding genes, 37 tRNA, eight rRNA, and two pseudogenes. Phylogenetic analysis revealed that this Peruvian chili pepper landrace is closely related to the undomesticated species C. galapagoense; all belong to the Capsiceae tribe.

Influence of arboreal components on the physical-chemical characteristics of the soil under four silvopastoral systems in northeastern Peru.

Silvopastoral systems (SPS) are presented as an alternative for the protection and recovery of soils; however, the relationship between the tree component and the physical-chemical characteristics of the soil is unknown. Thus, the objective of this study was to evaluate the physical-chemical characteristics of the soil under four silvopastoral systems (SPS), alder (Alnus acuminata), pine (Pinus patula), cypress (Cupressus macrocarpa), and pona (Ceroxylon quindiuense), and a treeless system (TS) in the Amazonas region. A completely randomized design (CRD) with five treatments and three replicates was used. The experimental units were sampled at two depths, 0-15 and 15-30 cm. The parameters evaluated were pH, electrical conductivity (dS/m), organic matter (%), phosphorus (ppm), potassium (ppm), cation exchange capacity (meq/100 g), porosity (%), mechanical resistance (kg/cm2), bulk density (gr/cm3), moisture (%) and total carbon (t/ha). The results were analyzed by analysis of variance (α = 0.05 %) and Tukey's test of means (p ≤ 0.05). The systems presented strong acidic pH values (4.11-5.61), which resulted in high organic matter contents in all systems (6.74-9.99 %). The highest phosphorus content was in the SPS with alder (12.64 ppm), and the highest potassium content was in the SPS with cypress (382.33 ppm). Porosity in all systems was higher than 60 %. The highest bulk density was between 15 and 30 cm, and the highest percentage of moisture was in the surface layer (0-15 cm). The mechanical strength was higher in the SPS with cypress (2.62 kg/cm2). For all the systems evaluated, the highest carbon stock was found in the first 15 cm. The SPS with pine had the best soil characteristics and carbon sequestration (149.05 t/ha).

The complete chloroplast genome of the national tree of Peru, quina (Cinchona officinalis L., Rubiaceae).

Here, we report the first complete chloroplast (cp) genome of Cinchona officinalis. This cp genome has a 156,984 bp in length with typical quadripartite structure, containing a large single copy (LSC) region (83,929 bp) and an 18,051 bp small single-copy (SSC) region, separated by two inverted repeat (IR) regions (27,502 bp). The total GC content was 37.75%. Quina tree chloroplast genome possesses 135 genes that consisted of 89 protein-coding genes, 37 tRNA, eight rRNA, and one pseudogene. Phylogenetic analysis showed that C. officinalis is sister to C. pubescens and sister to them is Isertia laevis; all belong to the Cinchonoideae sub-family.

Congruence between morphological and molecular markers inferred from the analysis of the intra-morphotype genetic diversity and the spatial structure of Oxalis tuberosa Mol.

Oxalis tuberosa is an important crop cultivated in the highest Andean zones. A germplasm collection is maintained ex situ by CIP, which has developed a morphological markers system to classify the accessions into morphotypes, i.e. groups of morphologically identical accessions. However, their genetic uniformity is currently unknown. The ISSR technique was used in two experiments to determine the relationships between both morphological and molecular markers systems. The intra-morphotype genetic diversity, the spatial structures of the diversity and the congruence between both markers systems were determined. In the first experience, 44 accessions representing five morphotypes, clearly distinct from each other, were analyzed. At the molecular level, the accessions exactly clustered according to their morphotypes. However, a genetic variability was observed inside each morphotype. In the second experiment, 34 accessions gradually differing from each other on morphological base were analyzed. The morphological clustering showed no geographical structure. On the opposite, the molecular analysis showed that the genetic structure was slightly related to the collection site. The correlation between both markers systems was weak but significant. The lack of perfect congruence between morphological and molecular data suggests that the morphological system may be useful for the morphotypes management but is not appropriate to study the genetic structure of the oca. The spatial structure of the genetic diversity can be related to the evolution of the species and the discordance between the morphological and molecular structures may result from similar selection pressures at different places leading to similar forms with a different genetic background.