Whole-genome Resequencing Data Research Articles

Evolution of chromosomal inversions across an avian radiation.

Chromosomal inversions are structural mutations that can play a prominent role in adaptation and speciation. Inversions segregating across species boundaries (trans-species inversions) are often taken as evidence for ancient balancing selection or adaptive introgression but can also be due to incomplete lineage sorting (ILS). Using whole-genome resequencing data from 18 populations of 11 recognized munia species in the genus Lonchura (N = 176 individuals), we identify four large para- and pericentric inversions ranging in size from 4 to 20 Mb. All four inversions co-segregate across multiple species and predate the numerous speciation events associated with the rapid radiation of this clade across the prehistoric Sahul (Australia, New Guinea) and Bismarck Archipelago. Using coalescent theory, we infer that trans-specificity is improbable for neutrally segregating variation despite substantial ILS characterizing this young radiation. Instead, maintenance of all three autosomal inversions (chr1, chr5, chr6) is best explained by selection acting along eco-geographic clines not observed for the collinear parts of the genome. In addition, the large sex chromosome inversion largely aligns with species boundaries and shows signatures of repeated positive selection for both alleles. This study provides evidence for trans-species inversion polymorphisms involved in both adaptation and speciation. It further highlights the importance of informing selection inference using a null model of neutral evolution derived from the collinear part of the genome.

Sex chromosome turnover in hybridizing stickleback lineages

Abstract Recent discoveries of sex chromosome diversity across the tree of life have challenged the canonical model of conserved sex chromosome evolution and evoked new theories on labile sex chromosomes that maintain less differentiation and undergo frequent turnover. However, theories of labile sex chromosome evolution lack direct empirical support due to the paucity of case studies demonstrating ongoing sex chromosome turnover in nature. Two divergent lineages (viz. WL & EL) of nine-spined sticklebacks (Pungitius pungitius) with different sex chromosomes (linkage group [LG] 12 in the EL, unknown in the WL) hybridize in a natural secondary contact zone in the Baltic Sea, providing an opportunity to study ongoing turnover between coexisting sex chromosomes. In this study, we first identify an 80 kbp genomic region on LG3 as the sex-determining region (SDR) using whole-genome resequencing data of family crosses of a WL population. We then verify this region as the SDR in most other WL populations and demonstrate a potentially ongoing sex chromosome turnover in admixed marine populations where the evolutionarily younger and homomorphic LG3 sex chromosome replaces the older and heteromorphic LG12 sex chromosome. The results provide a rare glimpse of sex chromosome turnover in the wild and indicate the possible existence of additional yet undiscovered sex chromosome diversity in Pungitius sticklebacks.

Whole-genome resequencing revealed the population structure and selection signal of 4 indigenous Chinese laying ducks

The assessment of animal genetic structure had significant importance for the preservation and breeding of animal germplasm resources. Selection signals are genotype markers generated during the process of biological evolution, and the detection of selection signals could reveal the direction of species evolution. The aim of this study was to generate a whole-genome resequencing data from Jinding duck, Shanma duck, Youxian Partridge duck, and Taiwan Brown tsaiya duck to reveal their population structure and selection signals. The population structure analysis revealed significant genetic differences among the 4 indigenous laying ducks, indicating their independent lineage. Specifically, Shanma duck and Youxian partridge duck were closely and likely originated from a common ancestor. In addition, selection sweep analysis was performed using the population genetic differentiation coefficient (Fst) and nucleotide diversity ratio (π ratio). The top 5% was used as the threshold for the Fst and π ratio, and the 2 thresholds were combined to identify selected genomic regions. In the selected regions of the 3 comparison groups, 136, 143, and 268 candidate genes were detected. Further screening of all candidate genes revealed that 35 candidate genes appeared simultaneously in 3 comparative groups, with 16 genes annotated. The 16 genes were analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. The results revealed 5 functional genes (AQP3, PIK3C3, NOL6, RPP25, and DCTN3) that may be related to important economic traits in laying ducks and involved mainly invasopressin-regulated water reabsorption, ribosome biogenesis, and the PI3K signaling pathway. The results provide insights into the protection and exploitation of genetic resources of Chinese indigenous laying ducks.

Assessing the population genetic structure of yellow croaker in China: Insights into the ecological and genetic consequences of artificial breeding on natural populations

The large yellow croaker, a crucial marine economic fish in China, experienced a drastic decline in wild resource in the early 1980s due to overfishing and habitat degradation. Since the year 2000, the government has implemented habitat protection measures and initiated long-term, large-scale releases aimed at restoring this ecologically and economically important fish resource. Despite these efforts, there has been no definitive success in increasing the biomass of this overexploited species, and the underlying reasons for this lack of success remain elusive. Leveraging whole-genome resequencing data from a large number of samples encompassing both cultured and wild populations, this study offers the most exhaustive analysis of the genetic structure of large yellow croaker to date. The results revealed no significant geographical structure among the wild populations along the coast of China, overturning the long-held belief that designated them into three genetic management units. Additionally, it highlighted the absence of genetic introgression from cultured to wild populations. Furthermore, this study identified adaptive changes in cultured populations due to artificial selection and aquaculture environments, potentially reducing their adaptability to natural environments. This reduced adaptability may impact the survival rate of released populations and their hybridization with wild populations, providing insights into the limited success of previous release efforts. In light of these findings, this study recommends implementing conservation and management strategies based on a single genetic management unit for the large yellow croaker. When selecting breeding populations for hatchery releases, emphasis should be placed on genetic background considerations rather than solely relying on the origin of individuals captured in the sea. These insights hold significant implications for conservation and restoration of the large yellow croaker resource, particularly in the management of germplasm resources.

Resequencing of Brazilian locally adapted cattle breeds revealed variants in candidate genes and transcription factors for meat fatty acid profile.

The beef cattle industry has experienced a shift driven by a market demand for healthier meat, cost efficiency and environmental sustainability in recent years. Consequently, there has been a growing focus on the fatty acids content and functions of meat in cattle breeding programmes. Besides, a deeper understanding of the biological mechanisms influencing the expression of different phenotypes related to fatty acid profiles is crucial. In this study, we aimed to identify Single-Nucleotide Variants (SNV) and Insertion/Deletion (InDels) DNA variants in candidate genes related to fatty acid profiles described in genomic, transcriptomic and proteomic studies conducted in beef cattle breeds. Utilizing whole-genome re-sequencing data from Brazilian locally adapted bovine breeds, namely Caracu and Pantaneiro, we identified SNVs and InDels associated with 23,947 genes. From these, we identified 318 candidate genes related to fatty acid profiles that contain variants. Subsequently, we select only genes with SNVs and InDels in their promoter, 5' UTR and coding region. Through the gene-biological process network, approximately 19 genes were highlighted. Furthermore, considering the studied trait and a literature review, we selected the main transcription factors (TF). Functional analysis via gene-TF network allowed us to identify the 30 most likely candidate genes for meat fatty acid profile in cattle. LIPE, MFSD2A and SREBF1 genes were highlighted in networks due to their biological importance. Further dissection of these genes revealed 15 new variants found in promoter regions of Caracu and Pantaneiro sequences. The gene networks facilitated a better functional understanding of genes and TF, enabling the identification of variants potentially related to the expression of candidate genes for meat fatty acid profiles in cattle.

Genomic Inbreeding and Runs of Homozygosity Analysis of Cashmere Goat.

Cashmere goats are valuable genetic resources which are famous worldwide for their high-quality fiber. Runs of homozygosity (ROHs) have been identified as an efficient tool to assess inbreeding level and identify related genes under selection. However, there is limited research on ROHs in cashmere goats. Therefore, we investigated the ROH pattern, assessed genomic inbreeding levels and examined the candidate genes associated with the cashmere trait using whole-genome resequencing data from 123 goats. Herein, the Inner Mongolia cashmere goat presented the lowest inbreeding coefficient of 0.0263. In total, we identified 57,224 ROHs. Seventy-four ROH islands containing 50 genes were detected. Certain identified genes were related to meat, fiber and milk production (FGF1, PTPRM, RERE, GRID2, RARA); fertility (BIRC6, ECE2, CDH23, PAK1); disease or cold resistance and adaptability (PDCD1LG2, SVIL, PRDM16, RFX4, SH3BP2); and body size and growth (TMEM63C, SYN3, SDC1, STRBP, SMG6). 135 consensus ROHs were identified, and we found candidate genes (FGF5, DVL3, NRAS, KIT) were associated with fiber length or color. These findings enhance our comprehension of inbreeding levels in cashmere goats and the genetic foundations of traits influenced by selective breeding. This research contributes significantly to the future breeding, reservation and use of cashmere goats and other goat breeds.

Dissecting Selective Signatures and Candidate Genes in Grandparent Lines Subject to High Selection Pressure for Broiler Production and in a Local Russian Chicken Breed of Ushanka.

Breeding improvements and quantitative trait genetics are essential to the advancement of broiler production. The impact of artificial selection on genomic architecture and the genetic markers sought remains a key area of research. Here, we used whole-genome resequencing data to analyze the genomic architecture, diversity, and selective sweeps in Cornish White (CRW) and Plymouth Rock White (PRW) transboundary breeds selected for meat production and, comparatively, in an aboriginal Russian breed of Ushanka (USH). Reads were aligned to the reference genome bGalGal1.mat.broiler.GRCg7b and filtered to remove PCR duplicates and low-quality reads using BWA-MEM2 and bcftools software; 12,563,892 SNPs were produced for subsequent analyses. Compared to CRW and PRW, USH had a lower diversity and a higher genetic distinctiveness. Selective sweep regions and corresponding candidate genes were examined based on ZFST, hapFLK, and ROH assessment procedures. Twenty-seven prioritized chicken genes and the functional projection from human homologs suggest their importance for selection signals in the studied breeds. These genes have a functional relationship with such trait categories as body weight, muscles, fat metabolism and deposition, reproduction, etc., mainly aligned with the QTLs in the sweep regions. This information is pivotal for further executing genomic selection to enhance phenotypic traits.

Influence of breed and environment on leukocyte telomere length in cattle.

High milk yield is associated with reduced longevity in high-producing dairy cattle breeds. Pre-term culling leads to high replacement heifer demand and economic losses for the dairy industry. Selection for this trait is limited because of low heritability and difficulties in phenotype measurement. Telomeres are elements found at the ends of chromosomes, consisting of repetitive DNA sequences, several thousand base pairs in length, coupled with nucleoprotein complexes. Eventually, in humans and most other animals, telomere length reduces with age. When telomeric DNA is truncated to a critical length, cell ageing, cell cycle arrest, and apoptosis are induced. As a result, telomere length can be considered as a predictor of health risks and an individual's lifespan. The leukocyte telomere length may be used as a proxy phenotype of productive lifespan to improve cattle selection. Our objectives were to assess the effects of breed and breed group (dairy vs. beef) on the leukocyte telomere length and to estimate the effect of cold climate on this trait in Kalmyk cattle populations from the South (Rostov Oblast) and Far North (Republic of Sakha) regions of Russia. The leukocyte telomere lengths were estimated computationally from whole-genome resequencing data. We leveraged data on leukocyte telomere length, sex, and age of 239 animals from 17 cattle breeds. The breed factor had a significant effect on leukocyte telomere length across our sample. There was no difference in leukocyte telomere length between dairy and beef groups. The population factor had a significant effect on leukocyte telomere length in Kalmyk animals. In conclusion, we found that breed, but not breed group (dairy vs. beef), was significantly associated with leukocyte telomere length in cattle. Residence in colder climates was associated with longer leukocyte telomere length in Kalmyk breed cattle.

Study on the origin of the Baise horse based on whole-genome resequencing.

The Baise horse, an indigenous horse breed mainly distributed in the Baise region of Guangxi province in southwest China, has a long history as draft animal. However, there is a lack of research regarding the origin and ancestral composition of the Baise horse. In this study, whole-genome resequencing data from 236 horses of seven Chinese indigenous horse breeds, five foreign horse breeds, and four Przewalski's horses were used to investigate the relationships between the Baise horse and other horse breeds. The results showed that foreign horse breeds had no significant impact on the formation of the Baise horse. The two southwestern horse populations, the Debao pony and the Jinjiang horse, exhibit the closest genetic affinity with the Baise horse. This is consistent with their adjacent geographical distribution. Analysis of the migration route revealed a gene flow from the Chakouyi horse into the Baise horse. In summary, our results confirm that the formation of the Baise horse did not involve participation from foreign breeds. Geographical distance emerges as a crucial factor in determining the genetic relationships with the Baise horse. Gene flows of indigenous horse breeds along ancient routes of trade activities had played a role in the formation of the Baise horse.

Whole-mitogenome analysis unveils previously undescribed genetic diversity in cane toads across their invasion trajectory.

Invasive species offer insights into rapid adaptation to novel environments. The iconic cane toad (Rhinella marina) is an excellent model for studying rapid adaptation during invasion. Previous research using the mitochondrial NADH dehydrogenase 3 (ND3) gene in Hawai'ian and Australian invasive populations found a single haplotype, indicating an extreme genetic bottleneck following introduction. Nuclear genetic diversity also exhibited reductions across the genome in these two populations. Here, we investigated the mitochondrial genomics of cane toads across this invasion trajectory. We created the first reference mitochondrial genome for this species using long-read sequence data. We combined whole-genome resequencing data of 15 toads with published transcriptomic data of 125 individuals to construct nearly complete mitochondrial genomes from the native (French Guiana) and introduced (Hawai'i and Australia) ranges for population genomic analyses. In agreement with previous investigations of these populations, we identified genetic bottlenecks in both Hawai'ian and Australian introduced populations, alongside evidence of population expansion in the invasive ranges. Although mitochondrial genetic diversity in introduced populations was reduced, our results revealed that it had been underestimated: we identified 45 mitochondrial haplotypes in Hawai'ian and Australian samples, none of which were found in the native range. Additionally, we identified two distinct groups of haplotypes from the native range, separated by a minimum of 110 base pairs (0.6%). These findings enhance our understanding of how invasion has shaped the genetic landscape of this species.

Genomic analyses elucidate S-locus evolution in response to intra-specific losses of distyly in Primula vulgaris.

Distyly, a floral dimorphism that promotes outcrossing, is controlled by a hemizygous genomic region known as the S-locus. Disruptions of genes within the S-locus are responsible for the loss of distyly and the emergence of homostyly, a floral monomorphism that favors selfing. Using whole-genome resequencing data of distylous and homostylous individuals from populations of Primula vulgaris and leveraging high-quality reference genomes of Primula we tested, for the first time, predictions about the evolutionary consequences of transitions to selfing on S-genes. Our results reveal a previously undetected structural rearrangement in CYPᵀ associated with the shift to homostyly and confirm previously reported, homostyle-specific, loss-of-function mutations in the exons of the S-gene CYPᵀ. We also discovered that the promoter and intronic regions of CYPᵀ in distylous and homostylous individuals are conserved, suggesting that down-regulation of CYPᵀ via mutations in its promoter and intronic regions is not a cause of the shift to homostyly. Furthermore, we found that hemizygosity is associated with reduced genetic diversity in S-genes compared with their paralogs outside the S-locus. Additionally, the shift to homostyly lowers genetic diversity in both the S-genes and their paralogs, as expected in primarily selfing plants. Finally, we tested, for the first time, long-standing theoretical models of changes in S-locus genotypes during early stages of the transition to homostyly, supporting the assumption that two copies of the S-locus might reduce homostyle fitness.

Detection of two homologous major QTLs and development of diagnostic molecular markers for sucrose content in peanut.

We identified two stable and homologous major QTLs for sucrose content in peanut, and developed breeder-friendly molecular markers for marker-assisted selection breeding. Sucrose content is a crucial quality trait for edible peanuts, and increasing sucrose content is a key breeding objective. However, the genetic basis of sucrose content in peanut remains unclear, and major quantitative trait loci (QTLs) for sucrose content have yet to be identified. In this study, a high-density genetic map was constructed based on whole-genome re-sequencing data from a peanut RIL population. This map consisted of 2,042 bins and 24,142 SNP markers, making it one of the most comprehensive maps to date in terms of marker density. Two major QTLs (qSCA06.2 and qSCB06.2) were identified, explaining 31.41% and 24.13% of the phenotypic variance, respectively. Notably, these two QTLs were located in homologous genomic regions between the A and B subgenomes. The elite allele of qSCA06.2 was exclusive to Valencia-type, while the elite allele of qSCB06.2 existed in other peanut types. Importantly, the distribution of alleles from two homologous QTLs in the RIL population and diverse germplasm accessions consistently demonstrated that only the combination of elite allelic genotypes from both QTLs/genes resulted in a significantly dominant phenotype, accompanied by a substantial increase in sucrose content. The newly developed diagnostic markers for these QTLs were confirmed to be reliable and could facilitate future breeding efforts to enhance sucrose content using marker-assisted selection techniques. Overall, this study highlights the co-regulation of sucrose content by two major homologous QTLs/genes and provides valuable insights into the genetic basis of sucrose in peanuts.

Whole-Genome Scanning for Selection Signatures Reveals Candidate Genes Associated with Growth and Tail Length in Sheep.

Compared to Chinese indigenous sheep, Western sheep have rapid growth rate, larger physique, and higher meat yield. These excellent Western sheep were introduced into China for crossbreeding to expedite the enhancement of production performance and mutton quality in local breeds. Here, we investigated population genetic structure and genome-wide selection signatures among the Chinese indigenous sheep and the introduced sheep based on whole-genome resequencing data. The PCA, N-J tree and ADMIXTURE results showed significant genetic difference between Chinese indigenous sheep and introduced sheep. The nucleotide diversity (π) and linkage disequilibrium (LD) decay results indicated that the genomic diversity of introduced breeds were lower. Then, Fst & π ratio, XP-EHH, and de-correlated composite of multiple signals (DCMS) methods were used to detect the selection signals. The results showed that we identified important candidate genes related to growth rate and body size in the introduced breeds. Selected genes with stronger selection signatures are associated with growth rate (CRADD), embryonic development (BVES, LIN28B, and WNT11), body size (HMGA2, MSRB3, and PTCH1), muscle development and fat metabolism (MSTN, PDE3A, LGALS12, GGPS1, and SAR1B), wool color (ASIP), and hair development (KRT71, KRT74, and IRF2BP2). Thus, these genes have the potential to serve as candidate genes for enhancing the growth traits of Chinese indigenous sheep. We also identified tail-length trait-related candidate genes (HOXB13, LIN28A, PAX3, and VEGFA) in Chinese long-tailed breeds. Among these genes, HOXB13 is the main candidate gene for sheep tail length phenotype. LIN28A, PAX3, and VEGFA are related to embryonic development and angiogenesis, so these genes may be candidate genes for sheep tail type traits. This study will serve as a foundation for further genetic improvement of Chinese indigenous sheep and as a reference for studies related to growth and development of sheep.

Genome sequences and population genomics reveal climatic adaptation and genomic divergence between two closely related sweetgum species.

Understanding the genetic basis of population divergence and adaptation is an important goal in population genetics and evolutionary biology. However, the relative roles of demographic history, gene flow, and/or selective regime in driving genomic divergence, climatic adaptation, and speciation in non-model tree species are not yet fully understood. To address this issue, we generated whole-genome resequencing data of Liquidambar formosana and L. acalycina, which are broadly sympatric but altitudinally segregated in the Tertiary relict forests of subtropical China. We integrated genomic and environmental data to investigate the demographic history, genomic divergence, and climatic adaptation of these two sister species. We inferred a scenario of allopatric species divergence during the late Miocene, followed by secondary contact during the Holocene. We identified multiple genomic islands of elevated divergence that mainly evolved through divergence hitchhiking and recombination rate variation, likely fostered by long-term refugial isolation and recent differential introgression in low-recombination genomic regions. We also found some candidate genes with divergent selection signatures potentially involved in climatic adaptation and reproductive isolation. Our results contribute to a better understanding of how late Tertiary/Quaternary climatic change influenced speciation, genomic divergence, climatic adaptation, and introgressive hybridization in East Asia's Tertiary relict flora. In addition, they should facilitate future evolutionary, conservation genomics, and molecular breeding studies in Liquidambar, a genus of important medicinal and ornamental values.

Conservation genomics analysis reveals recent population decline and possible causes in bumblebee Bombus opulentus.

Bumblebees are a genus of pollinators (Bombus) that play important roles in natural ecosystem and agricultural production. Several bumblebee species have been recorded as under population decline, and the proportion of species experiencing population decline within subgenus Thoracobombus is higher than average. Bombus opulentus is 1 species in Thoracobombus, but little is known about its recent population dynamics. Here, we employed conservation genomics methods to investigate the population dynamics of B. opulentus during the recent past and identify the likely environmental factors that may cause population decline. Firstly, we placed the scaffold-level of B. opulentus reference genome sequence onto chromosome-level using Hi-C technique. Then, based on this reference genome and whole-genome resequencing data for 51 B. opulentus samples, we reconstructed the population structure and effective population size (Ne ) trajectories of B. opulentus and identified genes that were under positive selection. Our results revealed that the collected B. opulentus samples could be divided into 2 populations, and 1 of them experienced a recent population decline; the declining population also exhibited lower genetic diversity and higher inbreeding levels. Genes related to high-temperature tolerance, immune response, and detoxication showed signals of positive selection in the declining population, suggesting that climate warming and pathogen/pesticide exposures may contribute to the decline of this B. opulentus population. Taken together, our study provided insights into the demography of B. opulentus populations and highlighted that populations of the same bumblebee species could have contrasting Ne trajectories and population decline could be caused by a combination of various stressors.

Morpho-physiological, Genomic, and Transcriptional Diversities in Response to Potassium Deficiency in Rapeseed (Brassica napus L.) Genotypes.

Variations in the resistance to potassium (K) deficiency among rapeseed genotypes emphasize complicated regulatory mechanisms. In this study, a low-K-sensitivity accession (L49) responded to K deficiency with smaller biomasses, severe leaf chlorosis, weaker photosynthesis ability, and deformed stomata morphology compared to a low-K resistant accession (H280). H280 accumulated more K+ than L49 under low K. Whole-genome resequencing (WGS) revealed a total of 5,538,622 single nucleotide polymorphisms (SNPs) and 859,184 insertions/deletions (InDels) between H280 and L49. RNA-seq identified more differentially expressed K+ transporter genes with higher expression in H280 than in L49 under K deficiency. Based on the K+ profiles, differential expression profiling, weighted gene coexpression network analysis, and WGS data between H280 and L49, BnaC4.AKT1 was proposed to be mainly responsible for root K absorption-mediated low K resistance. BnaC4.AKT1 was expressed preferentially in the roots and localized on the plasma membrane. An SNP and an InDel found in the promoter region of BnaC4.AKT1 were proposed to be responsible for its differential expression between rapeseed genotypes. This study identified a gene resource for improving low-K resistance. It also facilitates an integrated knowledge of the differential physiological and transcriptional responses to K deficiency in rapeseed genotypes.

Divergent Selection in Low Recombination Regions Shapes the Genomic Islands in Two Incipient Shorebird Species.

Speciation in the face of gene flow is usually associated with a heterogeneous genomic landscape of divergence in nascent species pairs. However, multiple factors, such as divergent selection and local recombination rate variation, can influence the formation of these genomic islands. Examination of the genomic landscapes of species pairs that are still in the early stages of speciation provides an insight into this conundrum. In this study, population genomic analyses were undertaken using a wide range of sampling and whole-genome resequencing data from 96 unrelated individuals of Kentish plover (Charadrius alexandrinus) and white-faced plover (Charadrius dealbatus). We suggest that the two species exhibit varying levels of population admixture along the Chinese coast and on the Taiwan Island. Genome-wide analyses for introgression indicate that ancient introgression had occurred in Taiwan population, and gene flow is still ongoing in mainland coastal populations. Furthermore, we identified a few genomic regions with significant levels of interspecific differentiation and local recombination suppression, which contain several genes potentially associated with disease resistance, coloration, and regulation of plumage molting and thus may be relevant to the phenotypic and ecological divergence of the two nascent species. Overall, our findings suggest that divergent selection in low recombination regions may be a main force in shaping the genomic islands in two incipient shorebird species.

Ecological Diversification in an Adaptive Radiation of Plants: The Role of De Novo Mutation and Introgression.

Adaptive radiations are characterized by rapid ecological diversification and speciation events, leading to fuzzy species boundaries between ecologically differentiated species. Adaptive radiations are therefore key systems for understanding how species are formed and maintained, including the role of de novo mutations versus preexisting variation in ecological adaptation and the genome-wide consequences of hybridization events. For example, adaptive introgression, where beneficial alleles are transferred between lineages through hybridization, may fuel diversification in adaptive radiations and facilitate adaptation to new environments. In this study, we employed whole-genome resequencing data to investigate the evolutionary origin of hummingbird-pollinated flowers and to characterize genome-wide patterns of phylogenetic discordance and introgression in Penstemon subgenus Dasanthera, a small and diverse adaptive radiation of plants. We found that magenta hummingbird-adapted flowers have apparently evolved twice from ancestral blue-violet bee-pollinated flowers within this radiation. These shifts in flower color are accompanied by a variety of inactivating mutations to a key anthocyanin pathway enzyme, suggesting that independent de novo loss-of-function mutations underlie the parallel evolution of this trait. Although patterns of introgression and phylogenetic discordance were heterogenous across the genome, a strong effect of gene density suggests that, in general, natural selection opposes introgression and maintains genetic differentiation in gene-rich genomic regions. Our results highlight the importance of both de novo mutation and introgression as sources of evolutionary change and indicate a role for de novo mutation in driving parallel evolution in adaptive radiations.

Genome-wide discovery of InDels and validation of PCR-Based InDel markers for earliness in a RIL population and genotypes of lentil (Lens culinaris Medik.).

The systematic identification of insertion/deletion (InDel) length polymorphisms from the entire lentil genome can be used to map the quantitative trait loci (QTL) and also for the marker-assisted selection (MAS) for various linked traits. The InDels were identified by comparing the whole-genome resequencing (WGRS) data of two extreme bulks (early- and late-flowering bulk) and a parental genotype (Globe Mutant) of lentil. The bulks were made by pooling 20 extreme recombinant inbred lines (RILs) each, derived by crossing Globe Mutant (late flowering parent) with L4775 (early flowering parent). Finally, 734,716 novel InDels were identified, which is nearly one InDel per 5,096 bp of lentil genome. Furthermore, 74.94% of InDels were within the intergenic region and 99.45% displayed modifier effects. Of these, 15,732 had insertions or deletions of 20 bp or more, making them amenable to the development of PCR-based markers. An InDel marker I-SP-356.6 (chr. 3; position 356,687,623; positioned 174.5 Kb from the LcFRI gene) was identified as having a phenotypic variance explained (PVE) value of 47.7% for earliness when validated in a RIL population. Thus, I-SP-356.6 marker can be deployed in MAS to facilitate the transfer of the earliness trait to other elite late-maturing cultivars. Two InDel markers viz., I-SP-356.6 and I-SP-383.9 (chr. 3; linked to LcELF3a gene) when tested in 9 lentil genotypes differing for maturity duration, clearly distinguished three early (L4775, ILL7663, Precoz) and four late genotypes (Globe Mutant, MFX, L4602, L830). However, these InDels could not be validated in two genotypes (L4717, L4727), suggesting either absence of polymorphism and/or presence of other loci causing earliness. The identified InDel markers can act as valuable tools for MAS for the development of early maturing lentil varieties.

Cuticular proteins in codling moth (Cydia pomonella) respond to insecticide and temperature stress

The insect cuticle consists of chitin and cuticular proteins (CPs), which stabilize the body shape and provide an effective physical barrier against the external environment. They are also potential target sites for developing environmentally friendly insect management through the utilization of physiology-based methods. The codling moth, Cydia pomonella, is a pest afflicting fruit orchards worldwide. This study used a comparative genomic approach, whole-genome resequencing, and transcriptome data to understand the role that CPs played in the environmental adaptation of the codling moth. A total of 182 putative CPs were identified in the codling moth genome, which were classified into 12 CP families. 119 CPR genes, including 54 RR-1, 60 RR-2, and 5 RR-3 genes were identified and accounted for 65.4% of the total CPs. Eight and seven gene clusters are formed in RR1 and RR2 subfamily and the ancestor-descendant relationship was explained. Five CPAP genes were highly expressed during the egg stage and exposed to high temperature, which indicated their potential role in aiding codling moth eggs in acclimating to varying external heat conditions. Moreover, six CPs belonging to the CPR and CPLCP families were identified in association with insecticide resistance by population resequencing. Their expression levels increased after exposure to insecticides, suggesting they might be involved in codling moth resistance to the insecticides azinphos-methyl or deltamethrin. Our results provide insight into the evolution of codling moth CPs and their association with high temperature adaptation and insecticide resistance, and provide an additional information required for further analysis of CPs in environmental adaptation.